Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 1,693 records · Page 94
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Hydrogen doping in HfO{sub 2} resistance change random access memory

The structures and energies of hydrogen-doped monoclinic hafnium dioxide were calculated using density-functional theory. The electronic interactions are described within the LDA + U formalism, where on-site Coulomb corrections are applied to the 5d orbital electrons of Hf atoms and 2p orbital electrons of the O atoms. The effects of charge state, defect-defect interactions, and hydrogenation are investigated and compared with experiment. It is found that hydrogenation of HfO{sub 2} resistance-change random access memory devices energetically stabilizes the formation of oxygen vacancies and conductive vacancy filaments through multiple mechanisms, leading to improved switching characteristic and device yield.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Parallel-in-time quantum simulation via Page and Wootters quantum time

In the past few decades, researchers have created a veritable zoo of quantum algorithms by drawing inspiration from classical computing, information theory, and even from physical phenomena. Here, we present quantum algorithms for parallel-in-time simulations that are inspired by the Page and Wootters formalism. In this framework, and thus in our algorithms, the classical time variable of quantum mechanics is promoted to the quantum realm by introducing a Hilbert space of “clock” qubits that are then entangled with the “system” qubits. We show that our algorithms can compute temporal properties over 𝑁 different times of many-body systems by only using log⁡(𝑁) clock qubits. As such, we achieve an exponential trade-off between time and spatial complexities. In addition, we rigorously prove that the entanglement created between the system qubits and the clock qubits has operational meaning, as it encodes valuable information about the system’s dynamics. We also provide a circuit depth estimation of all the protocols, showing a running time advantage in computation times over traditional sequential-in-time algorithms. In particular, for the case when the dynamics are determined by the Aubry-Andre model, we present a hybrid method for which our algorithms have a depth that only scales as 𝒪⁡(log⁡(𝑁)⁢𝑛). As a by-product, we can relate the previous schemes to the problem of equilibration of an isolated quantum system, thus indicating that our framework enables a new dimension for studying dynamical properties of many-body systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2012

On the variational noncommutative poisson geometry

We outline the notions and concepts of the calculus of variational multivectors within the Poisson formalism over the spaces of infinite jets of mappings from commutative (non)graded smooth manifolds to the factors of noncommutative associative algebras over the equivalence under cyclic permutations of the letters in the associative words. We state the basic properties of the variational Schouten bracket and derive an interesting criterion for (non)commutative differential operators to be Hamiltonian (and thus determine the (non)commutative Poisson structures). We place the noncommutative jet-bundle construction at hand in the context of the quantum string theory.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2013

Determination of the maximum energy loss for electron stopping power calculations and its effect on backscattering electron yield in Monte-Carlo simulations applying continuous slowing-down approximation

The maximum energy loss for electron stopping power calculations by the full Penn algorithm within the dielectric formalism is determined with taking into account the contribution of electron and plasmon excitations. Use of these calculated electron stopping powers in Monte Carlo simulations applying continuous slowing down approximation gives the backscattering electron yields in much better agreement with experimental data than previous other theoretical results. The muffin-tin model is used to describe the electron elastic scattering by atom bound in solids with taking into account the exchange correlation and polarization effect.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Exploring nonmultiplicativity in the geometric measure of entanglement

The geometric measure of entanglement (GME) quantifies how close a multipartite quantum state is to the set of separable states under the Hilbert-Schmidt inner product. The GME can be nonmultiplicative, meaning that the closest product state to two states is entangled across subsystems. In this work, we explore the GME in two families of states: those that are invariant under bilateral orthogonal (𝑂⊗𝑂) transformations, and mixtures of singlet states. In both cases, a region of GME nonmultiplicativity is identified around the antisymmetric projector state. We employ state-of-the-art numerical optimization methods and models to quantitatively analyze nonmultiplicativity in these states for 𝑑=3. Here, we also investigate a constrained form of GME that measures closeness to the set of real product states and show that this measure can be nonmultiplicative even for real separable states.

Quantum correlations in quantum information↗
DOE OSTI2024

Stochastic modeling of x-ray superfluorescence

An approach to modeling the dynamics of x-ray amplified spontaneous emission and superfluorescence, the phenomenon of collective x-ray emission initiated by intense pulses of x-ray free-electron lasers, is developed based on stochastic partial differential equations. The equations are derived from first principles, and the relevant approximations, derivation steps, and extensions specific to stimulated x-ray emission are presented. The resulting equations take the form of three-dimensional generalized Maxwell-Bloch equations augmented with noise terms for both field and atomic variables. The derived noise terms possess specific correlation properties that enable the correct reconstruction of spontaneous emission. Consequently, the developed theoretical formalism is universally suitable for describing all stages of stimulated x-ray emission: spontaneous emission, amplified spontaneous emission, and superfluorescence. Here, we present numerical examples that illustrate various properties of the emitted field, including spatiotemporal coherence and spectral-angular and polarization characteristics. We anticipate that the proposed theoretical framework will establish a robust foundation for interpreting measurements in stimulated x-ray emission spectroscopy, modeling x-ray laser oscillators, and describing other experiments leveraging x-ray superfluorescence.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Applied nonrelativistic conformal field theory: Scattering-length and effective-range corrections to rate of production of three neutrons at low relative momenta

Due to an accidentally large s-wave scattering length, in a relatively wide range of energy, neutrons are approximately described by the nonrelativistic conformal field theory of unitarity fermions, perturbed by one relevant and an infinite number of irrelevant operators. We develop a formalism which provides a nonperturbative definition of local operators in that nonrelativistic conformal field theory. We compute the scattering-length and effective-range corrections to the two-point functions of primary charge-three operators using the technique of conformal perturbation theory. These calculations allow us to find the first corrections to the scale-invariant behavior of the rate of nuclear reactions with three neutrons in the final state in the regime when the neutrons have small relative momenta.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Superfluid Weight Bounds from Symmetry and Quantum Geometry in Flat Bands

Flat-band superconductivity has theoretically demonstrated the importance of band topology to correlated phases. In two dimensions, the superfluid weight, which determines the critical temperature through the Berezinksii-Kosterlitz-Thouless criteria, is bounded by the Fubini-Study metric at zero temperature. We show this bound is nonzero within flat bands whose Wannier centers are obstructed from the atoms—even when they have identically zero Berry curvature. Next, we derive general lower bounds for the superfluid weight in terms of momentum space irreps in all 2D space groups, extending the reach of topological quantum chemistry to superconducting states. We find that the bounds can be naturally expressed using the formalism of real space invariants (RSIs) that highlight the separation between electronic and atomic degrees of freedom. Finally, using exact Monte Carlo simulations on a model with perfectly flat bands and strictly local obstructed Wannier functions, we find that an attractive Hubbard interaction results in superconductivity as predicted by the RSI bound beyond mean field. Hence, obstructed bands are distinguished from trivial bands in the presence of interactions by the nonzero lower bound imposed on their superfluid weight.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Elevating density functional theory to chemical accuracy for water simulations through a density-corrected many-body formalism

Abstract Density functional theory (DFT) has been extensively used to model the properties of water. Albeit maintaining a good balance between accuracy and efficiency, no density functional has so far achieved the degree of accuracy necessary to correctly predict the properties of water across the entire phase diagram. Here, we present density-corrected SCAN (DC-SCAN) calculations for water which, minimizing density-driven errors, elevate the accuracy of the SCAN functional to that of “gold standard” coupled-cluster theory. Building upon the accuracy of DC-SCAN within a many-body formalism, we introduce a data-driven many-body potential energy function, MB-SCAN(DC), that quantitatively reproduces coupled cluster reference values for interaction, binding, and individual many-body energies of water clusters. Importantly, molecular dynamics simulations carried out with MB-SCAN(DC) also reproduce the properties of liquid water, which thus demonstrates that MB-SCAN(DC) is effectively the first DFT-based model that correctly describes water from the gas to the liquid phase.

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

Spin-filtering, giant magnetoresistance, rectifying and negative differential resistance effects in planar four-coordinate Fe complex with graphene nanoribbon electrodes

By using the nonequilibrium Green's function formalism combined with the density functional theory, we have investigated the spin-polarized transport properties of a planar four-coordinate Fe complex sandwiched between two zigzag-edge graphene nanoribbon (ZGNR) electrodes, where the ZGNRs are modulated by external magnetic field. The results show that the system can exhibit perfect dual spin-filtering and spin-rectifying effects at a wide bias range, giant magnetoresistance effect with large magnetoresistance ratio at small bias, and obvious negative differential resistance behavior. The mechanisms are proposed for these phenomena.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

MO-F-16A-04: Case Study: Estimation of Peak Skin Dose Following a Physician Reported “High Dose” Case and Sentinel Event Considerations

Purpose: This work offers as a teaching example a reported high dose fluoroscopy case and the workflow the institution followed to self-report a radiation overdose sentinel event to the Joint Commission. Methods: Following the completion of a clinical case in a hybrid OR room with a reported air kerma of >18 Gy at the Interventional Reference Point (IRP) the physicians involved in the case referred study to the institution's Radiation Safety Committee (RSC) for review. The RSC assigned a Diagnostic Medical Physicist (DMP) to estimate the patient's Peak Skin Dose (PSD) and analyze the case. Following the DMP's analysis and estimate of a PSD of >15 Gy the institution's adverse event committee was convened to discuss the case and to self-report the case as a radiation overdose sentinel event to the Joint Commission. The committee assigned a subgroup to perform the root cause analysis and develop institutional responses to the event. Results: The self-reporting of the sentinel event and the associated root cause analysis resulted in several institutional action items that are designed to improve process and safety. A formal reporting and analysis mechanism was adopted to review fluoroscopy cases with air kerma greater than 6 Gy at the IRP. An improved and formalized radiation safety training program for physicians using fluoroscopy equipment was implemented. Additionally efforts already under way to monitor radiation exposure in the Radiology department were expanded to include all fluoroscopy equipment capable of automated dose reporting. Conclusion: The adverse event review process and the root cause analysis following the self-reporting of the sentinel event resulted in policies and procedures that are expected to improve the quality and safe usage of fluoroscopy throughout the institution.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Verification and Validation of the Stakeholder Tool for Assessing Radioactive Transportation (START) – 22323

The U.S. Department of Energy (DOE) Office of Integrated Waste Management is planning for the eventual transportation, storage, and disposal of spent nuclear fuel (SNF) and high-level radioactive waste (HLW) from nuclear power plant and DOE sites. The Stakeholder Tool for Assessing Radioactive Transportation (START) is a web-based, geospatial decision-support tool developed for evaluating routing options and other aspects of transporting SNF and HLW, covering rail, truck, barge, and intermodal infrastructure and operations in the continental United States. The verification and validation (V&V) process is intended to independently assess START to provide confidence in the ability of START to accurately provide intended results. The V&V process checks the START tool using a variety of methods, ranging from independent hand calculations to comparison of START performance and results to those of other codes. The V&V activity was conducted independently from the START development team with opportunities to provide feedback and collaborate throughout the process. The V&V analyzed attributes of transportation routes produced by START, including route distance and both population and population density captured within buffer zones around routes. Population in the buffer zone, population density in the buffer zone, and route distance were all identified as crucial outputs of the START code and were subject to V&V tasks. Some of the improvements identified through the V&V process were standardizing the underlying population data in START, changing the projection of the population raster data, and changes to the methodology used for population density to improve its applicability for expected users. This collaboration also led to suggested improvements to some of the underlying shape file segments within START. A formal series of test routes went through the V&V process from START version 3.2.1 following the implementation of the recommended improvements to START. The V&V of the buffer zone population and distance reported in START were achieved by utilizing a custom workflow developed in QGIS. This platform was chosen because it is fully independent of START. The buffer zone population and route distances in START showed excellent agreement to independent V&V test results. Over 200 route test cases were run; in all cases, the percent difference in the population within the buffer zone was less than +/- 5%, and the majority of cases were below +/- 1%. These cases were also duplicated in ArcMap, which is Environmental Systems Research Institute’s (ESRI) GIS desktop application. The ArcMap results also show good agreement; however, they cannot be considered independent of START because START also uses ESRI software. The V&V tasks are currently addressing population density and estimated radiation dose in START as well as implementing further automation of the V&V process to allow larger suites of test routes to be run more easily and quickly by less experienced GIS users. These improvements to the V&V process will make it a valuable tool for the development team and the stakeholders so that each new iteration of START can easily be checked for quality and consistency. This paper and presentation will describe the V&V process, findings and implications, and ongoing V&V activities.

START, DOE, Transporation, dose, GIS↗
DOE OSTI2016

Second order gyrokinetic theory for particle-in-cell codes

The main idea of the gyrokinetic dynamical reduction consists in a systematical removal of the fast scale motion (the gyromotion) from the dynamics of the plasma, resulting in a considerable simplification and a significant gain of computational time. The gyrokinetic Maxwell–Vlasov equations are nowadays implemented in for modeling (both laboratory and astrophysical) strongly magnetized plasmas. Different versions of the reduced set of equations exist, depending on the construction of the gyrokinetic reduction procedure and the approximations performed in the derivation. The purpose of this article is to explicitly show the connection between the general second order gyrokinetic Maxwell–Vlasov system issued from the modern gyrokinetic theory and the model currently implemented in the global electromagnetic Particle-in-Cell code ORB5. Necessary information about the modern gyrokinetic formalism is given together with the consistent derivation of the gyrokinetic Maxwell–Vlasov equations from first principles. The variational formulation of the dynamics is used to obtain the corresponding energy conservation law, which in turn is used for the verification of energy conservation diagnostics currently implemented in ORB5. This work fits within the context of the code verification project VeriGyro currently run at IPP Max-Planck Institut in collaboration with others European institutions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Non-gaussian signatures of general inflationary trajectories

We carry out a numerical calculation of the bispectrum in generalised trajectories of canonical, single-field inflation. The trajectories are generated in the Hamilton-Jacobi (HJ) formalism based on Hubble Slow Roll (HSR) parameters. The calculation allows generally shape and scale dependent bispectra, or dimensionless f{sub NL}, in the out-of-slow-roll regime. The distributions of f{sub NL} for various shapes and HSR proposals are shown as an example of how this procedure can be used within the context of Monte Carlo exploration of inflationary trajectories. We also show how allowing out-of-slow-roll behaviour can lead to a bispectrum that is relatively large for equilateral shapes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Solving gyrokinetic systems with higher-order time dependence

We discuss theoretical and numerical aspects of gyrokinetics as a Lagrangian field theory when the field perturbation is introduced into the symplectic part. A consequence is that the field equations and particle equations of motion in general depend on the time derivatives of the field. The most well-known example is when the parallel vector potential is introduced as a perturbation, where a time derivative of the field arises only in the equations of motion, so an explicit equation for the fields may still be written. We will consider the conceptually more problematic case where the time-dependent fields appear in both the field equations and equations of motion, but where the additional term in the field equations is formally small. The conceptual issues were described by Burby (J. Plasma Phys., vol. 82 (3), 2016, 905820304): these terms lead to apparent additional degrees of freedom to the problem, so that the electric field now requires an initial condition, which is not required in low-frequency (Darwin) Vlasov–Maxwell equations. Also, the small terms in the Euler–Lagrange equations are a singular perturbation, and these two issues are interlinked. For well-behaved problems the apparent additional degrees of freedom are spurious, and the physically relevant solution may be directly identified. Because we needed to assume that the system is well behaved for small perturbations when deriving gyrokinetic theory, we must continue to assume that when solving it, and the physical solutions are thus the regular ones. The spurious nature of the singular degrees of freedom may also be seen by changing coordinate systems so the varying field appears only in the Hamiltonian. We then describe how methods appropriate for singular perturbation theory may be used to solve these asymptotic equations numerically. We then describe a proof-of-principle implementation of these methods for an electrostatic strong-flow gyrokinetic system; two basic test cases are presented to illustrate code functionality.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Multiparticle integral and differential correlation functions

This paper formalizes the use of integral and differential cumulants for measurements of multiparticle event-by-event transverse momentum fluctuations, rapidity fluctuations, as well as net-charge fluctuations. This enables the introduction of multiparticle balance functions, defined based on differential correlation functions (factorial cumulants), that suppress two- and three-prong resonance decays effects and enable measurements of underlying long-range correlations obeying quantum number conservation constraints. These multiparticle balance functions satisfy simple sum rules determined by quantum number conservation. It is additionally shown that these multiparticle balance functions arise as an intrinsic component of high-order net-charge cumulants. This implies that the magnitude of these cumulants, measured in a specific experimental acceptance, is strictly constrained by charge conservation and primarily determined by the rapidity and momentum width of these balance functions. The paper also presents techniques to reduce the computation time of differential correlation functions up to order n = 10 based on the methods of moments. Published by the American Physical Society 2024

Physics↗
DOE OSTI2026

Multi-Artifact Analysis of Self-Admitted Technical Debt in Scientific Software

Context: Self-admitted technical debt (SATD) occurs when developers acknowledge shortcuts in code. In scientific software (SSW), such debt poses unique risks to the validity and reproducibility of results. Objective: This study aims to identify, categorize, and evaluate scientific debt, a specialized form of SATD in SSW, and assess the extent to which traditional SATD categories capture these domain-specific issues. Method: We conduct a multi-artifact analysis across code comments, commit messages, pull requests, and issue trackers from 23 open-source SSW projects. We construct and validate a curated dataset of scientific debt, develop a multi-source SATD classifier to guide SATD management, and conduct a practitioner validation to assess the practical relevance of scientific debt. Results: Our classifier performs strongly across 900,358 artifacts from 23 SSW projects. SATD is most prevalent in pull requests and issue trackers, underscoring the value of multi-artifact analysis. Models trained on traditional SATD often miss scientific debt, emphasizing the need for its explicit detection in SSW. Practitioner validation confirmed that scientific debt is both recognizable and useful in practice. Conclusions: Scientific debt represents a unique form of SATD in SSW that that is not adequately captured by traditional categories and requires specialized identification and management. Our dataset, classification analysis, and practitioner validation results provide the first formal multi-artifact perspective on scientific debt, highlighting the need for tailored SATD detection approaches in SSW.

Melin, Eric [Boise State University]↗
DOE OSTI2015

Fermi orbital derivatives in self-interaction corrected density functional theory: Applications to closed shell atoms

A recent modification of the Perdew-Zunger self-interaction-correction to the density-functional formalism has provided a framework for explicitly restoring unitary invariance to the expression for the total energy. The formalism depends upon construction of Löwdin orthonormalized Fermi-orbitals which parametrically depend on variational quasi-classical electronic positions. Derivatives of these quasi-classical electronic positions, required for efficient minimization of the self-interaction corrected energy, are derived and tested, here, on atoms. Total energies and ionization energies in closed-shell singlet atoms, where correlation is less important, using the Perdew-Wang 1992 Local Density Approximation (PW92) functional, are in good agreement with experiment and non-relativistic quantum-Monte-Carlo results albeit slightly too low.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗