A bridge between trust and control: computational workflows meet automated battery cycling

Compliance with good research data management practices means trust in the integrity of the data, and it is achievable by full control of the data gathering process. In this work, we demonstrate tooling which bridges these two aspects, and illustrate its use in a case study of automated battery cycling. We successfully interface off-the-shelf battery cycling hardware with the computational workflow management software AiiDA, allowing us to control experiments, while ensuring trust in the data by tracking its provenance. We design user interfaces compatible with this tooling, which span the inventory, experiment design, and result analysis stages. Other features, including monitoring of workflows and import of externally generated and legacy data are also implemented. Finally, the full software stack required for this work is made available in a set of open-source packages.

The rule of four: anomalous distributions in the stoichiometries of inorganic compounds

Why are materials with specific characteristics more abundant than others? This is a fundamental question in materials science and one that is traditionally difficult to tackle, given the vastness of compositional and configurational space. We highlight here the anomalous abundance of inorganic compounds whose primitive unit cell contains a number of atoms that is a multiple of four. This occurrence-named here the rule of four-has to our knowledge not previously been reported or studied. Here, we first highlight the rule's existence, especially notable when restricting oneself to experimentally known compounds, and explore its possible relationship with established descriptors of crystal structures, from symmetries to energies. We then investigate this relative abundance by looking at structural descriptors, both of global (packing configurations) and local (the smooth overlap of atomic positions) nature. Contrary to intuition, the overabundance does not correlate with low-energy or high-symmetry structures; in fact, structures which obey the rule of four are characterized by low symmetries and loosely packed arrangements maximizing the free volume. We are able to correlate this abundance with local structural symmetries, and visualize the results using a hybrid supervised-unsupervised machine learning method.

Magnetostriction-Driven Muon Localization in an Antiferromagnetic Oxide

Magnetostriction results from the coupling between magnetic and elastic degrees of freedom. Though it is associated with a relatively small energy, we show that it plays an important role in determining the site of an implanted muon, so that the energetically favorable site can switch on crossing a magnetic phase transition. This surprising effect is demonstrated in the cubic rocksalt antiferromagnet MnO which undergoes a magnetostriction-driven rhombohedral distortion at the Néel temperature T_{N}=118  K. Above T_{N}, the muon becomes delocalized around a network of equivalent sites, but below T_{N} the distortion lifts the degeneracy between these equivalent sites. Our first-principles simulations based on Hubbard-corrected density-functional theory and molecular dynamics are consistent with the experimental data and help to resolve a long-standing puzzle regarding muon data on MnO, as well as having wider applicability to other magnetic oxides.

Mechanistic Insights on Functionalization of Graphene with Ozone

The exposure of graphene to O3 results in functionalization of its lattice with epoxy, even at room temperature. This reaction is of fundamental interest for precise lattice patterning, however, is not well understood. Herein, using van der Waals density functional theory (vdW-DFT) incorporating spin-polarized calculations, we find that O3 strongly physisorbs on graphene with a binding energy of -0.46 eV. It configures in a tilted position with the two terminal O atoms centered above the neighboring graphene honeycombs. A dissociative chemisorption follows by surpassing an energy barrier of 0.75 eV and grafting an epoxy group on graphene reducing the energy of the system by 0.14 eV from the physisorbed state. Subsequent O3 chemisorption is preferred on the same honeycomb, yielding two epoxy groups separated by a single C-C bridge. We show that capturing the onset of spin in oxygen during chemisorption is crucial. We verify this finding with experiments where an exponential increase in the density of epoxy groups as a function of reaction temperature yields an energy barrier of 0.66 eV, in agreement with the DFT prediction. These insights will help efforts to obtain precise patterning of the graphene lattice.

Topological magnons driven by the Dzyaloshinskii-Moriya interaction in the centrosymmetric ferromagnet Mn5Ge3

The phase of the quantum-mechanical wave function can encode a topological structure with wide-ranging physical consequences, such as anomalous transport effects and the existence of edge states robust against perturbations. While this has been exhaustively demonstrated for electrons, properties associated with the elementary quasiparticles in magnetic materials are still underexplored. Here, we show theoretically and via inelastic neutron scattering experiments that the bulk ferromagnet Mn5Ge3 hosts gapped topological Dirac magnons. Although inversion symmetry prohibits a net Dzyaloshinskii-Moriya interaction in the unit cell, it is locally allowed and is responsible for the gap opening in the magnon spectrum. This gap is predicted and experimentally verified to close by rotating the magnetization away from the c-axis with an applied magnetic field. Hence, Mn5Ge3 realizes a gapped Dirac magnon material in three dimensions. Its tunability by chemical doping or by thin film nanostructuring defines an exciting new platform to explore and design topological magnons. More generally, our experimental route to verify and control the topological character of the magnons is applicable to bulk centrosymmetric hexagonal materials, which calls for systematic investigation.

Projectability disentanglement for accurate and automated electronic-structure Hamiltonians

Maximally-localized Wannier functions (MLWFs) are broadly used to characterize the electronic structure of materials. Generally, one can construct MLWFs describing isolated bands (e.g. valence bands of insulators) or entangled bands (e.g. valence and conduction bands of insulators, or metals). Obtaining accurate and compact MLWFs often requires chemical intuition and trial and error, a challenging step even for experienced researchers and a roadblock for high-throughput calculations. Here, we present an automated approach, projectability-disentangled Wannier functions (PDWFs), that constructs MLWFs spanning the occupied bands and their complement for the empty states, providing a tight-binding picture of optimized atomic orbitals in crystals. Key to the algorithm is a projectability measure for each Bloch state onto atomic orbitals, determining if that state should be kept identically, discarded, or mixed into the disentanglement. We showcase the accuracy on a test set of 200 materials, and the reliability by constructing 21,737 Wannier Hamiltonians.

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.

Unraveling the Oxidation of a Graphitic Lattice: Structure Determination of Oxygen Clusters

Unraveling the oxidation of graphitic lattice is of great interest for atomic-scale lattice manipulation. Herein, we build epoxy cluster, atom by atom, using Van der Waals' density-functional theory aided by Clar's aromatic π-sextet rule. We predict the formation of cyclic epoxy trimers and its linear chains propagating along the armchair direction of the lattice to minimize the system's energy. Using low-temperature scanning tunneling microscopy on oxidized graphitic lattice, we identify linear chains as bright features that have a threefold symmetry, and which exclusively run along the armchair direction of the lattice confirming the theoretical predictions.

Infrared-active phonons in one-dimensional materials and their spectroscopic signatures

Dimensionality provides a clear fingerprint on the dispersion of infrared-active, polar-optical phonons. For these phonons, the local dipoles parametrized by the Born effective charges drive the LO-TO splitting of bulk materials; this splitting actually breaks down in two-dimensional materials. Here, we develop the theory for one-dimensional (1D) systems-nanowires, nanotubes, and atomic and polymeric chains. Combining an analytical model with the implementation of density-functional perturbation theory in 1D boundary conditions, we show that the dielectric splitting in the dispersion relations collapses as x 2 log ( x ) at the zone center. The dielectric properties and the radius of the 1D materials are linked by the present work to these red shifts, opening infrared and Raman characterization avenues.

Electronic Structure of Few-Layer Black Phosphorus from μ-ARPES

Black phosphorus (BP) stands out among two-dimensional (2D) semiconductors because of its high mobility and thickness dependent direct band gap. However, the quasiparticle band structure of ultrathin BP has remained inaccessible to experiment thus far. Here we use a recently developed laser-based microfocus angle resolved photoemission (μ-ARPES) system to establish the electronic structure of 2-9 layer BP from experiment. Our measurements unveil ladders of anisotropic, quantized subbands at energies that deviate from the scaling observed in conventional semiconductor quantum wells. We quantify the anisotropy of the effective masses and determine universal tight-binding parameters, which provide an accurate description of the electronic structure for all thicknesses.

Thermal conductivity of glasses: first-principles theory and applications

Predicting the thermal conductivity of glasses from first principles has hitherto been a very complex problem. The established Allen-Feldman and Green-Kubo approaches employ approximations with limited validity-the former neglects anharmonicity, the latter misses the quantum Bose-Einstein statistics of vibrations-and require atomistic models that are very challenging for first-principles methods. Here, we present a protocol to determine from first principles the thermal conductivity κ(T) of glasses above the plateau (i.e., above the temperature-independent region appearing almost without exceptions in the κ(T) of all glasses at cryogenic temperatures). The protocol combines the Wigner formulation of thermal transport with convergence-acceleration techniques, and accounts comprehensively for the effects of structural disorder, anharmonicity, and Bose-Einstein statistics. We validate this approach in vitreous silica, showing that models containing less than 200 atoms can already reproduce κ(T) in the macroscopic limit. We discuss the effects of anharmonicity and the mechanisms determining the trend of κ(T) at high temperature, reproducing experiments at temperatures where radiative effects remain negligible.

Expansion of the Materials Cloud 2D Database

Two-dimensional (2D) materials are among the most promising candidates for beyond-silicon electronic, optoelectronic, and quantum computing applications. Recently, their recognized importance sparked a push to discover and characterize novel 2D materials. Within a few years, the number of experimentally exfoliated or synthesized 2D materials went from a few to more than a hundred, with the number of theoretically predicted compounds reaching a few thousand. In 2018 we first contributed to this effort with the identification of 1825 compounds that are either easily (1036) or potentially (789) exfoliable from experimentally known 3D compounds. Here, we report on a major expansion of this 2D portfolio thanks to the extension of the screening protocol to an additional experimental database (MPDS) as well as the updated versions of the two databases (ICSD and COD) used in our previous work. This expansion leads to the discovery of an additional 1252 monolayers, bringing the total to 3077 compounds and, notably, almost doubling the number of easily exfoliable materials to 2004. We optimize the structural properties of all these monolayers and explore their electronic structure with a particular emphasis on those rare large-bandgap 2D materials that could be precious in isolating 2D field-effect-transistor channels. Finally, for each material containing up to 6 atoms per unit cell, we identify the best candidates to form commensurate heterostructures, balancing requirements on supercell size and minimal strain.

Accurate Prediction of Hall Mobilities in Two-Dimensional Materials through Gauge-Covariant Quadrupolar Contributions

Despite considerable efforts, accurate computations of electron-phonon and carrier transport properties of low-dimensional materials from first principles have remained elusive. By building on recent advances in the description of long-range electrostatics, we develop a general approach to the calculation of electron-phonon couplings in two-dimensional materials. We show that the nonanalytic behavior of the electron-phonon matrix elements depends on the Wannier gauge, but that a missing Berry connection restores invariance to quadrupolar order. We showcase these contributions in a MoS_{2} monolayer, calculating intrinsic drift and Hall mobilities with precise Wannier interpolations. We also find that the contributions of dynamical quadrupoles to the scattering potential are essential, and that their neglect leads to errors of 23% and 76% in the room-temperature electron and hole Hall mobilities, respectively.

Testing Koopmans spectral functionals on the analytically solvable Hooke's atom

Koopmans spectral functionals are a class of orbital-density-dependent functionals designed to accurately predict spectroscopic properties. They do so markedly better than their Kohn-Sham density-functional theory counterparts, as demonstrated in earlier works on benchmarks of molecules and bulk systems. This work is a complementary study where-instead of comparing against real, many-electron systems-we test Koopmans spectral functionals on Hooke's atom, a toy two-electron system that has analytical solutions for particular strengths of its harmonic confining potential. As these calculations clearly illustrate, Koopmans spectral functionals do an excellent job of describing Hooke's atom across a range of confining potential strengths. This work also provides broader insights into the features and capabilities of Koopmans spectral functionals more generally.

Unraveling the effects of inter-site Hubbard interactions in spinel Li-ion cathode materials

Accurate first-principles predictions of the structural, electronic, magnetic, and electrochemical properties of cathode materials can be key in the design of novel efficient Li-ion batteries. Spinel-type cathode materials Li_xMn₂O₄ and Li_xMn_1.5Ni_0.5O₄ are promising candidates for Li-ion battery technologies, but they present serious challenges when it comes to their first-principles modeling. Here, we use density-functional theory with extended Hubbard functionals-DFT+U+V with on-site U and inter-site V Hubbard interactions-to study the properties of these transition-metal oxides. The Hubbard parameters are computed from first-principles using density-functional perturbation theory. We show that while U is crucial to obtain the right trends in properties of these materials, V is essential for a quantitative description of the structural and electronic properties, as well as the Li-intercalation voltages. This work paves the way for reliable first-principles studies of other families of cathode materials without relying on empirical fitting or calibration procedures.

First-Principles Thermodynamics of CsSnI3

CsSnI₃ is a promising ecofriendly solution for energy harvesting technologies. It exists at room temperature in either a black perovskite polymorph or a yellow 1D double-chain, which irreversibly deteriorates in the air. In this work, we unveil the relative thermodynamic stability between the two structures with a first-principles sampling of the CsSnI₃ finite-temperature phase diagram, discovering how it is driven by anomalously large quantum and anharmonic ionic fluctuations. Thanks to a comprehensive treatment of anharmonicity, the simulations deliver a remarkable agreement with known experimental data for the transition temperatures of the orthorhombic, rhombohedral, and cubic perovskite structures and the thermal expansion coefficient. We disclose how the perovskite polymorphs are the ground state above 270 K and discover an abnormal decrease in heat capacity upon heating in the cubic black perovskite. Our results also significantly downplay the Cs⁺ rattling modes' contribution to mechanical instability. The remarkable agreement with experiments validates our methodology, which can be systematically applied to all metal halides.

Solids that are also liquids: elastic tensors of superionic materials

Superionics are fascinating materials displaying both solid- and liquid-like characteristics: as solids, they respond elastically to shear stress; as liquids, they display fast-ion diffusion at normal conditions. In addition to such scientific interest, superionics are technologically relevant for energy, electronics, and sensing applications. Characterizing and understanding their elastic properties is, e.g., urgently needed to address their feasibility as solid-state electrolytes in all-solid-state batteries. However, static approaches to elasticity assume well-defined reference positions around which atoms vibrate, in contrast with the quasi-liquid motion of the mobile ions in fast ionic conductors. Here, we derive the elastic tensors of superionics from ensemble fluctuations in the isobaric-isothermal ensemble, exploiting extensive Car-Parrinello simulations. We apply this approach to paradigmatic Li-ion conductors, and complement with a block analysis to compute statistical errors. Static approaches sampled over the trajectories often overestimate the response, highlighting the importance of a dynamical treatment in determining elastic tensors in superionics.

In Situ Nucleation-Decoupled and Site-Specific Incorporation of Å-Scale Pores in Graphene Via Epoxidation

Generating pores in graphene by decoupled nucleation and expansion is desired to achieve a fine control over the porosity, and is desired to advance several applications. Herein, epoxidation is introduced, which is the formation of nanosized epoxy clusters on the graphitic lattice as nucleation sites without forming pores. In situ gasification of clusters inside a transmission electron microscope shows that pores are generated precisely at the site of the clusters by surpassing an energy barrier of 1.3 eV. Binding energy predictions using ab initio calculations combined with the cluster nucleation theory reveal the structure of the epoxy clusters and indicate that the critical cluster is an epoxy dimer. Finally, it is shown that the cluster gasification can be manipulated to form Å-scale pores which then effectively sieve gas molecules based on their size. This decoupled cluster nucleation and pore formation will likely pave the way for an independent control of pore size and density.

Koopmans Spectral Functionals in Periodic Boundary Conditions

Koopmans spectral functionals aim to describe simultaneously ground-state properties and charged excitations of atoms, molecules, nanostructures, and periodic crystals. This is achieved by augmenting standard density functionals with simple but physically motivated orbital-density-dependent corrections. These corrections act on a set of localized orbitals that, in periodic systems, resemble maximally localized Wannier functions. At variance with the original, direct supercell implementation (Phys. Rev. X 2018, 8, 021051), we discuss here (i) the complex but efficient formalism required for a periodic boundary code using explicit Brillouin zone sampling and (ii) the calculation of the screened Koopmans corrections with density functional perturbation theory. In addition to delivering improved scaling with system size, the present development makes the calculation of band structures with Koopmans functionals straightforward. The implementation in the open-source Quantum ESPRESSO distribution and the application to prototypical insulating and semiconducting systems are presented and discussed.

Phonon-Assisted Luminescence in Defect Centers from Many-Body Perturbation Theory

Phonon-assisted luminescence is a key property of defect centers in semiconductors, and can be measured to perform the readout of the information stored in a quantum bit, or to detect temperature variations. The investigation of phonon-assisted luminescence usually employs phenomenological models, such as that of Huang and Rhys, with restrictive assumptions that can fail to be predictive. In this work, we predict luminescence and study exciton-phonon couplings within a rigorous many-body perturbation theory framework, an analysis that has never been performed for defect centers. In particular, we study the optical emission of the negatively charged boron vacancy in 2D hexagonal boron nitride, which currently stands out among defect centers in 2D materials thanks to its promise for applications in quantum information and quantum sensing. We show that phonons are responsible for the observed luminescence, which otherwise would be dark due to symmetry. We also show that the symmetry breaking induced by the static Jahn-Teller effect is not able to describe the presence of the experimentally observed peak at 1.5 eV.

Structure Evolution of Graphitic Surface upon Oxidation: Insights by Scanning Tunneling Microscopy

Oxidation of graphitic materials has been studied for more than a century to synthesize materials such as graphene oxide, nanoporous graphene, and to cut or unzip carbon nanotubes. However, the understanding of the early stages of oxidation is limited to theoretical studies, and experimental validation has been elusive. This is due to (i) challenging sample preparation for characterization because of the presence of highly mobile and reactive epoxy groups formed during oxidation, and (ii) gasification of the functional groups during imaging with atomic resolution, e.g., by transmission electron microscopy. Herein, we utilize a low-temperature scanning tunneling microscope (LT-STM) operating at 4 K to solve the structure of epoxy clusters form upon oxidation. Three distinct nanostructures corresponding to three stages of evolution of vacancy defects are found by quantitatively verifying the experimental data by the van der Waals density functional theory. The smallest cluster is a cyclic epoxy trimer. Their observation validates the theoretical prediction that epoxy trimers minimize the energy in the cyclic structure. The trimers grow into honeycomb superstructures to form larger clusters (1-3 nm). Vacancy defects evolve only in the larger clusters (2-3 nm) in the middle of the cluster, highlighting the role of lattice strain in the generation of vacancies. Semiquinone groups are also present and are assigned at the carbon edge in the vacancy defects. Upon heating to 800 °C, we observe cluster-free vacancy defects resulting from the loss of the entire epoxy population, indicating a reversible functionalization of epoxy groups.

Nanoscale-Resolved Surface-to-Bulk Electron Transport in CsPbBr3 Perovskite

Describing the nanoscale charge carrier transport at surfaces and interfaces is fundamental for designing high-performance optoelectronic devices. To achieve this, we employ time- and angle-resolved photoelectron spectroscopy with ultraviolet pump and extreme ultraviolet probe pulses. The resulting high surface sensitivity reveals an ultrafast carrier population decay associated with surface-to-bulk transport, which was tracked with a sub-nanometer spatial resolution normal to the surface, and on a femtosecond time scale, in the case of the inorganic CsPbBr₃ lead halide perovskite. The decay time exhibits a pronounced carrier density dependence, which is attributed via modeling to enhanced diffusive transport and concurrent recombination. The transport is found to approach an ordinary diffusive regime, limited by electron-hole scattering, at the highest excitation fluences. This approach constitutes an important milestone in our capability to probe hot-carrier transport at solid interfaces with sub-nanometer resolution in a theoretically and experimentally challenging, yet technologically relevant, high-carrier-density regime.

Energy Level Alignment at the Cobalt Phosphate/Electrolyte Interface: Intrinsic Stability vs Interfacial Chemical Reactions in 5 V Lithium Ion Batteries

The intrinsic stability of the 5 V LiCoPO₄-LiCo₂P₃O₁₀ thin-film (carbon-free) cathode material coated with MoO₃ thin layer is studied using a comprehensive synchrotron electron spectroscopy in situ approach combined with first-principle calculations. The atomic-molecular level study demonstrates fully reversible electronic properties of the cathode after the first electrochemical cycle. The polyanionic oxide is not involved in chemical reactions with the fluoroethylene-containing liquid electrolyte even when charged to 5.1 V vs Li⁺/Li. The high stability of the cathode is explained on the basis of the developed energy level model. In contrast, the chemical composition of the cathode-electrolyte interface evolves continuously by involving MoO₃ in the decomposition reaction with consequent leaching of oxide from the surface. The proposed mechanisms of chemical reactions are attributed to external electrolyte oxidation via charge transfer from the relevant electron level to the MoO₃ valence band state and internal electrolyte oxidation via proton transfer to the solvents. This study provides a deeper insight into the development of both a doping strategy to enhance the electronic conductivity of high-voltage cathode materials and an efficient surface coating against unfavorable interfacial chemical reactions.

Gate Control of Spin-Layer-Locking FETs and Application to Monolayer LuIO

A recent 2D spinFET concept proposes to switch electrostatically between two separate sublayers with strong and opposite intrinsic Rashba effects, exploiting the spin-layer-locking mechanism in centrosymmetric materials with local dipole fields. Here, we propose a novel monolayer material within this family, lutetium oxide iodide (LuIO). It displays one of the largest Rashba effects among 2D materials (up to k_R = 0.08 Å^-1), leading to a π/2 rotation of the spins over just 1 nm. The monolayer was predicted to be exfoliable from its experimentally known 3D bulk counterpart, with a binding energy lower than graphene. We characterize and simulate the interplay of the two gate-controlled parameters for such devices: doping and spin channel selection. We show that the ability to split the spin channels in energy diminishes with doping, leading to specific gate-operation guidelines that can apply to all devices based on spin-layer locking.

Ab Initio Electron-Phonon Interactions in Correlated Electron Systems

Electron-phonon (e-ph) interactions are pervasive in condensed matter, governing phenomena such as transport, superconductivity, charge-density waves, polarons, and metal-insulator transitions. First-principles approaches enable accurate calculations of e-ph interactions in a wide range of solids. However, they remain an open challenge in correlated electron systems (CES), where density functional theory often fails to describe the ground state. Therefore reliable e-ph calculations remain out of reach for many transition metal oxides, high-temperature superconductors, Mott insulators, planetary materials, and multiferroics. Here we show first-principles calculations of e-ph interactions in CES, using the framework of Hubbard-corrected density functional theory (DFT+U) and its linear response extension (DFPT+U), which can describe the electronic structure and lattice dynamics of many CES. We showcase the accuracy of this approach for a prototypical Mott system, CoO, carrying out a detailed investigation of its e-ph interactions and electron spectral functions. While standard DFPT gives unphysically divergent and short-ranged e-ph interactions, DFPT+U is shown to remove the divergences and properly account for the long-range Fröhlich interaction, allowing us to model polaron effects in a Mott insulator. Our work establishes a broadly applicable and affordable approach for quantitative studies of e-ph interactions in CES, a novel theoretical tool to interpret experiments in this broad class of materials.

Shear and Breathing Modes of Layered Materials

Layered materials (LMs), such as graphite, hexagonal boron nitride, and transition-metal dichalcogenides, are at the center of an ever-increasing research effort, due to their scientific and technological relevance. Raman and infrared spectroscopies are accurate, non-destructive approaches to determine a wide range of properties, including the number of layers, N, and the strength of the interlayer interactions. We present a general approach to predict the complete spectroscopic fan diagrams, i.e., the relations between frequencies and N for the optically active shear and layer-breathing modes of any multilayer comprising N ≥ 2 identical layers. In order to achieve this, we combine a description of the normal modes in terms of a one-dimensional mechanical model, with symmetry arguments that describe the evolution of the point group as a function of N. Group theory is then used to identify which modes are Raman- and/or infrared-active, and to provide diagrams of the optically active modes for any stack composed of identical layers. We implement the method and algorithms in an open-source tool to assist researchers in the prediction and interpretation of such diagrams. Our work will underpin future efforts on Raman and infrared characterization of known, and yet not investigated, LMs.

OPTIMADE, an API for exchanging materials data

The Open Databases Integration for Materials Design (OPTIMADE) consortium has designed a universal application programming interface (API) to make materials databases accessible and interoperable. We outline the first stable release of the specification, v1.0, which is already supported by many leading databases and several software packages. We illustrate the advantages of the OPTIMADE API through worked examples on each of the public materials databases that support the full API specification.

Electronic Structure of Water from Koopmans-Compliant Functionals

Obtaining a precise theoretical description of the spectral properties of liquid water poses challenges for both molecular dynamics (MD) and electronic structure methods. The lower computational cost of the Koopmans-compliant functionals with respect to Green's function methods allows the simulations of many MD trajectories, with a description close to the state-of-art quasi-particle self-consistent GW plus vertex corrections method (QSGW + f_xc). Thus, we explore water spectral properties when different MD approaches are used, ranging from classical MD to first-principles MD, and including nuclear quantum effects. We have observed that different MD approaches lead to up to 1 eV change in the average band gap; thus, we focused on the band gap dependence with the geometrical properties of a system to explain such spread. We have evaluated the changes in the band gap due to variations in the intramolecular O-H bond distance and HOH angle, as well as the intermolecular hydrogen bond O···O distance and the OHO angles. We have observed that the dominant contribution comes from the O-H bond length; the O···O distance plays a secondary role, and the other geometrical properties do not significantly influence the gap. Furthermore, we analyze the electronic density of states (DOS), where the KIPZ functional shows good agreement with the DOS obtained with state-of-art approaches employing quasi-particle self-consistent GW plus vertex corrections. The O-H bond length also significantly influences the DOS. When nuclear quantum effects are considered, broadening of the peaks driven by the broader distribution of the O-H bond lengths is observed, leading to a closer agreement with the experimental photoemission spectra.

Electronic-structure methods for materials design

The accuracy and efficiency of electronic-structure methods to understand, predict and design the properties of materials has driven a new paradigm in research. Simulations can greatly accelerate the identification, characterization and optimization of materials, with this acceleration driven by continuous progress in theory, algorithms and hardware, and by adaptation of concepts and tools from computer science. Nevertheless, the capability to identify and characterize materials relies on the predictive accuracy of the underlying physical descriptions, and on the ability to capture the complexity of realistic systems. We provide here an overview of electronic-structure methods, of their application to the prediction of materials properties, and of the different strategies employed towards the broader goals of materials design and discovery.

Prediction of Phonon-Mediated Superconductivity with High Critical Temperature in the Two-Dimensional Topological Semimetal W2N3

Two-dimensional superconductors attract great interest both for their fundamental physics and for their potential applications, especially in the rapidly growing field of quantum computing. Despite intense theoretical and experimental efforts, materials with a reasonably high transition temperature are still rare. Even more rare are those that combine superconductivity with a nontrivial band topology that could potentially give rise to exotic states of matter. Here, we predict a remarkably high superconducting critical temperature of 21 K in the easily exfoliable, topologically nontrivial 2D semimetal W₂N₃. By studying its electronic and superconducting properties as a function of doping and strain, we also find large changes in the electron-phonon interactions that make this material a unique platform to study different coupling regimes and test the limits of current theories of superconductivity. Last, we discuss the possibility of tuning the material to achieve coexistence of superconductivity and topologically nontrivial edge states.

Theoretical Prediction of Two-Dimensional Materials, Behavior, and Properties

Predictive modeling of two-dimensional (2D) materials is at the crossroad of two current rapidly growing interests: 2D materials per se, massively sought after and explored in experimental laboratories, and materials theoretical-computational models in general, flourishing on a fertile mix of condensed-matter physics and chemistry with advancing computational technology. Here the general methods and specific techniques of modeling are briefly overviewed, along with a somewhat philosophical assessment of what "prediction" is, followed by selected practical examples for 2D materials, from structures and properties, to device functionalities and synthetic routes for their making. We conclude with a brief sketch-outlook of future developments.

Role of OH Intermediates during the Au Oxide Electro-Reduction at Low pH Elucidated by Electrochemical Surface-Enhanced Raman Spectroscopy and Implicit Solvent Density Functional Theory

Molecular understanding of the electrochemical oxidation of metals and the electro-reduction of metal oxides is of pivotal importance for the rational design of catalyst-based devices where metal(oxide) electrodes play a crucial role. Operando monitoring and reliable identification of reacting species, however, are challenging tasks because they require surface-molecular sensitive and specific experiments under reaction conditions and sophisticated theoretical calculations. The lack of molecular insight under operating conditions is largely due to the limited availability of operando tools and to date still hinders a quick technological advancement of electrocatalytic devices. Here, we present a combination of advanced density functional theory (DFT) calculations considering implicit solvent contributions and time-resolved electrochemical surface-enhanced Raman spectroscopy (EC-SERS) to identify short-lived reaction intermediates during the showcase electro-reduction of Au oxide (AuOx) in sulfuric acid over several tens of seconds. The EC-SER spectra provide evidence for temporary Au-OH formation and for the asynchronous adsorption of (bi)sulfate ions at the surface during the reduction process. Spectral intensity fluctuations indicate an OH/(bi)sulfate turnover period of 4 s. As such, the presented EC-SERS potential jump approach combined with implicit solvent DFT simulations allows us to propose a reaction mechanism and prove that short-lived Au-OH intermediates also play an active role during the AuOx electro-reduction in acidic media, implying their potential relevance also for other electrocatalytic systems operating at low pH, like metal corrosion, the oxidation of CO, HCOOH, and other small organic molecules, and the oxygen evolution reaction.

Materials Cloud, a platform for open computational science

Materials Cloud is a platform designed to enable open and seamless sharing of resources for computational science, driven by applications in materials modelling. It hosts (1) archival and dissemination services for raw and curated data, together with their provenance graph, (2) modelling services and virtual machines, (3) tools for data analytics, and pre-/post-processing, and (4) educational materials. Data is citable and archived persistently, providing a comprehensive embodiment of entire simulation pipelines (calculations performed, codes used, data generated) in the form of graphs that allow retracing and reproducing any computed result. When an AiiDA database is shared on Materials Cloud, peers can browse the interconnected record of simulations, download individual files or the full database, and start their research from the results of the original authors. The infrastructure is agnostic to the specific simulation codes used and can support diverse applications in computational science that transcend its initial materials domain.

AiiDA 1.0, a scalable computational infrastructure for automated reproducible workflows and data provenance

The ever-growing availability of computing power and the sustained development of advanced computational methods have contributed much to recent scientific progress. These developments present new challenges driven by the sheer amount of calculations and data to manage. Next-generation exascale supercomputers will harden these challenges, such that automated and scalable solutions become crucial. In recent years, we have been developing AiiDA (aiida.net), a robust open-source high-throughput infrastructure addressing the challenges arising from the needs of automated workflow management and data provenance recording. Here, we introduce developments and capabilities required to reach sustained performance, with AiiDA supporting throughputs of tens of thousands processes/hour, while automatically preserving and storing the full data provenance in a relational database making it queryable and traversable, thus enabling high-performance data analytics. AiiDA's workflow language provides advanced automation, error handling features and a flexible plugin model to allow interfacing with external simulation software. The associated plugin registry enables seamless sharing of extensions, empowering a vibrant user community dedicated to making simulations more robust, user-friendly and reproducible.

2-D Materials for Ultrascaled Field-Effect Transistors: One Hundred Candidates under the Ab Initio Microscope

Due to their remarkable properties, single-layer 2-D materials appear as excellent candidates to extend Moore's scaling law beyond the currently manufactured silicon FinFETs. However, the known 2-D semiconducting components, essentially transition metal dichalcogenides, are still far from delivering the expected performance. Based on a recent theoretical study that predicts the existence of more than 1800 exfoliable 2-D materials, we investigate here the 100 most promising contenders for logic applications. Their current versus voltage characteristics are simulated from first-principles, combining density functional theory and advanced quantum transport calculations. Both n- and p-type configurations are considered, with gate lengths ranging from 15 down to 5 nm. From this large collection of electronic materials, we identify 13 compounds with electron and hole currents potentially much higher than those in future Si FinFETs. The resulting database widely expands the design space of 2-D transistors and provides original guidelines to the materials and device engineering community.

Evidence of Large Polarons in Photoemission Band Mapping of the Perovskite Semiconductor CsPbBr_{3}

Lead-halide perovskite (LHP) semiconductors are emergent optoelectronic materials with outstanding transport properties which are not yet fully understood. We find signatures of large polaron formation in the electronic structure of the inorganic LHP CsPbBr_{3} by means of angle-resolved photoelectron spectroscopy. The experimental valence band dispersion shows a hole effective mass of 0.26±0.02 m_{e}, 50% heavier than the bare mass m_{0}=0.17 m_{e} predicted by density functional theory. Calculations of the electron-phonon coupling indicate that phonon dressing of the carriers mainly occurs via distortions of the Pb-Br bond with a Fröhlich coupling parameter α=1.81. A good agreement with our experimental data is obtained within the Feynman polaron model, validating a viable theoretical method to predict the carrier effective mass of LHPs ab initio.

Quantum ESPRESSO toward the exascale

Quantum ESPRESSO is an open-source distribution of computer codes for quantum-mechanical materials modeling, based on density-functional theory, pseudopotentials, and plane waves, and renowned for its performance on a wide range of hardware architectures, from laptops to massively parallel computers, as well as for the breadth of its applications. In this paper, we present a motivation and brief review of the ongoing effort to port Quantum ESPRESSO onto heterogeneous architectures based on hardware accelerators, which will overcome the energy constraints that are currently hindering the way toward exascale computing.

Wannier90 as a community code: new features and applications

Wannier90 is an open-source computer program for calculating maximally-localised Wannier functions (MLWFs) from a set of Bloch states. It is interfaced to many widely used electronic-structure codes thanks to its independence from the basis sets representing these Bloch states. In the past few years the development of Wannier90 has transitioned to a community-driven model; this has resulted in a number of new developments that have been recently released in Wannier90 v3.0. In this article we describe these new functionalities, that include the implementation of new features for wannierisation and disentanglement (symmetry-adapted Wannier functions, selectively-localised Wannier functions, selected columns of the density matrix) and the ability to calculate new properties (shift currents and Berry-curvature dipole, and a new interface to many-body perturbation theory); performance improvements, including parallelisation of the core code; enhancements in functionality (support for spinor-valued Wannier functions, more accurate methods to interpolate quantities in the Brillouin zone); improved usability (improved plotting routines, integration with high-throughput automation frameworks), as well as the implementation of modern software engineering practices (unit testing, continuous integration, and automatic source-code documentation). These new features, capabilities, and code development model aim to further sustain and expand the community uptake and range of applicability, that nowadays spans complex and accurate dielectric, electronic, magnetic, optical, topological and transport properties of materials.

Bulk and Surface Electronic Structure of the Dual-Topology Semimetal Pt_{2}HgSe_{3}

We report high-resolution angle-resolved photoemission measurements on single crystals of Pt_{2}HgSe_{3} grown by high-pressure synthesis. Our data reveal a gapped Dirac nodal line whose (001) projection separates the surface Brillouin zone in topological and trivial areas. In the nontrivial k-space range, we find surface states with multiple saddle points in the dispersion, resulting in two van Hove singularities in the surface density of states. Based on density-functional theory calculations, we identify these surface states as signatures of a topological crystalline state, which coexists with a weak topological phase.

Large-scale synthesis of crystalline g-C3N4 nanosheets and high-temperature H2 sieving from assembled films

Poly(triazine imide) (PTI), a crystalline g-C₃N₄, hosting two-dimensional nanoporous structure with an electron density gap of 0.34 nm, is highly promising for high-temperature hydrogen sieving because of its high chemical and thermal robustness. Currently, layered PTI is synthesized in potentially unsafe vacuum ampules in milligram quantities. Here, we demonstrate a scalable and safe ambient pressure synthesis route leading to several grams of layered PTI platelets in a single batch with 70% yield with respect to the precursor. Solvent exfoliation under anhydrous conditions led to single-layer PTI nanosheets evidenced by the observation of triangular g-C₃N₄ nanopores. Gas permeation studies confirm that PTI nanopores can sieve He and H₂ from larger molecules. Last, high-temperature H₂ sieving from PTI nanosheet-based membranes, prepared by the scalable filter coating technique, is demonstrated with H₂ permeance reaching 1500 gas permeation units, with H₂/CO₂, H₂/N₂, and H₂/CH₄ selectivities reaching 10, 50, and 60, respectively, at 250°C.

Operando XANES from first-principles and its application to iridium oxide

Density-functional theory calculations augmented with a continuum description of the electrochemical environment are implemented to simulated X-ray absorption spectra as a function of the applied potential.

Relative Abundance of [Formula: see text] Topological Order in Exfoliable Two-Dimensional Insulators

Quantum spin Hall insulators make up a class of two-dimensional materials with a finite electronic band gap in the bulk and gapless helical edge states. In the presence of time-reversal symmetry, [Formula: see text] topological order distinguishes the topological phase from the ordinary insulating one. Some of the phenomena that can be hosted in these materials, from one-dimensional low-dissipation electronic transport to spin filtering, could be promising for many technological applications in the fields of electronics, spintronics, and topological quantum computing. Nevertheless, the rarity of two-dimensional materials that can exhibit nontrivial [Formula: see text] topological order at room temperature hinders development. Here, we screen a comprehensive database we recently identified of 1825 monolayers that can be exfoliated from experimentally known compounds to search for novel quantum spin Hall insulators. Using density-functional and many-body perturbation theory simulations, we identify 13 monolayers that are candidates for quantum spin Hall insulators including high-performing materials such as AsCuLi₂ and (platinum) jacutingaite (Pt₂HgSe₃). We also identify monolayer Pd₂HgSe₃ (palladium jacutingaite) as a novel Kane-Mele quantum spin Hall insulator and compare it with platinum jacutingaite. A handful of promising materials are mechanically stable and exhibit [Formula: see text] topological order, either unperturbed or driven by small amounts of strain. Such screening highlights a relative abundance of [Formula: see text] topological order of around 1% and provides an optimal set of candidates for experimental efforts.

Equipartition of Energy Defines the Size-Thickness Relationship in Liquid-Exfoliated Nanosheets

Liquid phase exfoliation is a commonly used method to produce 2D nanosheets from a range of layered crystals. However, such nanosheets display broad size and thickness distributions and correlations between area and thickness, issues that limit nanosheet application potential. To understand the factors controlling the exfoliation process, we have liquid-exfoliated 11 different layered materials, size-selecting each into fractions before using AFM to measure the nanosheet length, width, and thickness distributions for each fraction. The resultant data show a clear power-law scaling of nanosheet area with thickness for each material. We have developed a simple nonequilibrium thermodynamics-based model predicting that the power-law prefactor is proportional to both the ratios of in-plane-tearing/out-of-plane-peeling energies and in-plane/out-of-plane moduli. By comparing the experimental data with the modulus ratio calculated from first-principles, we find close agreement between experiment and theory. This supports our hypothesis that energy equipartition holds between nanosheet tearing and peeling during sonication-assisted exfoliation.

Valley-Engineering Mobilities in Two-Dimensional Materials

Two-dimensional materials are emerging as a promising platform for ultrathin channels in field-effect transistors. To this aim, novel high-mobility semiconductors need to be found or engineered. Although extrinsic mechanisms can in general be minimized by improving fabrication processes, the suppression of intrinsic scattering (driven, for example, by electron-phonon interactions) requires modification of the electronic or vibrational properties of the material. Because intervalley scattering critically affects mobilities, a powerful approach to enhance transport performance relies on engineering the valley structure. We show here the power of this strategy using uniaxial strain to lift degeneracies and suppress scattering into entire valleys, dramatically improving performance. This is shown in detail for arsenene, where a 2% strain stops scattering into four of the six valleys and leads to a 600% increase in mobility. The mechanism is general and can be applied to many other materials, including in particular the isostructural antimonene and blue phosphorene.

Koopmans Meets Bethe-Salpeter: Excitonic Optical Spectra without GW

The Bethe-Salpeter equation (BSE) can be applied to compute from first-principles optical spectra that include the effects of screened electron-hole interactions. As input, BSE calculations require single-particle states, quasiparticle energy levels, and the screened Coulomb interaction, which are typically obtained with many-body perturbation theory, whose cost limits the scope of possible applications. This work tries to address this practical limitation, instead deriving spectral energies from Koopmans-compliant functionals and introducing a new methodology for handling the screened Coulomb interaction. The explicit calculation of the W matrix is bypassed via a direct minimization scheme applied on top of a maximally localized Wannier function basis. We validate and benchmark this approach by computing the low-lying excited states of the molecules in Thiel's set and the optical absorption spectrum of a C₆₀ fullerene. The results show the same trends as quantum chemical methods and are in excellent agreement with previous simulations carried out at the time-dependent density functional theory or G₀ W₀-BSE level. Conveniently, the new framework reduces the parameter space controlling the accuracy of the calculation, thereby simplifying the simulation of charge-neutral excitations, offering the potential to expand the applicability of first-principles spectroscopies to larger systems of applied interest.

Continuum models of the electrochemical diffuse layer in electronic-structure calculations

Continuum electrolyte models represent a practical tool to account for the presence of the diffuse layer at electrochemical interfaces. However, despite the increasing popularity of these in the field of materials science, it remains unclear which features are necessary in order to accurately describe interface-related observables such as the differential capacitance (DC) of metal electrode surfaces. We present here a critical comparison of continuum diffuse-layer models that can be coupled to an atomistic first-principles description of the charged metal surface in order to account for the electrolyte screening at electrified interfaces. By comparing computed DC values for the prototypical Ag(100) surface in an aqueous solution to experimental data, we validate the accuracy of the models considered. Results suggest that a size-modified Poisson-Boltzmann description of the electrolyte solution is sufficient to qualitatively reproduce the main experimental trends. Our findings also highlight the large effect that the dielectric cavity parameterization has on the computed DC values.

Koopmans-Compliant Functionals and Potentials and Their Application to the GW100 Test Set

Koopmans-compliant (KC) functionals have been shown to provide accurate spectral properties through a generalized condition of piecewise linearity of the total energy as a function of the fractional addition/removal of an electron to/from any orbital. We analyze the performance of different KC functionals on a large and standardized set of 100 molecules, the GW100 test set, comparing vertical ionization potentials (taken as opposite of the orbital energies) to those obtained from accurate quantum chemistry methods, and to experimental results. We find excellent agreement, with a mean absolute error of 0.20 eV for the KIPZ functional on the first ionization potential, which is state-of-the-art for both density functional theory (DFT)-based calculations and many-body perturbation theory. We highlight similarities and differences between KC functionals and other electronic-structure approaches, such as dielectric-dependent hybrid functionals and Green's function methods, both from a theoretical and from a practical point of view, arguing that KC potentials can be considered as local and orbital-dependent approximations to the electronic self-energy, already including approximate vertex corrections.