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,549 records · Page 86

Exchange-correlation contributions to thermodynamic properties of the homogeneous electron gas from a cumulant Green's function approach

Thermodynamic properties of the interacting homogeneous electron gas are calculated using a finite-temperature cumulant Green's function approach over a broad range of densities and temperatures up to the warm dense matter regime $T ~ T_F$, where $T_F$ is the Fermi degeneracy temperature. These properties can be separated into independent particle and exchange-correlation contributions, and our focus here is on the latter. In this work our approach is based on the Galitskii-Migdal-Koltun and electron number sum rules from the finite temperature many-body Green's function formalism, together with an extension of the cumulant Green's function to finite temperature. Previously this approach yielded exchange-correlation energies and potentials in good agreement with quantum Monte-Carlo calculations. Here the method is extended for various thermodynamic quantities including the chemical potential, total energy, Helmholtz free-energy, electronic equation of state, specific heat, and isothermal compressibility, which optionally include spin dependence. We find that the exchange-correlation contributions are weakly varying at low temperature but exhibit significant temperature dependence in the WDM regime, as well as a crossover from exchange- to correlation-dominated behavior. In contrast to the $T = 0^+$ limit, we also find that renormalization effects are largely but not completely suppressed at finite temperature. Comparisons with other approaches at various levels of approximation are also discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Learning Conservation Laws in Unknown Quantum Dynamics

We present a learning algorithm for discovering conservation laws given as sums of geometrically local observables in quantum dynamics. This includes conserved quantities that arise from local and global symmetries in closed and open quantum many-body systems. The algorithm combines the classical shadow formalism for estimating expectation values of observable and data analysis techniques based on singular value decompositions and robust polynomial interpolation to discover all such conservation laws in unknown quantum dynamics with rigorous performance guarantees. Our method can be directly realized in quantum experiments, which we illustrate with numerical simulations, using closed and open quantum system dynamics in a ℤ 2 gauge theory and in many-body localized spin chains.

Machine learning↗

SU-F-J-130: Margin Determination for Hypofractionated Partial Breast Irradiation

Purpose: To determine the Planning Target Volume (PTV) margin for Hypofractionated Partial Breast Irradiation (HPBI) using the van Herk formalism (M=2.5∑+0.7σ). HPBI is a novel technique intended to provide local control in breast cancer patients not eligible for surgical resection, using 40 Gy in 5 fractions prescribed to the gross disease. Methods: Setup uncertainties were quantified through retrospective analysis of cone-beam computed tomography (CBCT) data sets, collected prior to (prefraction) and after (postfraction) treatment delivery. During simulation and treatment, patients were immobilized using a wing board and an evacuated bag. Prefraction CBCT was rigidly registered to planning 4-dimensional computed tomography (4DCT) using the chest wall and tumor, and translational couch shifts were applied as needed. This clinical workflow was faithfully reproduced in Pinnacle (Philips Medical Systems) to yield residual setup and intrafractional error through translational shifts and rigid registrations (ribs and sternum) of prefraction CBCT to 4DCT and postfraction CBCT to prefraction CBCT, respectively. All ten patients included in this investigation were medically inoperable; the median age was 84 (range, 52–100) years. Results: Systematic (and random) setup uncertainties (in mm) detected for the left-right, craniocaudal and anteroposterior directions were 0.4 (1.5), 0.8 (1.8) and 0.4 (1.0); net uncertainty was determined to be 0.7 (1.5). Rotations >2° in any axis occurred on 8/72 (11.1%) registrations. Conclusion: Preliminary results suggest a non-uniform setup margin (in mm) of 2.2, 3.3 and 1.7 for the left-right, craniocaudal and anteroposterior directions is required for HPBI, given its immobilization techniques and online setup verification protocol. This investigation is ongoing, though published results from similar studies are consistent with the above findings. Determination of margins in breast radiotherapy is a paradigm shift, but a necessary step in moving towards hypofractionated regiments, which may ultimately redefine the standard of care for this select patient population.

60 APPLIED LIFE SCIENCES↗

The Effect of the Polarization Characteristics of Probe Light on the Signal of Optically Detected Magnetic Resonance in Magnetometric and Gyroscopic Quantum Sensors

We consider the effect of the polarization characteristics of probe light on the signal of optically detected magnetic resonance in quantum sensors, including quantum magnetometers based on the phenomenon of electron paramagnetic resonance and quantum gyroscopes employing both the electron and nuclear magnetic resonance. Relationships between the magnetic resonance signal magnitude and parameters of the optical system elements, which are based on the Stokes–Mueller formalism, are derived and verified. It is found that the main destructive influence in the signal in a standard two-beam scheme is produced by phase delays introduced by both metallic and dielectric mirrors. Methods for compensation of this destructive influence are proposed and verified.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Multi-phenomenology Yield Characterization

This report serves as the first delivery of a four-year applied science effort to transform and advance the error bounds for the yield estimate of an explosion. Each year’s delivery will be in this form, culminating in the submission of this work for peer review to a scientific journal. Importantly, the yearly progress reports can then also be viewed as expanding drafts working towards a formal journal article submission. For the first tranche of funding, we collaborated with Air Force Technical Applications Center (AFTAC) scientists to identify unclassified real-world data that demonstrate and validate our advanced error propagation methods. Collaboration includes visits to AFTAC and telecons. For this development, we illustrate the fusion of seismic, acoustic, optical, and surface effect signatures from an explosion. The mathematics and code being adapted to this specific application (Williams et al., 2021) involves physics models of multiple sensor signatures. We have also identified related physics models and have integrated them into code. Current methods of underground explosion yield estimation for the Threshold Test Ban Treaty (TTBT) have served the US treaty monitoring mission well for decades. A research objective of the Defense Nuclear Nonproliferation Research and Development (DNN R&D) office of the National Nuclear Security Administration (NNSA) has always been to provide new technical capabilities for monitoring lower thresholds. The general model and error propagation code to be developed in this project is based on significant advances in error modeling and propagation needed to analyze data at lower detection thresholds. The second tranche of funding for this project began on May 1, 2022, and planned work for the second tranche includes: i) completing the integration of physical model code into the general error model framework; this code accommodates a wide range of linear/nonlinear source models, fixed/ random effects, and frequentist/Bayesian analyses (the purpose of which is not to dictate to users how to analyze data, but instead to allow users the maximum flexibility in their work); ii) illustrative application of code to identified data, and; iii) initial planning with AFTAC researchers on delivery of code to the Common Development Environment at AFTAC, and continued writing of the planned final journal article submission (year two of this progress report), with particular emphasis on descriptions of data identified for this effort.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Semi-classical Kinetic Theory for Massive Spin-half Fermions with Leading-order Spin Effects

We consider the quantum kinetic-theory description for interacting massive spin-half fermions using the Wigner function formalism. We derive a general kinetic theory description assuming that the spin effects appear at the classical and quantum level. To track the effect of such different contributions we use the semi-classical expansion method to obtain the generalized dynamical equations including spin, analogous to classical Boltzmann equation. This approach can be used to obtain a collision kernel involving local as well as non-local collisions among the microscopic constituent of the system and eventually, a framework of spin hydrodynamics ensuring the conservation of the energy-momentum tensor and total angular momentum tensor.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

RF surface resistance of YBa2Cu3O(7-x) thin films

The excitement engendered by the discovery of the new T sub c oxide superconductors has led to much speculation about practical applications of thin films of these materials in digital and analog electronic devices. Most of these envisioned applications involve high frequency signals for which a detailed knowledge of the surface impedance of the novel superconductors is very important. We have measured the surface resistance of thin films of YBaCuO in the frequency range 0.5 less than f less than 17 GHz using a stripline-resonator method. The stripline procedure also was used to measure the surface resistance of high quality gold and aluminum films; the resistance values obtained agree with values predicted from the measured dc resistance using the Pippard formalism for the anomalous skin effect. The YBaCuO were produced by a multilayer deposition process. The films are formed by e-beam evaporation of 24 layers of Y, Ba, and Cu. Films with the highest transition temperature were obtained using yttria-stabilized zirconia (YSZ) substrates. After deposition, the films are transferred to a furnace where they are annealed in flowing O2 at 850 C for 2h. The furnace then is turned off and allowed to cool to 100 C in about 16 h. Auger profiling of the films made by this process shows that the concentrations of Y, Ba, Cu, and O are uniform to within 1 percent throughout the thickness of the film.

Source record↗

Comparison of Several Numerical Methods for Simulation of Compressible Shear Layers

An investigation is conducted on several numerical schemes for use in the computation of two-dimensional, spatially evolving, laminar variable-density compressible shear layers. Schemes with various temporal accuracies and arbitrary spatial accuracy for both inviscid and viscous terms are presented and analyzed. All integration schemes use explicit or compact finite-difference derivative operators. Three classes of schemes are considered: an extension of MacCormack's original second-order temporally accurate method, a new third-order variant of the schemes proposed by Rusanov and by Kutier, Lomax, and Warming (RKLW), and third- and fourth-order Runge-Kutta schemes. In each scheme, stability and formal accuracy are considered for the interior operators on the convection-diffusion equation U(sub t) + aU(sub x) = alpha U(sub xx). Accuracy is also verified on the nonlinear problem, U(sub t) + F(sub x) = 0. Numerical treatments of various orders of accuracy are chosen and evaluated for asymptotic stability. Formally accurate boundary conditions are derived for several sixth- and eighth-order central-difference schemes. Damping of high wave-number data is accomplished with explicit filters of arbitrary order. Several schemes are used to compute variable-density compressible shear layers, where regions of large gradients exist.

Kennedy, Christopher A.↗

On the Use of Kronecker Operators for the Solution of Generalized Stochastic Petri Nets

We discuss how to describe the Markov chain underlying a generalized stochastic Petri net using Kronecker operators on smaller matrices. We extend previous approaches by allowing both an extensive type of marking-dependent behavior for the transitions and the presence of immediate synchronizations. The derivation of the results is thoroughly formalized, including the use of Kronecker operators in the treatment of the vanishing markings and the computation of impulse-based reward measures. We use our techniques to analyze a model whose solution using conventional methods would fail because of the state-space explosion. In the conclusion, we point out ideas to parallelize our approach.

Ciardo, Gianfranco↗

GPU Lossless Hyperspectral Data Compression System for Space Applications

On-board lossless hyperspectral data compression reduces data volume in order to meet NASA and DoD limited downlink capabilities. At JPL, a novel, adaptive and predictive technique for lossless compression of hyperspectral data, named the Fast Lossless (FL) algorithm, was recently developed. This technique uses an adaptive filtering method and achieves state-of-the-art performance in both compression effectiveness and low complexity. Because of its outstanding performance and suitability for real-time onboard hardware implementation, the FL compressor is being formalized as the emerging CCSDS Standard for Lossless Multispectral & Hyperspectral image compression. The FL compressor is well-suited for parallel hardware implementation. A GPU hardware implementation was developed for FL targeting the current state-of-the-art GPUs from NVIDIA(Trademark). The GPU implementation on a NVIDIA(Trademark) GeForce(Trademark) GTX 580 achieves a throughput performance of 583.08 Mbits/sec (44.85 MSamples/sec) and an acceleration of at least 6 times a software implementation running on a 3.47 GHz single core Intel(Trademark) Xeon(Trademark) processor. This paper describes the design and implementation of the FL algorithm on the GPU. The massively parallel implementation will provide in the future a fast and practical real-time solution for airborne and space applications.

Graphic Processor Units↗

Sufficient Conditions for the Equivalent Convergence of Sequences of Different Approximants for Two-Dimensional Continued Fractions

By the method of fundamental inequalities, we study the equivalent convergence of two figured approximants for two-dimensional continued fractions, namely, the approximants obtained from the problem of correspondence of two-dimensional continued fractions to some formal double power series and the problem of equivalence of two-dimensional continued fractions. We also establish sufficient conditions under which the indicated two approximants of two-dimensional continued fractions converge to the same limit.

97 MATHEMATICS AND COMPUTING↗

Assessing the impact of substrate-level enzyme regulations limiting ethanol titer in Clostridium thermocellum using a core kinetic model

Clostridium thermocellum is a promising candidate for consolidated bioprocessing because it can directly ferment cellulose to ethanol. Despite significant efforts, achieved yields and titers fall below industrially relevant targets. This implies that there still exist unknown enzymatic, regulatory, and/or possibly thermodynamic bottlenecks that can throttle back metabolic flow. By (i) elucidating internal metabolic fluxes in wild-type C. thermocellum grown on cellobiose via 13 C-metabolic flux analysis ( 13 C-MFA), (ii) parameterizing a core kinetic model, and (iii) subsequently deploying an ensemble-docking workflow for discovering substrate-level regulations, this paper aims to reveal some of these factors and expand our knowledgebase governing C. thermocellum metabolism. Generated 13 C labeling data were used with 13 C-MFA to generate a wild-type flux distribution for the metabolic network. Notably, flux elucidation through MFA alluded to serine generation via the mercaptopyruvate pathway. Using the elucidated flux distributions in conjunction with batch fermentation process yield data for various mutant strains, we constructed a kinetic model of C. thermocellum core metabolism (i.e. k-ctherm138). Subsequently, we used the parameterized kinetic model to explore the effect of removing substrate-level regulations on ethanol yield and titer. Upon exploring all possible simultaneous (up to four) regulation removals we identified combinations that lead to many-fold model predicted improvement in ethanol titer. In addition, by coupling a systematic method for identifying putative competitive inhibitory mechanisms using K-FIT kinetic parameterization with the ensemble-docking workflow, we flagged 67 putative substrate-level inhibition mechanisms across central carbon metabolism supported by both kinetic formalism and docking analysis.

59 BASIC BIOLOGICAL SCIENCES↗

Polarizable Frozen Density Embedding with External Orthogonalization

In this work, we report a polarizable subsystem density functional theory to describe electronic properties of molecules embedded on a metal cluster. Interaction between the molecule and metal cluster is described using frozen density embedding (FDE). Substituting the nonadditive kinetic potential (NAKP) by approximate functionals is circumvented by enforcing external orthogonality (EO) through a projection operator. The computationally expensive freeze/thaw (FT) cycles are bypassed by including a polarization term in the embedding operator. Furthermore, the combination of polarization and EO permits supermolecular basis set calculations, which was not possible for strongly interacting systems with existing kinetic energy functionals. To test the method, we described the ground state density of pyridine, water, and benzene on a silver cluster. Performing FT on top of EO results in exact density embedding for this category of systems and is thus used for benchmarking the method. We find that the density is reproduced to within 0.15e, and the dipole and quadrupole moments are within 18% of the reference points for subsystem separations ranging from bonding to noninteracting distances. Additionally, our formalism allows the flexibility of incorporating different density functionals to the molecular and the metallic subsystems reducing the overall computational cost.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Scrutinizing formally Ni IV centers through the lenses of core spectroscopy, molecular orbital theory, and valence bond theory

Nickel K- and L 2,3 -edge X-ray absorption spectra (XAS) are discussed for 16 complexes and complex ions with nickel centers spanning a range of formal oxidation states from II to IV. K-edge XAS alone is shown to be an ambiguous metric of physical oxidation state for these Ni complexes. Meanwhile, L 2,3 -edge XAS reveals that the physical d-counts of the formally Ni IV compounds measured lie well above the d 6 count implied by the oxidation state formalism. The generality of this phenomenon is explored computationally by scrutinizing 8 additional complexes. The extreme case of NiF 6 2- is considered using high-level molecular orbital approaches as well as advanced valence bond methods. The emergent electronic structure picture reveals that even highly electronegative F-donors are incapable of supporting a physical d 6 Ni IV center. The reactivity of Ni IV complexes is then discussed, highlighting the dominant role of the ligands in this chemistry over that of the metal centers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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 (↗

Entanglement entropy of electronic excitations

A new perspective into correlation effects in electronically excited states is provided through quantum information theory. The entanglement between the electron and hole quasiparticles is examined, and it is shown that the related entanglement entropy can be computed from the eigenvalue spectrum of the well-known natural transition orbital (NTO) decomposition. Non-vanishing entanglement is obtained whenever more than one NTO pair is involved, i.e., in the case of a multiconfigurational or collective excitation. An important implication is that in the case of entanglement it is not possible to gain a complete description of the state character from the orbitals alone, but more specific analysis methods are required to decode the mutual information between the electron and hole. Moreover, the newly introduced number of entangled states is an important property by itself giving information about excitonic structure. The utility of the formalism is illustrated in the cases of the excited states of two interacting ethylene molecules, the conjugated polymer para-phenylene vinylene, and the naphthalene molecule.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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↗