Search NASA⌕ Search

SEARCH · Search NASA

Results for “optimized energy density functionals”

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 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Supramolecular Assembly of One-Dimensional Coordination Polymers for Efficient Separation of Xenon and Krypton

Efficient separation and purification of xenon (Xe) from krypton (Kr) represent an industrially crucial but challenging process. While the adsorption-based separation of these atomic gases represents an energy-efficient process, achieving highly selective adsorbents remains a difficult task. Here, we demonstrate a supramolecular assembly of coordination polymers, termed as M(II)-dhbq (M = Mg, Mn, Co, and Zn; dhbq = 2,5-dihydroxy-1,4-benzoquinone), with high-density open metal sites (5.3 nm –3 ) and optimal pore size (5.5 Å), which are able to selectively capture Xe among other chemically inert gases including Kr, Ar, N 2 , and O 2 . Among M(II)-dhbq materials, Mn-dhbq exhibits the highest Xe uptake capacity of 3.1 mmol/g and a Xe/Kr selectivity of 11.2 at 298 K and 1.0 bar, outperforming many state-of-the-art adsorbents reported so far. Remarkably, the adsorption selectivity of Mn-dhbq for Xe/O 2 , Xe/N 2 , and Xe/Ar at ambient conditions reaches as high as 70.0, 139.3, and 64.0, respectively. Direct breakthrough experiments further confirm that all M(II)-dhbq materials can efficiently discriminate Xe atoms from other inert gases. Furthermore, it is revealed from the density functional theory calculations that the strong affinity between Xe and the coordination polymer is mainly attributed to the polarization by open metal sites.

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

Final Report: Hot Carrier Collection in Thin Film Silicon with Tailored Nanocrystalline/Amorphous Structure

This project developed, characterized, and perfected a new type of highly tunable nanocrystalline silicon (nc-Si:H) incorporating quantum confined silicon nanoparticles (SiNPs). A dual zone deposition process and system were developed and demonstrated. The depositions of SiNPs, the amorphous phase, and co-deposited material were characterized and optimized. Material design and interpretation of results were guided by new theoretical tools that examined both the electronic structure and carrier dynamics of this hybrid material. Heterojunction and p-i-n solar cells were demonstrated and characterized. Photo-thin-film-transistors allowed mobility to be studied as a function SiNP density in the films. Rapid (hot) transfer of carriers from the amorphous matrix to the quantum confined SiNPs was observed and connected to reduced photo-degradation. The results carry quantum confined Si dots from a novelty to materials that can be harnessed for PV and optoelectronic applications. The growth process is broadly extendable with alternative amorphous matrices, novel layered structures, and alternative NPs easily accessible. The hot carrier effects hold the potential for third generation photovoltaics.

14 SOLAR ENERGY↗
DOE OSTI2017

Development of a Multiphase Adjoint Capability in OpenFOAM

Thermal hydraulic predictions of flow regimes within nuclear reactor cores require significant computational resources and accuracy to ensure that the core design does not violate the thermal limits of materials. Multiphase computational fluid dynamics (CFD) has become a valuable computational capability for predicting flow characteristics as they relate to applications such as nuclear reactor safety and mixing phenomena around spacer grids within fuel assemblies. Multiphase CFD requires closure relationships for both the effects of turbulence as well as the transfer of mass, momentum, and energy between phases. Understanding the sensitivity of the parameters used in closure relationships is important for improving accuracy, but performing rigorous sensitivity analysis is computationally expensive. Adjoint mathematics have been utilized for reactor physics and fluid dynamics applications with respect to verification, sensitivity analysis, adaptive mesh refinement, and optimization. The usage of adjoint methods to predict the sensitivities of a variety of closure relationship parameters within a given problem has been shown to return sensitivity responses with significantly fewer computations than normal 'large-scale' methods. As the complexity of multiphase CFD methods increases, the demand for accurate verification tools and improvements in computational efficiency also increases, making the development of adjoint capabilities an attractive option. Closure relationships that describe interphase momentum physics in multiphase CFD modelling rely heavily on prescribed constants determined from experimental data. Due to the computational efficiency of an adjoint approach, the time to derive these closure relationships could be dramatically reduced, and data assimilation techniques could be used to return probability density functions of closure models derived using experimental data. This work develops a multiphase adjoint capability for the Eulerian-Eulerian two field, two phase equations. Specifically, the adjoint capability has been developed for the sub-cooled nucleate boiling solver boilEulerFoam developed by Dr. Nam Dinh of NC State University within the open-source CFD library OpenFOAM.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
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

Dynamic Metal–Support Interaction Dictates Cu Nanoparticle Sintering on Al 2 O 3 Surfaces

Nanoparticle sintering remains a critical challenge in heterogeneous catalysis. In this work, we present a unified deep potential (DP) model based on the Perdew–Burke–Ernzerhof approximation of density functional theory for Cu nanoparticles on three Al 2 O 3 surfaces (γ-Al 2 O 3 (100), γ-Al 2 O 3 (110), and α-Al 2 O 3 (0001)). Using DP-accelerated simulations, we reveal that the nanoparticle size-mobility relationship strongly depends on the supporting surface. The diffusion of nanoparticles on the two γ-Al 2 O 3 surfaces is almost independent of the size of the nanoparticle, while the diffusion on α-Al 2 O 3 (0001) decreases rapidly with increasing size. Interestingly, nanoparticles with fewer than 55 atoms diffuse several times faster on α-Al 2 O 3 (0001) than on γ-Al 2 O 3 (100) at 800 K while expected to be more sluggish based on their larger binding energy at 0 K. The diffusion on α-Al 2 O 3 (0001) is facilitated by dynamic metal–support interaction (MSI), where Al atoms move out of the surface plane to optimize contact with the nanoparticle and relax back to the plane as the nanoparticle moves away. In contrast, the MSI on γ-Al 2 O 3 (100) and on γ-Al 2 O 3 (110) is dominated by more stable and directional Cu–O bonds, consistent with the limited diffusion observed on these surfaces. Our extended MD simulations provide insight into the sintering processes, showing that the dispersity of the nanoparticles strongly influences the coalescence driven by nanoparticle diffusion. We observed that the coalescence of Cu 13 nanoparticles on α-Al 2 O 3 (0001) can occur in a short time (10 ns) at 800 K even with an initial internanoparticle distance increased to 3 nm, while the coalescence on the two γ-Al 2 O 3 surfaces are inhibited significantly by increasing the initial internanoparticle distance. These findings demonstrate that the dynamics of the supporting surface is crucial to understanding the sintering mechanism and offer guidance for designing sinter-resistant catalysts by engineering the support morphology.

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

Measurement of the dynamic charge response of materials using low-energy, momentum-resolved electron energy-loss spectroscopy (M-EELS)

One of the most fundamental properties of an interacting electron system is its frequency- and wave-vector-dependent density response function, \chi({\bf q},\omega) . The imaginary part, \chi''({\bf q},\omega) , defines the fundamental bosonic charge excitations of the system, exhibiting peaks wherever collective modes are present. \chi χ quantifies the electronic compressibility of a material, its response to external fields, its ability to screen charge, and its tendency to form charge density waves. Unfortunately, there has never been a fully momentum-resolved means to measure \chi({\bf q},\omega) at the meV energy scale relevant to modern electronic materials. Here, we demonstrate a way to measure \chi χ with quantitative momentum resolution by applying alignment techniques from x-ray and neutron scattering to surface high-resolution electron energy-loss spectroscopy (HR-EELS). This approach, which we refer to here as “M-EELS”, allows direct measurement of \chi''({\bf q},\omega) with meV resolution while controlling the momentum with an accuracy better than a percent of a typical Brillouin zone. We apply this technique to finite-q excitations in the optimally-doped high temperature superconductor, Bi _2 2 Sr _2 2 CaCu _2 2 O _{8+x} 8 + x (Bi2212), which exhibits several phonons potentially relevant to dispersion anomalies observed in ARPES and STM experiments. Our study defines a path to studying the long-sought collective charge modes in quantum materials at the meV scale and with full momentum control.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Ion-transport engineering of alkaline-earth hydrides for hydride electrolyte applications

The heavier alkaline-earth hydrides (AeH 2 ; Ae = Ca, Sr, Ba) are promising materials for hydrogen energy applications, due to their excellent ionic conductivity and thermal stability. We use first-principles calculations to elucidate their defect chemistry and ion-transport mechanisms. Based on calculated formation energies, we find that hydrogen vacancies ( V H ) are present in large concentrations, and hence, transport is vacancy-mediated. The vacancy migration barrier depends strongly on the charge state of the vacancy, with V H + yielding the lowest barriers. The activation energy for ionic transport, which is the sum of the migration barrier and formation energy, can be further reduced by lowering the formation energy of the vacancies. This goal can be achieved by doping with acceptors, such as alkali-metal impurities, which induces larger concentrations of V H + . We show that, with optimized conditions for doping, ionic conductivities can be improved by several orders of magnitude; for BaH 2 , the conductivity is on par with proton conductivity in the best proton-conducting oxides. For BaH 2 , we find K to be the optimal dopant for enhancing conductivity, and for CaH 2 and SrH 2 , we find Na to be optimal. Here, the resulting large ionic conductivity of hydrogen makes the heavier alkaline-earth hydrides strong candidates for electrolytes in hydrogen fuel cells.

08 HYDROGEN↗
DOE OSTI2019

Atlas of putative minima and low-lying energy networks of water clusters n = 3 - 25

We report a database consisting of the putative minima and ~3.2 million local minima lying within 5 kcal/mol from the putative minima for water clusters of sizes n = 3 – 25 using an improved version of the Monte Carlo Temperature Basin Paving (MCTBP) global optimization procedure in conjunction with the ab-initio based, flexible, polarizable Thole-Type Model (TTM2.1-F, version 2.1) interaction potential for water. Several of the low-lying structures, as well as lowlying penta-coordinated water networks obtained with the TTM2.1-F potential, were further refined at the Møller-Plesset second order perturbation MP2/aug-cc-pVTZ level of theory. In total, we have identified 3,138,303 networks corresponding to local minima of the clusters n = 3 – 25, whose Cartesian coordinates and relative energies can be obtained from the webpage https://sites.uw.edu/wdbase/. Networks containing penta-coordinated water molecules start to appear at n = 11 and, quite surprisingly, are energetically close (within 1 – 3 kcal/mol) to the putative minima, a fact that has been confirmed from the MP2 calculations. This large database of water cluster minima spanning quite dissimilar hydrogen bonding networks is expected to influence the development and assessment of the accuracy of interaction potentials for water as well as lower scaling electronic structure methods (such as different density functionals). Furthermore, it can also be used in conjunction with data science approaches (including but no limited to neural networks, machine and deep learning) to understand the properties of water, nature’s most important substance.

Rakshit, Avijit↗
DOE OSTI2016

Vibrational spectroscopic and structural investigations on fullerene: A DFT approach

The molecular structure of fullerene (C{sub 60}) molecule was optimized by the DFT/B3LYP method with 6-31G and 6-31G(d,p) basis sets using Gaussian 09 program. The vibrational frequencies were calculated for the optimized molecular structure of the molecule. The calculated vibrational frequencies confirm that the molecular structure of the molecule was located at the minimum energy potential energy surface. The calculated vibrational frequencies were assigned on the basis of functional group analysis and also confirmed using the GaussView 05 software. The frontier molecular orbitals analysis was carried out. The FMOs related molecular properties were predicted. The higher ionization potential, higher electron affinity, higher softness, lower band gap energy and lower hardness values were obtained, which confirm that the fullerene molecule has a higher molecular reactivity. The Mulliken atomic charge distribution of the molecule was also calculated. Hence, these results play an important role due to its potential applications as drug delivery devices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Toward Long-Life, Ultrahigh-Nickel Layered Oxide Cathodes for Lithium-Ion Batteries: Optimizing the Interphase Chemistry with a Dual-Functional Polymer

Boosting the Ni content in LiMO 2 (M = Ni, Co, Mn, etc.) layered oxides is a promising way to establish high-energydensity, low-cost cathodes, but the poor cathode surface stability is a daunting challenge for their practical viability. Herein, by constructing a dual-functional binder framework with a conductive polymer - polyaniline (PANI), the ultrahigh-Ni layered oxide cathode (LiNi 0.94 Co 0.06 O 2 ) exhibits significantly improved cyclability, with a capacity retention greatly increased from 47% to 81% over 1,000 cycles in full cells. It is demonstrated that the acidic species (e.g. HF) in the electrolyte can be efficiently scavenged through a protonation process of PANI, hence the cathode surface reactivity is greatly suppressed and the rock-salt phase propagation into bulk structure is considerably alleviated. Furthermore, the PANI binder system effectively prevents both the cathode-electrolyte interphase (CEI) and anode-electrolyte interphase (AEI) from degrading to a thick “triple-layer” architecture upon extensive cycling, resulting in more robust, thinner CEI and AEI with regulated interphasial chemistry. Moreover, the delocalized π-conjugated electrons along the backbone of PANI facilitate fast electron transfer and promote rate capability even at low temperatures (-20 °C). Finally, this work sheds light on rational binder engineering for developing high-energy-density lithium-ion batteries with acceptable cycle life.

25 ENERGY STORAGE↗
DOE OSTI2021

Substituted ALPO-5 Zeolites as Promising O 2 Sorption Pump Materials: A Density Functional Theory Study

Many chemical processes depend on having an environment that is low in oxygen partial pressure (P O2 < 100 Pa); sorption pumps are a promising route to establishing that environment by either oxygen pumping or oxygen separation from an inert gas. Near ambient sorption-based processes rely on either pressure or thermal swings, requiring no moving parts, no electricity, and neither very high nor very low temperatures. In this work, we use ab initio calculations to explore zeolites as a class of materials for sorption based oxygen pumping/separation. Our calculations indicate that while the neat ALPO-5 zeolite does not adsorb O 2 , O 2 does adsorb on zeolites selectively substituted with transition metals and metalloids and hence, can enable separation and pumping. ALPO-5 substituted with Si, Ge, Sn, Pd, Pt, Ti, V, Cr, Mn, Zr, Mo, Hf, W, Ce, and Pr provides adsorption energies ranging from -0.19 to -3.92 eV, (-) indicates exothermic process. Additionally, we provide a comprehensive understanding of what controls the adsorption energy: 1) the substi-tutions must be able to adopt an oxidation state that is more positive than the cation it replaces, and 2) the size of the pore into which the O 2 adsorbs to the wall. Additionally, we conduct a ther-modynamic analysis of a thermal swing cycle to approximate the optimal O 2 binding energy for 2 low energy O 2 pumping/separation. We find that the minimum energy cost likely occurs when the adsorption energy is in the range of 0.75 – 1.00 eV (72 – 97 kJ·mol -1 ), which corresponds to Ge, V, Pt, or Ce substituted ALPO-5.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Electrospun Ti–Zr Oxide Heterostructures Enable Strongly Anchored Ultralow-Ir Anodes for Durable Acidic Oxygen Evolution

Proton-exchange-membrane water electrolysis (PEMWE) requires acidic oxygen-evolution-reaction (OER) anodes that combine high activity, high durability, and low Ir loading. Here, we report a Ti-Zr composite electrospun oxide (ESO) nanorod support that enables ultralow-Ir anodes for high-performance PEMWE. Zr-containing Ti oxide heterostructures stabilize anatase-rich TiO2, tune the local oxygen-coordination environment, and strengthen interfacial anchoring of IrOx under acidic anodic conditions. The electrospun nanorod network further creates an open, mechanically coherent catalyst layer that improves Ir utilization, ionomer penetration, and mass transport. At an anode loading of 0.2 mgIr cm-2, the optimized Ir/TiZr20-ESO anode delivers a PEMWE mass activity of 0.99 A mgIr-1 at 1.45 V, 28.3 and 43.0 times higher than commercial Ir black and commercial IrO2/TiO2, respectively. The same anode reaches 3.0 and 4.0 A cm-2 at 1.75 and 1.83 V, respectively, and sustains 2000 h operation at 2.0 A cm-2. Also, accelerated stress tests up to 525 hours over 31,500 cycles confirm promising long-term durability, with an insignificant performance decay of 0.4 μV per cycle. Density functional theory indicates that the Ti-Zr oxide heterostructure suppresses Ti demetallation and strengthens IrO2 interfacial binding, rationalizing the improved high-current-density stability.

25 ENERGY STORAGE↗
DOE OSTI2022

Enhancing glycerol electrooxidation from synergistic interactions of platinum and transition metal carbides

To effectively utilize glycerol as a fuel for electrochemical fuel cells, it is necessary to optimize catalysts for effective C-C bond cleavage and complete oxidation of reaction intermediates to achieve maximum efficiency. Here, the current work showed that the synergistic interactions of platinum (Pt) with transition metal carbide (TMC) substrates, such as tungsten carbide (WC) and tantalum carbide (TaC), fulfilled these criteria. The TMC-supported Pt catalysts showed higher activity and selectivity for complete glycerol oxidation than commercial 10 wt% Pt/C. In-situ FTIR analysis revealed that 5 wt% Pt/WC was the most effective catalyst among those tested for complete glycerol oxidation at 0.9 V vs RHE. In-situ X-ray absorption fine structure characterization and density functional theory calculations provided additional insight into the synergistic interactions for glycerol oxidation over Pt/TMC catalysts.

25 ENERGY STORAGE↗
DOE OSTI2014

Integrated Nano-Domains of Disordered and Ordered Spinel Phases in LiNi 0.5 Mn 1.5 O 4 for Li-Ion Batteries

Recent calculations and experimental data suggest that understanding the local ordering behavior of Ni/Mn will be critical to optimize the electrochemical properties of LiNi 0.5 Mn 1.5 O 4 (LNMO) high voltage spinel. In this study, we systematically controlled the evolution of Ni and Mn ordering in LNMO samples by annealing them at 700 °C in air for different dwelling times, followed by quenching to room temperature. The long- and short-range ordering behavior of Ni and Mn were analyzed by combining neutron powder diffraction, X-ray powder diffraction (XRD), transmission electron microscopy (TEM), and Fourier transform infrared spectroscopy (FT-IR) data. The results show that the fraction of ordered phase increases rapidly during initial annealing at 700 °C for 6 h, and accompanied by decreasing amounts of secondary phases. Annealing longer than 6 h led to the growth in size of ordered domains (i.e., increased segregation of ordered and disordered domains) along with a slow increase in the fraction of ordered phase. The dependence of open circuit voltages (OCVs) of the LNMO on the degree of ordering agrees well with recent calculations using the density functional theory. The increase in the degree of ordering increases the open circuit voltage (OCV) and the initial capacity but reduces cycle life and rate capability. The LNMO delivered optimal battery performances (capacity, cycle life, and rate capability) after annealing at 700 °C for 2 h. This partially ordered sample showed the respective advantages from both disordered and ordered spinels: better spinel-phase purity (thus, higher initial capacity) from the ordered LNMO and better cycle life and rate capability from the disordered LNMO.

25 ENERGY STORAGE↗
DOE OSTI2024

An integrated high-throughput robotic platform and active learning approach for accelerated discovery of optimal electrolyte formulations

Solubility of redox-active molecules is an important determining factor of the energy density in redox flow batteries. However, the advancement of electrolyte materials discovery has been constrained by the absence of extensive experimental solubility datasets, which are crucial for leveraging data-driven methodologies. In this study, we design and investigate a highly automated workflow that synergizes a high-throughput experimentation platform with a state-of-the-art active learning algorithm to significantly enhance the solubility of redox-active molecules in organic solvents. Our platform identifies multiple solvents that achieve a remarkable solubility threshold exceeding 6.20 M for the archetype redox-active molecule, 2,1,3-benzothiadiazole, from a comprehensive library of more than 2000 potential solvents. Significantly, our integrated strategy necessitates solubility assessments for fewer than 10% of these candidates, underscoring the efficiency of our approach. Our results also show that binary solvent mixtures, particularly those incorporating 1,4-dioxane, are instrumental in boosting the solubility of 2,1,3-benzothiadiazole. Beyond designing an efficient workflow for developing high-performance redox flow batteries, our machine learning-guided high-throughput robotic platform presents a robust and general approach for expedited discovery of functional materials.

25 ENERGY STORAGE↗
DOE OSTI2023

Record high room temperature resistance switching in ferroelectric-gated Mott transistors unlocked by interfacial charge engineering

Abstract The superior size and power scaling potential of ferroelectric-gated Mott transistors makes them promising building blocks for developing energy-efficient memory and logic applications in the post-Moore’s Law era. The close to metallic carrier density in the Mott channel, however, imposes the bottleneck for achieving substantial field effect modulation via a solid-state gate. Previous studies have focused on optimizing the thickness, charge mobility, and carrier density of single-layer correlated channels, which have only led to moderate resistance switching at room temperature. Here, we report a record high nonvolatile resistance switching ratio of 38,440% at 300 K in a prototype Mott transistor consisting of a ferroelectric PbZr 0.2 Ti 0.8 O 3 gate and an R NiO 3 ( R : rare earth)/La 0.67 Sr 0.33 MnO 3 composite channel. The ultrathin La 0.67 Sr 0.33 MnO 3 buffer layer not only tailors the carrier density profile in R NiO 3 through interfacial charge transfer, as corroborated by first-principles calculations, but also provides an extended screening layer that reduces the depolarization effect in the ferroelectric gate. Our study points to an effective material strategy for the functional design of complex oxide heterointerfaces that harnesses the competing roles of charge in field effect screening and ferroelectric depolarization effects.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

An improved guess for the variational calculation of charge-transfer excitations in large systems

Ab initio quantum-chemical methods that perform well for computing the electronic ground state are not straightforwardly transferable to electronically excited states, particularly in large molecular systems. Wave function theory offers high accuracy, but is often prohibitively expensive. Methods based on time-dependent density functional theory (TD-DFT) are crucially sensitive to the chosen exchange-correlation functional (XCF) parameterization, and system-specific tuning protocols were therefore proposed to address the method's robustness. Methods based on the variational relaxation of the excited-state electron density showcased promising results for the calculation of charge-transfer excitations, but the complex shape of the electronic hypersurface makes convergence to a specific excited state much more difficult than for the ground state when standard variational techniques are applied. We address the latter aspect by providing suitable initial guesses, which we obtain by two separate constrained algorithms. Combined with the squared-gradient minimization algorithm for all-electrons relaxation in a freeze-and-release scheme (FRZ-SGM), we demonstrate that orbital-optimized density functional theory (OO-DFT) calculations can reliably converge to the charge-transfer states of interest even for large molecular systems. We test the FRZ-SGM method on a phenothiazine-anthraquinone CT excitation in a supramolecular Pd(II) coordination cage complex as a function of the cage conformation. This compound has been studied experimentally prior to our work. We compare this freeze-and-release scheme to two XCF reparameterizations, which were recently proposed as low-cost TD-DFT-based alternatives to variational methods. Two dye-semiconductor complexes, which were previously investigated in the context of photovoltaic applications, serve as a second example to investigate the convergence and stability of the FRZ-SGM approach. Our results demonstrate that FRZ-SGM provides reliable convergence for charge-transfer excited states and avoids variational collapse to lower-lying electronic states, whereas time-dependent DFT calculations with an adequate tuning procedure for the range-separation parameter provide a computationally efficient initial estimate of the corresponding energies, with a computational cost comparable to that of configuration-interaction singles (CIS) calculations.

Bogo, Nicola↗
DOE OSTI2020

A Generalized Variational Principle with Applications to Excited State Mean Field Theory

We present a generalization of the variational principle that is compatible with any Hamiltonian eigenstate that can be specified uniquely by a list of properties. This variational principle appears to be compatible with a wide range of electronic structure methods, including mean field theory, density functional theory, multireference theory, and quantum Monte Carlo. Like the standard variational principle, this generalized variational principle amounts to the optimization of a nonlinear function that, in the limit of an arbitrarily flexible wave function, has the desired Hamiltonian eigenstate as its global minimum. Unlike the standard variational principle, it can target excited states and select individual states in cases of degeneracy or near-degeneracy. As an initial demonstration of how this approach can be useful in practice, we employ it to improve the optimization efficiency of excited state mean field theory by an order of magnitude. With this improved optimization, we are able to demonstrate that the accuracy of the corresponding second-order perturbation theory rivals that of singles-and-doubles equation-of-motion coupled cluster in a substantially broader set of molecules than could be explored by our previous optimization methodology.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗