Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “State functions”

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,171 records · Page 65
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Fast model-based scenario optimization in NSTX-U enabled by analytic gradient computation

Model-based optimization offers a systematic approach to advanced scenario planning. In this case, the feedforward-control inputs (actuator trajectories) that are needed to attain and sustain a desired scenario are obtained by solving a nonlinear constrained optimization problem. This class of problems generally minimize a cost function that measures the difference between desired and actual plasma states. Several numerical optimization algorithms, such as sequential quadratic programming, require repeated calculation of the cost function gradients with respect to the input trajectories. Calculating these gradients numerically can be computationally intensive, increasing the time needed to solve the feedforward-control optimization problem. Here, this work introduces a method to analytically calculate these cost function gradients from the current profile evolution model. This can significantly reduce the computational time and allow for fast feedforward-control optimization, which would eventually enable optimal scenario planning between discharges. The performance of the feedforward optimizer with analytical gradients is compared to a traditional optimization algorithm based on numerical gradients for different NSTX-U scenarios. The plasma dynamics in the optimization algorithm are simulated using the Control Oriented Transport SIMulator (COTSIM). Results of the work show that analytical gradients consistently reduce the computation time while achieving trajectories that are comparable to those obtained by traditional optimization algorithms based on numerical gradients.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Crystal Structure of Colloidally Prepared Metastable Ag 2 Se Nanocrystals

Structural polymorphism is known for many bulk materials; however, on the nanoscale metastable polymorphs tend to form more readily than in the bulk, and with more structural variety. One such metastable polymorph observed for colloidal Ag 2 Se nanocrystals has traditionally been referred to as the “tetragonal” phase. While there are reports on the chemistry and properties of this metastable polymorph, its crystal structure, and therefore electronic structure, has yet to be determined. We report that an anti-PbCl 2 -like structure type (space group P2 1 /n) more accurately describes the powder X-ray diffraction and X-ray total scattering patterns of colloidal Ag 2 Se nanocrystals prepared by several different methods. Density functional theory (DFT) calculations indicate that this anti-PbCl 2 -like Ag 2 Se polymorph is a dynamically stable, narrow-band-gap semiconductor. Here, the anti-PbCl 2 -like structure of Ag 2 Se is a low-lying metastable polymorph at 5–25 meV/atom above the ground state, depending on the exchange-correlation functional used.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2024

Electronic Structure of Excitons in Hematite Fe 2 O 3

We present a computational study of the electronic structure of excitons in hematite Fe 2 O 3 using density functional theory (DFT) and time-dependent DFT theory. Excitons resulting from solar light absorption are precursors to photogenerated charge carriers, which are the active species in solar-to-fuel redox conversions. Upon vertical excitation, the lowest energy exciton is found to have a multideterminantal wave function, describing an electron excited from O 2p states to Fe 3d states of the nearest basal plane, possibly supporting a recent XUV-derived exciton structure with a radius of a single Fe–O bond. Subsequently, the exciton self-traps into stable electron–hole pair structures whose wave functions are well-approximated by single excitations conveniently labeled “HOMO–LUMO”. Self-trapping yields electron–hole pair structures whereby the electron state (electron polaron-like) is separated from the hole state (hole polaron-like) by 3, 5, 7, 9, ... basal planes in structures referred to as Exc-3, Exc-5, Exc-7, Exc-9, ... The natural transition orbitals (electron–hole states) exhibit a strong “localized” character. The hole state is best described as a (FeO 6 ) + octahedral moiety carrying ~0.65 h + charge, and overall ~70% of the hole is assigned to O 2p atomic states, mostly on the four (4) equatorial O atoms of the moiety. The electron state is best described as a (FeFe) – moiety carrying ~0.80 e – charge in Fe 3d atomic states. The lattice distortions around the hole site exhibit Fe–O bond shortening due to the removal of electron density from an orbital state with O–Fe–O antibonding character. The electron site shows Fe–Fe distance shortening due to excess electron density added to an orbital state with in-phase Fe–Fe interactions. Furthermore, these excitonic structures with increasing electron–hole separation can be viewed as the onset of exciton separation into isolated electrons and hole polarons.

36 MATERIALS SCIENCE↗
DOE OSTI2024

In Situ Alloying Enabled by Active Liquid Metal Filler for Self‐Healing Composite Polymer Electrolytes

Abstract Solid inorganics, known for kinetically inhibiting polymer crystallization and enhancing ionic conductivity, have attracted significant attention in solid polymer electrolytes. However, current composite polymer electrolytes (CPEs) are still facing challenges in Li metal batteries, falling short of inhibiting severe dendritic growth and resulting in very limited cycling life. This study introduces Ga 62.5 In 21.5 Sn 16 (Galinstan) liquid metal (LM) as an active liquid alternative to conventional passive solid fillers, aiming at realizing self‐healing protection against dendrite problems. Compared to solid inorganics, for example silica, LM droplets could more significantly reduce polymer crystallinity and enhance Li‐ion conductivity due to their liquid nature, especially at temperatures below the polymer melting point. More importantly, LMs are unraveled as dynamic chemical traps, which are capable of blocking and consuming lithium dendrites upon contact via in situ alloying during battery operation and further inhibiting dendritic growth due to the lower deposition energy barrier of the formed Li‐LM alloy. As a proof of concept, by strategically designing an asymmetric CPE with the active LM filling, a solid‐state Li/LiFePO 4 battery achieves promising full‐cell functionality with notable rate performance and stable cycle life. This active filler‐mediated self‐healing approach could bring new insights into the battery design in versatile solid‐state systems.

Wu, Kai↗
DOE OSTI2025

Achieving designed texture and flows in bulk active nematics using optimal control theory

Being intrinsically nonequilibrium, active materials can potentially perform functions that would be thermodynamically forbidden in passive materials. However, active systems have diverse local attractors that correspond to distinct dynamical states, many of which exhibit chaotic turbulent-like dynamics and thus cannot perform work or useful functions. Designing such a system to choose a specific dynamical state is a formidable challenge. Motivated by recent advances enabling optogenetic control of experimental active materials, we describe an optimal control theory framework that identifies a spatiotemporal sequence of light-generated activity that drives an active nematic system toward a prescribed dynamical steady state. Active nematics are unstable to spontaneous defect proliferation and chaotic streaming dynamics in the absence of control. We demonstrate that optimal control theory can compute activity fields that redirect the dynamics into a variety of alternative dynamical programs and functions. This includes dynamically reconfiguring between states, selecting and stabilizing emergent behaviors that do not correspond to attractors, and are hence unstable in the uncontrolled system. Furthermore, our results provide a roadmap to leverage optical control methods to rationally design structure, dynamics, and function in a wide variety of active materials.

Complex systems theory↗
DOE OSTI2023

State Preparation of Antisymmetrized Geminal Power on a Quantum Computer without Number Projection

The antisymmetrized geminal power (AGP) is equivalent to the number projected Bardeen–Cooper–Schrieffer (PBCS) wave function. It is also an elementary symmetric polynomial (ESP) state. We generalize previous research on deterministically implementing the Dicke state to a state preparation algorithm for an ESP state, or equivalently AGP, on a quantum computer. Our method is deterministic and has polynomial cost, and it does not rely on number symmetry breaking and restoration. We also show that our circuit is equivalent to a disentangled unitary paired coupled cluster operator and a layer of unitary Jastrow operator acting on a single Slater determinant. Here, the method presented herein highlights the ability of disentangled unitary coupled cluster to capture nontrivial entanglement properties that are hardly accessible with traditional Hartree–Fock based electronic structure methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Bayesian refinement of covariant energy density functionals

The last five years have seen remarkable progress in our quest to determine the equation of state of neutron rich matter. Here, recent advances across the theoretical, experimental, and observational landscape have been incorporated in a Bayesian framework to refine existing covariant energy density functionals previously calibrated by the properties of finite nuclei. In particular, constraints on the maximum neutron star mass from pulsar timing, on stellar radii from the NICER mission, on tidal deformabilities from the LIGO-Virgo collaboration, and on the dynamics of pure neutron matter as predicted from chiral effective field theories have resulted in significant refinements to the models, particularly to those predicting a stiff symmetry energy. Still, even after these improvements, we find it challenging to reproduce simultaneously the neutron skin thickness of both 208 Pb and 48 Ca recently reported by the PREX/CREX collaboration.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Photocatalytic Transfer Hydrogenation in Water: Insight into Mechanism and Catalyst Speciation

Organometallic catalysts with strong metal-to-ligand charge-transfer (MLCT) optical transitions and reactive metal-hydride bonds provide opportunities for new excited-state reactivity. Toward this end, we have investigated photoinduced transfer hydrogenation reactions for a series of Cp*IrH complexes with functionalized bipyridyl ligands. Overcoming the baseline reactivity of this class of catalyst toward hydrogen evolution in aqueous media presents a major challenge for focusing this photochemical reactivity into desired transformations of more complex organic substrates. In this work, a photochemically induced aqueous ketone reduction catalyzed by [Cp*Ir(bpy)H][Cl] is reported and benchmarked to the competing hydrogen evolution reaction (HER). A detailed mechanistic study is conducted of the reduction of cyclohexanone with sodium formate delivering the reducing equivalent. Examining the kinetic profile of the cyclohexanone reduction revealed an unexpected induction period, suggesting a more complex underlying reaction mechanism where product selectivity was not simply a function of direct competition from a common excited-state reactive intermediate. The reaction kinetics were analyzed by variable time normalization analysis (VTNA), and the catalyst resting state was interrogated through operando light-emitting diode (LED) illumination H-1 NMR spectroscopy. Overall, these and other mechanistic studies provide insight into a photochemical reaction pathway to ketone reduction that is distinct from that previously reported for HER.

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

Accurate prediction of core-level spectra of radicals at density functional theory cost via square gradient minimization and recoupling of mixed configurations

State-specific orbital optimized approaches are more accurate at predicting core-level spectra than traditional linear-response protocols, but their utility had been restricted due to the risk of "variational collapse" down to the ground state. We employ the recently developed square gradient minimization algorithm to reliably avoid variational collapse and study the effectiveness of orbital optimized density functional theory (DFT) at predicting second period element 1s core-level spectra of open-shell systems. Several density functionals (including SCAN, B3LYP, and ωB97X-D3) are found to predict excitation energies from the core to singly occupied levels with high accuracy (≤0.3 eV RMS error) against available experimental data. Higher excited states are, however, more challenging by virtue of being intrinsically multiconfigurational. We thus present a configuration interaction inspired route to self-consistently recouple single determinant mixed configurations obtained from DFT, in order to obtain approximate doublet states. This recoupling scheme is used to predict the C K-edge spectra of the allyl radical, the O K-edge spectra of CO + , and the N K-edge of NO 2 with high accuracy relative to experiment, indicating substantial promise in using this approach for the computation of core-level spectra for doublet species [vs more traditional time dependent DFT, equation of motion coupled cluster singles and doubles (EOM-CCSD), or using unrecoupled mixed configurations]. Finally, we also present general guidelines for computing core-excited states from orbital optimized DFT.

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

Visualization of defect induced in-gap states in monolayer MoS 2

Atomic-scale intrinsic defects play a key role in controlling functional electronic properties of two-dimensional (2D) materials. Here, we present a low-temperature scanning–tunneling microscopy and spectroscopy investigation of a common point-defect in monolayer molybdenum disulfide (MoS 2 ). We employ a sample preparation method in which the film surface is never exposed to air so that the native dangling bonds surrounding the defects in the film are preserved. Molybdenum vacancies are identified by their three characteristic in-gap resonances by combining scanning–tunneling measurements with parallel Green’s function-based theoretical modeling. The relative energy shifts between the various in-gap states allow us to identify a relative charge difference between two of the observed vacancies. The role of the substrate on the band structure of the defective MoS 2 monolayer is unveiled. Our study highlights the effects of the substrate on the in-gap states of common defects found in MoS 2 providing a pathway in designing and optimizing 2D materials for electronic applications.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Unconventional surface state pairs in a high-symmetry lattice with anti-ferromagnetic band-folding

Many complex magnetic structures in a high-symmetry lattice can arise from a superposition of well-defined magnetic wave vectors. These “multi-q” structures have garnered much attention because of interesting real-space spin textures such as skyrmions. However, the role multi-q structures play in the topology of electronic bands in momentum space has remained rather elusive. Here we show that the type-I anti-ferromagnetic 1q, 2q and 3q structures in an face-centered cubic sublattice with band inversion, such as NdBi, can induce unconventional surface state pairs inside the band-folding hybridization bulk gap. Our density functional theory calculations match well with the recent experimental observation of unconventional surface states with hole Fermi arc-like features and electron pockets below the Neel temperature. We further show that these multi-q structures have Dirac and Weyl nodes. Our work reveals the special role that band-folding from anti-ferromagnetism and multi-q structures can play in developing new types of surface states.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

On the Language of Reliability: A System Engineer Perspective

In its classical definition, risk is defined by three elements: what can go wrong, what are its consequences, and how likely is it to occur? While this definition makes sense in a regulatory-based framework where for the current fleet of operating light water reactors (LWRs), the risks associated with nuclear power plants typically are characterized in terms of core damage and large early release frequency (LERF), this approach does not provide a useful snapshot of the health of the plant from a broader perspective. This is due to the very narrow context in which the term “risk” typically is defined as nuclear safety aspects that have the potential to impact public health. In this paper, we take the viewpoint of nuclear safety that is reflective of the current fleet of operating LWRs for which core damage frequency and LERF are appropriate metrics. For other advanced reactor designs, other more applicable technology neutral metrics of reactor safety metrics would be specified. A possible alternate path would start by redefining the word risk with a broader meaning that better reflects the needs of a system health and asset management decision-making process. Rather than asking how likely an event could occur (in probabilistic terms), we can ask how far this event is from occurring. Our approach starts by defining and quantifying component and system health in terms of a “distance” between its actual and limiting conditions, i.e., determination of the margin that exists between the current state/condition and the state where the component/system is no longer capable of achieving its intended function. A margin is a measure that is more reflective of the current state or performance of a component, and therefore more closely tied to decisions that are made on an ongoing basis. Further, we will show how, given the data available from plant equipment reliability and monitoring (e.g., pump vibration data) and prognostic (e.g., component remaining useful life estimation) data, a margin can be described and determined for all types of maintenance approaches (e.g., corrective or predictive maintenance). We show how classical reliability models (e.g., fault trees) can be used to quantify the system margin provided component margin values. In the approach described in this paper, the propagation of margin values through classical reliability models are not performed using classical probabilistic calculations applied to sets (as performed in a typical plant probabilistic risk assessment). Instead, we show how it is possible to propagate margin values through Boolean logic gates (i.e., AND and OR operators) through distance-based operations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Interpreting Tafel behavior of consecutive electrochemical reactions through combined thermodynamic and steady state microkinetic approaches

Assessing the reaction pathway of multi-electron-transfer reactions is an essential yet difficult task for the rational design of electrocatalysts. In this work, we develop a heuristic approach that combines thermodynamic adsorption energetics calculated through density functional theory with microkinetic modeling using the steady state approximation to interpret the potential-dependent Tafel behavior of consecutive electrochemical reactions. In doing so, we introduce a kinetic framework for ab initio calculations that ensures self-consistent adsorption energetics based on kinetically limited adsorbate coverages. The approach is applied to experimental results on CoO x (OH) 2-x single crystal electrocatalyst particles yielding coverage dependent mechanistic information and identification of the rate-limiting step with standard rate constants for the oxygen evolution reaction on the ($11\bar{2}0$) surfaces of the β-Co(OH) 2 , β-CoOOH, and CoO 2 bulk phases. Finally, this generalizable method enables catalyst benchmarking based on determining the active species involved and associated intrinsic reaction rate constants in consecutive multi-electron-transfer reactions.

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

Three-state harmonic models for photoinduced charge transfer

A widely used strategy for simulating the charge transfer between donor and acceptor electronic states in an all-atom anharmonic condensedphase system is based on invoking linear response theory to describe the system in terms of an effective spin-boson model Hamiltonian. Extending this strategy to photoinduced charge transfer processes requires also taking into consideration the ground electronic state in addition to the excited donor and acceptor electronic states. In this paper, we revisit the problem of describing such nonequilibrium processes in terms of an effective three-state harmonic model. We do so within the framework of nonequilibrium Fermi’s golden rule (NE-FGR) in the context of photoinduced charge transfer in the carotenoid–porphyrin–C 60 (CPC 60 ) molecular triad dissolved in explicit tetrahydrofuran (THF). To this end, we consider different ways for obtaining a three-state harmonic model from the equilibrium autocorrelation functions of the donor–acceptor, donor–ground, and acceptor–ground energy gaps, as obtained from all-atom molecular dynamics simulations of the CPC 60 /THF system. The quantum-mechanically exact time-dependent NE-FGR rate coefficients for two different charge transfer processes in two different triad conformations are then calculated using the effective three-state model Hamiltonians as well as a hierarchy of more approximate expressions that lead to the instantaneous Marcus theory limit. Our results show that the photoinduced charge transfer in CPC 60 /THF can be described accurately by the effective harmonic three-state models and that nuclear quantum effects are small in this system.

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

Constraining the particle-scale diversity of black carbon light absorption using a unified framework

Abstract. Atmospheric black carbon (BC), the strongest absorber of visible solar radiation in the atmosphere, manifests across a wide spectrum of morphologies and compositional heterogeneity. Phenomenologically, the distribution of BC among diverse particles of varied composition gives rise to enhancement of its light absorption capabilities by over twofold in comparison to that of nascent or unmixed homogeneous BC. This situation has challenged the modeling community to consider the full complexity and diversity of BC on a per-particle basis for accurate estimation of its light absorption. The conventionally adopted core–shell approximation, although computationally inexpensive, is inadequate not only in estimating but also capturing absorption trends for ambient BC. Here we develop a unified framework that encompasses the complex diversity in BC morphology and composition using a single metric, the phase shift parameter (ρBC), which quantifies how much phase shift the incoming light waves encounter across a particle compared to that in its absence. We systematically investigate variations in ρBC across the multi-space distribution of BC morphology, mixing state, mass, and composition as reported by field and laboratory observations. We find that ρBC>1 leads to decreased absorption by BC, which explains the weaker absorption enhancements observed in certain regional BC compared to laboratory results of similar mixing state. We formulate universal scaling laws centered on ρBC and provide physics-based insights regarding core–shell approximation overestimating BC light absorption. We conclude by packaging our framework in an open-source Python application to facilitate community-level use in future BC-related research. The package has two main functionalities. The first functionality is for forward problems, wherein experimentally measured BC mixing state and assumed BC morphology are input, and the aerosol absorption properties are output. The second functionality is for inverse problems, wherein experimentally measured BC mixing state and absorption are input, and the morphology of BC is returned. Further, if absorption is measured at multiple wavelengths, the package facilitates the estimation of the imaginary refractive index of coating materials by combining the forward and inverse procedures. Our framework thus provides a computationally inexpensive source for calculation of absorption by BC and can be used to constrain light absorption throughout the atmospheric lifetime of BC.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Materials Learning Algorithms (MALA): Scalable machine learning for electronic structure calculations in large-scale atomistic simulations

We present the Materials Learning Algorithms (MALA) package, a scalable machine learning framework designed to accelerate density functional theory (DFT) calculations suitable for large-scale atomistic simulations. Using local descriptors of the atomic environment, MALA models efficiently predict key electronic observables, including local density of states, electronic density, density of states, and total energy. The package integrates data sampling, model training and scalable inference into a unified library, while ensuring compatibility with standard DFT and molecular dynamics codes. We demonstrate MALA's capabilities with examples including boron clusters, aluminum across its solid-liquid phase boundary, and predicting the electronic structure of a stacking fault in a large beryllium slab. Scaling analyses reveal MALA's computational efficiency and identify bottlenecks for future optimization. With its ability to model electronic structures at scales far beyond standard DFT, MALA is well suited for modeling complex material systems, making it a versatile tool for advanced materials research.

Density functional theory↗
DOE OSTI2026

Low-Lying Excited States of Linear All- Trans Polyenes: Insights from Analytic Gradient and Nonadiabatic Coupling Calculations Based on Multireference Configuration Interaction

Polyenes serve as a rigorous test for theoretical models and electronic structure methods, playing a key role in advancing computational and theoretical chemistry. Here, we present a high-level theoretical investigation of linear, all-trans polyenes using energy gradients and nonadiabatic coupling vectors based on an MR-CISD wave function to describe electronic transitions involving the ground state (1 1 A g – ) and three low-lying excited states (2 1 A g – , 1 1 B u + , and 2 1 B u – ) of hexatriene, octatetraene, and decapentaene. This approach enables accurate evaluation of both adiabatic and vertical excitation and emission energies, yielding results in excellent agreement with experiment, as well as locating minima on the crossing seam between adiabatic states. Our results show that vertical excitation energies to the 1 1 B u + state are blue-shifted by 0.2–0.3 eV relative to the experimental absorption maximum, whereas the vertical emission energy from the 2 1 A g – state is red-shifted by ∼0.2 eV relative to the experimental emission maximum. Upon relaxation from the Franck–Condon geometry, the 2 1 A g – state stabilizes by around 1 eV, compared to 0.2–0.3 eV for the 1 1 B u + state. An analysis of the S 1 /S 0 crossing seam in hexatriene shows that its minimum involves asymmetric backbone deformations and provides an efficient channel for ultrafast internal conversion to the ground state, consistent with the absence of detectable fluorescence in this molecule. These results demonstrate the power of analytic gradients and nonadiabatic coupling vectors based on an MR-CISD wave function for accurately characterizing the electronic structure and photophysics of polyenes.

Excited states↗
DOE OSTI2024

Anomalous Hall crystals in rhombohedral multilayer graphene. II. General mechanism and a minimal model

Here we propose a minimal "three-patch model"for the anomalous Hall crystal (AHC), a topological electronic state that spontaneously breaks both time-reversal symmetry and continuous translation symmetry. The proposal for this state is inspired by the recently observed integer and fractional quantum Hall states in rhombohedral multilayer graphene at zero magnetic field. There, interaction effects appear to amplify the effects of a weak moiré potential, leading to the formation of stable, isolated Chern bands. It has been further shown that Chern bands are stabilized in mean-field calculations even without a moiré potential, enabling a realization of the AHC state. Our model is built on the dissection of the Brillouin zone into patches centered around high-symmetry points. Within this model, the wave functions at high-symmetry points fully determine the topology and energetics of the state. We extract two quantum geometrical phases of the noninteracting wave functions that control the stability of the topologically nontrivial AHC state. The model predicts that the AHC state wins over the topological trivial Wigner crystal in a wide range of parameters, and agrees very well with the results of full self-consistent Hartree-Fock calculations of the rhombohedral multilayer graphene Hamiltonian.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗