High-entropy engineering of the crystal and electronic structures in a Dirac material.

Dirac and Weyl semimetals are a central topic of contemporary condensed matter physics, and the discovery of new compounds with Dirac/Weyl electronic states is crucial to the advancement of topological materials and quantum technologies. Here we show a widely applicable strategy that uses high configuration entropy to engineer relativistic electronic states. We take the AMnSb2 (A = Ba, Sr, Ca, Eu, and Yb) Dirac material family as an example and demonstrate that mixing of Ba, Sr, Ca, Eu and Yb at the A site generates the compound (Ba0.38Sr0.14Ca0.16Eu0.16Yb0.16)MnSb2 (denoted as A5MnSb2), giving access to a polar structure with a space group that is not present in any of the parent compounds. A5MnSb2 is an entropy-stabilized phase that preserves its linear band dispersion despite considerable lattice disorder. Although both A5MnSb2 and AMnSb2 have quasi-two-dimensional crystal structures, the two-dimensional Dirac states in the pristine AMnSb2 evolve into a highly anisotropic quasi-three-dimensional Dirac state triggered by local structure distortions in the high-entropy phase, which is revealed by Shubnikov-de Haas oscillations measurements.

Simulation of Electrochemical Oxidation in Aqueous Environments under Applied Voltage Using Classical Molecular Dynamics.

Reactive molecular dynamics (MD) simulations of metal electrodes under an applied voltage in an explicit water environment were performed and compared to predictions from both other calculations and simulations and experimental measurements and observations. MD simulations using the third-generation charge-optimized many body (COMB3) potentials and the electrode COMB (eCOMB) approach allow for the simulation of an externally applied voltage by modifying the equations of motion during the charge equilibration step (QEq) of the MD simulation. Unlike previous work, which prevented charge transfer between the water and metal electrodes, this work coupled the water and metal together through the QEq, which leads to an accumulation of a negative charge on the water and a positive charge on the metal before any voltage is applied. The extent of this charge accumulation is characterized, and strategies to mitigate it are explored. Root mean square deviation plots of the metal surfaces are created for each simulation to compare the extent of oxidation and dissolution. Good agreement is found between the oxidation behavior of the simulated metal surfaces and the experimental observations.

Strong electron-phonon coupling driven pseudogap modulation and density-wave fluctuations in a correlated polar metal.

There is tremendous interest in employing collective excitations of the lattice, spin, charge, and orbitals to tune strongly correlated electronic phenomena. We report such an effect in a ruthenate, Ca3Ru2O7, where two phonons with strong electron-phonon coupling modulate the electronic pseudogap as well as mediate charge and spin density wave fluctuations. Combining temperature-dependent Raman spectroscopy with density functional theory reveals two phonons, B2P and B2M, that are strongly coupled to electrons and whose scattering intensities respectively dominate in the pseudogap versus the metallic phases. The B2P squeezes the octahedra along the out of plane c-axis, while the B2M elongates it, thus modulating the Ru 4d orbital splitting and the bandwidth of the in-plane electron hopping; Thus, B2P opens the pseudogap, while B2M closes it. Moreover, the B2 phonons mediate incoherent charge and spin density wave fluctuations, as evidenced by changes in the background electronic Raman scattering that exhibit unique symmetry signatures. The polar order breaks inversion symmetry, enabling infrared activity of these phonons, paving the way for coherent light-driven control of electronic transport.

koopmans: An Open-Source Package for Accurately and Efficiently Predicting Spectral Properties with Koopmans Functionals.

Over the past decade we have developed Koopmans functionals, a computationally efficient approach for predicting spectral properties with an orbital-density-dependent functional framework. These functionals impose a generalized piecewise linearity condition to the entire electronic manifold, ensuring that orbital energies match the corresponding electron removal/addition energy differences (in contrast to semilocal DFT, where a mismatch between the two lies at the heart of the band gap problem and, more generally, the unreliability of Kohn-Sham orbital energies). This strategy has proven to be very powerful, yielding molecular orbital energies and solid-state band structures with comparable accuracy to many-body perturbation theory but at greatly reduced computational cost while preserving a functional formulation. This paper reviews the theory of Koopmans functionals, discusses the algorithms necessary for their implementation, and introduces koopmans, an open-source package that contains all of the code and workflows needed to perform Koopmans functional calculations and obtain reliable spectral properties of molecules and materials.

Atomic-scale polarization switching in wurtzite ferroelectrics.

Ferroelectric wurtzites have the potential to revolutionize modern microelectronics because they are easily integrated with multiple mainstream semiconductor platforms. However, the electric fields required to reverse their polarization direction and unlock electronic and optical functions need substantial reduction for operational compatibility with complementary metal-oxide semiconductor (CMOS) electronics. To understand this process, we observed and quantified real-time polarization switching of a representative ferroelectric wurtzite (Al0.94B0.06N) at the atomic scale with scanning transmission electron microscopy. The analysis revealed a polarization reversal model in which puckered aluminum/boron nitride rings in the wurtzite basal planes gradually flatten and adopt a transient nonpolar geometry. Independent first-principles simulations reveal the details and energetics of the reversal process through an antipolar phase. This model and local mechanistic understanding are a critical initial step for property engineering efforts in this emerging material class.

Band Gap Narrowing in a High-Entropy Spinel Oxide Semiconductor for Enhanced Oxygen Evolution Catalysis.

High-entropy oxides (HEOs), which contain five or more metal cations that are generally thought to be randomly mixed in a crystalline oxide lattice, can exhibit unique and enhanced properties, including improved catalytic performance, due to synergistic effects. Here, we show that band gap narrowing emerges in a high-entropy aluminate spinel oxide, (Fe0.2Co0.2Ni0.2Cu0.2Zn0.2)Al2O4 (A5Al2O4). The 0.9 eV band gap of A5Al2O4 is narrower than the band gaps of all parent spinel oxides. First-principles calculations for multicomponent AAl2O4 spinels indicate that the band gap narrowing arises from the broadening of the energy distribution of the 3d states due to variations in the electronegativities and crystal field splitting across the 3d transition-metal series. As a catalyst for the oxygen evolution reaction in an alkaline electrolyte, A5Al2O4 reaches a current density of 10 mA/cm2 at an overpotential of 400 mV, outperforming all of the single-metal end members at an applied potential of 1.7 V vs RHE. Catalyst deactivation occurs after 5 h at 10 mA/cm2 and is attributed, based on elemental analysis and grazing-incidence X-ray diffraction, to the formation of a passivating layer that blocks the high-entropy oxide surface. This result helps to validate that the HEO is the active catalyst. The observation of band gap narrowing in A5Al2O4 expands the scope of synergistic properties exhibited by high-entropy materials and offers insight into the question of how the electronic structure of multicomponent oxide materials can be engineered via a high-entropy approach to achieve enhanced catalytic properties.

Using C-DFT to develop an e-ReaxFF force field for acetophenone radical anion.

Increased electricity usage over the past several decades has accelerated the need for efficient high-voltage power transmission with reliable insulating materials. Cross-linked polyethylene (XLPE) prepared via dicumyl peroxide (DCP) cross-linking has emerged as the insulator of choice for modern power cables. Although DCP cross-linking generates the desired XLPE product in high yield, other by-products are also produced. One such by-product, acetophenone, is particularly intriguing due to its aromaticity and positive electron affinity. In this work, constrained density functional theory (C-DFT) was utilized to develop an e-ReaxFF force field suitable for describing the acetophenone radical anion. Initial parameters were taken from the 2021 Akbarian e-ReaxFF force field, which was developed to describe XLPE chemistry. Then, C-DFT geometry optimizations were performed wherein an excess electron was constrained to each atom of acetophenone. The resulting C-DFT energy values for the various electronic positions were added to the e-ReaxFF training set. Next, an analogous set of structures was energy-minimized using e-ReaxFF, and equilibrium mixture compositions for the two methods were compared at multiple temperatures. Iterative fitting against C-DFT energy data was performed until satisfactory agreement was achieved. To test force field performance, molecular dynamics simulations were performed in e-ReaxFF and the resulting electronic distributions were qualitatively compared to unconstrained-DFT spin density data. By expanding our e-ReaxFF force field for XLPE, namely, adding the capability to describe acetophenone and its interactions with an excess electron, we move one step closer to a comprehensive molecular understanding of XLPE chemistry in a high-voltage power cable.

Environmental impact of amino acids on selenate-bearing hydrocalumite: Experimental and DFT studies.

Selenium (Se) radioactive wastes can be disposed through stabilization/solidification (S/S) based on the cementitious matrix on hydration products, where hydrocalumite (Ca2Al-LDH) is expected to play an important role in the retention of SeO42-. Natural organic matters (NOMs) are known to be a risk to affect the transportation and mobility of undesirable chemical species in the pedosphere which receives the low level radioactive wastes (LLW). In the present work, five amino acids were selected as the simplified models of NOMs in the pedosphere to explore their effects on the stability of Ca2Al-LDH after immobilized SeO42- under alkaline conditions. As the loading amount of amino acids on Ca2Al-LDH increasing, release of SeO42- was enhanced in HGly, H2Asp, and H2Cys series, while no enhancement was observed in HPhe and HTrp series. Density functional theory (DFT) calculation predicted ion-exchange of amino acids and CO32- with SeO42- in a unit cell of LDH model. The intercalation of Asp2- and CO32- caused 003 peaks in XRD sharper and d003 decreased from 8.15 Å to 7.70 Å which is assigned to Ca2Al-LDH(Asp, CO3). In H2Cys series, the 003 peaks were kept broad and SeO42- was still relatively maintained in LDH which was caused by the lower amounts of intercalated CO32- in the presence of H2Cys. Amino acids in the interlayer of Ca2Al-LDH have several possible configurations, where the most stable one is prone to be in a horizontal direction through hydrogen bonds and Ca-O chemical bonds. This provides an insight on the stability of selenate immobilized in hydrocalumite, which can be produced in cement disposing in the pedosphere for a long term of burying. Not only carbonate but also small molecular organic matters like amino acids possibly give environmental impact on the mobility of low level anionic radionuclides in LDH.

Single-Step Direct Laser Writing of Multimetal Oxygen Evolution Catalysts from Liquid Precursors.

We investigate a laser direct-write method to synthesize and deposit metastable, mixed transition metal oxides and evaluate their performance as oxygen evolution reaction catalysts. This laser processing method enabled the rapid synthesis of diverse heterogeneous alloy and oxide catalysts directly from cost-effective solution precursors, including catalysts with a high density of nanocrystalline metal alloy inclusions within an amorphous oxide matrix. The nanoscale heterogeneous structures of the synthesized catalysts were consistent with reactive force-field Monte Carlo calculations. By evaluating the impact of varying transition metal oxide composition ratios, we created a stable Fe0.63Co0.19Ni0.18Ox/C catalyst with a Tafel slope of 38.23 mV dec-1 and overpotential of 247 mV, a performance similar to that of IrO2. Synthesized Fe0.63Co0.19Ni0.18Ox/C and Fe0.14Co0.46Ni0.40Ox/C catalysts were experimentally compared in terms of catalytic performance and structural characteristics to determine that higher iron content and a less crystalline structure in the secondary matrix decrease the charge transfer resistance and thus is beneficial for electrocatalytic activity. This conclusion is supported by density-functional theory calculations showing distorted active sites in ternary metal catalysts are key for lowering overpotentials for the oxygen evolution reaction.

Environmental impact of amino acids on the release of selenate immobilized in hydrotalcite: Integrated interpretation of experimental and density-functional theory study.

The environmental impact of amino acids on the release of SeO42- immobilized into hydrotalcite (Mg2Al-LDH) which belongs to the layered double hydroxides (LDHs) family was investigated by experimental study and the observed layer structure of hydrotalcite was verified through density-functional theory (DFT) calculations. Glycine, l-cysteine, and l-aspartic acid, which have smaller molecular sizes, can release SeO42- largely due to intercalation, unstabilization of Mg2Al-LDH and simple dissolution, while l-tryptophan and l-phenylalanine caused limited SeO42- release due to their larger sizes and aromaticity. XRD patterns for the solid residues after intercalation of amino acids revealed that the layer distance of Mg2Al-LDH was partially expanded. The main peaks and shoulder features corresponding to d003 diffraction were well explained by DFT simulations using glycine as a model: the layer spacing of the main peak is responsible for the remaining SeO42- and singly stacked glycine molecule and the layer spacing of the shoulder peak was well explained by doubly stacked glycine molecules. Hydrogen bonds between amino acids and hydroxyl ions in the metallic layers of Mg2Al-LDH were responsible for the stable configuration of the intercalated Mg2Al-LDH. This study indicates potential limitations to the stability of low-level radioactive wastes of 79Se in repositories which are affected by smaller molecules of amino acids released through degradation of organic matters in the pedosphere.

Data-driven analysis of the electronic-structure factors controlling the work functions of perovskite oxides.

Tuning the work functions of materials is of practical interest for maximizing the performance of microelectronic and (photo)electrochemical devices, as the efficiency of these systems depends on the ability to control electronic levels at surfaces and across interfaces. Perovskites are promising compounds to achieve such control. In this work, we examine the work functions of more than 1000 perovskite oxide surfaces (ABO3) using data-driven (machine-learning) analysis and identify the factors that determine their magnitude. While the work functions of the BO2-terminated surfaces are sensitive to the energy of the hybridized oxygen p bands, the work functions of the AO-terminated surfaces exhibit a much less trivial dependence with respect to the filling of the d bands of the B-site atom and of its electronic affinity. This study shows the utility of interpretable data-driven models in analyzing the work functions of cubic perovskites from a limited number of electronic-structure descriptors.

Photophysics and Electronic Structure of Lateral Graphene/MoS2 and Metal/MoS2 Junctions.

Integration of semiconducting transition metal dichalcogenides (TMDs) into functional optoelectronic circuitries requires an understanding of the charge transfer across the interface between the TMD and the contacting material. Here, we use spatially resolved photocurrent microscopy to demonstrate electronic uniformity at the epitaxial graphene/molybdenum disulfide (EG/MoS2) interface. A 10× larger photocurrent is extracted at the EG/MoS2 interface when compared to the metal (Ti/Au)/MoS2 interface. This is supported by semi-local density functional theory (DFT), which predicts the Schottky barrier at the EG/MoS2 interface to be ∼2× lower than that at Ti/MoS2. We provide a direct visualization of a 2D material Schottky barrier through combination of angle-resolved photoemission spectroscopy with spatial resolution selected to be ∼300 nm (nano-ARPES) and DFT calculations. A bending of ∼500 meV over a length scale of ∼2-3 µm in the valence band maximum of MoS2 is observed via nano-ARPES. We explicate a correlation between experimental demonstration and theoretical predictions of barriers at graphene/TMD interfaces. Spatially resolved photocurrent mapping allows for directly visualizing the uniformity of built-in electric fields at heterostructure interfaces, providing a guide for microscopic engineering of charge transport across heterointerfaces. This simple probe-based technique also speaks directly to the 2D synthesis community to elucidate electronic uniformity at domain boundaries alongside morphological uniformity over large areas.

Phase-Selective Solution Synthesis of Perovskite-Related Cesium Cadmium Chloride Nanoparticles.

All-inorganic metal halide perovskite-related phases are semiconducting materials that are of significant interest for a wide range of applications. Nanoparticles of these materials are particularly useful because they permit solution processing while offering unique and tunable properties. Of the many metal halide systems that have been studied extensively, cesium cadmium chlorides remain underexplored, and synthetic routes to access them as nanoscale materials have not been established. Here we demonstrate that a simple solution-phase reaction involving the injection of a cesium oleate solution into a cadmium chloride solution produces three distinct cesium cadmium chlorides: hexagonal CsCdCl3 and the Ruddlesden-Popper layered perovskites Cs2CdCl4 and Cs3Cd2Cl7. The phase-selective synthesis emerges from differences in reagent concentrations, temperature, and injection rates. A key variable is the rate at which the cesium oleate solution is injected into the cadmium chloride solution, which is believed to influence the local Cs:Cd concentration during precipitation, leading to control over the phase that forms. Band structure calculations indicate that hexagonal CsCdCl3 is a direct band gap semiconductor while Cs2CdCl4 and Cs3Cd2Cl7 have indirect band gaps. The experimentally determined band gap values for CsCdCl3, Cs2CdCl4, and Cs3Cd2Cl7 are 5.13, 4.91, and 4.70 eV, respectively, which places them in a rare category of ultrawide-band-gap semiconductors.

Optimized utilization of COMB3 reactive potentials in LAMMPS.

An investigation to optimize the application of the third-generation charge optimized many-body (COMB3) interatomic potential and associated input parameters was carried out through the study of solid-liquid interactions in classical molecular dynamics simulations. The rates of these molecular interactions are understood through the wetting rates of water nano-droplets on a bare copper (111) surface. Implementing the Langevin thermostat, the influence of simulation time step, the number of atoms in the system, the frequency at which charge equilibration is performed, and the temperature relaxation rate are all examined. The results indicate that time steps of 0.4 fs are possible when using longer relaxation times for the system temperature, which is almost double the typical time step used for reactive potentials. The use of the charge equilibration allows for a fewer atomic layers to be used in the Cu slab. In addition, charge equilibrium schemes do not need to be performed every time step to ensure accurate charge transfer. Interestingly, the rate of wetting for the nanodroplets is dominantly dependent on the temperature relaxation time, which is predicted to significantly change the viscosity of the water droplets. This work provides a pathway for optimizing simulations using the COMB3 reactive interatomic potential.

Achieving Minimal Heat Conductivity by Ballistic Confinement in Phononic Metalattices.

Controlling the thermal conductivity of semiconductors is of practical interest in optimizing the performance of thermoelectric and phononic devices. The insertion of inclusions of nanometer size in a semiconductor is an effective means of achieving such control; it has been proposed that the thermal conductivity of silicon could be reduced to 1 W/m/K using this approach and that a minimum in the heat conductivity would be reached for some optimal size of the inclusions. Yet the experimental verification of this design rule has been limited. In this work, we address this question by studying the thermal properties of silicon metalattices that consist of a periodic distribution of spherical inclusions with radii from 7 to 30 nm, embedded into silicon. Experimental measurements confirm that the thermal conductivity of silicon metalattices is as low as 1 W/m/K for silica inclusions and that this value can be further reduced to 0.16 W/m/K for silicon metalattices with empty pores. A detailed model of ballistic phonon transport suggests that this thermal conductivity is close to the lowest achievable by tuning the radius and spacing of the periodic inhomogeneities. This study is a significant step in elucidating the scaling laws that dictate ballistic heat transport at the nanoscale in silicon and other semiconductors.

Energy-resolved distribution of electron traps for O/S-doped carbon nitrides by reversed double-beam photoacoustic spectroscopy and the photocatalytic reduction of Cr(vi).

We report for the first time to our knowledge the identification of heteroatom-doped and undoped C3N4 with the energy-resolved distribution of electron traps (ERDT) near the conduction band bottom position (CBB) using reversed double-beam photoacoustic spectroscopy. The ERDT/CBB pattern is used to classify the type of elemental doping in C3N4, related to photocatalytic efficiency.

Effects of surface charge and cluster size on the electrochemical dissolution of platinum nanoparticles using COMB3 and continuum electrolyte models.

We study the site-dependent dissolution of platinum nanoparticles under electrochemical conditions to assess their thermodynamic stability as a function of shape and size using empirical molecular dynamics and electronic-structure models. The third-generation charge optimized many-body potential is employed to determine the validity of uniform spherical representations of the nanoparticles in predicting dissolution potentials (the Kelvin model). To understand the early stages of catalyst dissolution, implicit solvation techniques based on the self-consistent continuum solvation method are applied. It is demonstrated that interfacial charge and polarization can shift the dissolution energies by amounts on the order of 0.74 eV depending on the surface site and nanoparticle shape, leading to the unexpected preferential removal of platinum cations from highly coordinated sites in some cases.

Spectroscopic and first-principles investigations of iodine species incorporation into ettringite: Implications for iodine migration in cement waste forms.

Low-level radioactive wastes are commonly immobilized in cementitious materials, where cement-based material can incorporate radionuclides into their crystal structure. Specifically, ettringite (Ca6Al2(OH)12(SO4)3∙26H2O) is known to stabilize anionic species, which is appealing for waste streams with radioactive iodine (129I) that persists as iodide (I-) and iodate (IO3-) in the cementitious nuclear waste repository. However, the structural information and immobilization mechanisms of iodine species in ettringite remain unclear. The present results suggested minimal I- incorporation into ettringite (0.05 %), whereas IO3- exhibited a high affinity for ettringite via anion substitution for SO42- (96 %). The combined iodine K-edge extended X-ray absorption fine structure (EXAFS) spectra and first-principles calculations using density functional theory (DFT) suggested that IO3- was stabilized in ettringite by hydrogen bonding and electrostatic forces. Substituting IO3- for SO42- was energetically favorable by -0.41 eV, whereas unfavorable substitution energy of 4.21 eV was observed for I- substitution. Moreover, the bonding charge density analysis of the substituted IO3- and I- anions into the ettringite structure revealed the interaction between intercalated ions with the structural water molecules. These results provided valuable insight into the long-term stabilization of anionic iodine species and their migration in cementitious nuclear waste repository or alkaline environments.

Vibrational probe of the origin of singlet exciton fission in TIPS-pentacene solutions.

We use native vibrational modes of the model singlet fission chromophore 6,13-bis(triisopropylsilylethynyl)pentacene (TIPS-Pn) to examine the origins of singlet fission in solution between molecules that are not tethered by a covalent linkage. We use the C-H stretch modes of TIPS side groups of TIPS-Pn to demonstrate that singlet fission does not occur by diffusive encounter of independent molecules in solution. Instead, TIPS-Pn molecules aggregate in solution through their TIPS side groups. This aggregation breaks the symmetry of the TIPS-Pn molecules and enables the formation of triplets to be probed through the formally symmetry forbidden symmetric alkyne stretch mode of the TIPS side groups. The alkyne stretch modes of TIPS-Pn are sensitive to the electronic excited states present during the singlet fission reaction and provide unique signatures of the formation of triplets following the initial separation of triplet pair intermediates. These findings highlight the opportunity to leverage structural information from vibrational modes to better understand intermolecular interactions that lead to singlet fission.

Understanding the influence of defects and surface chemistry on ferroelectric switching: a ReaxFF investigation of BaTiO3.

Ferroelectric materials such as barium titanate (BaTiO3) have a wide range of applications in nano scale electronic devices due to their outstanding properties. In this study, we developed an easily extendable atomistic ReaxFF reactive force field for BaTiO3 that can capture both its field- and temperature-induced ferroelectric hysteresis and corresponding changes due to surface chemistry and bulk defects. Using our force field, we were able to reproduce and explain a number of experimental observations: (1) the existence of a critical thickness of 4.8 nm below which ferroelectricity vanishes in BaTiO3; (2) migration and clustering of oxygen vacancies (OVs) in BaTiO3 and a reduction in the polarization and the Curie temperature due to the OVs; (3) domain wall interaction with the surface chemistry to influence the ferroelectric switching and polarization magnitude. This new computational tool opens up a wide range of possibilities for making predictions for realistic ferroelectric interfaces in energy-conversion, electronic and neuromorphic systems.