Search NASA⌕ Search

SEARCH · Search NASA

Results for “molecular cluster structure”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Molecular Examination of Ion-Pair Competition in Alkaline Aluminate Solutions Using In Situ Liquid SIMS

Understanding the structure and composition of aluminate complexes in extremely alkaline systems such as Bayer liquors has received enormous attentions due to its fundamental and industrial importance. However, obtaining direct molecular information of the underlying ion-ion interactions using traditional approaches such as NMR or Raman spectroscopy is challenging due to the weakness of these interactions and/or their complex overlapping spectral signatures. Here, we exploit in situ liquid secondary ion mass spectrometry (SIMS) as a new approach and show how it enables new insights. In contrast with traditional techniques, using SIMS we succeeded in acquiring information on dominant ion clusters in these alkaline systems. In Na + /K + mixed alkaline aluminate solutions we clearly observe preferential formation of Na + -anion clusters over K + -anion clusters. Further, evaluation of these clusters by density functional theory (DFT) calculations shows that these structures are stable and that their relative bond energies are consistent with their observed SIMS signal intensity differences. This demonstrates a key advantage of in situ liquid SIMS for overcoming ambiguities obscuring important information in these systems on constituent molecular clusters defined by relatively weak ion pair competition and ion-solvent interactions.

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

Materials Data on ScH16I3O8 by Materials Project

ScH10O7IH2OHI1IH3 crystallizes in the monoclinic P2_1/c space group. The structure is zero-dimensional and consists of four hydriodic acid molecules, four molecular hydrogen;hydroiodide molecules, four water molecules, and four ScH10O7I clusters. In each ScH10O7I cluster, Sc3+ is bonded in a 7-coordinate geometry to seven O2- atoms. There are a spread of Sc–O bond distances ranging from 1.97–2.35 Å. There are ten inequivalent H1+ sites. In the first H1+ site, H1+ is bonded in a single-bond geometry to one O2- atom. The H–O bond length is 0.98 Å. In the second H1+ site, H1+ is bonded in a single-bond geometry to one O2- atom. The H–O bond length is 0.98 Å. In the third H1+ site, H1+ is bonded in a single-bond geometry to one O2- atom. The H–O bond length is 0.99 Å. In the fourth H1+ site, H1+ is bonded in a single-bond geometry to one O2- atom. The H–O bond length is 0.99 Å. In the fifth H1+ site, H1+ is bonded in a single-bond geometry to one O2- atom. The H–O bond length is 0.98 Å. In the sixth H1+ site, H1+ is bonded in a single-bond geometry to one O2- atom. The H–O bond length is 0.99 Å. In the seventh H1+ site, H1+ is bonded in a single-bond geometry to one O2- and one I1- atom. The H–O bond length is 1.00 Å. The H–I bond length is 2.42 Å. In the eighth H1+ site, H1+ is bonded in a single-bond geometry to one O2- atom. The H–O bond length is 0.99 Å. In the ninth H1+ site, H1+ is bonded in a single-bond geometry to one O2- atom. The H–O bond length is 0.97 Å. In the tenth H1+ site, H1+ is bonded in a single-bond geometry to one O2- atom. The H–O bond length is 0.98 Å. There are seven inequivalent O2- sites. In the first O2- site, O2- is bonded in a distorted water-like geometry to one Sc3+ and two H1+ atoms. In the second O2- site, O2- is bonded in a single-bond geometry to one Sc3+, one H1+, and one O2- atom. The O–O bond length is 1.48 Å. In the third O2- site, O2- is bonded in a distorted water-like geometry to one Sc3+ and two H1+ atoms. In the fourth O2- site, O2- is bonded in a distorted water-like geometry to one Sc3+ and two H1+ atoms. In the fifth O2- site, O2- is bonded in a 1-coordinate geometry to one Sc3+ and one O2- atom. In the sixth O2- site, O2- is bonded in a distorted water-like geometry to one Sc3+ and one H1+ atom. In the seventh O2- site, O2- is bonded in a distorted water-like geometry to one Sc3+ and two H1+ atoms. I1- is bonded in a 1-coordinate geometry to one H1+ atom.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Materials Data on MgH20C4Br2(N4O3)2 by Materials Project

MgC4H4(N2O3)2(NH2)4(H2)2(H2Br)2 crystallizes in the monoclinic P2_1/c space group. The structure is zero-dimensional and consists of eight ammonia molecules, four hydrogen molecules, four molecular hydrogen;hydrobromide molecules, and two MgC4H4(N2O3)2 clusters. In each MgC4H4(N2O3)2 cluster, Mg2+ is bonded in an octahedral geometry to six O2- atoms. There are a spread of Mg–O bond distances ranging from 1.97–2.23 Å. There are two inequivalent C4+ sites. In the first C4+ site, C4+ is bonded in a 1-coordinate geometry to two N3- and one O2- atom. Both C–N bond lengths are 1.33 Å. The C–O bond length is 1.26 Å. In the second C4+ site, C4+ is bonded in a single-bond geometry to one O2- atom. The C–O bond length is 1.17 Å. There are two inequivalent N3- sites. In the first N3- site, N3- is bonded in a distorted single-bond geometry to one C4+ atom. In the second N3- site, N3- is bonded in a distorted single-bond geometry to one C4+ atom. There are two inequivalent H1+ sites. In the first H1+ site, H1+ is bonded in a single-bond geometry to one O2- atom. The H–O bond length is 0.98 Å. In the second H1+ site, H1+ is bonded in a single-bond geometry to one O2- atom. The H–O bond length is 0.98 Å. There are three inequivalent O2- sites. In the first O2- site, O2- is bonded in a bent 120 degrees geometry to one Mg2+ and one C4+ atom. In the second O2- site, O2- is bonded in a bent 150 degrees geometry to one Mg2+ and one C4+ atom. In the third O2- site, O2- is bonded in a distorted trigonal planar geometry to one Mg2+ and two H1+ atoms.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Catalyst design for small molecule activation of energy consequence

This project targets the conversion of ubiquitous small molecules (e.g. NO, CO, H 2 O) into viable precursors to synthetic fuels. Current state of the art catalyst design has not directly targeted transition metal complexes capable of mediating the multi-electron redox processes necessary to reduce the overpotential (energy loss) required achieve efficient activation of small molecule substrates. In this vein, a new strategy has been developed for the assembly of polynuclear architectures; allowing for the construction of tunable polymetallic centers that assemble easily within a pre-organized template (conferring stability, selectivity and tunability) that can effect multi-electron redox processes for reactions. Catalyst development has commenced with the following target design elements: (1) catalysts featuring multiple transition metal ions in the same reaction space to greatly expand accessible molecular redox capabilities; (2) catalysts are assembled in a polynucleating ligand framework that permits control over the cluster morphology as well as the local steric and electronic environment of the transition metal ions within the cluster. The high-tunability of the catalyst composition (metal content) and geometric flexibility has permitted a rigorous assessment of electronic-structure-to-function relationship to be developed, further guiding synthetic efforts to realize more potent catalysts. The numerous permutations possible showcase the high degree of generality to this approach with many synthetic handles to tune redox and reaction chemistry. Trinuclear complexes have been synthesized featuring homo- and hetero-trinuclear cores featuring a variety of first row transition metal ions (Cr→Ni). The molecular clusters have been shown to successfully mediate multi-electron redox processes in a cooperative fashion without requiring strong chemical reductants or oxidants. The reactive molecular complexes are being utilized to activate and breakdown the robust bonds within typical waste stream small molecules (e.g., greenhouse gases) and convert them into value-added commodity chemicals. Ultimately, the catalysts developed by this approach will be required to convert energy acquired via renewable resources (e.g., solar or wind) into synthetic fuels as an energy storage mechanism.

10 SYNTHETIC FUELS↗
DOE OSTI2021

Combining artificial intelligence and physics-based modeling to directly assess atomic site stabilities: from sub-nanometer clusters to extended surfaces

The performance of functional materials is dictated by chemical and structural properties of individual atomic sites. In catalysts, for instance, the thermodynamic stability of constituting atomic sites is a key descriptor from which more complex properties, such as molecular adsorption energies and reaction rates, can be derived. In this study, we present a widely applicable machine learning (ML) approach to instantaneously compute the stability of individual atomic sites in structurally and electronically complex nano-materials. Conventionally, we determine such site stabilities using computationally intensive first-principles calculations. With our approach, we predict the stability of atomic sites in sub-nanometer metal clusters of 3–55 atoms with mean absolute errors in the range of 0.11–0.14 eV. To extract physical insights from the ML model, we introduce a genetic algorithm (GA) for feature selection. This algorithm distills the key structural and chemical properties governing the stability of atomic sites in size-selected nanoparticles, allowing for physical interpretability of the models and revealing structure–property relationships. The results of the GA are generally model and materials specific. In the limit of large nanoparticles, the GA identifies features consistent with physics-based models for metal–metal interactions. By combining the ML model with the physics-based model, we predict atomic site stabilities in real time for structures ranging from sub-nanometer metal clusters (3–55 atom) to larger nanoparticles (147 to 309 atoms) to extended surfaces using a physically interpretable framework. Finally, we present a proof of principle showcasing how our approach can determine stable and active nanocatalysts across a generic materials space of structure and composition.

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

Microstructure‐Dependent Sodium Storage Mechanisms in Hard Carbon Anodes

Sustainable energy storage is essential to support the transition to renewables and meet the increasing demand for energy. Sodium‐ion batteries (NIBs) are attractive for grid‐scale energy storage due to the abundance and low cost of sodium, sustainability of other battery components, and electrochemical performance. Hard carbon (HC) is a leading anode material for NIBs, but its complex microstructure complicates the understanding of sodium storage mechanisms. Using X‐ray total scattering and density functional theory calculations, this study clarifies how HC's microstructural variations influence sodium storage across the slope (high potential) and plateau (low potential) regions of the potential capacity curve. In the slope region, sodium initially adsorbs at high‐binding energy defect sites and subsequently intercalates between graphene layers, adsorbing at low‐binding energy defect sites, correlating with different slopes observed during initial sodiation. Initial irreversibility arises from sodium trapping at surface defects and solid electrolyte interface formation. In the plateau region, sodium simultaneously intercalates and fills pores, influenced by pore size, interlayer spacing, and defect concentration. HCs with larger pore sizes form larger sodium clusters. In conclusion, the proposed mechanism underscores the role of microstructure engineering in enhancing HC performance and advancing NIBs for grid‐scale energy storage.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Fate of multiparticle resonances: From Q-balls to 3 He droplets

We study a system of N nonrelativistic particles which form a near-threshold resonance. Assuming no subset of these particles can form a bound state, the resonance can only decay through an “explosion” into N particles. We find that the decay width of the resonance scales as E Δ –5/2 in the limit when the energy E of the resonance goes to zero, where Δ is the ground-state energy of a system of N particles in a spherical harmonic trap with unit frequency. Here, the formula remains valid when some pairs of final particles have zero-energy s-wave resonance, but the Efimov effect is not present. In the limit of large N, we show that the final particles follow a Maxwell-Boltzmann distribution if they are bosons and a semicirclelike law if they are fermions. We expect our general result to be applicable to various systems that exist in nature. In particular, we argue that metastable 3 He droplets exist with the lifetime varying over many orders of magnitude ranging from a fraction of a nanosecond to values greatly exceeding the age of the Universe.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2026

Molecular Design Principles for Photosystem I-Based Biohybrid Solar Fuel Catalysts

Direct solar-to-chemical conversion offers a compelling route to clean, dispatchable energy. Photosystem I (PSI), an evolutionarily optimized light-driven oxidoreductase, can be repurposed for solar-fuel production by coupling its photochemistry to catalytic interfaces. However, the molecular determinants that govern productive electron transfer to abiotic catalysts remain poorly understood. Here, we present molecular structures of active PSI-Pt nanoparticle (PtNP) biohybrids that reveal how protein architecture controls catalyst access, binding geometry, and photocatalytic efficiency. Removal of stromal subunits exposes the electron transfer chain and enables PtNP binding proximal to the F X cluster, demonstrating that steric occlusion limits access to native acceptor regions in PSI. In contrast, in trimeric PSI, PtNPs bind at multiple sites per monomer, but only a subset are positioned within electron transfer distance of terminal cofactors, resulting in a heterogeneous population of productive and nonproductive configurations. Structural analyses and molecular dynamics simulations define the interface topology, electrostatics, and cofactor-to-nanoparticle distances that govern catalyst binding and electron transfer. These results establish that catalytic inefficiency arises not only from intrinsic electron transfer constraints but also from the distribution of binding geometries imposed by the protein scaffold. Together, these findings provide a molecular framework linking protein structure to biohybrid function and define design principles for engineering PSI-based solar fuel systems and protein-nanomaterial interfaces for light-driven catalysis.

biohybrid↗
DOE OSTI2024

Diffusion quantum Monte Carlo approach to the polaritonic ground state

Making and using polaritonic states (i.e., hybrid electron-photon states) for chemical applications has recently become one of the most prominent and active fields that connects the communities of chemistry and quantum optics. Modeling of such polaritonic phenomena using ab initio approaches calls for new methodologies, leading to the reinvention of many commonly used electronic structure methods, such as Hartree-Fock, density functional, and coupled cluster theories. Here in this work, we explore the formally exact diffusion quantum Monte Carlo approach to obtain numerical solutions to the polaritonic ground state during the dissociation of the H 2 molecular system. Here, we examine various electron-nuclear-photon properties throughout the dissociation, such as changes to the minimum of the cavity Born-Oppenheimer surface, the localization of the electronic wave function, and the average mode occupation. Finally, we directly compare our results to that obtained with state-of-the-art, yet approximate, polaritonic coupled cluster approaches.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Effects of interaction strength of associating groups on linear and star polymer dynamics

A small number of associating groups incorporated onto a polymer backbone have dramatic effects on the mobility and viscoelastic response of the macromolecules in melts. These associating groups assemble, driving the formation of clusters, whose lifetime affects the properties of the polymers. In this study, we probe the effects of the interaction strength on the structure and dynamics of two topologies, linear and star polymer melts, and further investigate blends of associative and non-associating polymers using molecular dynamics simulations. Polymer chains of approximately one entanglement length are described by a bead–spring model, and the associating groups are incorporated in the form of interacting beads with an interaction strength between them that is varied from 1 to 20 k B T. We find that, for all melts and blends, interaction of a few k B T between the associating groups drives cluster formation, where the size of the clusters increases with increasing interaction strength. These clusters act as physical crosslinkers, which slow the chain mobility. Blends of chains with and without associating groups macroscopically phase separate for interaction strength between the associating groups of a few k B T and above. For weakly interacting associating groups, the static structure function S(q) is well fit by functional form predicted by the random phase approximation where a clear deviation occurs as phase segregation takes place, providing a quantitative assessment of phase segregation.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Rare Earth and Transition Metal Containing Glasses

Transition metal (TM) and rare earth (RE) ions have been incorporated into many glass systems such as silicate, phosphate, and borosilicate-based oxide glasses, as well as in halide and chalcogenide glasses, that find applications ranging from optical, photonic, and magnetic devices, solid-state battery, to nuclear waste disposal. Understanding the structural role of RE and TM in these glasses can help to develop glass compositions for targeted applications with either high-optical emission efficiency, electrical conductivity, or chemical durability. In this chapter, we first provide a general introduction of the applications and structural features of RE and TM in glasses, then the critical aspects of molecular dynamics (MD) simulations of these glasses such as interatomic potentials, structural analysis tools to study RE and TM ions in glasses and their clustering behaviors, Quantitative Structure–Property Analysis (QSPR), diffusion and dynamic property calculations, and electronic structure calculations to understand electronic defects such as charge trapping and radiation effects are introduced. Three representative case studies are presented: the first one is on MD simulations of erbium- and europium-doped silica and silicate glasses, as well as cerium doped aluminophosphate glasses, that revealed the effect of glass composition on RE ion local structure and clustering behavior. Electronic structure calculations of cerium-doped glass show how the existence of multioxidation states help to mediate radiation-induced damages caused by excited electron–hole pairs was also discussed. The second one focuses on alkali vanadophosphate glasses where the existence of two vanadium oxidation states help to provide electronic conduction in the glasses while alkali ions provide ionic conduction. MD simulations were used to understand vanadium environments and other structural aspects in the phosphate glasses, as well as the ionic transport behaviors of alkali ions. The third case study is on zirconium-containing borosilicate and aluminosilicate glasses that find wide applications in nuclear waste disposal. MD simulations help to provide structural details of zirconium ions that are validated by diffraction and EXAFS spectra. The structural information was used to interpret changes of mechanical properties and chemical durability by using QSPR and other analyses-based MD-generated structure models.

Du, Jincheng↗
DOE OSTI2021

Critical Role of the Interfacial Layer in Associating Polymers with Microphase Separation

Polymers with dynamic (transient) bonds, often called associating polymers, have been attracting significant attention in recent years due to their unique viscoelastic properties, self-healing ability, and recyclability. Nevertheless, understanding the mechanisms and the factors controlling their macroscopic properties remains limited due to the higher complexity introduced by the dynamic bonds. Here in this study, small-angle X-ray scattering (SAXS), broadband dielectric spectroscopy (BDS), and rheology were applied to unravel the structure and dynamics of telechelic associating polymers with different molecular weights. SAXS measurements revealed phase separation of the functional end groups with an average cluster size of ~2–3 nm and the distance between clusters controlled by the chain length. Borrowing the interfacial layer model analysis of BDS data from the polymer nanocomposite field, we demonstrated the presence of an interfacial polymer layer with a thickness of ~0.7–0.9 nm surrounding these clusters. Rheological measurements showed quantitatively that the presence of the interfacial layer significantly alters the viscoelastic behavior of these materials, indicating the crucial role of the interfacial layer in defining the macroscopic mechanical properties of the studied telechelic materials. The presented results emphasize that phase separation of the functional groups in associating polymers leads to very significant changes of the viscoelastic properties, opening a promising avenue in the design of novel functional materials.

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

Determination of Molecular Structure and Dynamics of Molten Salts by Advanced Neutron and X-ray Scattering Measurements and Computer Modeling

The design and development of fully functional Molten-Salt Reactors (MSR) require detailed knowledge of the molten salt properties in order to understand and predict the salt’s behavior. Fundamental properties of interest include molecular structure, speciation, and dynamics (such as diffusion coefficients) of salt components and dissolved corrosion and fission products. Computer modeling is necessary to predict changes in physical and chemical properties due to irradiation, burning of dissolved fuel, and corrosion. The modeling requires experimental data, and advanced neutron and x-ray scattering and spectroscopy provide the most reliable and direct determination of the structure (Pair-Distribution Functions, PDF), and dynamics of ions in the melt. This project dealt with both fluoride and chloride salts. The PDFs have been measured by a combination of neutron and x-ray diffraction. We utilized the techniques of isotope substitutions, a very powerful tool available for neutron-scattering, to extract the details of the liquid structure. Although similar measurements have been done before, modern advanced neutron and x-ray-scattering techniques allow collecting the data at much higher resolution and in a wider range of temperatures. Importantly, we were among the first to study fluoride salts by neutron scattering. The importance of impurities and their effects on salt properties have become apparent recently and so new methods of salt purification were developed. We took advantage of these developments to produce reliable data, which have been used for computer simulations of both clean salts and those with added fission and corrosion products most relevant for MSRs. Ab initio molecular dynamics simulations have been performed to understand the multi-component liquid solution, in particular solubility of impurities and thermodynamic interactions in relation to the ionic-cluster structure of the fluid. We applied machine learning to regress from the simulation and experimental data in order to develop a fast-acting model that can handle molten salt with an arbitrary (≥ 10) number of chemical elements and be able to predict chemical potential as a function of composition and temperature. This project resulted in a number of experimental and computer-simulation publications, a patent application, and numerous conference presentations (American Physical Society, American Chemical Society, and The Electrochemical Society among others). Multiple students and postdocs participated and collaborated on aspects of this project. This project seeded new collaborations between MIT and other institutions, such as the University of Massachusetts Lowell, the University of Illinois Urbana-Champaign, the University of California Berkeley, and Oak Ridge and Los Alamos National Labs. As such, this project has had a broad and lasting impact beyond its original scientific scope.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Structural and electronic characteristics of amorphous Ge 8 Sb 2 Te 11

GeTe-rich phase-change materials can be utilized in rewriteable optical memory due to the large contrast in reflectivity between amorphous and crystalline phases. Here we explored the structure and electronic properties of amorphous Ge 8 Sb 2 Te 11 using ab initio molecular dynamics simulations. The results indicate that amorphous Ge 8 Sb 2 Te 11 is dominantly composed of Ge-, Sb- and Te-centered octahedrons with distortions, while 30.4 % of Ge-centered clusters are in the form of tetrahedrons which are randomly distributed. The 5-fold rings possess a large proportion, and the Gesingle bondTe and Sbsingle bondTe bonds present larger formation energies than other bonds, leading to the ABAB bonding sequence (A: Ge and Sb, B: Te). The lone pair electrons locating at the opposite direction of bonds possess a large fraction of 14.8 %, which may enhance the distortions of local clusters. Importantly, these structural properties lead to the high stability of amorphous Ge 8 Sb 2 Te 11 and thus long data retention in the information storage.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Molecular dynamics simulation of Fe-Si alloys using a neural network machine learning potential

Interatomic potential development using machine learning (ML) approaches has attracted a lot of attention in recent years because these potentials can effectively describe the structural and dynamical properties of complex materials at the atomistic level. Here, in this work, we present the development of a neural network (NN) deep ML interatomic potential for Fe-Si alloys, and we demonstrate the effectiveness of the NN-ML potential in predicting the structures and energies of liquid and crystalline phases of Fe-Si alloys in comparison with the results from ab initio molecular dynamics simulations or experimental data. The developed NN-ML potential is also used to perform molecular dynamics simulations to study the structures of Fe-Si alloys with various compositions under rapid solidification conditions. The short-ranged orders in the rapidly solidified Fe-Si alloys are also analyzed by a cluster alignment method.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Determining Partial Atomic Charges for Liquid Water: Assessing Electronic Structure and Charge Models

Partial atomic charges provide an intuitive and efficient way to describe the charge distribution and the resulting intermolecular electrostatic interactions in liquid water. Many charge models exist and it is unclear which model provides the best assignment of partial atomic charges in response to the local molecular environment. In this work, we systematically scrutinize various electronic structure methods and charge models (Mulliken, natural population analysis, CHelpG, RESP, Hirshfeld, Iterative Hirshfeld, and Bader) by evaluating their performance in predicting the dipole moments of isolated water, water clusters, and liquid water as well as charge transfer in the water dimer and liquid water. Although none of the seven charge models is capable of fully capturing the dipole moment increase from isolated water (1.85 D) to liquid water (about 2.9 D), the Iterative Hirshfeld method performs best for liquid water, reproducing its experimental average molecular dipole moment, yielding a reasonable amount of intermolecular charge transfer, and showing modest sensitivity to the local water environment. The performance of the charge model is dependent on the choice of the density functional and the quantum treatment of the environment. The computed molecular dipole moment of water generally increases with the percentage of the exact Hartree–Fock exchange in the functional, whereas the amount of charge transfer between molecules decreases. For liquid water, including two full solvation shells of surrounding water molecules (within about 5.5 Å of the central water) in the quantum chemical calculation converges the charges of the central water molecule. Furthermore, our final pragmatic quantum chemical charge-assigning protocol for liquid water is the Iterative Hirshfeld method with M06-HF/aug-cc-pVDZ and a quantum region cutoff radius of 5.5 Å.

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

Fermi-Löwdin-orbital self-interaction correction using the optimized-effective-potential method within the Krieger-Li-Iafrate approximation

Perdew-Zunger self-interaction correction (PZ-SIC) offers a route to remove self interaction errors on an orbital-by-orbital basis. A recent formulation of PZ-SIC by Pederson, Ruzsinszky and Perdew proposes restricting the unitary transformation to localized orbitals called Fermi-Löwdin orbitals. This formulation, called the FLOSIC method, simplifies PZ-SIC calculations and was implemented self-consistently using a Jacobi-like (FLOSIC-Jacobi) iteration scheme. In this work we implement the FLOSIC approach using the Krieger-Li-Iafrate (KLI) approximation to the optimized effective potential (OEP). We compare the results of present FLOSIC-KLI approach with FLOSIC-Jacobi scheme for atomic energies, atomization energies, ionization energies, barrier heights, polarizability of chains of hydrogen molecules etc. to validate the FLOSIC-KLI approach. The FLOSICKLI approach, which is within the realm of Kohn-Sham theory, predicts smaller energy gaps between frontier orbitals due to the lowering of eigenvalues of the lowest unoccupied orbitals. Results show that atomic energies, atomization energies, ionization energy as an absolute of highest occupied orbital eigenvalue, and polarizability of chains of hydrogen molecules between the two methods agree within 2%. Lastly the FLOSIC-KLI approach is used to determine the vertical ionization energies of water clusters.

74 ATOMIC AND MOLECULAR PHYSICS↗