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 199 records · Page 11

A combined experimental and computational thermodynamic investigation of the U-Th-O system

Here, the thermodynamics of the U–Th–O system have been assessed using the Calphad method. The compound energy formalism (CEF) and a partially ionic two-sublattice liquid model (TSLM) were used for the fluorite U 1–y Th y O 2±x , γ-(U,Th) 4 O 9 , and the U–Th–O melt. The O 2 activity of fluorite U 1–y Th y O 2±x with temperature and composition was determined by thermogravimetric analysis. Thermodynamic studies for the Th–O binary and U–Th–O ternary available in the open literature were critically reviewed. A self-consistent data set was selected and compiled with the equilibrium oxygen pressures determined by thermogravimetry in order to optimize the adjustable parameters of models selected to represent the phases in the Th–O and U–Th–O systems.

36 MATERIALS SCIENCE↗

On the Enthalpy and Entropy of Point Defect Formation in Crystals

A standard way to determine the formation enthalpy H and entropy S of point defect formation in crystals consists in the application of the Arrhenius equation for the defect concentration. In this work, we show that a formal use of this method actually gives the effective (apparent) values of these quantities, which appear to be significantly overestimated. The underlying physical reason lies in temperature-dependent formation enthalpy of the defects, which is controlled by temperature dependence of the elastic moduli. We present an evaluation of the “true” H- and S-values for aluminum, which are derived on the basis of experimental data by taking into account temperature dependence of the formation enthalpy related to temperature dependence of the elastic moduli. The knowledge of the “true” activation parameters is needed for a correct calculation of the defect concentration constituting thus an issue of major importance for different fundamental and application issues of condensed matter physics and chemistry.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Ultrabroadband infrared near-field spectroscopy and imaging of local resonators in percolative gold films

Percolation processes are ubiquitous in nature and are responsible for many critical phenomena such as first-order phase transitions and infectious epidemic networks. The optical properties of a percolative medium can generally be captured by the effective medium approximation (EMA) when the degree of percolation and the properties of the constituent materials are properly addressed. However, the important local collective responses of nanoclusters in the deep subwavelength regime are often only phenomenologically addressed in the standard EMA formalism. Additionally, a comprehensive method that measures local light–matter interactions and registers how the local responses influence global optical properties has yet to be established on a firm basis. In this paper, we use infrared nanoimaging/spectroscopy to investigate percolative gold films in the vicinity of the critical percolation threshold. We demonstrate experimentally and theoretically that the near-field spectra yield quantitative information of the characteristic length scale of the local gold clusters and their relative oscillator strengths. As a result, EMA analysis can be augmented with near-field nano-spectroscopy to yield better predictability of the far-field reflection spectrum at the corresponding spectral range.

36 MATERIALS SCIENCE↗

Atomic and Molecular Physics in Strong Fields

In this project, we have investigated a broad range of atomic, molecular, and optical (AMO) processes of current significance in strong-field and ultrafast science and technology. The main goals are to develop new theoretical formalisms and accurate computational methods for ab initio and in-depth non-perturbative study of time-resolved multiphoton quantum dynamics and spectroscopy, very-high-order harmonic generation, optimal control of ultrafast attosecond pulse generation, and resonance-enhanced below threshold harmonics in intense ultrashort laser fields. We have worked on the following projects and we have published 20 top 15 % SCI journal articles during the last 3.5 years (Jan. 1, 2016 to June 30, 2019), including one article in top raking Physical Review Letters, and 3 articles in top ranking (Nature) Scientific Reports, and one Physical Review A article was chosen as Editors’ suggestion paper, etc. (a) Development of new time-dependent methods for accurate and efficient solution of multi-dimensional time-dependent Schrödinger equation and time-dependent density functional theory. (b) Non-perturbative treatment of very-high order (1000th – 5000th order) harmonic generation and broadband supercontinuum spectra from XUV to soft x-ray regime. (c) Development of synchrosqueezing transform (SST) for precision time-frequency analysis of electron multi-rescattering dynamics. (d) Exploration of the dynamical origin of resonance-enhanced near- and below-threshold harmonics. (e) Probing the high-order harmonics generation in disordered semiconductor. (f) Multiphoton ionization of one-electron relativistic diatomic quasi-molecules in strong laser fields. and, (g) New development of non-adiabatic density functional theory. The theoretical and simulation projects we have pursued are closely related to the ongoing experimental researches in several leading laboratories in the US and Europe, etc. It is expected that the outcome of this project will provide timely and significant new information and insight for the advancement of strong-field atomic and molecular physics and ultrafast science and technology.

74 ATOMIC AND MOLECULAR PHYSICS↗

Predicting magnetic properties of single-molecule magnets from self-interaction-free density-functional theory (Final Report)

In this project we investigated electronic structure of an intermediate-sized copper-based molecule and magnetic and hyperfine properties of several small non-magnetic and magnetic molecules including transition-metal elements by applying self-interaction corrections to density functional theory (DFT). In a sufficient number of cases, DFT-calculated electronic structure and magnetic properties of single-molecule magnets qualitatively differ from corresponding experimental data. This is partly due to self-interacting electrons within the DFT formalism. Recently, an efficient method to correct the self-interactions was proposed, i.e., Fermi-Lowdin prbital (FLO) based self-interaction corrected (SIC) methodology, within DFT. Henceforth, this method is referred to as FLO-SIC method which exists in FLOSIC code. We used this FLO-SIC method for our studies of electronic structure and magnetic and hyperfine properties of small magnetic molecules and non-magnetic molecules. Our study will provide insight into predictions of magnetic properties of single-molecule magnets where self-interaction corrections play a critical role. There are two components of this project. In the first work, we studied the electronic structure of a planar mononuclear Cu-based molecule in two oxidation states, using DFT with the FLO-SIC method. We chose this system because it is small enough and it includes a transition metal element. We found that the standard FLO-SIC method takes too much compute time even for the small transition-metal molecule and so we slightly modified the method in order to expedite the process. In the dianionic state, we found that the FLO-SIC spin density agrees quantitatively with accurate quantum chemistry methods, while DFT spin density without self-interaction corrections are severely deviated from the quantum chemistry methods. We also showed that the energy gap between the highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) of the dianionic state is larger than that of the monoanionic state. This result is consistent with experimental data. In the second work, we investigated how the interaction between the electron spin and nuclear spin is affected by electron self-interactions within small non-magnetic and magnetic molecules. Such an interaction is called hyperfine interaction. For molecules without significant orbital angular momentum, the hyperfine interaction consists of Fermi contact and dipolar interaction terms. Since the Fermi contact term depends on electron spin density at the nuclear site, it would be highly affected by self-interaction corrections. Therefore, we calculated the hyperfine interaction for the small molecules using DFT with the slightly modified expedited FLO-SIC method which was obtained in the first work, and compared the results to experimental data and DFT calculations without self-interaction corrections. We found significant improvement of the Fermi contact term computed using the FLO-SIC method for small magnetic molecules. Overall, the first and second work provided positive outlook of application of the FLO-SIC method to magnetic molecules and systems including transition-metal elements.

36 MATERIALS SCIENCE↗

Discrete Mathematical Approaches to Graph-Based Traffic Analysis

Modern cyber defense and anlaytics requires general, formal models of cyber systems. Multi-scale network models are prime candidates for such formalisms, using discrete mathematical methods based in hierarchically-structured directed multigraphs which also include rich sets of labels. An exemplar of an application of such an approach is traffic analysis, that is, observing and analyzing connections between clients, servers, hosts, and actors within IP networks, over time, to identify characteristic or suspicious patterns. Towards that end, NetFlow (or more generically, IPFLOW) data are available from routers and servers which summarize coherent groups of IP packets flowing through the network. In this paper, we consider traffic analysis of Netflow using both basic graph statistics and two new mathematical measures involving labeled degree distributions and time interval overlap measures. We do all of this over the VAST test data set of 96M synthetic Netflow graph edges, against which we can identify characteristic patterns of simulated ground-truth network attacks.

Joslyn, Cliff A.↗

Methods for analysing the quality of the element base of quantum information technologies

Methods for analysing the quality of the element base of quantum information technologies are considered. Methods for simulating quantum operations with allowance for quantum noise, based on the formalism of quantum operations (operator sum method and Choi – Jamiolkowski isomorphism), are described. Results of simulating some one- and two-qubit quantum transformations with allowance for the amplitude and phase relaxation and the presence of depolarising noise are reported. The quantum state tomography, based on the root approach and maximum likelihood method, is considered. The method of guaranteed estimation of quantum tomography accuracy, which implies quantitative estimation of measurement results, is described. Tomography of a single-qubit transformation on a superconducting quantum processor (IBM) is performed. The practical purpose of the study is to provide high quality and efficiency of quantum information technologies. The results obtained can be used to improve the quality of logic gates of quantum computers and debug the procedures of quantum state control in quantum cryptography problems. The bibliography includes 28 references. (quantum technology)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Nuclear-electronic orbital reduced explicitly correlated Hartree-Fock approach: Restricted basis sets and open-shell systems

The nuclear electronic orbital (NEO) reduced explicitly correlated Hartree-Fock (RXCHF) approach couples select electronic orbitals to the nuclear orbital via Gaussian-type geminal functions. This approach is extended to enable the use of a restricted basis set for the explicitly correlated electronic orbitals and an open-shell treatment for the other electronic orbitals. The working equations are derived and the implementation is discussed for both extensions. The RXCHF method with a restricted basis set is applied to HCN and FHF{sup −} and is shown to agree quantitatively with results from RXCHF calculations with a full basis set. The number of many-particle integrals that must be calculated for these two molecules is reduced by over an order of magnitude with essentially no loss in accuracy, and the reduction factor will increase substantially for larger systems. Typically, the computational cost of RXCHF calculations with restricted basis sets will scale in terms of the number of basis functions centered on the quantum nucleus and the covalently bonded neighbor(s). In addition, the RXCHF method with an odd number of electrons that are not explicitly correlated to the nuclear orbital is implemented using a restricted open-shell formalism for these electrons. This method is applied to HCN{sup +}, and the nuclear densities are in qualitative agreement with grid-based calculations. Future work will focus on the significance of nonadiabatic effects in molecular systems and the further enhancement of the NEO-RXCHF approach to accurately describe such effects.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Tensor optimized Fermi sphere method for nuclear matter—Power series correlated wave function and a cluster expansion

A new formalism, called “tensor optimized Fermi sphere (TOFS) method”, is developed to treat the nuclear matter using a bare interaction among nucleons. In this method, the correlated nuclear matter wave function is taken to be a power series type, Ψ{sub N}=[∑{sub n=0}{sup N}(1∕n!)F{sup n}]Φ{sub 0} and an exponential type, Ψ{sub ex}=exp(F)Φ{sub 0}, with the uncorrelated Fermi-gas wave function Φ{sub 0}, where the correlation operator F can induce central, spin–isospin, tensor, etc. correlations, and Ψ{sub ex} corresponds to a limiting case of Ψ{sub N} (N→∞). In the TOFS formalism based on Hermitian form, it is shown that the energy per particle in nuclear matter with Ψ{sub ex} can be expressed in terms of a linked-cluster expansion. On the basis of these results, we present the formula of the energy per particle in nuclear matter with Ψ{sub N}. We call the Nth-order TOFS calculation for evaluating the energy with Ψ{sub N}, where the correlation functions are optimally determined in the variation of the energy. The TOFS theory is applied for the study of symmetric nuclear matter using a central NN potential with short-range repulsion. The calculated results are fairly consistent to those of other theories such as the Brueckner–Hartree–Fock approach etc.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Mixed quantum–classical approach to model non-adiabatic electron–nuclear dynamics: Detailed balance and improved surface hopping method

We develop a density matrix formalism to describe coupled electron–nuclear dynamics. To this end, we introduce an effective Hamiltonian formalism that describes electronic transitions and small (quantum) nuclear fluctuations along a classical trajectory of the nuclei. Using this Hamiltonian, we derive equations of motion for the electronic occupation numbers and for the nuclear coordinates and momenta. We show that, in the limit, when the number of nuclear degrees of freedom coupled to a given electronic transition is sufficiently high (i.e., the strong decoherence limit), the equations of motion for the electronic occupation numbers become Markovian. Furthermore, the transition rates in these (rate) equations are asymmetric with respect to the lower-to-higher energy transitions and vice versa. In thermal equilibrium, such asymmetry corresponds to the detailed balance condition. We also study the equations for the electronic occupations in the non-Markovian regime and develop a surface hopping algorithm based on our formalism. To treat the decoherence effects, we introduce additional “virtual” nuclear wave packets whose interference with the “real” (physical) wave packets leads to the reduction in coupling between the electronic states (i.e., decoherence) as well as to the phase shifts that improve the accuracy of the numerical approach. Remarkably, the same phase shifts lead to the detailed balance condition in the strong decoherence limit.

74 ATOMIC AND MOLECULAR PHYSICS↗

Conservation laws in coupled cluster dynamics at finite temperature

We extend the finite-temperature Keldysh non-equilibrium coupled cluster theory (Keldysh-CC) [A. F. White and G. K.-L. Chan, J. Chem. Theory Comput. 15, 6137–6253 (2019)] to include a time-dependent orbital basis. When chosen to minimize the action, such a basis restores local and global conservation laws (Ehrenfest’s theorem) for all one-particle properties while remaining energy conserving for time independent Hamiltonians. We present the time-dependent Keldysh orbital-optimized coupled cluster doubles method in analogy with the formalism for zero-temperature dynamics, extended to finite temperatures through the time-dependent action on the Keldysh contour. To demonstrate the conservation property and understand the numerical performance of the method, we apply it to several problems of nonequilibrium finite-temperature dynamics: a 1D Hubbard model with a time-dependent Peierls phase, laser driving of molecular H2, driven dynamics in warm-dense silicon, and transport in the single impurity Anderson model.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Scattering Observables from One- and Two-body Densities: Formalism and Application to $\pmb \gamma $ ${}^3\hbox {He}$ Scattering

We introduce the transition-density formalism, an efficient and general method for calculating the interaction of external probes with light nuclei. One- and two-body transition densities that encode the nuclear structure of the target are evaluated once and stored. They are then convoluted with an interaction kernel to produce amplitudes, and hence observables. Here, by choosing different kernels, the same densities can be used for any reaction in which a probe interacts perturbatively with the target. The method therefore exploits the factorisation between nuclear structure and interaction kernel that occurs in such processes. We study in detail the convergence in the number of partial waves for matrix elements relevant in elastic Compton scattering on 3 He. The results are fully consistent with our previous calculations in Chiral Effective Field Theory. But the new approach is markedly more computationally efficient, which facilitates the inclusion of more partial-wave channels in the calculation. We also discuss the usefulness of the transition-density method for other nuclei and reactions. Calculations of elastic Compton scattering on heavier targets like 4 He are straightforward extensions of this study, since the same interaction kernels are used. And the generality of the formalism means that our 3 He densities can be used to evaluate any 3 He elastic-scattering observable with contributions from one- and two-body operators.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Random Phase Approximation Correlation Energy Using Real-Space Density Functional Perturbation Theory

We present a real-space method for computing the random phase approximation (RPA) correlation energy within Kohn–Sham density functional theory, leveraging the low-rank nature of the frequency-dependent density response operator. In particular, we employ a cubic-scaling formalism based on density functional perturbation theory that circumvents the calculation of the response function matrix, instead relying on the ability to compute its product with a vector through the solution of the associated Sternheimer linear systems. We develop a large-scale parallel implementation of this formalism using the subspace iteration method in conjunction with the spectral quadrature method while employing the Kronecker product-based method for the application of the Coulomb operator and the conjugate orthogonal conjugate gradient method for the solution of the linear systems. We demonstrate convergence with respect to key parameters and verify the method’s accuracy by comparing with plane-wave results. We show that the framework achieves good strong scaling to many thousands of processors, reducing the time to solution for a lithium hydride system with 128 electrons to around 150 s on 4608 processors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

TU-F-BRE-06: Flattening-Filter-Free Beam Quality Correction Factor Determination Using Experimental and Monte Carlo Methods

Purpose: To investigate beam quality correction factors for the flattening-filter-free (FFF) energies of the TrueBeam™ accelerator based on a dosimetry formalism for small and nonstandard fields. Methods: Three detectors - an Exradin W1 scintillator, Sun Nuclear EDGE diode, and LiF(Mg,Tl) TLD-100 chips - were investigated to determine their applicability as tools to measure quality correction factors for ionization chambers in the small and nonstandard fields of the TrueBeam™. Volume-averaging effects and energy dependence were observed in fields ranging from 1×1 to 40×40 cm{sup 2} for 6 MV and 10 MV beam energies using both FFF and flattened beam modes. Correction factors were determined for three ionization chambers: an Exradin A12 Farmer-type chamber, an Exradin A1SL scanning chamber, and an Exradin A26 reference-class microchamber. Beam quality corrections were also obtained using a benchmarked model of the TrueBeam™ created with the BEAMnrc user code of EGSnrc. Results: All three detectors demonstrated measureable energy dependence in the megavoltage range. The EDGE diode was deemed the most appropriate tool for beam quality correction factor measurements due to its low energy dependence and small size; however, alanine will be used in the future to reduce energy dependent effects even further. Measured k{sub Qmsr,Q} corrections of up to 4% were found for the 6MV FFF and 10 MV FFF beams, corresponding to a discrepancy of up to 3% compared to TG-51-determined dose. Up to a 10% k{sub Qclin,Qmsr} correction was measured for small fields referenced to a 10×10 cm{sup 2} field of the same energy. Much larger corrections were determined using the Monte Carlo model, and these discrepancies require further investigation. Conclusion: Progress has been made toward determining beam quality correction factors for the small and nonstandard fields of the TrueBeam™ accelerator. Further work must be done to ensure greater accuracy in patient treatments with this new modality.

07 ISOTOPE AND RADIATION SOURCES↗

Toward a resolution of the NN controversy

Lattice QCD calculations of two-nucleon interactions have been underway for about a decade, but still haven't reached the pion mass regime necessary for matching onto effective field theories and extrapolating to the physical point. Furthermore, results from different methods, including the use of the Lüscher formalism with different types of operators, as well as the HALQCD potential method, do not agree even qualitatively at very heavy pion mass. We investigate the role that different operators employed in the literature may play on the extraction of spectra for use within the Lüscher method. We first explore expectations from Effective Field Theory solved within a finite volume, for which the exact spectrum may be computed given different physical scenarios. We then present preliminary lattice QCD results for two-nucleon spectra calculated using different operators on a common lattice ensemble.

Nicholson, A↗

Stochastic generation of explicit pore structures by thresholding Gaussian random fields

We provide a description and computational investigation of an efficient method to stochastically generate realistic pore structures. Smolarkiewicz and Winter introduced this specific method in pores resolving simulation of Darcy flows (Smolarkiewicz and Winter, 2010 [1]) without giving a complete formal description or analysis of the method, or indicating how to control the parameterization of the ensemble. We address both issues in this paper. The method consists of two steps. First, a realization of a correlated Gaussian field, or topography, is produced by convolving a prescribed kernel with an initial field of independent, identically distributed random variables. The intrinsic length scales of the kernel determine the correlation structure of the topography. Next, a sample pore space is generated by applying a level threshold to the Gaussian field realization: points are assigned to the void phase or the solid phase depending on whether the topography over them is above or below the threshold. Hence, the topology and geometry of the pore space depend on the form of the kernel and the level threshold. Manipulating these two user prescribed quantities allows good control of pore space observables, in particular the Minkowski functionals. Extensions of the method to generate media with multiple pore structures and preferential flow directions are also discussed. To demonstrate its usefulness, the method is used to generate a pore space with physical and hydrological properties similar to a sample of Berea sandstone. -- Graphical abstract: -- Highlights: •An efficient method to stochastically generate realistic pore structures is provided. •Samples are generated by applying a level threshold to a Gaussian field realization. •Two user prescribed quantities determine the topology and geometry of the pore space. •Multiple pore structures and preferential flow directions can be produced. •A pore space based on Berea sandstone is generated.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Analytical gradient-based optimization of CALPHAD model parameters

The calibration of CALPHAD (CALculation of PHAse Diagrams) models involves the solution of a very challenging high-dimensional multiobjective optimization problem. Traditional approaches to parameter fitting predominantly rely on gradient-free methods, which while robust, are computationally inefficient and often scale poorly with model complexity. In this work, we introduce and demonstrate a generalizable framework for analytic gradient-based optimization of the parameters of the CALPHAD model enabled by the recently formalized Jansson derivative technique. This method allows for efficient evaluation of gradients of thermodynamic properties at equilibrium with respect to model parameters, even in the presence of arbitrarily complex internal degrees of freedom. Leveraging these semi-analytic gradients, we employ the conjugate gradient (CG) method to optimize thermodynamic model parameters for four binary alloy systems: Cu-Mg, Fe-Ni, Cr-Ni, and Cr-Fe. Across all systems, CG achieves comparable or superior optimality relative to Bayesian ensemble Markov Chain Monte Carlo (MCMC) with improvements in computational efficiency ranging from one to three orders of magnitude. Furthermore, our results establish a new paradigm for CALPHAD assessments in which high fidelity data-rich model calibration becomes tractable using deterministic gradient-informed algorithms.

CALPHAD↗