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,153 records · Page 64

Learning intermolecular forces at liquid–vapor interfaces

By adopting a perspective informed by contemporary liquid-state theory, we consider how to train an artificial neural network potential to describe inhomogeneous, disordered systems. Here, we find that neural network potentials based on local representations of atomic environments are capable of describing some properties of liquid-vapor interfaces but typically fail for properties that depend on unbalanced long-ranged interactions that build up in the presence of broken translation symmetry. These same interactions cancel in the translationally invariant bulk, allowing local neural network potentials to describe bulk properties correctly. By incorporating explicit models of the slowly varying long-ranged interactions and training neural networks only on the short-ranged components, we can arrive at potentials that robustly recover interfacial properties. We find that local neural network models can sometimes approximate a local molecular field potential to correct for the truncated interactions, but this behavior is variable and hard to learn. Generally, we find that models with explicit electrostatics are easier to train and have higher accuracy. We demonstrate this perspective in a simple model of an asymmetric dipolar fluid, where the exact long-ranged interaction is known, and in an ab initio water model, where it is approximated.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Hydrogen bond arrangements in (H 2 O) 20, 24, 28 clathrate hydrate cages: Optimization and many-body analysis

Here we provide a detailed study of hydrogen bonding arrangements, relative stability, residual entropy, and an analysis of the many-body effects in the (H 2 O) 20 (D-cage), (H 2 O) 24 (T-cage), and (H 2 O) 28 (H-cage) hollow cages making up structures I (sI) and II (sII) of clathrate hydrate lattices. Based on the enumeration of the possible hydrogen bonding networks for a fixed oxygen atom scaffold, the residual entropy (S 0 ) of these three gas phase cages was estimated at 0.754 82, 0.754 44, and 0.754 17 · Nk b , where N is the number of molecules and k b is Boltzmann’s constant. A previously identified descriptor of enhanced stability based on the relative arrangement and connectivity of nearest-neighbor fragments on the polyhedral water cluster [strong-weak-effective-bond model] also applies to the larger hollow cages. The three cages contain a maximum of 7, 9, and 11 such preferable arrangements of trans nearest dimer pairs with one “free” OH bond on the donor molecule (t1d dimers). The Many-Body Expansion (MBE) up to the 4-body suggests that the many-body terms vary nearly linearly with the cluster binding energy. Using a hierarchical approach of screening the relative stability of networks starting from optimizations with the TIP4P, TTM2.1-F, and MB-pol classical potentials, subsequently refining at more accurate levels of electronic structure theory (DFT and MP2), and finally correcting for zero-point energy, we were able to identify a group of four low-lying isomers of the (H 2 O) 24 T-cage, two of which are antisymmetric and the other two form a pair of antipode configurations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Understanding electronic structure tunability by metal dopants for promoting MgB 2 hydrogenation

Hydrogen is a promising energy carrier, but its onboard application is limited by the need for compact, low-pressure storage solutions. Solid-state complex metal hydride systems, such as MgB 2 /Mg(BH 4 ) 2 , offer high storage capacities but suffer from sluggish kinetics and poor reversibility. One avenue for improving reactivity is to introduce metal dopants to alter electronic and atomic properties, but the role of these chemical additives remains poorly understood, particularly for the hydrogenation reaction. In this work, we used density functional theory calculations on model MgB 2 systems to rationalize the potential role of metal dopants in destabilizing B–B bonding within the MgB 2 lattice. We carried out detailed electronic structure analyses for 28 different metal dopant adatoms to identify properties that contribute to a dopant’s efficacy. Based on the simulation results, we propose that an intermediate ionic and covalent character of the bonds between adatoms and B atoms is desirable for facilitating charge redistribution, disrupting the B–B bond network, and promoting H 2 dissociation and H atom chemisorption on MgB 2 .

08 HYDROGEN↗

Staircases of passive and active scalar concentration in cellular flow

This paper develops a unified model for staircase formation in both passive and active scalar systems, building upon prior numerical studies by offering new heuristic and physical insights. While prior studies primarily reported numerical results, they did not explore the underlying unifying physics that governs both types of scalar transport; this work addresses that gap by identifying shared mechanisms across both cases. Results of studies of passive and active scalar staircase formation in cellular flows are presented. Staircase formation in cellular flows occurs due to the interplay of fast mixing within cells and slow transport across the inter-cell boundary. The cell boundary emerges as a de facto transport barrier. Special attention is focused on the effects of cellular fluctuations and noise upon staircase structure. A forced, fluctuating vortex array model is used to drive the underlying flow structure. Cellular Peclet number and staircase profile curvature are identified as figures-of-merit to quantify the resiliency of layering. These are related to simple, multi-scatterer scalar random walk models. Results for Peclet number and curvature scaling with flow excitation are presented. We also study staircases of magnetic potential evolving in two-dimensional magnetohydrodynamics as examples of layering of active scalar concentration. Formation of magnetic potential staircases is indeed observed. Flux expulsion inhibits the intercellular transport of magnetic potential and strengthens staircase barriers. Magnetic staircases can be supported against resistive decay by magnetic potential noise forcing. Implications for staircase formation in magnetic confinement experiments are discussed.

Control theory↗

Effect of net direct current on the properties of radio frequency sheaths: simulation and cross-code comparison

In order to understand, predict and control ion cyclotron range of frequency (ICRF) interactions with tokamak scrape-off layer plasmas, computational tools which can model radio frequency (RF) sheaths are needed. In particular, models for the effective surface impedance and DC rectified sheath potentials may be coupled with full wave RF simulation codes to predict self-consistent wave fields near surfaces and the resulting power dissipation and plasma-material interactions from ion sputtering. In this study, previous work assuming zero net DC current flow through the sheath is generalized to allow the surface to collect net positive or negative current, as is often observed in experiments. The waveforms, DC potential and RF admittance are investigated by means of analytical theory, nonlinear fluid and particle-in-cell (PIC) codes. Cross-code comparisons provide detailed model verification and elucidate the roles of ion and electron kinetics. When the sheath draws negative (positive) DC current, the voltage rectification is reduced (increased) compared with the zero-current case, and the magnitude of both the real and imaginary parts of the admittance are increased (reduced). In conclusion, a previous four-input parametrization of the sheath rectification and admittance properties is generalized to include a fifth parameter describing the DC sheath current.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Neutron scattering study of magnetic anisotropy in the tetragonal antiferromagnet Bi 2 CuO 4

In this work, we present a comprehensive study of magnon excitations in the tetragonal easy-plane antiferromagnet Bi 2 CuO 4 using inelastic neutron scattering and spin-wave analyses. The nature of low-energy magnons, and hence the anisotropy in this material, has been controversial. We show unambiguously that the low-energy magnon spectrum consists of a gapped and a gapless mode, which we attribute to out-of-plane and in-plane spin fluctuations, respectively. We modeled the observed magnon spectrum using linear spin-wave analysis of a minimal anisotropic spin model motivated by the lattice symmetry. By studying the magnetic field dependence of the (1,0,0) Bragg peak intensity and the in-plane magnon intensity, we observed a spin-flop transition in the ab plane at ~0.4 T, which directly indicates the existence of a small in-plane anisotropy that is classically forbidden. It is only by taking into account magnon zero-point fluctuations beyond the linear spin-wave approximation that we could explain this in-plane anisotropy and its magnitude, the latter of which is deduced from the critical field of the spin-flop transition. The microscopic origins of the observed anisotropic interactions are also discussed. We found that our data are inconsistent with a large Dzyaloshinskii-Moriya interaction, which suggests a potential departure of Bi 2 CuO 4 from the conventional theories of magnetic anisotropy for other cuprates.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Role of the likelihood for elastic scattering uncertainty quantification

In the last decade, uncertainty quantification (UQ) for optical model potentials (OMPs) has become a focal point for nuclear reaction theory, and several competing approaches for OMP UQ have recently been developed. Here, we clarify recent efforts to compare frequentist and Bayesian approaches in the context of OMP UQ [G. B. King et al., Phys. Rev. Lett. 122, 232502 (2019)]. We replicate a portion of that OMP UQ study but use independent statistical tools. Specifically, we compare two methods for OMP parameter inference from elastic scattering data: the Levenberg-Marquardt algorithm for χ 2 minimization on one hand and Markov chain Monte Carlo (MCMC) sampling on the other. Separately, we assess the common practice of using a renormalized likelihood (χ 2 /N), N being the number of data points, instead of the canonical weighted-least-squares likelihood (χ 2 ), as a way of accounting for unknown data correlations. Here, we show that for a generic linear model and for a five-parameter OMP analysis, frequentist and uniform-prior Bayesian approaches recover the same optimum and uncertainty estimates—not systematically larger uncertainties for the Bayesian approach, as was concluded in G. B. King et al., Phys. Rev. Lett. 122, 232502 (2019). Further, we show that if an additional, near-degenerate parameter is introduced into the same OMP analysis such that the parameter posterior becomes non-Gaussian, then covariance-based estimates of uncertainty become unreliable. Finally, we show that regardless of optimization approach, if χ 2 /N is used for the likelihood, the resulting parametric uncertainties increase by $\sqrt{N}$, and that this is responsible for the conclusions drawn in the revisited study. Based on our replication results, we find that a fortuitous cancellation of unreported errors and the renormalization factor can lead to improvement in empirical coverages, as was the case in the original comparative study. We emphasize that developing and applying a realistic likelihood function is an essential task in a UQ analysis, and that several recent UQ studies that employed a renormalized likelihood (i.e., including a 1/N factor) may have yielded unrealistically large uncertainties for elastic-scattering observables. If the parameter posterior deviates from multivariate-normal, a sampling-based approach like MCMC has a clear advantage over methods that assume the Laplace approximation holds. We note that empirical coverage can serve as an important internal check for the analyst whose model or data may have additional, unaccounted-for uncertainties.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Out-of-plane polarization and topological magnetic vortices in multiferroic CrPSe 3

Two-dimensional multiferroic materials are ideal systems for exploring new coupling mechanisms between different ferroic orders and producing novel quantum phenomena with potential applications. Here we employ first-principles density functional theory calculations to discover intrinsic ferroelectric and antiferroelectric phases of Cr⁢P⁢Se 3 , which show ferromagnetic order and compete with the centrosymmetric phase with an antiferromagnetic order. Our analysis shows that the electrical dipoles of such type-I multiferroic phases come from the out-of-plane displacements of phosphorus ions due to the stereochemically active lone pairs. The coupling between polar and magnetic orders creates the opportunity for tuning the magnetic ground state by switching from the centrosymmetric to the ferroelectric phase using an out-of-plane electric field. In ferroelectric and antiferroelectric phases, the combination of easy-plane anisotropy and Dzyaloshinskii-Moriya interactions indicates that they can host topological magnetic vortices such as meron pairs.

36 MATERIALS SCIENCE↗

Theoretical investigation of surface chemical reactivity of compositionally disordered multi-metallic alloys

Alloys are widely used in both academic and industrial research as heterogeneous catalysts for accelerating technologically important reactions ranging from hydrocarbon reforming, selective hydrogenation, to electrochemical reactions involving hydrogen and oxygen that are fundamental and indispensable to a sustainable future energy infrastructure. Most basic catalysis research has been limited to combinations of two elements that are usually ordered in some manner at the nanoscale because doing so allows observable results to be easily rationalized. It however leaves a huge material space consisting of multi metallic, compositionally disordered alloys unexplored. This research demonstrates a viable theoretical approach combining different types of first principles calculations (including Green’s function-based Korringa–Kohn–Rostoker coherent potential approximation (KKR–CPA) and planewave-based density functional theory) to enable the surface chemical reactivity of compositionally complex alloy surfaces to be investigated. We demonstrate the approach by investigating hydrogen adsorption and evolution on two types of compositionally disordered alloys: a quaternary CoCrFeNi high-entropy alloy, and binary Pt-3d base metal alloys.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Computational Imaging for Intelligence in Highly Scattering Aerosols (Final Report)

Natural and man-made degraded visual environments pose major threats to national security. The random scattering and absorption of light by tiny particles suspended in the air reduces situational awareness and causes unacceptable down-time for critical systems and operations. To improve the situation, we have developed several approaches to interpret the information contained within scattered light to enhance sensing and imaging in scattering media. These approaches were tested at the Sandia National Laboratory Fog Chamber facility and with tabletop fog chambers. Computationally efficient light transport models were developed and leveraged for computational sensing. The models are based on a weak angular dependence approximation to the Boltzmann or radiative transfer equation that appears to be applicable in both the moderate and highly scattering regimes. After the new model was experimentally validated, statistical approaches for detection, localization, and imaging of objects hidden in fog were developed and demonstrated. A binary hypothesis test and the Neyman-Pearson lemma provided the highest theoretically possible probability of detection for a specified false alarm rate and signal-to-noise ratio. Maximum likelihood estimation allowed estimation of the fog optical properties as well as the position, size, and reflection coefficient of an object in fog. A computational dehazing approach was implemented to reduce the effects of scatter on images, making object features more readily discernible. We have developed, characterized, and deployed a new Tabletop Fog Chamber capable of repeatably generating multiple unique fog-analogues for optical testing in degraded visual environments. We characterized this chamber using both optical and microphysical techniques. In doing so we have explored the ability of droplet nucleation theory to describe the aerosols generated within the chamber, as well as Mie scattering theory to describe the attenuation of light by said aerosols, and correlated the aerosol microphysics to optical properties such as transmission and meteorological optical range (MOR). This chamber has proved highly valuable and has supported multiple efforts inclusive to and exclusive of this LDRD project to test optics in degraded visual environments. Circularly polarized light has been found to maintain its polarization state better than linearly polarized light when propagating through fog. This was demonstrated experimentally in both the visible and short-wave infrared (SWIR) by imaging targets made of different commercially available retroreflective films. It was found that active circularly polarized imaging can increase contrast and range compared to linearly polarized imaging. We have completed an initial investigation of the capability for machine learning methods to reduce the effects of light scattering when imaging through fog. Previously acquired experimental long-wave images were used to train an autoencoder denoising architecture. Overfitting was found to be a problem because of lack of variability in the object type in this data set. The lessons learned were used to collect a well labeled dataset with much more variability using the Tabletop Fog Chamber that will be available for future studies. We have developed several new sensing methods using speckle intensity correlations. First, the ability to image moving objects in fog was shown, establishing that our unique speckle imaging method can be implemented in dynamic scattering media. Second, the speckle decorrelation over time was found to be sensitive to fog composition, implying extensions to fog characterization. Third, the ability to distinguish macroscopically identical objects on a far-subwavelength scale was demonstrated, suggesting numerous applications ranging from nanoscale defect detection to security. Fourth, we have shown the capability to simultaneously image and localize hidden objects, allowing the speckle imaging method to be effective without prior object positional information. Finally, an interferometric effect was presented that illustrates a new approach for analyzing speckle intensity correlations that may lead to more effective ways to localize and image moving objects. All of these results represent significant developments that challenge the limits of the application of speckle imaging and open important application spaces. A theory was developed and simulations were performed to assess the potential transverse resolution benefit of relative motion in structured illumination for radar systems. Results for a simplified radar system model indicate that significant resolution benefits are possible using data from scanning a structured beam over the target, with the use of appropriate signal processing.

58 GEOSCIENCES↗

Tritium Trapping Thermodynamics by Point Defects at Interfacial Fe-Al Phases of the Aluminide Coating

Density functional theory (DFT) simulations have been carried out to evaluate the potential for tritium trapping by metal vacancies in four phases (i.e., (Fe, Cr, Ni), AlFe 3 , NiAl, and AlFe) identified near the interface between the Al coating and 316 stainless steel (316 SS) cladding. In addition, an ab initio thermodynamics approach has been employed to predict the temperature and T 2 -partial pressure dependence on the thermodynamics of tritiated defects. Key results in this work suggest that metal vacancies in the four phases have the potential to favorably trap tritium species. This thermodynamic trend can be correlated to the high energy cost of having an interstitial tritium in the lattice, which has been calculated to be at least 0.48 eV. By combining the results of this study with previous theoretical works investigating tritium behavior in other Fe-Al phases identified in the aluminide coating, suggest that metal vacancies are generally able to trap tritium in various Fe-Al aluminide phases. Especially, it was found that Al vacancies are the most efficient to trap tritium, followed by Fe vacancies, then Ni and Cr vacancies. Strong interactions between tritium and the metal vacancies are occurring by the formation of Fe—T bonds. The formation of Cr—T, Ni—T, or Al—T bonds are found less energetically favorable than Fe—T bonds.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Thermodynamics of Tritium Trapping by Point Defects at Interfacial Cr-based Phases of the Al Coating

Density functional theory (DFT) simulations have been carried out to evaluate the potential for tritium trapping by metal vacancies in four phases Cr-containing phases (i.e., Cr3Si, Al0.3Cr0.7, Al8Cr5-HT, and Al8Cr5) identified near the interface between the Al coating and 316 stainless steel (316 SS) cladding. In addition, an ab initio thermodynamics approach has been employed to predict the temperature and T2-partial pressure dependence on the thermodynamics of singly tritiated defects. Key results in this work suggest that metal vacancies in the four phases have the potential to favorably trap tritium species, especially Si vacancies in Cr3Si phase. This overall thermodynamic trend can be correlated to the energy cost of having an interstitial tritium in the lattice, which has been calculated to range from 0.17 eV in Al8Cr5-HT to 0.89 eV in Cr3Si. A comparison of the results from this study with previous theoretical works investigating tritium behavior in other Fe-Al phases identified in the aluminide coating, suggests that metal vacancies are generally able to trap tritium in various Fe-Al aluminide phases. Especially, it was found that Si and Al vacancies would be the most efficient to trap for tritium, followed by Cr vacancies, then Fe and Ni vacancies. In the four Cr-based material phases investigated in this work, it is interesting to note that, in the absences of Fe or Ni species, strong interactions between tritium and the metal vacancies are always occurring by the formation of preferential Cr—T bonds (i.e., no Si—T or Al—T bonds were formed).

Sassi, Michel J.↗

Muon capture on deuteron using local chiral potentials

The muon capture reaction μ – + d → n + n + ν μ in the doublet hyperfine state is studied using nuclear potentials and consistent currents derived in the chiral effective field theory, which are local and expressed in coordinate space (the so-called Norfolk models). Only the largest contribution due to the 1 S 0 nn scattering state is considered. Particular attention is given to the estimate of theoretical uncertainty, for which four sources have been identified: 1) the model dependence, 2) the chiral-order convergence for the weak nuclear current, 3) the uncertainty in the single-nucleon axial form factor, and 4) the numerical technique adopted to solve the bound and scattering A = 2 systems. This last source of uncertainty has turned out to be essentially negligible. For the 1 S 0 doublet muon capture rate Γ D ( 1 S 0 ), we obtain Γ D ( 1 S 0 ) = 255.8(0.6) (4.4)(2.9) s –1 , where the three errors come from the first three sources of uncertainty. The value for Γ D ( 1 S 0 ) obtained within this local chiral framework is compared with previous calculations and found in very good agreement.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Bilocal Field Theory for Composite Scalar Bosons

We give a bilocal field theory description of a composite scalar with an extended binding potential that reduces to the Nambu–Jona-Lasinio (NJL) model in the pointlike limit. This provides a description of the internal dynamics of the bound state and features a static internal wave function, $\phi(\vec{r})$, in the center-of-mass frame that satisfies a Schrödinger–Klein–Gordon equation with eigenvalues m 2 . We analyze the “coloron” model (single perturbative massive gluon exchange) which yields a UV completion of the NJL model. This has a BCS-like enhancement of its interaction, ∝N c the number of colors, and is classically critical with g critical remarkably close to the NJL quantum critical coupling. Negative eigenvalues for m 2 lead to spontaneous symmetry breaking, and the Yukawa coupling of the bound state to constituent fermions is emergent.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Heat Transfer and Thermodynamics: a Compilation

A compilation is presented for the dissemination of information on technological developments which have potential utility outside the aerospace and nuclear communities. Studies include theories and mechanical considerations in the transfer of heat and the thermodynamic properties of matter and the causes and effects of certain interactions.

Source record↗

Computer simulation of screw dislocation in aluminum

The atomic structure in a 110 screw dislocation core for aluminum is obtained by computer simulation. The lattice statics technique is employed since it entails no artificially imposed elastic boundary around the defect. The interatomic potential has no adjustable parameters and was derived from pseudopotential theory. The resulting atomic displacements were allowed to relax in all three dimensions.

Esterling, D. M.↗

Planetary tides during the Maunder Sunspot Minimum

In order to test the tidal theory of sunspots, sun-centered planetary conjunctions and tidal potentials are reconstructed for the period of the Maunder Minimum (1645 to 1715). These are found to be effectively indistinguishable from patterns of conjunctions and power spectra of tidal potential in the modern era of a well-established 11-yr sunspot cycle. The pattern of planetary tidal forces during the Maunder Minimum is then reconstructed to investigate the possibility that multiple-planet forces were somehow fortuituously cancelled at that time; i.e., the positions of the slower moving planets in the late 17th and early 18th centuries were such that conjunctions and tidal potentials were reduced in number and force. Calculations of daily positions for Mercury, Venus, earth, and Jupiter as well as daily values of the tidal potential for the period from 1450 to 2000 indicate no striking dissimilarities between the time of the Maunder Minimum and any other period considered.

Smythe, C. M.↗

Adaptation of the TH Epsilon Mu formalism for the analysis of the equivalence principle in the presence of the weak and electroweak interaction

The TH epsilon mu formalism, used in analyzing equivalence principle experiments of metric and nonmetric gravity theories, is adapted to the description of the electroweak interaction using the Weinberg-Salam unified SU(2) x U(1) model. The use of the TH epsilon mu formalism is thereby extended to the weak interactions, showing how the gravitational field affects W sub mu (+ or -1) and Z sub mu (0) boson propagation and the rates of interactions mediated by them. The possibility of a similar extension to the strong interactions via SU(5) grand unified theories is briefly discussed. Also, using the effects of the potentials on the baryon and lepton wave functions, the effects of gravity on transition mediated in high-A atoms which are electromagnetically forbidden. Three possible experiments to test the equivalence principle in the presence of the weak interactions, which are technologically feasible, are then briefly outline: (1) K-capture by the FE nucleus (counting the emitted X-ray); (2) forbidden absorption transitions in high-A atoms' vapor; and (3) counting the relative Beta-decay rates in a suitable alpha-beta decay chain, assuming the strong interactions obey the equivalence principle.

Fennelly, A. J.↗