Correction: Atomate2: modular workflows for materials science
[This corrects the article DOI: 10.1039/D5DD00019J.].
Ganose, Alex M↗
Engineering topics
Sundararaman, Ravishankar
Publications and source records attributed to Sundararaman, Ravishankar.
BEAST DB: Grand-Canonical Database of Electrocatalyst Properties
We present BEAST DB, an open-source database comprised of ab initio electrochemical data computed using grand-canonical density functional theory in implicit solvent at consistent calculation parameters. The database contains over 20,000 surface calculations and covers a broad set of heterogeneous catalyst materials and electrochemical reactions. Calculations were performed at self-consistent fixed potential as well as constant charge to facilitate comparisons to the computational hydrogen electrode. This article presents common use cases of the database to rationalize trends in catalyst activity, screen catalyst material spaces, understand elementary mechanistic steps, analyze the electronic structure, and train machine learning models to predict higher fidelity properties. Users can interact graphically with the database by querying for individual calculations to gain a granular understanding of reaction steps or by querying for an entire reaction pathway on a given material using an interactive reaction pathway tool. BEAST DB will be periodically updated, with planned future updates to include advanced electronic structure data, surface speciation studies, and greater reaction coverage.
database↗
Static Subspace Approximation for Random Phase Approximation Correlation Energies: Implementation and Performance
Developing theoretical understanding of complex reactions and processes at interfaces requires using methods that go beyond semilocal density functional theory to accurately describe the interactions between solvent, reactants and substrates. Methods based on many-body perturbation theory, such as the random phase approximation (RPA), have previously been limited due to their computational complexity. However, this is now a surmountable barrier due to the advances in computational power available, in particular through modern GPU-based supercomputers. In this work, we describe the implementation of RPA calculations within BerkeleyGW and show its favorable computational performance on large complex systems relevant for catalysis and electrochemistry applications. Our implementation builds off of the static subspace approximation which, by employing a compressed representation of the frequency dependent polarizability, enables the evaluation of the RPA correlation energy with significant acceleration and systematically controllable accuracy. We find that the computational cost of calculating the RPA correlation energy scales only linearly with system size for systems containing up to 50 thousand bands, and is expected to scale quadratically thereafter. We also show excellent strong scaling results across several supercomputers, demonstrating the performance and portability of this implementation.
algorithmic development↗
Improving the reliability of machine learned potentials for modeling inhomogeneous liquids
The atomic-scale response of inhomogeneous fluids at interfaces and surrounding solute particles plays a critical role in governing chemical, electrochemical, and biological processes. Classical molecular dynamics simulations have been applied extensively to simulate the response of fluids to inhomogeneities directly, but are limited by the accuracy of the underlying interatomic potentials. Here, we use neural network potentials (NNPs) trained to ab initio simulations to accurately predict the inhomogeneous responses of two distinct fluids: liquid water and molten NaCl. Although NNPs can be readily trained to model complex bulk systems across a range of state points, we show that to appropriately model a fluid's response at an interface, relevant inhomogeneous configurations must be included in the training data. In order to sufficiently sample appropriate configurations of such inhomogeneous fluids, we develop protocols based on molecular dynamics simulations in the presence of external potentials. We demonstrate that NNPs trained on inhomogeneous fluid configurations can more accurately predict several key properties of fluids—including the density response, surface tension and size-dependent cavitation free energies—for liquid water and molten NaCl, compared to both empirical interatomic potentials and NNPs that are not trained on such inhomogeneous configurations. This work therefore provides a first demonstration and framework to extract the response of inhomogeneous fluids from first principles for classical density-functional treatment of fluids free from empirical potentials.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Highly confined, low-loss plasmonics based on two-dimensional solid-state defect lattices
Not Available
Ghorashi, Ali (ORCID:0000000258841756)↗
How spin relaxes and dephases in bulk halide perovskites
Abstract Spintronics in halide perovskites has drawn significant attention in recent years, due to their highly tunable spin-orbit fields and intriguing interplay with lattice symmetry. Here, we perform first-principles calculations to determine the spin relaxation time ( T 1 ) and ensemble spin dephasing time ( $${T}_{2}^{*}$$ T 2 * ) in a prototype halide perovskite, CsPbBr 3 . To accurately capture spin dephasing in external magnetic fields we determine the Landé g -factor from first principles and take it into account in our calculations. These allow us to predict intrinsic spin lifetimes as an upper bound for experiments, identify the dominant spin relaxation pathways, and evaluate the dependence on temperature, external fields, carrier density, and impurities. We find that the Fröhlich interaction that dominates carrier relaxation contributes negligibly to spin relaxation, consistent with the spin-conserving nature of this interaction. Our theoretical approach may lead to new strategies to optimize spin and carrier transport properties.
74 ATOMIC AND MOLECULAR PHYSICS↗
2023 Roadmap on molecular modelling of electrochemical energy materials
New materials for electrochemical energy storage and conversion are the key to the electrification and sustainable development of our modern societies. Molecular modelling based on the principles of quantum mechanics and statistical mechanics as well as empowered by machine learning techniques can help us to understand, control and design electrochemical energy materials at atomistic precision. Therefore, this roadmap, which is a collection of authoritative opinions, serves as a gateway for both the experts and the beginners to have a quick overview of the current status and corresponding challenges in molecular modelling of electrochemical energy materials for batteries, supercapacitors, CO 2 reduction reaction, and fuel cell applications.
25 ENERGY STORAGE↗
Accurate Prediction of HSE06 Band Structures for a Diverse Set of Materials Using Δ-Learning
Here we used machine learning (ML) to accurately predict eigenvalues of the hybrid HSE06 functional using eigenvalues computed by the less computationally expensive PBE functional and associated electronic features based on the k-point resolved atomic band character. The ML model was trained by using eigenvalues from only one k-point for each of the 168 compounds in the training set. The HSE06 eigenvalues across all k-points were then predicted for a separate set of 169 compounds with a mean absolute error (MAE) of 0.13 eV, representing a significant improvement over the error of PBE-computed eigenvalues relative to that of HSE06 (MAE = 0.96 eV). These accurately predicted eigenvalues result in remarkably accurate predictions for the band structures, projected density of states, and band gaps, even though the model was not explicitly trained on these other properties. Finally, we demonstrate that our ML model has a similar accuracy for both ternary and quaternary compounds well outside the initial training set and on systems with 112 and 160 atoms, demonstrating its potential to rapidly predict HSE06-quality electronic structures of complex materials that are practically unfeasible for HSE06.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
First-principles molten salt phase diagrams through thermodynamic integration
Precise prediction of phase diagrams in molecular dynamics simulations is challenging due to the simultaneous need for long time and large length scales and accurate interatomic potentials. Here, we show that thermodynamic integration from low-cost force fields to neural network potentials trained using density-functional theory (DFT) enables rapid first-principles prediction of the solid–liquid phase boundary in the model salt NaCl. We use this technique to compare the accuracy of several DFT exchange–correlation functionals for predicting the NaCl phase boundary and find that the inclusion of dispersion interactions is critical to obtain good agreement with experiment. Importantly, our approach introduces a method to predict solid–liquid phase boundaries for any material at an ab initio level of accuracy, with the majority of the computational cost at the level of classical potentials.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Revised Nitrogen Reduction Scaling Relations from Potential-Dependent Modeling of Chemical and Electrochemical Steps
The electrochemical nitrogen reduction reaction (NRR) is a promising route to enable carbon-free ammonia production. However, this reaction is limited by the poor activity and selectivity of current catalysts. The rational design of superior NRR electrocatalysts requires a detailed mechanistic understanding of current material limitations to inform how these might be overcome. The current understanding of how scaling limits NRR on metal catalysts is predicated on a simplified reaction pathway that considers only proton-coupled electron transfer (PCET) steps. Here, we apply grand-canonical density functional theory to investigate a more comprehensive NRR mechanism that includes both electrochemical and chemical steps on 30 metal surfaces in solvent under an applied potential. We applied ..phi..max, a grand-canonical adaptation of the Gmax thermodynamic descriptor, to evaluate trends in catalyst activity. This approach produces a ..phi..max "volcano" diagram for NRR activity scaling on metals that qualitatively differs from the scaling relations identified when only PCET steps are considered. NH3* desorption was found to limit the NRR activity for materials at the top of the volcano and truncate the volcano's peak at increasingly reducing potentials. These revised scaling relations may inform the rational design of superior NRR electrocatalysts. This approach is transferable to study other materials and reaction chemistries where both electrochemical and chemical steps are modeled under an applied potential.
electrocatalysis↗
Impact of solvation on the GW quasiparticle spectra of molecules
First-principles calculations for electrochemistry require accurate treatment of both electronic structure and solvation. The perturbative GW approximation starting from density functional theory (DFT) calculations accurately models materials systems with varying dimensionality. Continuum solvation models enable efficient treatment of solvation effects in DFT calculations, but their applications with beyond-DFT electronic structure methods such as GW have been limited. Here, we introduce the frequency-dependent liquid polarizability from a nonlocal continuum solvation model in the screened Coulomb interaction of full-frequency GW calculations with a solvated DFT starting point. In this study, we show that the liquid screening contributions substantially reduce the HOMO–LUMO gap of molecules by 3–5 eV, while solvent effects on the DFT starting point negligibly impact the GW gap. The resulting framework facilitates the simultaneous electronic and solvation accuracy needed for first-principles electrochemistry.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Spin-selective photocatalysis and quantum transport using ab initio density-matrix dynamics
The overarching goal of this research program is to quantitatively predict out-of-equilibrium quantum dynamics and transport with electrons and spin degrees of freedom coupled to phonons, photons, defects and liquid environments, specifically targeting spin-selective chemistry. Spin polarized carriers, which can be generated at room temperature by circularly-polarized light, provide unique opportunities for creating new reaction pathways and increasing selectivity towards specific pathways. If realized, this can substantially increase the level of reaction control possible in the field of heterogeneous catalysis. Spin, a quantum mechanical entity without a classical counterpart, requires quantum dynamics. Capturing both coherent and incoherent dynamics and accurate correlation of spins, electrons and phonon environments are essential to quantitative description of generation, relaxation and transport, as well as chemical reaction of spin-polarized carriers in molecules and materials promising for spin-photocatalysis.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Simulating dielectric spectra: A demonstration of the direct electric field method and a new model for the nonlinear dielectric response
Here we demonstrate a method to compute the dielectric spectra of fluids in molecular dynamics (MD) by directly applying electric fields to the simulation. We obtain spectra from MD simulations with low magnitude electric fields (≈0.01 V/Å) in agreement with spectra from the fluctuation–dissipation method for water and acetonitrile. We examine this method’s trade-off between noise at low field magnitudes and the nonlinearity of the response at higher field magnitudes. We then apply the Booth equation to describe the nonlinear response of both fluids at low frequency (0.1 GHz) and high field magnitude (up to 0.5 V/Å). We develop a model of the frequency-dependent nonlinear response by combining the Booth description of the static nonlinear dielectric response of fluids with the frequency-dependent linear dielectric response of the Debye model. We find good agreement between our model and the MD simulations of the nonlinear dielectric response for both acetonitrile and water.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Solvent effects determine the sign of the charges of maximum entropy and capacitance at silver electrodes
Fully harnessing electrochemical interfaces for reactions requires a detailed understanding of solvent effects in the electrochemical double layer. Predicting the significant impact of solvents on entropic and electronic properties of electrochemical interfaces has remained an open challenge of computational electrochemistry. Using molecular dynamics simulations of silver–water and silver–acetonitrile interfaces, we show that switching the solvent changes the signs for both the charge of maximum capacitance (CMC) and charge of maximum entropy (CME). Contrasting the capacitance and CME behavior of these two interfaces, we demonstrate that the preferred orientation of the solvent molecule and the corresponding charge density determine the sign of the CMC and CME and, hence, the qualitatively different charge asymmetry of the electrochemical interface.
Chemistry↗
First-principles free energies by hybrid thermodynamic integration for phase equilibria and fission product solubility in molten salts
The overarching goal of this research project is quantitative microscopic modeling of phase diagrams and solubility of actinides and lanthanide-fission products in molten salts for molten salt reactors. Application of molten salts as coolants and liquid fuels for next-generation reactors requires quantitative understanding of phase diagrams, equations of state and transport properties of binary and ternary molten salt mixtures, as well as solubility of fuel, fission products and reactor-wall corrosion products, but these remain beyond the reach of current first-principles computational techniques. We will address this fundamental capability gap using a novel hybrid framework for thermodynamic integration to directly predict absolute free energies of solid and molten salts at any specified conditions of composition, temperature and pressure. This framework will combine a hierarchy of first-principles, molecular dynamics and continuum techniques to deliver first-principles-level accuracy at a fraction of the computational cost of conventional ab initio molecular dynamics approaches. In particular, this technique is expected to provide first-principles accuracy, while off-loading most of the computational work on much less expensive classical molecular dynamics and continuum techniques.
36 MATERIALS SCIENCE↗
Generalized design principles for hydrodynamic electron transport in anisotropic metals
Interactions of charge carriers with lattice vibrations, or phonons, play a critical role in unconventional electronic transport of metals and semimetals. Recent observations of phonon-mediated collective electron flow in bulk semimetals, termed electron hydrodynamics, present new opportunities in the search for strong electron-electron interactions in high carrier density materials. Here we present the general transport signatures of such a second-order scattering mechanism, along with analytical limits at the Eliashberg level of theory. We study electronic transport, using ab initio calculations, in finite-size channels of semimetallic ZrSiS and TaAs 2 with and without topological band crossings, respectively. The order of magnitude separation between momentum-relaxing and momentum-conserving scattering length scales across a wide temperature range make both of them promising candidates for further experimental observation of electron hydrodynamics. More generally, our calculations suggest that the hydrodynamic transport regime does not, to first order, rely on the topological nature of the bands. Finally, we discuss general design principles guiding future search for hydrodynamic candidates, based on the analytical formulation and our ab initio predictions. Further, we find that systems with strong electron-phonon interactions, reduced electronic phase space, and suppressed phonon-phonon scattering at temperatures of interest are likely to feature hydrodynamic electron transport. We predict that layered and/or anisotropic semimetals composed of half-filled d shells and light group V/VI elements with lower crystal symmetry are promising candidates to observe hydrodynamic phenomena in the future.
36 MATERIALS SCIENCE↗
Improving the Accuracy of Atomistic Simulations of the Electrochemical Interface
Atomistic simulation of the electrochemical double layer is an ambitious undertaking, requiring quantum mechanical description of electrons, phase space sampling of liquid electrolytes, and equilibration of electrolytes over nanosecond time scales. All models of electrochemistry make different trade-offs in the approximation of electrons and atomic configurations, from the extremes of classical molecular dynamics of a complete interface with point-charge atoms to correlated electronic structure methods of a single electrode configuration with no dynamics or electrolyte. Here, we review the spectrum of simulation techniques suitable for electrochemistry, focusing on the key approximations and accuracy considerations for each technique. Additionally, we discuss promising approaches, such as enhanced sampling techniques for atomic configurations and computationally efficient beyond density functional theory (DFT) electronic methods, that will push electrochemical simulations beyond the present frontier.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗