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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 865 records · Page 48

All sources NASA NTRS DOE OSTI

DOE OSTI2015

Many-body Green’s function theory for electron-phonon interactions: The Kadanoff-Baym approach to spectral properties of the Holstein dimer

We present a Kadanoff-Baym formalism to study time-dependent phenomena for systems of interacting electrons and phonons in the framework of many-body perturbation theory. The formalism takes correctly into account effects of the initial preparation of an equilibrium state and allows for an explicit time-dependence of both the electronic and phononic degrees of freedom. The method is applied to investigate the charge neutral and non-neutral excitation spectra of a homogeneous, two-site, two-electron Holstein model. This is an extension of a previous study of the ground state properties in the Hartree (H), partially self-consistent Born (Gd) and fully self-consistent Born (GD) approximations published in Säkkinen et al. [J. Chem. Phys. 143, 234101 (2015)]. Here, the homogeneous ground state solution is shown to become unstable for a sufficiently strong interaction while a symmetry-broken ground state solution is shown to be stable in the Hartree approximation. Signatures of this instability are observed for the partially self-consistent Born approximation but are not found for the fully self-consistent Born approximation. By understanding the stability properties, we are able to study the linear response regime by calculating the density-density response function by time-propagation. This amounts to a solution of the Bethe-Salpeter equation with a sophisticated kernel. The results indicate that none of the approximations is able to describe the response function during or beyond the bipolaronic crossover for the parameters investigated. Overall, we provide an extensive discussion on when the approximations are valid and how they fail to describe the studied exact properties of the chosen model system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

CORAL: A framework for rigorous self-validated data modeling and integrative, reproducible data analysis

Abstract Background Many organizations face challenges in managing and analyzing data, especially when relevant datasets arise from multiple sources and methods. Analyzing heterogeneous datasets and additional derived data requires rigorous tracking of their interrelationships and provenance. This task has long been a Grand Challenge of data science and has more recently been formalized in the FAIR principles: that all data objects be Findable, Accessible, Interoperable, and Reusable, both for machines and for people. Adherence to these principles is necessary for proper stewardship of information, for testing regulatory compliance, for measuring the efficiency of processes, and for facilitating reuse of data-analytical frameworks. Findings We present the Contextual Ontology-based Repository Analysis Library (CORAL), a platform that greatly facilitates adherence to all 4 of the FAIR principles, including the especially difficult challenge of making heterogeneous datasets Interoperable and Reusable across all parts of a large, long-lasting organization. To achieve this, CORAL's data model requires that data generators extensively document the context for all data, and our tools maintain that context throughout the entire analysis pipeline. CORAL also features a web interface for data generators to upload and explore data, as well as a Jupyter notebook interface for data analysts, both backed by a common API. Conclusions CORAL enables organizations to build FAIR data types on the fly as they are needed, avoiding the expense of bespoke data modeling. CORAL provides a uniquely powerful platform to enable integrative cross-dataset analyses, generating deeper insights than are possible using traditional analysis tools.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Imposing strong correlated energy spread on relativistic bunches with transverse deflecting cavities

We demonstrate a novel scheme for imposing or removing large energy chirp on relativistic electron bunches. This method relies on the transverse-to-longitudinal mixing accomplished with a set of transverse deflecting cavities. The working principles of the new concept are explained using the first-order matrix formalism in the approximation of the linear single-particle dynamics. The scheme demonstrates versatility regarding the average beam energy and, hence, positioning along the accelerator beam line. The performance of the scheme is numerically investigated with the elegant particle-tracking code. It is shown that the associated nonlinear effects causing emittance deterioration for the extreme quantities of the total energy spread can be effectively compensated by optimization of the Twiss parameters while relying on the eigenemittance analysis. The impact of the longitudinal space charge effects on the beam dynamics in the proposed scheme is also investigated.

43 PARTICLE ACCELERATORS↗

DOE OSTI2024

Convex optimization of contour deformations

We discuss various formal aspects of contour deformations used to alleviate sign problems; most importantly, relating these contour deformations to a certain convex optimization problem. As a consequence of this connection we describe a general method for proving upper bounds on the average phase achievable by the contour deformation method. Using this method we show that Abelian lattice Yang-Mills in two spacetime dimensions possesses, for many values of the complex coupling, an exponential sign problem that cannot be removed via any contour deformation. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2016

Poster - 11: Radiation barrier thickness calculations for the GammaPod

A consortium of radiotherapy centers in North America is in the process of evaluating a novel new {sup 60}Co teletherapy device, called the GammaPod™ (Xcision Medical Systems, Columbia Maryland), designed specifically for breast SBRT. The GammaPod consists of 36 collimated {sup 60}Co sources with a total activity of 4320 Ci. The sources are housed in a hemispherical source carrier that rotates during treatment to produce a cylindrically symmetric cone of primary beam spanning 16° – 54° degrees from the horizontal. This unique beam geometry presents challenges when designing or evaluating room shielding for the purposes of meeting regulatory requirements, and for ensuring the safety of staff and the public in surrounding areas. Conventional methods for calculating radiation barrier thicknesses have been adapted so that barrier transmission factors for the GammaPod can be determined from a few relevant distances and characteristics of the primary beam. Simple formalisms have been determined for estimating shielding requirements for primary radiation (with a rotating and non-rotating source carrier), patient-scattered radiation, and leakage radiation. When making worst case assumptions, it was found that conventional barrier thicknesses associated with linac treatment suites are sufficient for shielding all sources of radiation from the GammaPod.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2016

Meson Resonances from Lattice QCD

There has been recent, significant, advances in the determination of the meson spectrum of QCD. Current efforts have focused on the development and application of finite-volume formalisms that allow for the determination of scattering amplitudes as well as resonance behavior in coupled channel systems. I will review some of these recent developments, and demonstrate the viability of the method in meson systems.

Edwards, Robert G.↗

DOE OSTI2016

Meson resonances on the lattice

Edwards, Robert G.↗

DOE OSTI2017

A Comparison of Social and Economic Engagement Approaches: Lessons from Recent International Radioactive Waste Projects - 17478

The central challenge of any proposed radioactive waste management project is securing not only regulatory permissions, but also a social license to operate. 'Social license' generally refers to a community's acceptance or approval of a specific project. Oftentimes development of social license occurs outside of, or in parallel to, formal regulatory or permitting processes, and requires sustained commitment and proposed economic benefits to communities by proponents to acquire, build, and sustain social capital in the context of trust-based relationships. Historically, radioactive waste management projects have defaulted to demagoguery and 'siege warfare' modes of thinking and action by both proponents and opponents. While radioactive waste management projects are always inherently controversial, it is only recently that project owners and developers have experimented with alternative forms of social engagement outside of, or in addition to, formal regulatory processes. This paper compares and contrasts recent projects in Canada, Australia, the United Kingdom, and the United States, by examining distinct approaches to social engagement, with a focus on new methods that have been used to promote mutual understanding, foster rational debate, build confidence, promote economic development, and secure social license. The projects examined span the life-cycle of radioactive waste management development, from ones that are purely aspirational, to ones in early planning stages, to ones seeking a host community, to one that is focused on sustainment of existing social license. Several notable lessons and innovative practices are identified that will be of interest to radioactive waste management project owners and developers, non-governmental organizations, potential host communities, and political decision-makers for such projects. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2021

Large-momentum effective theory

Since the parton model was introduced by Feynman more than 50 years ago, much has been learned about the partonic structure of the proton through a large body of high-energy experimental data and dedicated global fits. However, limited progress has been made in calculating partonic observables such as the parton distribution function (PDFs) from the fundamental theory of strong interactions, quantum chromodynamics (QCD). Recently some advocated for a formalism, large-momentum effective theory (LaMET), through which one can extract parton physics from the properties of the proton traveling at a moderate boost factor such as γ ~ 2 – 5 . The key observation behind this approach is that Lorentz symmetry allows the standard formalism of partons in terms of light-front operators to be replaced by an equivalent one with large-momentum states and time-independent operators of a universality class. With LaMET, the PDFs, generalized PDFs or generalized parton distributions, transverse-momentum-dependent PDFs, and light-front wave functions can all be extracted in principle from lattice simulations of QCD (or other nonperturbative methods) through standard effective field theory matching and running. Future lattice QCD calculations with exascale computational facilities could help one to understand the experimental data related to the hadronic structure, including those from the upcoming electron-ion colliders dedicated to exploring the partonic landscape of the proton. In this work, the progress made in the past few years in the development of the LaMET formalism and its applications is reviewed, with an emphasis on a demonstration of its effectiveness from initial lattice QCD simulations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2014

Covariant formulation of spatially non-symmetric kinetic equilibria in magnetized astrophysical plasmas

Astrophysical plasmas in the surrounding of compact objects and subject to intense gravitational and electromagnetic fields are believed to give rise to relativistic regimes. Theoretical and observational evidences suggest that magnetized plasmas of this type are collisionless and can persist for long times (e.g., with respect to a distant observer, coordinate, time), while exhibiting geometrical structures characterized by the absence of well-defined spatial symmetries. In this paper, the problem is posed whether such configurations can correspond to some kind of kinetic equilibrium. The issue is addressed from a theoretical perspective in the framework of a covariant Vlasov statistical description, which relies on the method of invariants. For this purpose, a systematic covariant variational formulation of gyrokinetic theory is developed, which holds without requiring any symmetry condition on the background fields. As a result, an asymptotic representation of the relativistic particle magnetic moment is obtained from its formal exact solution, in terms of a suitably defined invariant series expansion parameter (perturbative representation). On such a basis, it is shown that spatially non-symmetric kinetic equilibria can actually be determined, an example being provided by Gaussian-like distributions. As an application, the physical mechanisms related to the occurrence of a non-vanishing equilibrium fluid 4-flow are investigated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2015

Coupling a nano-particle with isothermal fluctuating hydrodynamics: Coarse-graining from microscopic to mesoscopic dynamics

We derive a coarse-grained description of the dynamics of a nanoparticle immersed in an isothermal simple fluid by performing a systematic coarse graining of the underlying microscopic dynamics. As coarse-grained or relevant variables, we select the position of the nanoparticle and the total mass and momentum density field of the fluid, which are locally conserved slow variables because they are defined to include the contribution of the nanoparticle. The theory of coarse graining based on the Zwanzing projection operator leads us to a system of stochastic ordinary differential equations that are closed in the relevant variables. We demonstrate that our discrete coarse-grained equations are consistent with a Petrov-Galerkin finite-element discretization of a system of formal stochastic partial differential equations which resemble previously used phenomenological models based on fluctuating hydrodynamics. Key to this connection between our “bottom-up” and previous “top-down” approaches is the use of the same dual orthogonal set of linear basis functions familiar from finite element methods (FEMs), both as a way to coarse-grain the microscopic degrees of freedom and as a way to discretize the equations of fluctuating hydrodynamics. Another key ingredient is the use of a “linear for spiky” weak approximation which replaces microscopic “fields” with a linear FE interpolant inside expectation values. For the irreversible or dissipative dynamics, we approximate the constrained Green-Kubo expressions for the dissipation coefficients with their equilibrium averages. Under suitable approximations, we obtain closed approximations of the coarse-grained dynamics in a manner which gives them a clear physical interpretation and provides explicit microscopic expressions for all of the coefficients appearing in the closure. Our work leads to a model for dilute nanocolloidal suspensions that can be simulated effectively using feasibly short molecular dynamics simulations as input to a FEM fluctuating hydrodynamic solver.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Streaming Generalized Canonical Polyadic Tensor Decompositions

In this paper, we develop a method which we call OnlineGCP for computing the Generalized Canonical Polyadic (GCP) tensor decomposition of streaming data. GCP differs from traditional canonical polyadic (CP) tensor decompositions as it allows for arbitrary objective functions which the CP model attempts to minimize. This approach can provide better fits and more interpretable models when the observed tensor data is strongly non-Gaussian. In the streaming case, tensor data is gradually observed over time and the algorithm must incrementally update a GCP factorization with limited access to prior data. In this work, we extend the GCP formalism to the streaming context by deriving a GCP optimization problem to be solved as new tensor data is observed, formulate a tunable history term to balance reconstruction of recently observed data with data observed in the past, develop a scalable solution strategy based on segregated solves using stochastic gradient descent methods, describe a software implementation that provides performance and portability to contemporary CPU and GPU architectures and integrates with Matlab for enhanced usability, and demonstrate the utility and performance of the approach and software on several synthetic and real tensor data sets.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2017

Dispersion relation for hadronic light-by-light scattering: two-pion contributions

In our third paper of a series dedicated to a dispersive treatment of the hadronic light-by-light (HLbL) tensor, we derive a partial-wave formulation for two-pion intermediate states in the HLbL contribution to the anomalous magnetic moment of the muon (g - 2) μ , including a detailed discussion of the unitarity relation for arbitrary partial waves. We show that obtaining a final expression free from unphysical helicity partial waves is a subtle issue, which we thoroughly clarify. As a by-product, we obtain a set of sum rules that could be used to constrain future calculations of γ*γ* → ππ. We validate the formalism extensively using the pion-box contribution, defined by two-pion intermediate states with a pion-pole left-hand cut, and demonstrate how the full known result is reproduced when resumming the partial waves. Using dispersive fits to high-statistics data for the pion vector form factor, we provide an evaluation of the full pion box, a$π-box\atop{μ}$ =-15.9(2) × 10 -11 . As an application of the partial-wave formalism, we present a first calculation of ππ-rescattering effects in HLbL scattering, with γ*γ* → ππ helicity partial waves constructed dispersively using ππ phase shifts derived from the inverse-amplitude method. In this way, the isospin-0 part of our calculation can be interpreted as the contribution of the f0(500) to HLbL scattering in (g - 2) μ . We also argue that the contribution due to charged-pion rescattering implements corrections related to the corresponding pion polarizability and show that these are moderate. Our final result for the sum of pion-box contribution and its S-wave rescattering corrections reads a$π-box\atop{μ}$ + a$ππ, π-pole LHC\atop{μ, J=0}$ = -24(1) × 10 -11 .

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2021

Interactions of two and three mesons including higher partial waves from lattice QCD

We study two- and three-meson systems composed either of pions or kaons at maximal isospin using Monte Carlo simulations of lattice QCD. Utilizing the stochastic LapH method, we are able to determine hundreds of two- and three-particle energy levels, in nine different momentum frames, with high precision. We fit these levels using the relativistic finite-volume formalism based on a generic effective field theory in order to determine the parameters of the two- and three-particle K-matrices. We find that the statistical precision of our spectra is sufficient to probe not only the dominant s-wave interactions, but also those in d waves. In particular, we determine for the first time a term in the three-particle K-matrix that contains two-particle d waves. We use three N f = 2 + 1 CLS ensembles with pion masses of 200, 280, and 340 MeV. This allows us to study the chiral dependence of the scattering observables, and compare to the expectations of chiral perturbation theory.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2017

Testing quantum gravity through dumb holes

We propose a method to test the effects of quantum fluctuations on black holes by analyzing the effects of thermal fluctuations on dumb holes, the analogs for black holes. The proposal is based on the Jacobson formalism, where the Einstein field equations are viewed as thermodynamical relations, and so the quantum fluctuations are generated from the thermal fluctuations. It is well known that all approaches to quantum gravity generate logarithmic corrections to the entropy of a black hole and the coefficient of this term varies according to the different approaches to the quantum gravity. It is possible to demonstrate that such logarithmic terms are also generated from thermal fluctuations in dumb holes. In this paper, we claim that it is possible to experimentally test such corrections for dumb holes, and also obtain the correct coefficient for them. This fact can then be used to predict the effects of quantum fluctuations on realistic black holes, and so it can also be used, in principle, to experimentally test the different approaches to quantum gravity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2017

Monte Carlo chord length sampling for d-dimensional Markov binary mixtures

The Chord Length Sampling (CLS) algorithm is a powerful Monte Carlo method that models the effects of stochastic media on particle transport by generating on-the-fly the material interfaces seen by the random walkers during their trajectories. This annealed disorder approach, which formally consists of solving the approximate Levermore–Pomraning equations for linear particle transport, enables a considerable speed-up with respect to transport in quenched disorder, where ensemble-averaging of the Boltzmann equation with respect to all possible realizations is needed. However, CLS intrinsically neglects the correlations induced by the spatial disorder, so that the accuracy of the solutions obtained by using this algorithm must be carefully verified with respect to reference solutions based on quenched disorder realizations. When the disorder is described by Markov mixing statistics, such comparisons have been attempted so far only for one-dimensional geometries, of the rod or slab type. In this work we extend these results to Markov media in two-dimensional (extruded) and three-dimensional geometries, by revisiting the classical set of benchmark configurations originally proposed by Adams, Larsen and Pomraning and extended by Brantley. In particular, we examine the discrepancies between CLS and reference solutions for scalar particle flux and transmission/reflection coefficients as a function of the material properties of the benchmark specifications and of the system dimensionality.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

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↗