Search NASA⌕ Search

SEARCH · Search NASA

Results for “formal methods”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Mechanism of Catechol Oxidation by the Coupled Bicupric Active Site of Tyrosinase: Completion of the Oxygenase/Oxidase Reaction Cycle

Tyrosinase contains a coupled binuclear copper (CBC) active site, which in its bicuprous form (deoxy-Ty) binds O 2 to form a side-on peroxide [Cu(II) 2 O 2 ] 2+ intermediate (oxy-Ty) that performs the regioselective monooxygenation of monophenols to catechols and their subsequent 2e − oxidation to quinones. Previously, we used spectroscopic, kinetic, and computational methods to elucidate the mechanism of the initial steps in the monooxygenation reaction. Herein, we investigated the final step in catecholate formation and provide experimental and computational results elucidating the catechol oxidation reaction, formally a proton-coupled 2e − process. We employed single-turnover stoppedflow absorption to observe the elusive catecholate-Cu(II) 2 OH intermediate via the anaerobic reaction of resting Cu(II) 2 OH Ty (met-Ty) with the catecholate substrate. This intermediate was cryo-trapped and characterized by electron paramagnetic resonance and X-ray absorption spectroscopies. These experimental results were correlated to quantum mechanics/molecular mechanics (QM/MM) and QM calculations to describe the complete catalytic cycle of tyrosinase, revealing that the 2e − oxidation of catechol to quinone involves two steps: (i) coupled 1H + /1e − transfer from the bound monoanionic catecholate to the bridged hydroxide of the met-Ty active site, followed by (ii) the second 1e − transfer. The latter is calculated to be the rate-limiting step in catechol oxidation, confirmed by experimental solvent kinetic isotope effect studies. Our computational data suggest that quinone release from the protein provides the necessary driving force for this reaction. These results establish a detailed description of this oxidase cycle of the tyrosinase mechanism and more broadly provide molecular-level insights into the diverse reactivity of CBC sites in biology.

QM/MM modeling↗
DOE OSTI2024

The cluster decomposition of the configurational energy of multicomponent alloys

Abstract The cluster expansion method (CEM) is a widely used lattice-based technique in the study of multicomponent alloys. Despite its prevalent use, a clear understanding of expansion terms is lacking. We present a modern mathematical formalism of the CEM and introduce thecluster decomposition—a unique and basis-independent decomposition for functions of the atomic configuration in a crystal. We identify the cluster decomposition as an invariant ANOVA decomposition; and demonstrate how functional analysis of variance and sensitivity analysis can be used to interpret interactions among species. Furthermore, we show how the mathematical structure of the cluster decomposition enables numerical evaluation that scales with the number of clusters and is independent of the number of species. Overall, our work enables rigorous interpretations of interactions among species, provides opportunities to explore parameter estimation beyond linear regression, introduces a numerical efficient implementation, and enables analysis of cluster expansions based on established mathematical and statistical principles.

Chemistry↗
DOE OSTI2025

Fast methods for multisite charge transfer. Processes II. Analytic nuclear gradients and nonadiabatic dynamics for cCASSCF(1,n) and cCASSCF(2n-1,n) wavefunctions

In this work we derive and implement analytic nuclear gradients and derivative couplings for a constrained complete active space self-consistent field with a small active space designed to model electron or hole transfer. Using a Lagrangian formalism, we are able to differentiate both the CASSCF energy and the constraint (which is required for smooth surfaces over a wide range of parameter space), and the resulting efficient algorithm can be immediately applied to nonadiabatic dynamics simulations of charge transfer processes. Here, we run initial surface-hopping simulations of a proton coupled electron transfer event for a phenoxyl–phenol system.

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

X-ray Thomson scattering absolute intensity from the f-sum rule in the imaginary-time domain

We present a formally exact and simulation-free approach for the normalization of X-ray Thomson scattering (XRTS) spectra based on the f-sum rule of the imaginary-time correlation function (ITCF). Our method works for any degree of collectivity, over a broad range of temperatures, and is applicable even in nonequilibrium situations. In addition to giving us model-free access to electronic correlations, this new approach opens up the intriguing possibility to extract a plethora of physical properties from the ITCF based on XRTS experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

cclib 2.0: An updated architecture for interoperable computational chemistry

Interoperability in computational chemistry is elusive, impeded by the independent development of software packages and idiosyncratic nature of their output files. The cclib library was introduced in 2006 as an attempt to improve this situation by providing a consistent interface to the results of various quantum chemistry programs. The shared API across programs enabled by cclib has allowed users to focus on results as opposed to output and to combine data from multiple programs or develop generic downstream tools. Initial development, however, did not anticipate the rapid progress of computational capabilities, novel methods, and new programs; nor did it foresee the growing need for customizability. Here, we recount this history and present cclib 2, focused on extensibility and modularity. We also introduce recent design pivots—the formalization of cclib’s intermediate data representation as a tree-based structure, a new combinator-based parser organization, and parsed chemical properties as extensible objects.

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

Synthesis, Structure, and Properties of the Complex Zintl Phase Eu 9 Zn 4.5 As 9 : A Candidate Topological Insulator and Thermoelectric Material

Reported are the synthesis and detailed analysis of the crystal and electronic structure of the novel Zintl phase Eu 9 Zn 4.5 As 9 . This material was identified in the densely populated Eu–Zn–As compositional space. For structure determination and for property measurements, suitable single crystals of this compound were grown from either Sn- or Pb-flux. Single-crystal X-ray diffraction methods indicate that Eu 9 Zn 4.5 As 9 crystallizes in the orthorhombic crystal system with the space group Pnma (a = 12.1953(7) Å, b = 4.3730(2) Å, c = 42.674(2) Å) and is formally isostructural to Ca 9 Mn 4+x Sb 9 , the less common “9–4–9” type. The structure is heavily disordered, with multiple partially occupied sites, yet, according to the Zintl-Klemm formalism, a charge-balanced composition (Eu 2+ ) 9 (Zn 2+ ) 4.5 (As 3− ) 9 is attained. Electronic structure calculations for a model, disorder-free structure indicate no energy gap between the valence and the conduction bands and suggest (semi)metallic behavior. Preliminary susceptibility measurements confirm the expected divalent nature of Eu 2+ ([Xe] 4 f 7 ground state).

Zintl phases↗
DOE OSTI2025

Exact closed-form unitary transformations of fermionic operators

Unitary transformations play a fundamental role in many-body physics, and except for special cases, they are not expressible in closed form. We present closed-form expressions for unitary transformations generated by a single fermionic operator for Hermitian and anti-Hermitian generators. We demonstrate the usefulness of these expressions in formal analyses of unitary transformations and numerical applications to Hamiltonian downfolding in quantum computing and Heisenberg dynamics. Furthermore, this work paves the way for new analytical treatments of unitary transformations and numerical many-body methods for fermions.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2026

Reduced salen in a monomeric Fe complex

We report the first monomeric formally M I transition metal salen complex without ion pairing. A combination of spectroscopy, crystallography, and computations indicates that reduction of an Fe II salen complex gives substantial reduced ligand character. We also present a computational method for accurately predicting 57 Fe Mössbauer parameters of Fe salen complexes.

Feldman, Dana M. [Yale University, New Haven, CT (↗
DOE OSTI2024

First-order formalism for β functions in bosonic sigma models from supersymmetry breaking

We consider the renormalization group flow equation for the two-dimensional sigma models with the Kähler target space. The first-order formulation allows us to treat perturbations in these models as current-current deformations. We demonstrate, however, that the conventional first-order formalism misses certain anomalies in the measure, and should be amended. We reconcile beta functions obtained within the conformal perturbation theory for the current-current deformations with traditional “geometric” results obtained in the background field methods, in this way resolving the peculiarities pointed out in O. Gamayun et al. [Peculiarities of beta functions in sigma models, J. High Energy Phys. 10 (2023) 097]. The result is achieved by the supersymmetric completion of the first-order sigma model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Extracting scattering amplitudes for arbitrary two-particle systems with one-particle left-hand cuts via lattice QCD

We derive a general formalism that relates the spectrum of two-particle systems in a finite volume to physical scattering amplitudes, taking into account the presence of any left-hand branch cuts due to single-particle exchanges. The method first relates the finite-volume spectrum to an infinite-volume short-range quantity, denoted ${\mathcal{M}}_0$, and then relates the latter to the physical scattering amplitudes via known integral equations. The derivation of both relations is performed using all-orders perturbation theory and is exact up to neglected exponentially suppressed volume dependence. The relations hold for arbitrary two-particle systems with any number of coupled channels, non-identical and non-degenerate particles, and any intrinsic spin.

algorithms↗
DOE OSTI2025

Post-hoc reweighting of hadron production in the Lund string model

We present a method for reweighting flavor selection in the Lund string fragmentation model. This is the process of calculating and applying event weights enabling fast and exact variation of hadronization parameters on pre-generated event samples. The procedure is post hoc, requiring only a small amount of additional information stored per event, and allowing for efficient estimation of hadronization uncertainties without repeated simulation. Weight expressions are derived from the hadronization algorithm itself, and validated against direct simulation for a wide range of observables and parameter shifts. The hadronization algorithm can be viewed as a hierarchical Markov process with stochastic rejections, a structure common to many complex simulations outside of high-energy physics. This perspective makes the method modular, extensible, and potentially transferable to other domains. We demonstrate the approach in Pythia, including both coverage considerations and timing benefits. For the purpose of this paper, our goal is to develop and demonstrate the the formalism, and we therefore exclude several model variations for baryon production (popcorn model, junction production) needed for proton collisions. These will be the topic of a future paper.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Income Verification Strategies for Income-Based Solar Programs

The Inflation Reduction Act has created substantial new programs that support adoption of solar power by low-income households, including the $7 billion Solar For All program and the Low-Income Communities Bonus Credit Program, which increases the investment tax credit for certain types of deployment. In addition, a growing number of states are using solar programs to reduce energy burdens and create energy justice opportunities for low-income households and disadvantaged communities. Verifying the income of participating customers is an important component of these programs. Program managers are seeking strategies to verify a large number of subscribing customers in an accurate, timely, and cost-efficient manner. To help inform program managers, Berkeley Lab investigated how a number of energy and non-energy programs manage income verification. The most common approach is to require proof through tax documents, pay stubs, or other formal income documentation, which can pose an impediment to enrolling eligible customers and create a paperwork burden for administrators. In order to reduce the burden for both the applicant and the program manager, some programs use alternative methods. We identify three common alternative verification methods: -Categorical eligibility: Customers enrolled in other, similar income-verified assistance programs are automatically eligible for enrollment in other income-qualified programs. -Geographic eligibility: Eligibility is based on the customer’s location within a specified area, typically a low-income or disadvantaged community or census tract, and; -“Self-attestation”: The participant claims eligibility with or without further documentation. We describe these options, their pros and cons, give examples of how they are used, and explore how some low-income programs address administrative issues, audits, or other quality control measures. Finally, we explore the risk of mistaken verifications (finding a participant eligible when they are not) in the different strategies. While this memo was initiated by a request relating to income-based community solar programs, the methods are applicable to any program with income eligibility requirements in the energy or non-energy sector. Funding was provided for this research by the Solar Energy Technologies Office of the US Department of Energy, through the National Community Solar Partnership.

14 SOLAR ENERGY↗
DOE OSTI2024

Weather effects on the lifecycle of U.S. Department of Defense equipment replacement (WELDER)

Extreme weather has a direct and significant impact on buildings and infrastructure, resulting in billions of dollars of damage each year. This problem continues to grow as climate patterns change and buildings are exposed to new and different hazards than what they were designed to withstand. In order to better plan for the long-range sustainment, restoration, modernization, and eventual recapitalization of these buildings, organizations with large building portfolios, such as the U.S. Department of Defense (DoD), must have an awareness of the risks that these extreme weather events present. This research aimed to develop an approach to estimate condition loss and reduction in service life for the components of a building due to extreme weather hazards, to understand the risks that may be present in certain buildings and building systems. To achieve this objective, a damage association matrix was developed that categorizes climate hazards, the damage modes that they produce, and the individual component types impacted. This damage matrix formally links state-of-the-art climate model output, which provides projections of the probability of various climate hazards with a damage effects model that quantifies the consequence on component-level condition and service life. This method is applied to an actual portfolio of buildings in a particular geographic location and with a pre-defined component inventory that comprises the building. This approach can be aggregated to the system-, facility-, and site-level thus helping support billions of dollars in recapitalization decisions related to restoration/modernization of facilities.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Source size of x rays from self-modulated laser wakefield accelerators

A comparative study of x-ray sources generated with different mechanisms from self-modulated laser wakefield acceleration (SM-LWFA) electrons was performed to compare the source size or spatial resolution for use in high energy density science applications. We examine the source size of betatron, inverse Compton scattering, and bremsstrahlung radiation with a Fresnel diffraction based formalism and a modified x-ray ray tracing model. We observe the dependence of source size on the radiation generation process, laser parameters, and compare to what is possible in other regimes of LWFA, as well as current methods. This information is significant as we begin to explore the use of light sources driven by SM-LWFA for use as a diagnostic at large-scale laser facilities where blowout regime LWFA is not possible.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Shock Hugoniot calculations using on-the-fly machine learned force fields with ab initio accuracy

We present a framework for computing the shock Hugoniot using on-the-fly machine learned force field (MLFF) molecular dynamics simulations. In particular, we employ an MLFF model based on the kernel method and Bayesian linear regression to compute the free energy, atomic forces, and pressure, in conjunction with a linear regression model between the internal and free energies to compute the internal energy, with all training data generated from Kohn–Sham density functional theory (DFT). We verify the accuracy of the formalism by comparing the Hugoniot for carbon with recent Kohn–Sham DFT results in the literature. In so doing, we demonstrate that Kohn–Sham calculations for the Hugoniot can be accelerated by up to two orders of magnitude, while retaining ab initio accuracy. We apply this framework to calculate the Hugoniots of 14 materials in the FPEOS database, comprising 9 single elements and 5 compounds, between temperatures of 10 kK and 2 MK. We find good agreement with first principles results in the literature while providing tighter error bars. In addition, we confirm that the inter-element interaction in compounds decreases with temperature.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

High pressure structural and lattice dynamics study of α-In2Se3

Layered α-In2Se3 has been studied using a combined in situ synchrotron angle-dispersive powder x-ray diffraction and Raman spectroscopy study in a diamond anvil cell up to 60+ GPa, at room temperature. Helium, which remains fairly hydrostatic up to the highest pressure in this study, was used as the pressure-transmitting medium. The results from both experimental methods reveal a pressure-induced structural phase transition from α-In2Se3 to a monoclinic β'-In2Se3 structure at ≈1 GPa, in agreement with previous studies. Based on our detailed measurements using both experimental techniques and the F-f formalism, the β'-In2Se3 structure remains stable up to 45 GPa, without a clear indication of a phase transition toward the previously reported β-In2Se3 phase. Above this pressure, In2Se3 adopts a disordered solid-solution-like orthorhombic structure, phase IV. The results are discussed in comparison with the relevant previous studies of α-In2Se3 under pressure.

Feng, Shiyu↗
DOE OSTI2024

A linear response framework for quantum simulation of bosonic and fermionic correlation functions

Response functions are a fundamental aspect of physics; they represent the link between experimental observations and the underlying quantum many-body state. However, this link is often under-appreciated, as the Lehmann formalism for obtaining response functions in linear response has no direct link to experiment. Within the context of quantum computing, and via a linear response framework, we restore this link by making the experiment an inextricable part of the quantum simulation. This method can be frequency- and momentum-selective, avoids limitations on operators that can be directly measured, and can be more efficient than competing methods. As prototypical examples of response functions, we demonstrate that both bosonic and fermionic Green’s functions can be obtained, and apply these ideas to the study of a charge-density-wave material on the ibm_auckland superconducting quantum computer. The linear response method provides a robust framework for using quantum computers to study systems in physics and chemistry.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

A Review of the Lawrence Livermore Nuclear Accident Dosimeter 1980s-present

A Nuclear Accident Dosimetry program is a federal requirement for all facilities that have the potential to have a criticality accident. Personnel Nuclear Accident Dosimeter (PNAD) theory and analytical procedures are driven by various scientific needs and interacting regulations. A brief history of the status of USA Department of Energy (DOE) nuclear accident dosimetry regulations, recommendations, and performance testing criteria are given. Then, the history of the Lawrence Livermore National Laboratory (LLNL) PNAD is explored, including changes in the physical dosimeter and adjustments of the analysis method through the last four decades. Finally, the performance of LLNL’s PNAD at criticality accident intercomparison training exercises since 2009 is explored. In general, reported neutron doses have been within or close to DOE-STD-1098 performance criteria while reported gamma doses have been outside of DOE-STD-1098 performance criteria. Reported total absorbed doses have varied in meeting ANSI/HPS N13.3 and ANSI/HPS N13.3 (R2019) performance criteria. Dosimetry staff retirement and turnover have left historical knowledge gaps, yet provided opportunities within the NAD program at LLNL. This review paper serves as an overview of the history and status of the NAD program. Brief technical, procedural and programmatic recommendations to improve LLNL’s NAD program are given. Technical recommendations include investigating orientation factors through modeling or empirical experimentation, investigating gamma dosimetry methods for high-dose scenarios, and exploring other dosimetric methods for simpler, quicker NAD analysis. Procedural recommendations include better documentation of conversion factor (activity-to-fluence and fluence-to-dose) derivations and spectrum uses, and updated analysis spreadsheets or simple Graphic User Interfaces for dose calculations. In conclusion, programmatic recommendations include formalized training for NAD analysts, and having multiple SMEs trained on the NAD program.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗