Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density functional theory development”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,225 records · Page 68
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Size-Dependent Energy of Ni Nanoparticles on Graphene Films on Ni(111) and Adhesion Energetics by Adsorption Calorimetry

The use of carbon supports for late transition-metal nanoparticle catalysts has grown substantially in recent years due to efforts to develop electrocatalysts for clean energy applications and catalysts for new aqueous-phase biomass-related conversions and due to the evolution of new carbon materials with unique properties (e.g., graphene, carbon nanotubes, and so forth). However, much less is known about the bonding energetics of catalytic metal nanoparticles on carbon supports in comparison with oxide supports, which are more common for thermal catalysis. Here, we report the growth morphology and heats of adsorption of Ni vapor deposited onto graphene/Ni(111) at 300 K and 100 K using metal vapor single-crystal adsorption calorimetry and He + low-energy ion scattering (LEIS). These results provide the Ni chemical potential versus particle size, and the Ni/graphene adhesion energy. LEIS intensities suggest that Ni grows as flat-topped face-centered cubic islands with a nearly constant thickness of ~1.5 nm when deposited at 300 K. At 100 K, Ni grows as smaller nanoparticles, well modeled as hemispherical hexagonal close-packed nanoparticles with a density of ~2 × 10 16 particles/m 2 . The Ni chemical potential as a function of average particle diameter in the 0.5 to 4 nm range at 100 K was determined from the heats of Ni gas adsorption. Further, by fitting the measured chemical potential as a function of diameter, we determined an adhesion energy of 3.6 J/m 2 for large Ni particles on graphene/Ni(111). This adhesion energy is in good agreement with previous scanning tunneling microscopy and density functional theory investigations of Ni/graphene/Ni(111).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Noise and error analysis and optimization in particle-based kinetic plasma simulations

In this paper we analyze the noise in macro-particle methods used in plasma physics and fluid dynamics, leading to approaches for minimizing the total error, focusing on electrostatic models in one dimension. We begin by describing kernel density estimation for continuous values of the spatial variable x, expressing the kernel in a form in which its shape and width are represented separately. The covariance matrix of the noise in the density is computed, first for uniform true density. The bandwidth of the covariance matrix C(x,y) is related to the width of the kernel. A feature that stands out is the presence of constant negative terms in the elements of the covariance matrix both on and off-diagonal. These negative correlations are related to the fact that the total number of particles is fixed at each time step; they also lead to the property ∫C(x,y)dy = 0. We investigate the effect of these negative correlations on the electric field computed by Gauss's law, finding that the noise in the electric field is related to a process called the Ornstein-Uhlenbeck bridge, leading to a covariance matrix of the electric field with variance significantly reduced relative to that of a Brownian process. For non-constant density, p(x), still with continuous x, we analyze the total error in the density estimation and discuss it in terms of bias-variance optimization (BVO). For some characteristic length l, determined by the density and its second derivative, and kernel width h, having too few particles within h leads to too much variance; for h that is large relative to l, there is too much smoothing of the density. The optimum between these two limits is found by BVO. For kernels of the same width, it is shown that this optimum (minimum) is weakly sensitive to the kernel shape. Next, we repeat the analysis for x discretized on a grid. In this case the charge deposition rule is determined by a particle shape. An important property to be respected in the discrete system is the exact preservation of total charge on the grid; this property is necessary to ensure that the electric field is equal at both ends, consistent with periodic boundary conditions. We find that if the particle shapes satisfy a partition of unity property, the particle charge deposited on the grid is conserved exactly. Further, if the particle shape is expressed as the convolution of a kernel with another kernel that satisfies the partition of unity, then the particle shape obeys the partition of unity. This property holds for kernels of arbitrary width, including widths that are not integer multiples of the grid spacing. Furthermore, we show results relaxing the approximations used to do BVO optimization analytically, by doing numerical computations of the total error as a function of the kernel width, on a grid in x. The comparison between numerical and analytical results shows good agreement over a range of particle shapes. We discuss the practical implications of our results, including the criteria for design and implementation of computationally efficient particle shapes that take advantage of the developed theory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Graphene-supported single atom catalysts for high performance lithium-oxygen batteries

The optimal choice of d-block metals in single atom catalysts (SACs) is crucial for designing efficient electrocatalysts for activating the Oxygen reduction reaction (ORR)/ Oxygen evolution reaction (OER) in lithiumoxygen batteries (LOBs). Herein, we used the Quantum Mechanics methods to understand the origin of reactivity for a series of 16 d-block metals supported on nitrogen-doped graphene as SACs for ORR and OER in LOBs. Based on the Gibbs free energy calculations, we found that among the 16 SACs investigated, Zn-SAC exhibits the highest electrochemical activity with the lowest overpotential of 0.17 V. Here we then used machine learning (ML) to develop an intrinsic descriptor, phi, that correlates the catalytic activity with electronic and chemical properties of the catalytic centers at the M-N 4 active site on graphene surface. We established a linear relationship between phi and the catalytic activity that provides guidance for designing efficient SACs for electrocatalysis in LOBs. To validate these predictions, we report electrochemical measurements showing that Zn-SAC exhibits an ultra-stable cyclability with reduced overpotentials over Mo-SAC and nitrogen-doped graphene (NG), confirming our theoretical prediction. This fundamental work provides a deep understanding on the rational design of efficient SACs for OER/ ORR in LOBs.

25 ENERGY STORAGE↗
DOE OSTI2024

Efficient Mixed-Precision Matrix Factorization of the Inverse Overlap Matrix in Electronic Structure Calculations with AI-Hardware and GPUs

In recent years, a new kind of accelerated hardware has gained popularity in the artificial intelligence (AI) community which enables extremely high-performance tensor contractions in reduced precision for deep neural network calculations. In this article, we exploit Nvidia Tensor cores, a prototypical example of such AI-hardware, to develop a mixed precision approach for computing a dense matrix factorization of the inverse overlap matrix in electronic structure theory, S –1 . This factorization of S –1 , written as ZZT = S –1 , is used to transform the general matrix eigenvalue problem into a standard matrix eigenvalue problem. Here we present a mixed precision iterative refinement algorithm where Z is given recursively using matrix–matrix multiplications and can be computed with high performance on Tensor cores. To understand the performance and accuracy of Tensor cores, comparisons are made to GPU-only implementations in single and double precision. Additionally, we propose a nonparametric stopping criteria which is robust in the face of lower precision floating point operations. The algorithm is particularly useful when we have a good initial guess to Z, for example, from previous time steps in quantum-mechanical molecular dynamics simulations or from a previous iteration in a geometry optimization.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Magnetocaloric Effect in Lightly‐Doped Fe 5 Si 3 Single Crystals

Abstract Development of promising new materials for above room temperature magnetic cooling applications relies on careful balancing of structure and composition to maximize accessible metastable phases that can drive a strong magnetocaloric effect (MCE). However, the working temperatures of these materials may fall outside of desired application windows. In this work, it is shown that it is possible to control metastable phase stability temperatures of Fe 5 Si 3 through selection of appropriate spin and charge doping. Here, the parent material's desired structure appears only within a narrow temperature range from 1098 to 1303 K. Doping with Mn and P is shown to allow stabilization of the parent's high temperature phase and resulting MCE to room temperature. The structural and magnetic properties, and the magnetocaloric effect of single crystal Fe 4.83 Mn 0.16 Si 2.91 P 0.09 (FMSP) are investigated experimentally and theoretically. A first‐order magneto‐elastic transition is observed at 348 K, where magnetic onset is accompanied by a change in lattice volume without an apparent change in crystal symmetry. Although the trace Mn and P doping are found to decrease the T C , the maximum magnetic entropy change Δ S Max ( T ) and the relative cooling power (RCP) of FMSP are enhanced compared to polycrystalline Fe 5 Si 3 . As a result, an intrinsically broader entropy change over a larger temperature span is generated in the lightly doped single crystal of Fe 5 Si 3 . The magnetic moment of the system is also enhanced. Density functional theory (DFT) calculations are performed to gain microscopic insights into the experimental findings. The results suggest that the hexagonal Fe 5 Si 3 is a new giant room temperature MCE material that is on par with La–Fe–Si and Fe‐Mn‐P‐Si systems.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Soft and transferable pseudopotentials from multi-objective optimization

Ab initio pseudopotentials are a linchpin of modern molecular and condensed matter electronic structure calculations. In this work, we employ multi-objective optimization to maximize pseudopotential softness while maintaining high accuracy and transferability. To accomplish this, we develop a formulation in which softness and accuracy are simultaneously maximized, with accuracy determined by the ability to reproduce all-electron energy differences between Bravais lattice structures, whereupon the resulting Pareto frontier is scanned for the softest pseudopotential that provides the desired accuracy in established transferability tests. We employ an evolutionary algorithm to solve the multi-objective optimization problem and apply it to generate a comprehensive table of optimized norm-conserving Vanderbilt (ONCV) pseudopotentials (https://github.com/SPARC-X/SPMS-psps). Here, we show that the resulting table is softer than existing tables of comparable accuracy, while more accurate than tables of comparable softness. The potentials thus afford the possibility to speed up calculations in a broad range of applications areas while maintaining high accuracy.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

High entropy alloys as catalysts: A focused review

A brief literature review of recent experimental and computational efforts on the use of high entropy alloys (HEAs) as catalysts is presented, while also sharing some perspectives and future insights. To fully broach the vast compositional possibilities of HEA materials, integrating computational modeling with high-throughput experimental synthesis and validation is necessary to accelerate their design and development. Once identified, specific HEAs can be a class of materials for the next generation of catalysts when addressing global challenges related to energy independence, commodity chemical production, environmental remediation, abating emissions, and modernized domestic supply chain resilience.

42 - ENGINEERING↗
DOE OSTI2023

Exploring Ligand-to-Metal Charge-Transfer States in the Photo-Ferrioxalate System Using Excited-State Specific Optimization

The photo-ferrioxalate system (PFS), [Fe(III)(C 2 O 4 )] 3– , more than an exact chemical actinometer, has been extensively applied in wastewater and environment treatment. Despite many experimental efforts to improve clarity, important aspects of the mechanism of ferrioxalate photolysis are still under debate. In this paper, we employ the recently developed WΓ-CASSCF to investigate the ligand-to-metal charge-transfer states that are key to ferrioxalate photolysis. This investigation provides a qualitative picture of these states and key potential energy surface features related to the photolysis. Furthermore, our theoretical results are consistent with the prompt charge-transfer picture seen in recent experiments and clarify some features that are not visible in experiments. Two ligand-to-metal charge-transfer states contribute to the photolysis of ferrioxalate, and the avoided crossing barrier between them is low compared with the initial photoexcitation energy. Our data also clarify that one Fe–O bond cleaves first, followed by the C–C bond and the other Fe–O bond.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Classification of Adsorbed Hydrocarbons Based on Bonding Configurations of the Adsorbates and Surface Site Stabilities

The design of heterogeneous catalysts can be accelerated by identifying relevant descriptors that accurately and effectively link the binding and activation energies to reactivity. Herein, we investigated scaling relations between binding energies of various hydrocarbon-based adsorbates on three different Pt surfaces and metal binding energies estimated via the recently developed α-scheme model. In this study, we find that the scaling slopes are similar for certain groups of adsorbates, which then can be classified based on their spatial and electronic structure enabling fast description of binding strengths for each member of the class. Hence, our findings show that the binding energies of simple hydrocarbons CH x , x = {0,1,2,3,4}, and CHCH 2 can be used to identify the binding energies of more complex hydrocarbon-based adsorbates. We introduce this classification to establish a generalizable scheme in which complex hydrogenation/dehydrogenation processes of higher hydrocarbons can be predicted via the binding energies of simpler hydrocarbon-based species and ultimately through surface site stabilities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Gamow-Teller Strength in Ca 48 and Ni 78 with the Charge-Exchange Subtracted Second Random-Phase Approximation

Here, we develop a fully self-consistent subtracted second random-phase approximation for charge-exchange processes with Skyrme energy-density functionals. As a first application, we study Gamow-Teller excitations in the doubly magic nucleus 48 Ca, the lightest double-β emitter that could be used in an experiment, and in 78 Ni, the single-beta-decay rate of which is known. The amount of Gamow-Teller strength below 20 or 30 MeV is considerably smaller than in other energy-density-functional calculations and agrees better with experiment in 48 Ca, as does the beta-decay rate in 78 Ni. These important results, obtained without ad hoc quenching factors, are due to the presence of two-particle–two-hole configurations. Their density progressively increases with excitation energy, leading to a long high-energy tail in the spectrum, a fact that may have implications for the computation of nuclear matrix elements for neutrinoless double-β decay in the same framework.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Accelerating Quantum Materials Development with Advances in Transmission Electron Microscopy

Quantum materials are driving a technology revolution in sensing, communication, and computing, while simultaneously testing many core theories of the past century. Materials such as topological insulators, complex oxides, superconductors, quantum dots, color center-hosting semiconductors, and other types of strongly correlated materials can exhibit exotic properties such as edge conductivity, multiferroicity, magnetoresistance, superconductivity, single photon emission, and optical-spin locking. These emergent properties arise and depend strongly on the material’s detailed atomic-scale structure, including atomic defects, dopants, and lattice stacking. In this review, we describe how progress in the field of electron microscopy (EM), including in situ and in operando EM, can accelerate advances in quantum materials and quantum excitations. We begin by describing fundamental EM principles and operation modes. We then discuss various EM methods such as (i) EM spectroscopies, including electron energy loss spectroscopy (EELS), cathodoluminescence (CL), and electron energy gain spectroscopy (EEGS); (ii) four-dimensional scanning transmission electron microscopy (4D-STEM); (iii) dynamic and ultrafast EM (UEM); (iv) complementary ultrafast spectroscopies (UED, XFEL); and (v) atomic electron tomography (AET). We describe how these methods could inform structure–function relations in quantum materials down to the picometer scale and femtosecond time resolution, and how they enable precision positioning of atomic defects and high-resolution manipulation of quantum materials. For each method, we also describe existing limitations to solve open quantum mechanical questions, and how they might be addressed to accelerate progress. Among numerous notable results, our review highlights how EM is enabling identification of the 3D structure of quantum defects; measuring reversible and metastable dynamics of quantum excitations; mapping exciton states and single photon emission; measuring nanoscale thermal transport and coupled excitation dynamics; and measuring the internal electric field and charge density distribution of quantum heterointerfaces- all at the quantum materials’ intrinsic atomic and near atomic-length scale. Finally, we conclude by describing open challenges for the future, including achieving stable sample holders for ultralow temperature (below 10K) atomic-scale spatial resolution, stable spectrometers that enable meV energy resolution, and high-resolution, dynamic mapping of magnetic and spin fields. With atomic manipulation and ultrafast characterization enabled by EM, quantum materials will be poised to integrate into many of the sustainable and energy-efficient technologies needed for the 21st century.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Towards a self-consistent analysis of the anisotropic galaxy two- and three-point correlation functions on large scales: application to mock galaxy catalogues

We establish a practical method for the joint analysis of anisotropic galaxy two- and three-point correlation functions (2PCF and 3PCF, respectively) on the basis of the decomposition formalism of the 3PCF using tripolar spherical harmonics. We perform such an analysis with MultiDark-Patchy mock catalogues to demonstrate and understand the benefit of the anisotropic 3PCF. We focus on scales above $80\, h^{-1}\, {\rm Mpc}$, and use information from the shape and the baryon acoustic oscillation (BAO) signals of the 2PCF and 3PCF. We also apply density field reconstruction to increase the signal-to-noise ratio of BAO in the 2PCF measurement, but not in the 3PCF measurement. In particular, we study in detail the constraints on the angular diameter distance and the Hubble parameter. We build a model of the bispectrum or 3PCF that includes the non-linear damping of the BAO signal in redshift space. We carefully account for various uncertainties in our analysis including theoretical models of the 3PCF, window function corrections, biases in estimated parameters from the fiducial values, the number of mock realizations to estimate the covariance matrix, and bin size. The joint analysis of the 2PCF and 3PCF monopole and quadrupole components shows a $30{{\ \rm per\ cent}}$ and $20{{\ \rm per\ cent}}$ improvement in Hubble parameter constraints before and after reconstruction of the 2PCF measurements, respectively, compared to the 2PCF analysis alone. Furthermore, this study clearly shows that the anisotropic 3PCF increases cosmological information from galaxy surveys and encourages further development of the modelling of the 3PCF on smaller scales than we consider.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Spin-state gaps and self-interaction-corrected density functional approximations: Octahedral Fe(II) complexes as case study

Accurate prediction of a spin-state energy difference is crucial for understanding the spin crossover phenomena and is very challenging for density functional approximations, especially for local and semi-local approximations due to delocalization errors. Here, we investigate the effect of the self-interaction error removal from the local spin density approximation (LSDA) and Perdew–Burke–Ernzerhof generalized gradient approximation on the spin-state gaps of Fe(II) complexes with various ligands using recently developed locally scaled self-interaction correction (LSIC) by Zope et al. [J. Chem. Phys. 151, 214108 (2019)]. The LSIC method is exact for one-electron density, recovers the uniform electron gas limit of the underlying functional, and approaches the well-known Perdew–Zunger self-interaction correction (PZSIC) as a particular case when the scaling factor is set to unity. Our results, when compared with reference diffusion Monte Carlo results, show that the PZSIC method significantly overestimates spin-state gaps favoring low spin states for all ligands and does not improve upon density functional approximations. The perturbative LSIC-LSDA using PZSIC densities significantly improves the gaps with a mean absolute error of 0.51 eV but slightly overcorrects for the stronger CO ligands. Finally, the quasi-self-consistent LSIC-LSDA, such as coupled-cluster single double and perturbative triple [CCSD(T)], gives a correct sign of spin-state gaps for all ligands with a mean absolute error of 0.56 eV, comparable to that of CCSD(T) (0.49 eV).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Percolative phase transition in few-layered MoSe 2 field-effect transistors using Co and Cr contacts

The metal-to-insulator phase transition (MIT) in two-dimensional (2D) materials under the influence of a gating electric field has revealed interesting electronic behavior and the need for a deeper fundamental understanding of electron transport processes, while attracting much interest in the development of next-generation electronic and optoelectronic devices. Although the mechanism of the MIT in 2D semiconductors is a topic under debate in condensed matter physics, our work demonstrates the tunable percolative phase transition in few-layered MoSe 2 field-effect transistors (FETs) using different metallic contact materials. Here, we attempted to understand the MIT through temperature-dependent electronic transport measurements by tuning the carrier density in a MoSe 2 channel under the influence of an applied gate voltage. In particular, we have examined this phenomenon using the conventional chromium (Cr) and ferromagnetic cobalt (Co) as two metal contacts. For both Cr and Co, our devices demonstrated n-type behavior with a room-temperature field-effect mobility of 16 cm 2 V −1 s −1 for the device with Cr-contacts and 92 cm 2 V −1 s −1 for the device with Co-contacts, respectively. Further, with low temperature measurements at 50 K, the mobilities increased significantly to 65 cm 2 V −1 s −1 for the device with Cr and 394 cm 2 V −1 s −1 for the device with Co-contacts. By fitting our experimental data to the percolative phase transition theory, the temperature-dependent conductivity data show a transition from an insulating-to-metallic behavior at a bias of ∼28 V for Cr-contacts and ∼20 V for Co-contacts. This cross-over of the conductivity can be attributed to an increase in carrier density as a function of the gate bias in temperature-dependent transfer characteristics. By extracting the critical exponents, we find that the transport behavior in the device with Co-contacts aligns closely with the 2D percolation theory. In contrast, the devices with Cr-contacts deviate significantly from the 2D limit at low temperatures.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Toward a First-Principles Framework for Predicting Collective Properties of Electrolytes

Conspectus Given the universal importance of electrolyte solutions, it is natural to expect that we have a nearly complete understanding of the fundamental properties of these solutions (e.g., the chemical potential) and that we can therefore explain, predict, and control the phenomena occurring in them. In fact, reality falls short of these expectations. But, recent advances in the simulation and modeling of electrolyte solutions indicate that it should soon be possible to make progress toward these goals. In this Account, we will discuss the use of first-principles interaction potentials based in quantum mechanics (QM) to enhance our understanding of electrolyte solutions. Specifically, we will focus on the use of quantum density functional theory (DFT) combined with molecular dynamics simulation (DFT-MD) as the foundation for our approach. The overarching concept is to understand and accurately reproduce the balance between local or short-ranged (SR) structural details and long-range (LR) correlations, allowing the prediction of the thermodynamics of both single ions in solution as well as the collective interactions characterized by activity/osmotic coefficients. In doing so, relevant collective motions and driving forces characterized by chemical potentials can be determined. Here, in this Account, we will make the case that understanding electrolyte solutions requires a faithful QM representation of the SR nature of the ion–ion, ion–water, and water–water interactions. However, the number of molecules that is required for collective behavior makes the direct application of high-level QM methods that contain the best SR physics untenable, making methods that balance accuracy and efficiency a practical goal. Alternatives such as continuum solvent models (CSMs) and empirically based classical molecular dynamics have been extensively employed to resolve this problem but without yet overcoming the fundamental issue of SR accuracy. We will demonstrate that accurately describing the SR interaction is imperative for predicting both intrinsic properties, namely, at infinite dilution, and collective properties of electrolyte solutions. DFT has played an important role in our understanding of condensed phase systems, e.g., bulk liquid water, the air–water interface, ions in bulk, and at the air–water interface. This approach holds huge promise to provide benchmark calculations of electrolyte solution properties that will allow for the development and improvement of more efficient methods, as well as an enhanced understanding of fundamental phenomena. However, the standard protocol using the generalized gradient approximation with van der Waals (vdW) correction requires improvement in order to achieve a high level of quantitative accuracy. Simply simulating with higher level DFT functionals may not be the best route considering the significant computational cost. Alternative methods of incorporating information from higher levels of QM should be explored; e.g., using force matching techniques on small clusters, where high level benchmark calculations are possible, to develop ideal correction terms to the DFT functional is a promising possibility. We argue that DFT with statistical mechanics is becoming an increasingly useful framework enabling the prediction of collective electrolyte properties.

Duignan, Timothy T.↗
DOE OSTI2022

Study of self-interaction-errors in barrier heights using locally scaled and Perdew-Zunger self-interaction methods

We study the effect of self-interaction errors on the barrier heights of chemical reactions. For this purpose we use the well-known Perdew-Zunger [J. P. Perdew and A. Zunger, Phys. Rev. B,23, 5048 (1981)] self-interaction-correction(PZSIC), as well two variations of the recently developed, locally scaled self-interaction correction (LSIC) [Zope et al.,J. Chem. Phys.151, 214108 (2019)] to study the barrier heights of the BH76 benchmark dataset. Our results show that both PZSIC and especially the LSIC methods improve the barrier heights relative to the local density approximation(LDA). The version of LSIC that uses the iso-orbital indicator z as a scaling factor gives a more consistent improvement than an alternative version that uses an orbital-dependent factor based on the ratio of orbital densities to the total electron density. We show that LDA energies evaluated using the self-consistent and self-interaction-free PZSIC densities can be used to assess density-driven errors. Furthermore, the LDA reaction barrier errors for the BH76 set are found to contain significant density-driven errors for all types of reactions contained in the set, but the corrections due to adding SIC to the functional are much larger than those stemming from the density for the hydrogen transfer reactions and of roughy equal size for the non-hydrogen transfer reactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Multi-objective goal-directed optimization of de novo stable organic radicals for aqueous redox flow batteries

Abstract Advances in the field of goal-directed molecular optimization offer the promise of finding feasible candidates for even the most challenging molecular design applications. One example of a fundamental design challenge is the search for novel stable radical scaffolds for an aqueous redox flow battery that simultaneously satisfy redox requirements at the anode and cathode, as relatively few stable organic radicals are known to exist. To meet this challenge, we develop a new open-source molecular optimization framework based on AlphaZero coupled with a fast, machine-learning-derived surrogate objective trained with nearly 100,000 quantum chemistry simulations. The objective function comprises two graph neural networks: one that predicts adiabatic oxidation and reduction potentials and a second that predicts electron density and local three-dimensional environment, previously shown to be correlated with radical persistence and stability. With no hard-coded knowledge of organic chemistry, the reinforcement learning agent finds molecule candidates that satisfy a precise combination of redox, stability and synthesizability requirements defined at the quantum chemistry level, many of which have reasonable predicted retrosynthetic pathways. The optimized molecules show that alternative stable radical scaffolds may offer a unique profile of stability and redox potentials to enable low-cost symmetric aqueous redox flow batteries.

25 ENERGY STORAGE↗
DOE OSTI2021

Application of the locally self-consistent embedding approach to the Anderson model with non-uniform random distributions

Highlights: • Typical Medium Theory (TMT) for the Anderson Localization. • Locally Self-Consistent Multiple Scattering Method (LSMS) for Random Disordered Systems. • Linear Scaling Computational Method for Random Disordered Systems. We apply the recently developed embedding scheme for the locally self-consistent method to random disorder electrons systems. The method is based on the locally self-consistent multiple scattering theory and the typical medium theory. The locally self-consistent multiple scattering theory divides a system into many small designated local interaction zones. The subsystem within each local interaction zone is embedded in a self-consistent field from the typical medium theory. This approximation allows the study of random systems with large numbers of sites. We present results for the three dimensional Anderson model with different random disorder potential distributions. Using the typical density of states as an indicator of Anderson localization, we find that the method can capture the localization for commonly studied disorder potentials. These include the uniform distribution, the Gaussian distribution, and even the unbounded Cauchy distribution.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗