1
|
Woo HM, Allam O, Chen J, Jang SS, Yoon BJ. Optimal high-throughput virtual screening pipeline for efficient selection of redox-active organic materials. iScience 2023; 26:105735. [PMID: 36582827 PMCID: PMC9793274 DOI: 10.1016/j.isci.2022.105735] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/22/2022] [Revised: 11/16/2022] [Accepted: 12/01/2022] [Indexed: 12/12/2022] Open
Abstract
As global interest in renewable energy continues to increase, there has been a pressing need for developing novel energy storage devices based on organic electrode materials that can overcome the shortcomings of the current lithium-ion batteries. One critical challenge for this quest is to find materials whose redox potential (RP) meets specific design targets. In this study, we propose a computational framework for addressing this challenge through the effective design and optimal operation of a high-throughput virtual screening (HTVS) pipeline that enables rapid screening of organic materials that satisfy the desired criteria. Starting from a high-fidelity model for estimating the RP of a given material, we show how a set of surrogate models with different accuracy and complexity may be designed to construct a highly accurate and efficient HTVS pipeline. We demonstrate that the proposed HTVS pipeline construction and operation strategies substantially enhance the overall screening throughput.
Collapse
Affiliation(s)
- Hyun-Myung Woo
- Computational Science Initiative, Brookhaven National Laboratory, Upton, NY 11973, USA
| | - Omar Allam
- School of Materials Science and Engineering, Georgia Institute of Technology, Atlanta, GA 30332, USA
- George W. Woodruff School of Mechanical Engineering, Georgia Institute of Technology, Atlanta, GA 30332, USA
| | - Junhe Chen
- School of Materials Science and Engineering, Georgia Institute of Technology, Atlanta, GA 30332, USA
| | - Seung Soon Jang
- School of Materials Science and Engineering, Georgia Institute of Technology, Atlanta, GA 30332, USA
| | - Byung-Jun Yoon
- Computational Science Initiative, Brookhaven National Laboratory, Upton, NY 11973, USA
- Department of Electrical and Computer Engineering, Texas A&M University, College Station, TX 77843, USA
| |
Collapse
|
2
|
Bethkenhagen M, Sharma A, Suryanarayana P, Pask JE, Sadigh B, Hamel S. Properties of carbon up to 10 million kelvin from Kohn-Sham density functional theory molecular dynamics. Phys Rev E 2023; 107:015306. [PMID: 36797894 DOI: 10.1103/physreve.107.015306] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/03/2021] [Accepted: 12/22/2022] [Indexed: 06/18/2023]
Abstract
Accurately modeling dense plasmas over wide-ranging conditions of pressure and temperature is a grand challenge critically important to our understanding of stellar and planetary physics as well as inertial confinement fusion. In this work, we employ Kohn-Sham density functional theory (DFT) molecular dynamics (MD) to compute the properties of carbon at warm and hot dense matter conditions in the vicinity of the principal Hugoniot. In particular, we calculate the equation of state (EOS), Hugoniot, pair distribution functions, and diffusion coefficients for carbon at densities spanning 8 g/cm^{3} to 16 g/cm^{3} and temperatures ranging from 100 kK to 10 MK using the Spectral Quadrature method. We find that the computed EOS and Hugoniot are in good agreement with path integral Monte Carlo results and the sesame database. Additionally, we calculate the ion-ion structure factor and viscosity for selected points. All results presented are at the level of full Kohn-Sham DFT-MD, free of empirical parameters, average-atom, and orbital-free approximations employed previously at such conditions.
Collapse
Affiliation(s)
- Mandy Bethkenhagen
- Lawrence Livermore National Laboratory, Livermore, California 94550, USA
- École Normale Supérieure de Lyon, Université Lyon 1, Laboratoire de Géologie de Lyon, CNRS UMR 5276, 69364 Lyon, Cedex 07, France
| | - Abhiraj Sharma
- College of Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332, USA
| | - Phanish Suryanarayana
- College of Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332, USA
| | - John E Pask
- Lawrence Livermore National Laboratory, Livermore, California 94550, USA
| | - Babak Sadigh
- Lawrence Livermore National Laboratory, Livermore, California 94550, USA
| | - Sebastien Hamel
- Lawrence Livermore National Laboratory, Livermore, California 94550, USA
| |
Collapse
|
3
|
Clérouin J, Blanchet A, Blancard C, Faussurier G, Soubiran F, Bethkenhagen M. Equivalence between pressure- and structure-defined ionization in hot dense carbon. Phys Rev E 2022; 106:045204. [PMID: 36397512 DOI: 10.1103/physreve.106.045204] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/23/2022] [Accepted: 09/23/2022] [Indexed: 06/16/2023]
Abstract
The determination of the ionization of a system in the hot dense regime is a long standing issue. Recent studies have shown inconsistencies between standard predictions using average atom models and evaluations deduced from electronic transport properties computed with quantum molecular dynamics simulations [Bethkenhagen et al., Phys. Rev. Res. 2, 023260 (2020)]2643-156410.1103/PhysRevResearch.2.023260. Here, we propose a definition of the ionization based on its effect on the plasma structure as given by the pair distribution function (PDF), and on the concept of effective one-component plasma (eOCP). We also introduce a definition based on the total pressure and on a modelization of the electronic pressure. We show the equivalence of these definitions on two studies of carbon along the 100 eV isotherm and the 10 g/cm^{3} isochor. Simulations along the 100 eV isotherm are obtained with the newly implemented Ext. First principles molecular dynamics (Fpmd) method in Abinit for densities ranging from 1 to 500 g/cm^{3}and along the 10 g/cm^{3} isochor with the recently published Spectral quadrature DFT (Sqdft) simulations, between 8 and 860 eV. The resulting ionizations are compared to the predictions of the average-atom code Qaam which is based on the muffin-tin approximation. A disagreement between the eOCP and the actual PDFs (non-OCP behavior) is interpreted as the onset of bonding in the system.
Collapse
Affiliation(s)
- Jean Clérouin
- CEA-DAM-DIF, F-91297 Arpajon, France
- Université Paris-Saclay, CEA, Laboratoire Matière sous conditions extrêmes, 91680 Bruyères-le-Châtel, France
| | - Augustin Blanchet
- CEA-DAM-DIF, F-91297 Arpajon, France
- Université Paris-Saclay, CEA, Laboratoire Matière sous conditions extrêmes, 91680 Bruyères-le-Châtel, France
| | - Christophe Blancard
- CEA-DAM-DIF, F-91297 Arpajon, France
- Université Paris-Saclay, CEA, Laboratoire Matière sous conditions extrêmes, 91680 Bruyères-le-Châtel, France
| | - Gérald Faussurier
- CEA-DAM-DIF, F-91297 Arpajon, France
- Université Paris-Saclay, CEA, Laboratoire Matière sous conditions extrêmes, 91680 Bruyères-le-Châtel, France
| | - François Soubiran
- CEA-DAM-DIF, F-91297 Arpajon, France
- Université Paris-Saclay, CEA, Laboratoire Matière sous conditions extrêmes, 91680 Bruyères-le-Châtel, France
| | - Mandy Bethkenhagen
- CNRS, École Normale Supérieure de Lyon, Laboratoire de Géologie de Lyon LGLTPE UMR 5276, Centre Blaise Pascal, 46 allée d'Italie Lyon 69364, France
| |
Collapse
|
4
|
Wu CJ, Myint PC, Pask JE, Prisbrey CJ, Correa AA, Suryanarayana P, Varley JB. Development of a Multiphase Beryllium Equation of State and Physics-based Variations. J Phys Chem A 2021; 125:1610-1636. [PMID: 33587640 DOI: 10.1021/acs.jpca.0c09809] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We construct a family of beryllium (Be) multiphase equation of state (EOS) models that consists of a baseline ("optimal") EOS and variations on the baseline to account for physics-based uncertainties. The Be baseline EOS is constructed to reproduce a set of self-consistent data and theory including known phase boundaries, the principal Hugoniot, isobars, and isotherms from diamond-anvil cell experiments. Three phases are considered, including the known hexagonal closed-packed (hcp) phase, the liquid, and the theoretically predicted high-pressure body-centered cubic (bcc) phase. Since both the high-temperature liquid and high-pressure bcc phases lack any experimental data, we carry out ab initio density functional theory (DFT) calculations to obtain new information about the EOS properties for these two regions. At extremely high temperature conditions (>87 eV), DFT-based quantum molecular dynamics simulations are performed for multiple liquid densities using the state-of-the-art Spectral Quadrature methodology in order to validate our selected models for the ion- and electron-thermal free energies of the liquid. We have also performed DFT simulations of hcp and bcc with different exchange-correlation functionals to examine their impact on bcc compressibility, which bound the hcp-bcc transition pressure to within 4 ± 0.5 Mbar. Our baseline EOS predicts the first density maximum along the Hugoniot to be 4.4-fold in compression, while the hcp-bcc-liquid triple-point pressure is predicted to be at 2.25 Mbar. In addition to the baseline EOS, we have generated eight variations to accommodate multiple sources of potential uncertainties such as (1) the choice of free-energy models, (2) differences in theoretical treatments, (3) experimental uncertainties, and (4) lack of information. These variations are designed to provide a reasonable representation of nonstatistical uncertainties for the Be EOS and may be used to assess its sensitivity to different inertial-confinement fusion capsule designs.
Collapse
Affiliation(s)
- Christine J Wu
- Lawrence Livermore National Laboratory, Livermore, California 94588, United States
| | - Philip C Myint
- Lawrence Livermore National Laboratory, Livermore, California 94588, United States
| | - John E Pask
- Lawrence Livermore National Laboratory, Livermore, California 94588, United States
| | - Carrie J Prisbrey
- Lawrence Livermore National Laboratory, Livermore, California 94588, United States
| | - Alfredo A Correa
- Lawrence Livermore National Laboratory, Livermore, California 94588, United States
| | - Phanish Suryanarayana
- College of Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332, United States
| | - Joel B Varley
- Lawrence Livermore National Laboratory, Livermore, California 94588, United States
| |
Collapse
|
5
|
Luo Z, Qin X, Wan L, Hu W, Yang J. Parallel Implementation of Large-Scale Linear Scaling Density Functional Theory Calculations With Numerical Atomic Orbitals in HONPAS. Front Chem 2020; 8:589910. [PMID: 33324611 PMCID: PMC7726133 DOI: 10.3389/fchem.2020.589910] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/31/2020] [Accepted: 09/08/2020] [Indexed: 11/13/2022] Open
Abstract
Linear-scaling density functional theory (DFT) is an efficient method to describe the electronic structures of molecules, semiconductors, and insulators to avoid the high cubic-scaling cost in conventional DFT calculations. Here, we present a parallel implementation of linear-scaling density matrix trace correcting (TC) purification algorithm to solve the Kohn-Sham (KS) equations with the numerical atomic orbitals in the HONPAS package. Such a linear-scaling density matrix purification algorithm is based on the Kohn's nearsightedness principle, resulting in a sparse Hamiltonian matrix with localized basis sets in the DFT calculations. Therefore, sparse matrix multiplication is the most time-consuming step in the density matrix purification algorithm for linear-scaling DFT calculations. We propose to use the MPI_Allgather function for parallel programming to deal with the sparse matrix multiplication within the compressed sparse row (CSR) format, which can scale up to hundreds of processing cores on modern heterogeneous supercomputers. We demonstrate the computational accuracy and efficiency of this parallel density matrix purification algorithm by performing large-scale DFT calculations on boron nitrogen nanotubes containing tens of thousands of atoms.
Collapse
Affiliation(s)
| | - Xinming Qin
- Hefei National Laboratory for Physical Sciences at the Microscale, Department of Chemical Physics, and Synergetic Innovation Center of Quantum Information and Quantum Physics, University of Science and Technology of China, Hefei, China
| | | | - Wei Hu
- Hefei National Laboratory for Physical Sciences at the Microscale, Department of Chemical Physics, and Synergetic Innovation Center of Quantum Information and Quantum Physics, University of Science and Technology of China, Hefei, China
| | | |
Collapse
|
6
|
Sharma A, Hamel S, Bethkenhagen M, Pask JE, Suryanarayana P. Real-space formulation of the stress tensor for O(N) density functional theory: Application to high temperature calculations. J Chem Phys 2020; 153:034112. [PMID: 32716199 DOI: 10.1063/5.0016783] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present an accurate and efficient real-space formulation of the Hellmann-Feynman stress tensor for O(N) Kohn-Sham density functional theory (DFT). While applicable at any temperature, the formulation is most efficient at high temperature where the Fermi-Dirac distribution becomes smoother and the density matrix becomes correspondingly more localized. We first rewrite the orbital-dependent stress tensor for real-space DFT in terms of the density matrix, thereby making it amenable to O(N) methods. We then describe its evaluation within the O(N) infinite-cell Clenshaw-Curtis Spectral Quadrature (SQ) method, a technique that is applicable to metallic and insulating systems, is highly parallelizable, becomes increasingly efficient with increasing temperature, and provides results corresponding to the infinite crystal without the need of Brillouin zone integration. We demonstrate systematic convergence of the resulting formulation with respect to SQ parameters to exact diagonalization results and show convergence with respect to mesh size to the established plane wave results. We employ the new formulation to compute the viscosity of hydrogen at 106 K from Kohn-Sham quantum molecular dynamics, where we find agreement with previous more approximate orbital-free density functional methods.
Collapse
Affiliation(s)
- Abhiraj Sharma
- College of Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332, USA
| | - Sebastien Hamel
- Physics Division, Lawrence Livermore National Laboratory, Livermore, California 94550, USA
| | - Mandy Bethkenhagen
- College of Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332, USAPhysics Division, Lawrence Livermore National Laboratory, Livermore, California 94550, USACNRS, École Normale Supérieure de Lyon, Laboratoire de Géologie de Lyon LGLTPE UMR5276, Centre Blaise Pascal, 46 Allée D'Italie, Lyon 69364, France
| | - John E Pask
- Physics Division, Lawrence Livermore National Laboratory, Livermore, California 94550, USA
| | - Phanish Suryanarayana
- College of Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332, USA
| |
Collapse
|
7
|
Sharma A, Suryanarayana P. On the calculation of the stress tensor in real-space Kohn-Sham density functional theory. J Chem Phys 2018; 149:194104. [PMID: 30466280 DOI: 10.1063/1.5057355] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023] Open
Abstract
We present an accurate and efficient formulation of the stress tensor for real-space Kohn-Sham density functional theory calculations. Specifically, while employing a local formulation of the electrostatics, we derive a linear-scaling expression for the stress tensor that is applicable to simulations with unit cells of arbitrary symmetry, semilocal exchange-correlation functionals, and Brillouin zone integration. In particular, we rewrite the contributions arising from the self-energy and the nonlocal pseudopotential energy to make them amenable to the real-space finite-difference discretization, achieving up to three orders of magnitude improvement in the accuracy of the computed stresses. Using examples representative of static and dynamic calculations, we verify the accuracy and efficiency of the proposed formulation. In particular, we demonstrate high rates of convergence with spatial discretization, consistency between the computed energy and the stress tensor, and very good agreement with reference planewave results.
Collapse
Affiliation(s)
- Abhiraj Sharma
- College of Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332, USA
| | - Phanish Suryanarayana
- College of Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332, USA
| |
Collapse
|
8
|
Banerjee AS, Lin L, Suryanarayana P, Yang C, Pask JE. Two-Level Chebyshev Filter Based Complementary Subspace Method: Pushing the Envelope of Large-Scale Electronic Structure Calculations. J Chem Theory Comput 2018; 14:2930-2946. [PMID: 29660292 DOI: 10.1021/acs.jctc.7b01243] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We describe a novel iterative strategy for Kohn-Sham density functional theory calculations aimed at large systems (>1,000 electrons), applicable to metals and insulators alike. In lieu of explicit diagonalization of the Kohn-Sham Hamiltonian on every self-consistent field (SCF) iteration, we employ a two-level Chebyshev polynomial filter based complementary subspace strategy to (1) compute a set of vectors that span the occupied subspace of the Hamiltonian; (2) reduce subspace diagonalization to just partially occupied states; and (3) obtain those states in an efficient, scalable manner via an inner Chebyshev filter iteration. By reducing the necessary computation to just partially occupied states and obtaining these through an inner Chebyshev iteration, our approach reduces the cost of large metallic calculations significantly, while eliminating subspace diagonalization for insulating systems altogether. We describe the implementation of the method within the framework of the discontinuous Galerkin (DG) electronic structure method and show that this results in a computational scheme that can effectively tackle bulk and nano systems containing tens of thousands of electrons, with chemical accuracy, within a few minutes or less of wall clock time per SCF iteration on large-scale computing platforms. We anticipate that our method will be instrumental in pushing the envelope of large-scale ab initio molecular dynamics. As a demonstration of this, we simulate a bulk silicon system containing 8,000 atoms at finite temperature, and obtain an average SCF step wall time of 51 s on 34,560 processors; thus allowing us to carry out 1.0 ps of ab initio molecular dynamics in approximately 28 h (of wall time).
Collapse
Affiliation(s)
- Amartya S Banerjee
- Computational Research Division , Lawrence Berkeley National Laboratory , Berkeley , California 94720 , United States
| | - Lin Lin
- Computational Research Division , Lawrence Berkeley National Laboratory , Berkeley , California 94720 , United States.,Department of Mathematics , University of California , Berkeley , California 94720 , United States
| | - Phanish Suryanarayana
- College of Engineering , Georgia Institute of Technology , Atlanta , Georgia 30332 , United States
| | - Chao Yang
- Computational Research Division , Lawrence Berkeley National Laboratory , Berkeley , California 94720 , United States
| | - John E Pask
- Physics Division , Lawrence Livermore National Laboratory , Livermore , California 94550 , United States
| |
Collapse
|
9
|
Aarons J, Skylaris CK. Electronic annealing Fermi operator expansion for DFT calculations on metallic systems. J Chem Phys 2018; 148:074107. [PMID: 29471650 DOI: 10.1063/1.5001340] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Jolyon Aarons
- Department of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| | - Chris-Kriton Skylaris
- Department of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, United Kingdom
| |
Collapse
|
10
|
Abstract
Recent developments in the chemistry of the transuranic elements are surveyed, with particular emphasis on computational contributions. Examples are drawn from molecular coordination and organometallic chemistry, and from the study of extended solid systems. The role of the metal valence orbitals in covalent bonding is a particular focus, especially the consequences of the stabilization of the 5f orbitals as the actinide series is traversed. The fledgling chemistry of transuranic elements in the +II oxidation state is highlighted. Throughout, the symbiotic interplay of experimental and computational studies is emphasized; the extraordinary challenges of experimental transuranic chemistry afford computational chemistry a particularly valuable role at the frontier of the periodic table.
Collapse
Affiliation(s)
- Nikolas Kaltsoyannis
- School of Chemistry, The University of Manchester, Oxford Road, Manchester, M13 9PL, UK
| |
Collapse
|