Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “potential theory”

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 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Efficient Parameterization of Density Functional Tight-Binding for 5 f -Elements: A Th–O Case Study

Density functional tight binding (DFTB) models for f-element species are challenging to parametrize owing to the large number of adjustable parameters. The explicit optimization of the terms entering the semiempirical DFTB Hamiltonian related to f orbitals is crucial to generating a reliable parametrization for f-block elements, because they play import roles in bonding interactions. However, since the number of parameters grows quadratically with the number of orbitals, the computational cost for parameter optimization is much more expensive for the f-elements than for the main group elements. In this work we present a set of efficient approaches for mitigating the hurdle imposed by the large size of the parameter space. A novel group-by-orbital correction functions for two-center bond integrals was developed. With this approach the number of parameters is reduced, and it grows linearly with the number of elements, maintaining the accuracy and the number of parameters, in the case of f elements, by more than 40%. The parameter optimization step was accelerated by means of the mini-batch BFGS method. This method allows parameter optimizations with much larger training sets than other single batch methods. A stochastic optimizer was employed that helped overcome shallow local minima in the objective function. The proposed algorithm was used to parametrize the DFTB Hamiltonian for the Th–O system, which was subsequently applied to the study of ThO 2 nanoparticles. The training set consisted of 6322 unique structures, which is barely feasible with conventional optimization methods. The optimized parameter set, LANL-ThO, displays good agreement with DFT-calculated properties such as energies, forces, and structures for both clusters and bulk ThO 2 . Benefiting from the fewer number of parameters and lower computational costs for objective function evaluations, this new approach shows its potential applications in DFTB parametrization for elements with high angular momentum, which present a challenge to conventional methods.

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

Augmenting Molecular Graphs with Geometries via Machine Learning Interatomic Potentials

Accurate molecular property predictions require 3D geometries, which are typically obtained using expensive methods such as density functional theory (DFT). Here, we attempt to obtain molecular geometries by relying solely on machine learning interatomic potential (MLIP) models. To this end, we first curate a large-scale molecular relaxation dataset comprising 3.5 million molecules and 300 million snapshots. Then MLIP pre-trained models are trained with supervised learning to predict energy and forces given 3D molecular structures. Once trained, we show that the pre-trained models can be used in different ways to obtain geometries either explicitly or implicitly. First, it can be used to obtain approximate low-energy 3D geometries via geometry optimization. While these geometries do not consistently reach DFT-level chemical accuracy or convergence, they can still improve downstream performance compared to non-relaxed structures. To mitigate potential biases and enhance downstream predictions, we introduce geometry fine-tuning based on the relaxed 3D geometries. Second, the pre-trained models can be directly fine-tuned for property prediction when ground truth 3D geometries are available. Our results demonstrate that MLIP pre-trained models trained on relaxation data can learn transferable molecular representations to improve downstream molecular property prediction and can provide practically valuable but approximate molecular geometries that benefit property predictions. Our code is publicly available at: https://github.com/divelab/AIRS/.

Fu, Cong [Texas A & M Univ., College Station, TX (↗
DOE OSTI2020

Pilgrim : A thermal rate constant calculator and a chemical kinetics simulator

Pilgrim is a program written in Python and designed to use direct dynamics in the calculation of thermal rate constants of chemical reactions by the variational transition state theory (VTST), based on electronic structure calculations for the potential energy surface. Pilgrim can also simulate reaction mechanisms using kinetic Monte Carlo (KMC). For reaction processes with many elementary steps, the rate constant of each of these steps can be calculated by means of conventional transition state theory (TST) or by using VTST. In the current version, Pilgrim can evaluate thermal rates using the canonical version of reaction-path VTST, which requires the calculation of the minimum energy path (MEP) associated with each elementary step or transition structure. Multi-dimensional quantum effects can be incorporated through the small-curvature tunneling (SCT) approximation. These methodologies are available both for reactions involving a single structure of the reactants and the transition state and also for reactions involving flexible molecules with multiple conformations of the reactant and/or of the transition state. For systems with many conformers, the program can evaluate each of the elementary reaction rate constants by multipath canonical VTST or multi-structural VTST. Moreover, the reactant can be unimolecular or bimolecular. Torsional anharmonicity can be incorporated through either the MSTor or the Q2DTor programs. Dual-level calculations are also available in Pilgrim: automatic high-level single-point energies can be used to correct the energy of reactants, transition states, products, and MEP points using the interpolated single-point energies (ISPE) algorithm. When the rate constants of all the chemical processes of interest are known, by means of their calculation using Pilgrim or alternatively through analytical fits to the rate constants as functions of temperature, it is possible to simulate a multistep mechanism under specified laboratory conditions using KMC. Finally, this algorithm allows performing a kinetic simulation to monitor the evolution of each chemical species with time and obtain the product yields.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Emergent Recursive Multiscale Interaction in Complex Systems

This project studied the potential for multiscale group dynamics in complex social systems, including emergent recursive interaction. Current social theory on group formation and interaction focuses on a single scale (individuals forming groups) and is largely qualitative in its explanation of mechanisms. We combined theory, modeling, and data analysis to find evidence that these multiscale phenomena exist, and to investigate their potential consequences and develop predictive capabilities. In this report, we discuss the results of data analysis showing that some group dynamics theory holds at multiple scales. We introduce a new theory on communicative vibration that uses social network dynamics to predict group life cycle events. We discuss a model of behavioral responses to the COVID-19 pandemic that incorporates influence and social pressures. Finally, we discuss a set of modeling techniques that can be used to simulate multiscale group phenomena.

42 ENGINEERING↗
DOE OSTI2020

EFT triangles in the same-sign WW scattering process at the HL-LHC and HE-LHC

We investigate the Beyond Standard Model discovery potential in the framework of the Effective Field Theory (EFT) for the same-sign $WW$ scattering process in purely leptonic $W$ decay modes at the High-Luminosity and High-Energy phases of the Large Hadron Collider (LHC). The commonly applied procedure of studying the $WWWW$ quartic coupling by testing dimension-8 operators of the EFT one at a time is used in this work. In the considered process there is no experimental handle on the $WW$ invariant mass, and it has previously been shown that the discovery potential at 14 TeV is rather slim. In this paper we report the results calculated for a 27 TeV machine and compare them with the discovery potential obtained at 14 TeV. We find that while the respective discovery regions shift to lower values of the Wilson coefficients, the overall discovery potential of this procedure does not get significantly larger with a higher beam energy.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Low-energy physics in neutrino LArTPCs

Here, in this paper, we review scientific opportunities and challenges related to detection and reconstruction of low-energy (less than 100 MeV) signatures in liquid argon time-projection chamber (LArTPC) neutrino detectors. LArTPC neutrino detectors designed for performing precise long-baseline oscillation measurements with GeV-scale accelerator neutrino beams also have unique sensitivity to a range of physics and astrophysics signatures via detection of event features at and below the few tens of MeV range. In addition, low-energy signatures are an integral part of GeV-scale accelerator neutrino interaction final-states, and their reconstruction can enhance the oscillation physics sensitivities of LArTPC experiments. New physics signals from accelerator and natural sources also generate diverse signatures in the low-energy range, and reconstruction of these signatures can increase the breadth of Beyond the Standard Model scenarios accessible in LArTPC-based searches. A variety of experimental and theory-related challenges remain to realizing this full range of potential benefits. Neutrino interaction cross-sections and other nuclear physics processes in argon relevant to sub-hundred-MeV LArTPC signatures are poorly understood, and improved theory and experimental measurements are needed; pion decay-at-rest sources and charged particle and neutron test beams are ideal facilities for improving this understanding. There are specific calibration needs in the low-energy range, as well as specific needs for control and understanding of radiological and cosmogenic backgrounds. Low-energy signatures, whether steady-state or part of a supernova burst or larger GeV-scale event topology, have specific triggering, DAQ and reconstruction requirements that must be addressed outside the scope of conventional GeV-scale data collection and analysis pathways. Novel concepts for future LArTPC technology that enhance low-energy capabilities should also be explored to help address these challenges.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

A static quantum embedding scheme based on coupled cluster theory

Here, we develop a static quantum embedding scheme that utilizes different levels of approximations to coupled cluster (CC) theory for an active fragment region and its environment. To reduce the computational cost, we solve the local fragment problem using a high-level CC method and address the environment problem with a lower-level Møller–Plesset (MP) perturbative method. This embedding approach inherits many conceptual developments from the hybrid second-order Møller–Plesset (MP2) and CC works by Nooijen [J. Chem. Phys. 111, 10815 (1999)] and Bochevarov and Sherrill [J. Chem. Phys. 122, 234110 (2005)]. We go beyond those works here by primarily targeting a specific localized fragment of a molecule and also introducing an alternative mechanism to relax the environment within this framework. We will call this approach MP-CC. We demonstrate the effectiveness of MP-CC on several potential energy curves and a set of thermochemical reaction energies, using CC with singles and doubles as the fragment solver, and MP2-like treatments of the environment. The results are substantially improved by the inclusion of orbital relaxation in the environment. Using localized bonds as the active fragment, we also report results for N=N bond breaking in azomethane and for the central C–C bond torsion in butadiene. We find that when the fragment Hilbert space size remains fixed (e.g., when determined by an intrinsic atomic orbital approach), the method achieves comparable accuracy with both a small and a large basis set. Additionally, our results indicate that increasing the fragment Hilbert space size systematically enhances the accuracy of observables, approaching the precision of the full CC solver.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1982

Application of a transonic potential flow code to the static aeroelastic analysis of three-dimensional wings

Since the aerodynamic theory is nonlinear, the method requires the coupling of two iterative processes - an aerodynamic analysis and a structural analysis. A full potential analysis code, FLO22, is combined with a linear structural analysis to yield aerodynamic load distributions on and deflections of elastic wings. This method was used to analyze an aeroelastically-scaled wind tunnel model of a proposed executive-jet transport wing and an aeroelastic research wing. The results are compared with the corresponding rigid-wing analyses, and some effects of elasticity on the aerodynamic loading are noted.

Whitlow, W., Jr.↗
DOE OSTI2023

Hybrid interatomic potential for Sn

To design materials for extreme applications, it is important to understand and predict phase transitions and their influence on material properties under high pressures and temperatures. Atomistic modeling can be a useful tool to assess these behaviors. However, this can be difficult due to the lack of fidelity of the interatomic potentials in reproducing this high pressure and temperature extreme behavior. Here, in this work, a hybrid EAM-R—which is the combination of embedded atom method (EAM) and rapid artificial neural network potential—for Tin (Sn) is described which is capable of accurately modeling the complex sequence of phase transitions between different metallic polymorphs as a function of pressure. This hybrid approach ensures that a basic empirical potential like EAM is used as a lower energy bound. By using the final activation function, the neural network contribution to energy must be positive, assuring stability over the whole configuration space. This implementation has the capacity to reproduce density functional theory results at 6 orders of magnitude slower than a pair potential for molecular dynamics simulation, including elastic and plastic characteristics and relative energies of each phase. Using calculations of the Gibbs free energy, it is demonstrated that the potential precisely predicts the experimentally observed phase changes at temperatures and pressures across the whole phase diagram. At 10.2 GPa, the present potential predicts a first-order phase transition between body-centered tetragonal (BCT) β-Sn and another polymorph of BCT-Sn. This structure transforms into body-centered cubic near the experimentally reported value at 33 GPa. Thus, the Sn potential developed in this paper can be used to study complex deformation mechanisms under extreme conditions of high pressure and strain rates unlike existing potentials. Moreover, the framework developed in this paper can be extended for different material systems with complex phase diagrams.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A Multicontinuum-Theory-Based Approach to the Analysis of Fiber-Reinforced Polymer Composites with Degraded Stiffness and Strength Properties Due to Moisture Absorption

Marine energy generation technologies such as wave and tidal power have great potential in meeting the need for renewable energy in the years ahead. Yet, many challenges remain associated with marine-based systems because of the corrosive environment. Conventional materials like metals are subject to rapid corrosive breakdown, crippling the lifespan of structures in such environments. Fiber-reinforced polymer composites offer an appealing alternative in their strength and corrosion resistance, but can experience degradation of mechanical properties as a result of moisture absorption. An investigation is conducted to test the application of a technique for micromechanical analysis of composites, known as multicontinuum theory and demonstrated in past works, as a mechanism for predicting the effects of prolonged moisture absorption on the performance of fiber-reinforced composites. Experimental tensile tests are performed on composite coupons with and without prolonged exposure to a salt water solution to obtain stiffness and strength properties. Multicontinuum theory is applied in conjunction with micromechanical modeling to deduce the effects of moisture absorption on the behavior of constituent materials within the composites. The results are consistent with experimental observations when guided by known mechanisms and trends from previous studies, indicating multicontinuum theory as a potentially effective tool in predicting the long-term performance of composites in marine environments.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Linear and nonlinear optical responses in the chiral multifold semimetal RhSi

Chiral topological semimetals are materials that break both inversion and mirror symmetries. They host interesting phenomena such as the quantized circular photogalvanic effect (CPGE) and the chiral magnetic effect. In this work, we report a comprehensive theoretical and experimental analysis of the linear and nonlinear optical responses of the chiral topological semimetal RhSi, which is known to host multifold fermions. We show that the characteristic features of the optical conductivity, which display two distinct quasi-linear regimes above and below 0.4 eV, can be linked to excitations of different kinds of multifold fermions. The characteristic features of the CPGE, which displays a sign change at 0.4 eV and a large non-quantized response peak of around 160 μA/V 2 at 0.7 eV, are explained by assuming that the chemical potential crosses a flat hole band at the Brillouin zone center. Our theory predicts that, in order to observe a quantized CPGE in RhSi, it is necessary to increase the chemical potential as well as the quasiparticle lifetime. More broadly, our methodology, especially the development of the broadband terahertz emission spectroscopy, could be widely applied to study photogalvanic effects in noncentrosymmetric materials and in topological insulators in a contact-less way and accelerate the technological development of efficient infrared detectors based on topological semimetals.

Materials Science↗
DOE OSTI2024

Cluster expansion by transfer learning for phase stability predictions

Recent progress towards universal machine-learned interatomic potentials holds considerable promise for materials discovery. Yet the accuracy of these potentials for predicting phase stability may still be limited. In contrast, cluster expansions provide accurate phase stability predictions but are computationally demanding to parameterize from first principles, especially for structures of low dimension or with a large number of components, such as interfaces or multimetal catalysts. We overcome this trade-off via transfer learning. Using Bayesian inference, we incorporate prior statistical knowledge from machine-learned and physics-based potentials, enabling us to sample the most informative configurations and to efficiently fit first-principles cluster expansions. Furthermore, this algorithm is tested on Pt:Ni, showing robust convergence of the mixing energies as a function of sample size with reduced statistical fluctuations.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Mass of a weakly measured photon

Bohmian mechanics has garnered significant attention as an interpretation of quantum theory since the paradigmatic experiments which inferred the average trajectories of photons in the nonrelativistic regime. These experiments were largely motivated by Wiseman's formulation of Bohmian mechanics, which grounded these trajectories in weak measurements. Recently, Wiseman's framework was extended to the relativistic regime by expressing the velocity field of single photons in terms of weak values of the photon energy and momentum. Here, we propose an operational, weak value-based definition for the Bohmian “local mass” of relativistic single particles. For relativistic wave functions satisfying the scalar Klein-Gordon equation, this mass coincides with the effective mass defined by de Broglie in his relativistic pilot-wave theory, a quantity closely connected with the quantum potential that is responsible for Bohmian trajectory self-bending and the anomalous photoelectric effect. Here, we demonstrate the relationship between the photon trajectories and the mass in an interferometric setup.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1988

Tunable superlattice p-i-n photodetectors: Characteristics, theory, and applications

Extended measurements and theory on the recently developed monolithic wavelength demultiplexer consisting of voltage-tunable superlattice p-i-n photodetectors in a waveguide configuration are discussed. It is shown that the device is able to demultiplex and detect two optical signals with a wavelength separation of 20 nm directly into different electrical channels at a data rate of 1 Gb/s and with a crosstalk attenuation varying between 20 and 28 dB, depending on the polarization. The minimum acceptable crosstalk attenuation at a data rate of 100 Mb/s is determined to be 10 dB. The feasibility of using the device as a polarization angle sensor for linearly polarized light is also demonstrated. A theory for the emission of photogenerated carriers out of the quantum wells is included, since this is potentially a speed limiting mechanism in these detectors. It is shown that a theory of thermally assisted tunneling by polar optical phonon interaction is able to predict emission times consistent with the observed temporal response.

Larsson, Anders↗
NASA NTRS1991

On the thermodynamic framework of generalized coupled thermoelastic-viscoplastic-damage modeling

A complete potential based framework using internal state variables is put forth for the derivation of reversible and irreversible constitutive equations. In this framework, the existence of the total (integrated) form of either the (Helmholtz) free energy or the (Gibbs) complementary free energy are assumed a priori. Two options for describing the flow and evolutionary equations are described, wherein option one (the fully coupled form) is shown to be over restrictive while the second option (the decoupled form) provides significant flexibility. As a consequence of the decoupled form, a new operator, i.e., the Compliance operator, is defined which provides a link between the assumed Gibb's and complementary dissipation potential and ensures a number of desirable numerical features, for example the symmetry of the resulting consistent tangent stiffness matrix. An important conclusion reached, is that although many theories in the literature do not conform to the general potential framework outlined, it is still possible in some cases, by slight modifications of the used forms, to restore the complete potential structure.

Arnold, S. M.↗
NASA NTRS2004

Network Theory: A Primer and Questions for Air Transportation Systems Applications

A new understanding (with potential applications to air transportation systems) has emerged in the past five years in the scientific field of networks. This development emerges in large part because we now have a new laboratory for developing theories about complex networks: The Internet. The premise of this new understanding is that most complex networks of interest, both of nature and of human contrivance, exhibit a fundamentally different behavior than thought for over two hundred years under classical graph theory. Classical theory held that networks exhibited random behavior, characterized by normal, (e.g., Gaussian or Poisson) degree distributions of the connectivity between nodes by links. The new understanding turns this idea on its head: networks of interest exhibit scale-free (or small world) degree distributions of connectivity, characterized by power law distributions. The implications of scale-free behavior for air transportation systems include the potential that some behaviors of complex system architectures might be analyzed through relatively simple approximations of local elements of the system. For air transportation applications, this presentation proposes a framework for constructing topologies (architectures) that represent the relationships between mobility, flight operations, aircraft requirements, and airspace capacity, and the related externalities in airspace procedures and architectures. The proposed architectures or topologies may serve as a framework for posing comparative and combinative analyses of performance, cost, security, environmental, and related metrics.

Holmes, Bruce J.↗
DOE OSTI2026

Thermodynamics and its prediction and CALPHAD modeling: Review, state of the art, and perspectives

Thermodynamics is a science concerning the state of a system, whether it is stable, metastable, or unstable, when interacting with its surroundings. The combined law of thermodynamics derived by Gibbs about 150 years ago laid the foundation of thermodynamics. In Gibbs combined law, the entropy production due to internal processes was not included, and the 2nd law was thus practically removed from the Gibbs combined law, so it is only applicable to systems under equilibrium, thus commonly termed as equilibrium or Gibbs thermodynamics. Gibbs further derived the classical statistical thermodynamics in terms of the probability of configurations in a system in the later 1800's and early 1900's. With the quantum mechanics (QM) developed in 1920's, the QM-based statistical thermodynamics was established and connected to classical statistical thermodynamics at the classical limit as shown by Landau in the 1940's. In 1960's the development of density functional theory (DFT) by Kohn and co-workers enabled the QM prediction of properties of the ground state of a system. On the other hand, the entropy production due to internal processes in non-equilibrium systems was studied separately by Onsager in 1930's and Prigogine and co-workers in the 1950's. In 1960's to 1970's the digitization of thermodynamics was developed by Kaufman in the framework of the CALculation of PHAse Diagrams (CALPHAD) modeling of individual phases with internal degrees of freedom. CALPHAD modeling of thermodynamics and atomic transport properties has enabled computational design of complex materials in the last 50 years. Our recently termed zentropy theory integrates DFT and statistical mechanics through the replacement of the internal energy of each individual configuration by its DFT-predicted free energy. The zentropy theory is capable of accurately predicting the free energy of individual phases, transition temperatures and properties of magnetic and ferroelectric materials with free energies of individual configurations solely from DFT-based calculations and without fitting parameters, and is being tested for other phenomena including superconductivity, quantum criticality, and black holes. Those predictions include the singularity at critical points with divergence of physical properties, negative thermal expansion, and the strongly correlated physics. Furthermore, those individual configurations may thus be considered as the genomic building blocks of individual phases in the spirit of the materials genome®. This has the potential to shift the paradigm of CALPHAD modeling from being heavily dependent on experimental inputs to becoming fully predictive with inputs solely from DFT-based calculations and machine learning models built on those calculations and existing experimental data through newly developed and future open-source tools. Furthermore, through the combined law of thermodynamics including the internal entropy production, it is shown that the kinetic coefficient matrix of independent internal processes is diagonal with respect to the conjugate potentials in the combined law, and the cross phenomena that the phenomenological Onsager flux and reciprocal relationships are due to the dependence of the conjugate potential of a molar quantity on nonconjugate molar quantities and other potentials, which can be predicted by the zentropy theory and CALPHAD modeling.

42 ENGINEERING↗
DOE OSTI2025

Atomic Energy Accuracy of Neural Network Potentials: Harnessing Pretraining and Transfer Learning

Machine learning-based interatomic potentials (MLIPs) have transformed the prediction of potential energy surfaces (PESs), achieving accuracy comparable to ab initio calculations. However, atomic energy predictions, often assumed to lack physical meaning, remain underexplored. In this study, we demonstrate that inaccuracies in atomic energy predictions reduce the robustness and transferability of Neural Network Potentials (NNPs) and atomic energy error can be masked in total energy predictions due to error cancellation. Here, we validate this finding using challenging configurations involving deformation and failure under tensile loading. By pretraining atomic energy predictions using empirical potentials and applying transfer learning with density functional theory (DFT) data, we achieve notable improvements in the accuracy of total energy, forces, and stress predictions. Furthermore, this approach enhances the robustness and transferability of NNPs, emphasizing the importance of atomic energy predictions in developing high-quality and reliable MLIPs.

Active Learning↗