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 1,081 records · Page 60
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Securely Aggregated Coded Matrix Inversion

Coded computing is a method for mitigating straggling workers in a centralized computing network, by using erasure-coding techniques. Federated learning is a decentralized model for training data distributed across client devices. In this work we propose approximating the inverse of an aggregated data matrix, where the data is generated by clients; similar to the federated learning paradigm, while also being resilient to stragglers. To do so, we propose a coded computing method based on gradient coding. We modify this method so that the coordinator does not access the local data at any point; while the clients access the aggregated matrix in order to complete their tasks. Here, the network we consider is not centrally administrated, and the communications which take place are secure against potential eavesdroppers.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

An investigation into the potential for wind turbines to cause barotrauma in bats

The high rates of bat mortality caused by operating wind turbines is a concern for wind energy and wildlife stakeholders. One theory that explains the mortality is that bats are not only killed by impact trauma, but also by barotrauma that results from exposure to the pressure variations caused by rotating turbine blades. To date, no published research has calculated the pressure changes that bats may be exposed to when flying near wind turbines and then used these data to estimate the likelihood that turbines cause barotrauma in bats. To address this shortcoming, we performed computational fluid dynamics simulations of a wind turbine and analytical calculations of blade-tip vortices to estimate the characteristics of the sudden pressure changes bats may experience when flying near a utility-scale wind turbine. Because there are no data available that characterize the pressure changes that cause barotrauma in bats, we compared our results to changes in pressure levels that cause barotrauma and mortality in other mammals of similar size. This comparison shows that the magnitude of the low-pressures bats experience when flying near wind turbines is approximately 8 times smaller than the pressure that causes mortality in rats, the smallest mammal for which data are available. The magnitude of the high-pressures that bats may experience are approximately 80 times smaller than the exposure level that causes 50% mortality in mice, which have a body mass similar to several bat species that are killed by wind turbines. Further, our results show that for a bat to experience the largest possible magnitude of low- and high-pressures, they must take very specific and improbable flight paths that skim the surface of the blades. Even a small change in the flight path results in the bat being hit by the blade or experiencing a much smaller pressure change. Accordingly, if bats have a physiological response to rapid low- and high-pressure exposure that is similar to other mammals, we conclude that it is unlikely that barotrauma is responsible for a significant number of turbine-related bat fatalities, and that impact trauma is the likely cause of the majority of wind-turbine-related bat fatalities.

17 WIND ENERGY↗
DOE OSTI2020

Weak Transitions in Light Nuclei

Nuclei are used for high-precision tests of the Standard Model and for studies of physics beyond the Standard Model. Without a thorough understanding of nuclei, we will not be able to meaningfully interpret the growing body of experimental data nor will we be able to disentangle new physics signals from underlying nuclear effects. This calls for accurate calculations of nuclear structure and reactions. In this work, we focus on electroweak decays in nuclei with mass number A ≤ 10 and report on ab initio Quantum Monte Carlo calculations of reduced matrix elements entering beta decays and electron captures in nuclei with mass number A ≤ 10. The many-body wave functions are calculated using selected Norfolk two- and three-nucleon potential models and associated one- and two-body axial currents at tree-level obtained from a chiral effective field theory with pions, nucleons, and Δ. The agreement with the experimental data is satisfactory except for transitions in A = 8 nuclei. In this specific case, the theory significantly underpredicts the experimental data, which indicates the need of further improvements in the corresponding nuclear wave functions. In this study, emphasis is placed on the contributions of two-body axial currents that are carefully analyzed using two-body transition densities. This allow us to study the spatial distribution and short-range behavior of two-body dynamics. In particular, the transition densities when scaled to peak at 1.0 exhibit universal short-range behavior across the considered nuclei, while they differ in the long-range tails.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Star-exoplanet interactions: A growing interdisciplinary field in heliophysics

Traditionally, heliophysics is characterized as the study of the near-Earth space environment, where plasmas and neutral gases originating from the Earth, the Sun, and other solar system bodies interact in ways that are detectable only through in-situ or close-range (usually within ~10 AU) remote sensing. As a result, heliophysics has data from the space environment around a handful of solar system objects, in particular the Sun and Earth. Comparatively, astrophysics has data from an extensive array of objects, but is more limited in temporal, spatial, and wavelength information from any individual object. Thus, our understanding of planetary space environments as a complex, multi-dimensional network of specific interacting systems may in the past have seemed to have little to do with the highly diverse space environments detected through astrophysical methods. Recent technological advances have begun to bridge this divide. Exoplanetary studies are opening up avenues to study planetary environments beyond our solar system, with missions like Kepler, TESS, and JWST, along with increasing capabilities of ground-based observations. At the same time, heliophysics studies are pushing beyond the boundaries of our heliosphere with Voyager, IBEX, and the future IMAP mission. The interdisciplinary field of star-exoplanet interactions is a critical, growing area of study that enriches heliophysics. A multidisciplinary approach to heliophysics enables us to better understand universal processes that operate in diverse environments, as well as the evolution of our solar system and extreme space weather. The expertise, data, theory, and modeling tools developed by heliophysicists are crucial in understanding the space environments of exoplanets, their host stars, and their potential habitability. The mutual benefit that heliophysics and exoplanetary studies offer each other depends on strong, continuing solar system-focused and Earth-focused heliophysics studies. The heliophysics discipline requires new targeted funding to support inter-divisional opportunities, including small multi-disciplinary research projects, large collaborative research teams, and observations targeting the heliophysics of planetary and exoplanet systems. Here we discuss areas of heliophysics-relevant exoplanetary research, observational opportunities and challenges, and ways to promote the inclusion of heliophysics within the wider exoplanetary community.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Graph-based quantum response theory and shadow Born–Oppenheimer molecular dynamics

Graph-based linear scaling electronic structure theory for quantum-mechanical molecular dynamics simulations [A. M. N. Niklasson et al., J. Chem. Phys. 144, 234101 (2016)] is adapted to the most recent shadow potential formulations of extended Lagrangian Born–Oppenheimer molecular dynamics, including fractional molecular-orbital occupation numbers [A. M. N. Niklasson, J. Chem. Phys. 152, 104103 (2020) and A. M. N. Niklasson, Eur. Phys. J. B 94, 164 (2021)], which enables stable simulations of sensitive complex chemical systems with unsteady charge solutions. The proposed formulation includes a preconditioned Krylov subspace approximation for the integration of the extended electronic degrees of freedom, which requires quantum response calculations for electronic states with fractional occupation numbers. For the response calculations, we introduce a graph-based canonical quantum perturbation theory that can be performed with the same natural parallelism and linear scaling complexity as the graph-based electronic structure calculations for the unperturbed ground state. Further, the proposed techniques are particularly well-suited for semi-empirical electronic structure theory, and the methods are demonstrated using self-consistent charge density-functional tight-binding theory both for the acceleration of self-consistent field calculations and for quantum-mechanical molecular dynamics simulations. Graph-based techniques combined with the semi-empirical theory enable stable simulations of large, complex chemical systems, including tens-of-thousands of atoms.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

In-situ Synthesis of PuCl 3 and Corresponding Raman and Density Functional Theory Vibrational Modes

The experimental and calculated Raman spectrum of PuCl 3 has been reported for the first time. PuCl 3 is a primary species found in plutonium metal refinement, specifically in pyrochemical salt processes including multicycle direct oxide reduction, metal chlorination, and electrorefining. As such, Raman signatures of PuCl 3 could serve as potential forensic indicators of material process history. A novel technique for synthesizing PuCl 3 from the in situ chlorination of plutonium metal with HCl was developed to establish these signatures. Cerium metal surrogates were utilized to ensure optimization of the plutonium experiments and to minimize personnel exposure, and all experiments were carried out in a Raman reaction chamber designed for air-tight, high vacuum environments. In situ Raman spectroscopy was employed in conjunction with density functional theory (DFT) to investigate the vibrational modes of PuCl 3 . Associated mixed oxy and hydroxyl phases are also reported. The combined Raman and DFT results have eliminated inconsistencies in Raman mode assignments for the MX 3 family of metal chlorides having P6 3 /m symmetry, and IR modes derived from the DFT calculations are additionally presented. The data observed in this study are of potential interest to nuclear forensic analyses, nuclear sample aging, nuclear energy, plutonium processing, and stockpile stewardship.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Multifidelity Neural Network Formulations for Prediction of Reactive Molecular Potential Energy Surfaces

Here, this paper focuses on the development of multifidelity modeling approaches using neural network surrogates, where training data arising from multiple model forms and resolutions are integrated to predict high-fidelity response quantities of interest at lower cost. We focus on the context of quantum chemistry and the integration of information from multiple levels of theory. Important foundations include the use of symmetry function-based atomic energy vector constructions as feature vectors for representing structures across families of molecules and single-fidelity neural network training capabilities that learn the relationships needed to map feature vectors to potential energy predictions. These foundations are embedded within several multifidelity topologies that decompose the high-fidelity mapping into model-based components, including sequential formulations that admit a general nonlinear mapping across fidelities and discrepancy-based formulations that presume an additive decomposition. Methodologies are first explored and demonstrated on a pair of simple analytical test problems and then deployed for potential energy prediction for C 5 H 5 using B2PLYP-D3/6-311++G(d,p) for high-fidelity simulation data and Hartree–Fock 6-31G for low-fidelity data. For the common case of limited access to high-fidelity data, our computational results demonstrate that multifidelity neural network potential energy surface constructions achieve roughly an order of magnitude improvement, either in terms of test error reduction for equivalent total simulation cost or reduction in total cost for equivalent error.

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

Information Entropy as Quantifier of Potential Predictability in the Tropical Indo-Pacific Basin

Global warming is posed to modify the modes of variability that control much of the climate predictability at seasonal to interannual scales. The quantification of changes in climate predictability over any given amount of time, however, remains challenging. Here we build upon recent advances in non-linear dynamical systems theory and introduce the climate community to an information entropy quantifier based on recurrence. The entropy, or complexity of a system is associated with microstates that recur over time in the time-series that define the system, and therefore to its predictability potential. A computationally fast method to evaluate the entropy is applied to the investigation of the information entropy of sea surface temperature in the tropical Pacific and Indian Oceans, focusing on boreal fall. In this season the predictability of the basins is controlled by two regularly varying non-linear oscillations, the El Niño-Southern Oscillation and the Indian Ocean Dipole. We compute and compare the entropy in simulations from the CMIP5 catalog from the historical period and RCP8.5 scenario, and in reanalysis datasets. Discrepancies are found between the models and the reanalysis, and no robust changes in predictability can be identified in future projections. The Indian Ocean and the equatorial Pacific emerge as troublesome areas where the modeled entropy differs the most from that of the reanalysis in many models. A brief investigation of the source of the bias points to a poor representation of the ocean mean state and basins' connectivity at the Indonesian Throughflow.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Freely jointed chain models with extensible links

We report analytical relations for the mechanical response of single polymer chains are valuable for modeling purposes, on both the molecular and the continuum scale. These relations can be obtained using statistical thermodynamics and an idealized single-chain model, such as the freely jointed chain model. To include bond stretching, the rigid links in the freely jointed chain model can be made extensible, but this almost always renders the model analytically intractable. Here, an asymptotically correct statistical thermodynamic theory is used to develop analytic approximations for the single-chain mechanical response of this model. The accuracy of these approximations is demonstrated using several link potential energy functions. This approach can be applied to other single-chain models, and to molecular stretching in general.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Fundamentals of Spark-Plasma Sintering (Final Report)

This project was focused on the development of the general theory of spark-plasma sintering (SPS) taking into account the role of both thermal and non-thermal factors in the acceleration of SPS mass transport. SPS is a particular kind of field-assisted sintering, which provides potentially revolutionary capabilities to the processing of materials into configurations previously unattainable. This approach significantly improves the processing time- and quality-wise. It carries the potential of maintaining the nano and sub-micron structure in nano-powder-based materials after consolidation. SPS gains particular prominence in connection with its exceptional potential of rapid and ultra-rapid processing of very hard-to-deform materials, which would typically require lengthy consolidation times at significantly elevated temperatures under conditions of conventional powder pressing or sintering. The achievement of the ultimate goal of the formulation of the general theory of spark plasma sintering required the identification of the contributions of all the thermal and non-thermal factors in the enhancement of mass transport under SPS processing conditions. As our research indicated, this enhancement is mostly reduced to the creation of the conditions of controlled non-equilibrium, which is especially important in ultra-rapid field-assisted sintering techniques. Thus, the main project objective was the analysis of the SPS physical basis at multiple scales specifically exploring the role of electric current in the acceleration of mass transport and with an emphasis on the conditions of controlled non-equilibrium.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Exploring the strong-coupling region of SU( N ) Seiberg-Witten theory

We consider the Seiberg-Witten solution of pure N = 2 gauge theory in four dimensions, with gauge group SU(N). A simple exact series expansion for the dependence of the 2(N – 1) Seiberg-Witten periods a I (u), a DI (u) on the N – 1 Coulomb-branch moduli un is obtained around the Z 2N -symmetric point of the Coulomb branch, where all u n vanish. This generalizes earlier results for N = 2 in terms of hypergeometric functions, and for N = 3 in terms of Appell functions. Using these and other analytical results, combined with numerical computations, we explore the global structure of the Kähler potential K = 1/2Σ I Im(a¯ I a DI ), which is single valued on the Coulomb branch. Evidence is presented that K is a convex function, with a unique minimum at the Z 2N -symmetric point. Finally, we explore candidate walls of marginal stability in the vicinity of this point, and their relation to the surface of vanishing Kähler potential.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

A multi-dimensional Child–Langmuir law for any diode geometry

While prior theoretical studies of multi-dimensional space-charge limited current (SCLC) assumed emission from a small patch on infinite electrodes, none have considered emission from an entire finite electrode. In this paper, we apply variational calculus (VC) and conformal mapping, which have previously been used to derive analytic solutions for SCLC density (SCLCD) for nonplanar one-dimensional geometries, to obtain mathematical relationships for any multi-dimensional macroscopic diode with finite cathode and anode. We first derive a universal mathematical relationship between space-charge limited potential and vacuum potential for any diode and apply this technique to determine SCLCD for an eccentric spherical diode. We then apply VC and the Schwartz–Christoffel transformation to derive an exact equation for SCLCD in a general two-dimensional planar geometry with emission from a finite emitter. Particle-in-cell simulations using VSim agreed within 4%–13% for a range of ratios of emitter width to gap distance using the thinnest electrodes practical for the memory constraints of our hardware, with the difference partially attributed to the theory's assumption of infinitesimally thin electrodes. After generalizing this approach to determine SCLCD for any orthogonal diode as a function of only the vacuum capacitance and vacuum potential, we derive an analytical formulation of the three-dimensional Child–Langmuir law for finite parallel rectangular and disk geometries. These results demonstrate the utility for calculating SCLCD for any diode geometry using vacuum capacitance and vacuum potential, which are readily obtainable for many diode geometries, to guide experiment and simulation development.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Electrochemical Degradation of Pt 3 Co Nanoparticles Investigated by Off-Lattice Kinetic Monte Carlo Simulations with Machine-Learned Potentials

In fuel cell applications, the durability of catalysts is critical for large-scale industrial implementation. However, limited synthesis controllability and spectroscopic resolution impede a comprehensive understanding of degradation mechanisms at the atomic level. In this study, we develop a machine-learned potential (MLP) to simulate the degradation processes for Pt3Co nanoparticles. The precision of MLP is determined to be comparable to that of density functional theory calculations. Using off-lattice kinetic Monte Carlo simulations with MLP, we successfully replicate established experimental trends and offer a logical resolution to ongoing debates regarding atomic orderings. Based on the simulation results, we suggest design principles for Pt3Co nanoparticles that combine high activity and durability. Finally, we validate the wide applicability of our method by successfully applying it to Pt3Ni and Pt3Co0.5Ni0.5 nanoparticles. The research serves as a guideline for developing MLPs for alloy electrochemical catalysts and lays the foundation for designing more durable and active fuel-cell catalysts.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Redox Potentials for Hydrogen Species in Acetonitrile and Tetrahydrofuran

The properties of the fundamental hydrogen species, proton (H+), hydrogen atom (H•), and hydride anion (H-) are critical to a vast range of chemical processes, yet their thermodynamic properties in nonaqueous solvents are not well established. A hybrid supermolecule-continuum approach is used to predict the Gibbs free energies of solvation (?G°solv) and standard redox potentials (E°) for the 2H+/H2 and H•/H- couples in acetonitrile (MeCN) and tetrahydrofuran (THF) following the approach previously used for water. Gas phase geometry optimizations were done using density functional theory (DFT) with the ?B97XD functional, Møller-Plesset second-order perturbation theory (MP2) with augmented correlation-consistent basis sets and composite-correlated molecular orbital theory (G3MP2 and G3MP2B3). Additional single point calculations in gas phase were performed for H- affinity for MeCN and THF at the coupled cluster CCSD(T) level using MP2 geometries. Solvation was included using the self-consistent reaction field SMD model. For H+, the ?G°solv values are predicted to be -252.2 kcal/mol in MeCN and -261.0 kcal/mol in THF. The redox potentials for the 2H+/H2 couple are +0.49 V in MeCN and +0.11 V in THF relative to the aqueous SHE, consistent with the available experimental data. The ?G°solv for H- is predicted to be -80 kcal/mol in MeCN and -68 in THF, and for H•, ?G°solv is predicted to be -1.3 kcal/mol in THF and -2.0 kcal/mol in MeCN. These solvation energies yield respective calculated redox potentials for the H•/H- couple of -0.14 V in MeCN and -0.68 V in THF relative to the aqueous SHE. H2 solvation is slightly exergonic in MeCN (-3.3 kcal/mol) and near thermoneutral in THF (+1.3 kcal/mol). The pKa(H2) is predicted to be 45.4 in MeCN and 47.8 in THF. The proton-coupled electron transfer (PCET) accounting term CG was evaluated giving values on the ferrocene scale within approximately 2 kcal/mol of literature in both solvents. These results provide benchmarks for hydrogen redox chemistry in MeCN and THF and improve consistency across E°, ?G°solv, pKa, and CG. These results provide benchmarks for hydrogen redox chemistry in nonaqueous media challenging previous estimates based on indirect extrapolations and assumptions.

Duda, Damian P.↗
DOE OSTI2025

Embedded random phase approximation for magnetic systems: H 2 dissociative adsorption on Fe(110)

The random phase approximation (RPA), a method for treating electron correlation, has been shown to be superior to standard density functional theory (DFT) approximations in numerous cases. However, the RPA’s computational cost is substantially higher than that of DFT, particularly restricting its application to extended surfaces. The recently introduced embedded RPA (emb-RPA) approach [Wei et al., J. Chem. Phys. 159(19), 194108 (2023)] reduces this computational cost by approximately two orders of magnitude. While previous applications of emb-RPA focused on non-spin-polarized systems, here we extend the approach to ferromagnetic ones. Unlike other embedded correlated wavefunction methods, such as embedded complete active space self-consistent field theory, emb-RPA is advantageous for spin-polarized systems because the RPA is compatible with unrestricted DFT solutions, which are eigenfunctions of the spin angular momentum operator S z but not the total spin-squared operator S 2 . By applying emb-RPA with specific magnetization constraints, we achieved a speedup of two to three orders of magnitude (one order when accounting for the one-time embedding potential optimization cost) with only small errors (∼50 meV) compared to full periodic RPA. Moreover, emb-RPA significantly reduces the over-binding errors of DFT approximations. In conclusion, we anticipate that the acceleration enabled by the spin-polarized emb-RPA approach will broaden the applicability of RPA to magnetic materials.

Density functional theory↗
DOE OSTI2023

Dynamic Molecular Investigation of the Solid-Electrolyte Interphase of an Anode-Free Lithium Metal Battery Using In Situ Liquid SIMS and Cryo-TEM

Solid electrolyte interphase (SEI) has been widely perceived to play a critical role in the stable cycling of rechargeable batteries. However, associated with the fragile and air-sensitive nature of the SEI layer, delineation of the formation process and the nature of SEI remains a big challenge. Here, we use in situ liquid time-of-flight secondary ion mass spectroscopy (TOF-SIMS), cryo- transmission electron microscope (TEM) and density functional theory (DFT) calculation to delineate molecular process on the formation of SEI layer under the dynamic operating condition. We discover that the onset potential for SEI layer formation and the thickness of the SEI show dependence on the solvation shell structure. Using LiCoO 2 as a cathode and Cu film as an anode, the SEI is noticed to start to form at around 2.0 V and reach its final thickness (irreversible part, ~ 40-50 nm) at about 3.0 V in the 1 M LiPF 6 –EC/DMC electrolyte, while for the case of 1 M LiFSI–DME, the SEI starts to form at around 1.5 V and reaches its final thickness (~ 20 nm) at about 2.0 V. The in situ TOF-SIMS clearly indicates the outer SEI layer formation and dissipation upon charging and discharging, implying a continued evolution of electrolyte structure with extended cycling. In conclusion, the present work establishes a direct correlation between the molecular signature of SEI layer with solvation feature of electrolytes in lithium batteries, providing insights for tailoring SEI layer toward improved electrochemical properties of lithium batteries.

25 ENERGY STORAGE↗
DOE OSTI2020

Exploring the role of electronic structure on photo-catalytic behavior of carbon-nitride polymorphs

A fully self-consistent density-functional theory (DFT) with improved functionals is used herein to provide a comprehensive account of structural, electronic, and optical properties of C 3 N 4 polymorphs. Using our recently developed van Leeuwen-Baerends (vLB) corrected local-density approximation (LDA), we implemented LDA + vLB within full-potential N th -order muffin-tin orbital (FP-NMTO) method and show that it improves structural properties and band gaps compared to semi-local functionals (LDA/GGA). We demonstrate that the LDA + vLB predicts band-structure and work-function for well-studied 2D-graphene and bulk-Si in very good agreement with experiments, and more exact hybrid functional (HSE) calculations as implemented in the Quantum-Espresso (QE) package. The structural and electronic-structure (band gap) properties of C 3 N 4 polymorphs calculated using FP-NMTO-LDA + vLB is compared with more sophisticated hybrid-functional calculations. We also perform detailed investigation of photocatalytic behavior using QE-HSE method of C 3 N 4 polymorphs through work-function, band (valence and conduction) position with respect to water reduction and oxidation potential. Our results show $γ$-C 3 N 4 as the best candidate for photocatalysis among all the C 3 N 4 polymorphs but it is dynamically unstable at ‘zero’ pressure. We show that $γ$-C 3 N 4 can be stabilized under hydrostatic-pressure, which improves its photocatalytic behavior relative to water reduction and oxidation potentials.

36 MATERIALS SCIENCE↗