Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “phase space 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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Mapping valence electron distributions with multipole density formalism using 4D-STEM

Recent advancement in aberration correction and detector technology opened a door to various applications using 4D-STEM, which yields a diffraction pattern for each scanning position within a crystal unit-cell in scanning transmission electron microscopy (STEM) and generates incredible amounts of data in momentum space. Currently 4D-STEM analysis relies on the center-of-mass of the diffraction patterns in electric field and charge density mapping. It only derives the total projected charge density and is limited to phase objects, e.g. extremely thin samples. Here, we propose a new analytical method to accurately map aspherical valence electron distributions with atom-centered multipolar functions formalism using the whole 4D-STEM dataset. We demonstrate that, with the full dynamical calculations for various sample thicknesses, the method is sensitive not only to the miniscule charge transfer, but also to the atomic site symmetry and aspherical electron orbitals. The process of the refinement is much more robust and reliable than quantitative convergent beam electron diffraction.

4D-STEM↗
DOE OSTI2023

Accurate prediction of short-range order and its effect on thermodynamic, structural, and electronic properties of disordered alloys: exemplified in archetypical Cu 3 Au

Electronic-structure methods based on density-functional theory (DFT) were used to quantify the effect of chemical short-range order (SRO) on thermodynamic, structural, and electronic properties of archetypal face-centered-cubic (fcc) Cu3Au alloy. We showed that SRO can be tuned to alter bonding and lattice dynamics (i.e., phonons) and detail how these properties are changed with SRO. Thermodynamically favorable SRO significantly improved the phase stability of fcc Cu3Au from -0.0343 eV-atom -1 to –0.0682 eV-atom -1 . We used our DFT-based linear-response theory to predict SRO and its electronic origin, and accurately estimate the observed transition temperature, ordering instability (L1 2 ), and Warren-Cowley SRO parameters, in agreement with experiments. The accurate prediction of real-space SRO gives an edge over computationally and resource intensive approaches such as monte-carlo methods or experiments, which will enable large scale molecular dynamic simulations by providing supercells with optimized SRO. Here we also analyzed phonon dispersion and estimated the vibrational entropy change (from 9kB at 300 K to 6kB at 100 K) in fcc Cu3Au. We established from SRO analysis that exclusion of chemical interactions may lead to a skewed view of true properties in chemically complex alloys. The first-principles methods described in this work are generally applicable to any arbitrary solid-solution alloys, including multi-principal-element alloys, therefore, holds promise for designing technologically useful materials.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Ab Initio Uncertainty Quantification of Neutrinoless Double-Beta Decay in 76 Ge

The observation of neutrinoless double-beta (0⁢𝜈⁢𝛽⁢𝛽) decay would offer proof of lepton number violation, demonstrating that neutrinos are Majorana particles, while also helping us understand why there is more matter than antimatter in the Universe. If the decay is driven by the exchange of the three known light neutrinos, a discovery would, in addition, link the observed decay rate to the neutrino mass scale through a theoretical quantity known as the nuclear matrix element (NME). Accurate values of the NMEs for all nuclei considered for use in 0⁢𝜈⁢𝛽⁢𝛽 experiments are therefore crucial for designing and interpreting those experiments. Here, we report the first comprehensive ab initio uncertainty quantification of the 0⁢𝜈⁢𝛽⁢𝛽-decay NME, in the key nucleus 76 Ge. Here, our method employs nuclear strong and weak interactions derived within chiral effective field theory and recently developed many-body emulators. Our result, with a conservative treatment of uncertainty, is an NME of 2.60$^{+1.28}_{−1.36}$, which, together with the best-existing half-life sensitivity and phase-space factor, sets an upper limit for effective neutrino mass of 187$^{+205}_{−62}$ meV. The result is important for designing next generation germanium detectors aiming to cover the entire inverted hierarchy region of neutrino masses.

Ab initio calculations↗
DOE OSTI2024

Jacobian-scaled K-means clustering for physics-informed segmentation of reacting flows

This work introduces Jacobian-scaled K-means (JSK-means) clustering, which is a physicsinformed clustering strategy centered on the K-means framework. The method allows for the injection of underlying physical knowledge into the clustering procedure through a distance function modification: instead of leveraging conventional Euclidean distance vectors, the JSKmeans procedure operates on distance vectors scaled by matrices obtained from dynamical system Jacobians evaluated at the cluster centroids. The goal of this work is to show how the JSKmeans algorithm - without modifying the input dataset - produces clusters that capture regions of dynamical similarity, in that the clusters are redistributed towards high-sensitivity regions in phase space and are described by similarity in the source terms of samples instead of the samples themselves. The algorithm is demonstrated on a complex reacting flow simulation dataset (a channel detonation configuration), where the dynamics in the thermochemical composition space are known through the highly nonlinear and stiff Arrhenius-based chemical source terms. Interpretations of cluster partitions in both physical space and composition space reveal how JSK-means shifts clusters produced by standard K-means towards regions of high chemical sensitivity (e.g., towards regions of peak heat release rate near the detonation reaction zone). Furthermore, the findings presented here illustrate the benefits of utilizing Jacobian-scaled distances in clustering techniques, and the JSK-means method in particular displays promising potential for improving former partition-based modeling strategies in reacting flow (and other multi-physics) applications.

Clustering↗
DOE OSTI2023

On formulations for modeling pressurized cracks within phase-field methods for fracture

Over the past few decades, the phase-field method for fracture has seen widespread appeal due to the many benefits associated with its ability to regularize a sharp crack geometry. Along the way, several different models for including the effects of pressure loads on the crack faces have been developed. Further, this work investigates the performance of these models and compares them to a relatively new formulation for incorporating crack-face pressure loads. It is shown how the new formulation can be obtained either by modifying the trial space in the traditional variational principle or by postulating a new functional that is dependent on the rates of the primary variables. The key differences between the new formulation and existing models for pressurized cracks in a phase-field setting are highlighted. Model-based simulations developed with discretized versions of the new formulation and existing models are then used to illustrate the advantages and differences. In order to analyze the results, a domain form of the J-integral is developed for diffuse cracks subjected to pressure loads. Results are presented for a one-dimensional cohesive crack, steady crack growth, and crack nucleation from a pressurized enclosure.

42 ENGINEERING↗
DOE OSTI2022

Global calculation of two-neutrino double- β decay within the finite amplitude method in nuclear density functional theory

Two-neutrino double-beta (2νββ) decay has been used to constrain the neutron-proton part of effective interactions, which in turn is used to compute the nuclear matrix elements for neutrinoless double-beta decay, the observation of which would have important consequences for fundamental physics. We carefully examine 2νββ matrix elements within the proton-neutron quasiparticle random-phase approximation with nuclear energy density functionals. Here we work with functionals that are fit globally to single-beta-decay half-lives and charge-exchange giant-resonance energies, but not to 2νββ half-lives themselves, to evaluate the 2νββ nuclear matrix elements for all important nuclei, including those whose half-lives have not yet been measured. Such a comprehensive evaluation in large model spaces without configuration truncation requires an efficient computational scheme; we employ a double contour integration within the finite amplitude method. The results generally reproduce the nuclear matrix element extracted from half-lives well, without the use of any of those half-lives in the fitting procedure. We present predictions of the matrix elements in a total of 27 nuclei with half-lives that are still unmeasured.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Predictive Design of Hybrid Improper Ferroelectric Double Perovskite Oxides

The computational design of suitable multiferroic double perovskite oxides requires finding materials that exhibit sizable polarization, magnetization, and coupling between them. Oxides with the chemical formula of AA'BB'O 6 with building blocks of ABO 3 single perovskite oxides in centrosymmetric Pnma symmetry are strong candidates that have been reported to satisfy such criteria. The system lowers to noncentrosymmetric, polar P2 1 symmetry if A/A' layered and B/B' rocksalt cation orderings are imposed. A detailed compositional search over a variety of chemical spaces followed by evaluating their polarization may lead to the identification of more of these compounds with ferroelectric ordering. The standard density functional theory practices to estimate polarization within the Berry phase formalism require the systems to be perfectly insulating. The number of compounds that can be evaluated using this method is therefore limited. In this work, we introduce a predictive learning strategy based on importance sampling to build a series of machine learning models using results from first-principles simulations to predict polarization and the corresponding switching barrier. The geometry-driven features related to charge states and cationic radii play key roles in predicting the switching barrier with complementary contributions from the key structural mode-based order parameters. These modes become important to draw reasonable predictions of polarization components from machine learning models. In conclusion, our predictive models identify candidates with high polarizations and low switching barriers from a pool of double perovskite oxides, suitable for future investigation for their potential applications in spintronic devices.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Throughput Optimization of Molybdenum Carbide Nanoparticle Catalysts in a Continuous Flow Reactor Using Design of Experiments

Transition metal carbides (TMCs) have attracted significant attention because of their applications toward a wide range of catalytic transformations. However, the practicality of their synthesis is still limited because of the harsh conditions in which most TMCs are prepared. Recently, a solution-phase synthesis of phase-pure a-MoC1-x nanoparticles was presented. While this synthetic route yielded nanoparticles with exceptional catalytic performance, the reaction parameter space was not explored, and catalyst throughput was not optimized for scale-up. Continuous flow platforms coupled with statistical design of experiments (DoE) can provide a powerful method for understanding the reaction parameter space for optimizations. Here, we demonstrate the use of statistical DoE in tandem with response surface methodology for a parametric screening analysis to optimize the throughput of a MoC1-x nanoparticle synthesis utilizing a millifluidic flow reactor. A full factorial design was implemented to evaluate four input variables (reaction temperature, flow rate, solvent fraction of oleylamine, and precursor concentration) that carry statistically significant effects on three responses (throughput, residence time, and isolated yield). A Doehlert matrix was implemented to investigate each significant variable at a higher number of levels to optimize throughput. Our results give a nonintuitive set of experimental conditions that resulted in an optimized throughput of 2.2 g h-1. This translates to a 50-fold increase in throughput compared to the previously reported batch method. The catalytic performance of the MoC1-x nanoparticles produced under optimized throughput was demonstrated in the CO2 hydrogenation reaction. This DoE screening analysis and throughput optimization of MoC1-x synthesis open the door to an increased feasibility for scale-up.

design of experiments↗
DOE OSTI2021

Natural Charge-Transfer Analysis: Eliminating Spurious Charge-Transfer States in Time-Dependent Density Functional Theory via Diabatization, with Application to Projection-Based Embedding

For many types of vertical excitation energies, linear-response time-dependent density functional theory (LR-TDDFT) offers a useful degree of accuracy combined with unrivaled computational efficiency, although charge-transfer excitation energies are often systematically and dramatically underestimated, especially for large systems and those that contain explicit solvent. As a result, low energy electronic spectra of solution-phase chromophores often contain tens to hundreds of spurious charge-transfer states, making LR-TDDFT needlessly expensive in bulk solution. More nefariously, intensity borrowing by the low-energy charge-transfer states can affect intensities of the valence excitations even if those excitation energies are accurate. At higher excitation energies, it is difficult to distinguish spurious CT states from genuine charge-transfer-to-solvent (CTTS) excitations. In this work, we introduce an automated diabatization scheme that enables fast and effective screening of the CTTS acceptor space in bulk solution. Our procedure introduces the concept of “natural charge-transfer orbitals”, which provide a means to isolate characteristic pairs of orbitals that are most likely to participate in a CTTS excitation. The projection of these orbitals onto solvent-centered virtual orbitals provides a criterion for defining the most important solvent molecules in a given excitation. We apply this method to analyze an ab initio molecular dynamics (MD) trajectory of I-(aq) and report the lowest-energy CTTS band in the absorption spectrum. Our results are in excellent agreement with experimental measurements for bulk I-(aq), and only one-third of the water molecules in the I-(H2O)96 simulation cell need to be described with LR-TDDFT in order to obtain excitation energies that are converged to < 0.1 eV. The tools introduced herein will improve the accuracy, efficiency, and usability of LR-TDDFT in solution-phase environments

Carter-Fenk, Kevin D.↗
DOE OSTI2021

Crystal Structure of Colloidally Prepared Metastable Ag 2 Se Nanocrystals

Structural polymorphism is known for many bulk materials; however, on the nanoscale metastable polymorphs tend to form more readily than in the bulk, and with more structural variety. One such metastable polymorph observed for colloidal Ag 2 Se nanocrystals has traditionally been referred to as the “tetragonal” phase. While there are reports on the chemistry and properties of this metastable polymorph, its crystal structure, and therefore electronic structure, has yet to be determined. We report that an anti-PbCl 2 -like structure type (space group P2 1 /n) more accurately describes the powder X-ray diffraction and X-ray total scattering patterns of colloidal Ag 2 Se nanocrystals prepared by several different methods. Density functional theory (DFT) calculations indicate that this anti-PbCl 2 -like Ag 2 Se polymorph is a dynamically stable, narrow-band-gap semiconductor. Here, the anti-PbCl 2 -like structure of Ag 2 Se is a low-lying metastable polymorph at 5–25 meV/atom above the ground state, depending on the exchange-correlation functional used.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2022

Improving the Accuracy of Atomistic Simulations of the Electrochemical Interface

Atomistic simulation of the electrochemical double layer is an ambitious undertaking, requiring quantum mechanical description of electrons, phase space sampling of liquid electrolytes, and equilibration of electrolytes over nanosecond time scales. All models of electrochemistry make different trade-offs in the approximation of electrons and atomic configurations, from the extremes of classical molecular dynamics of a complete interface with point-charge atoms to correlated electronic structure methods of a single electrode configuration with no dynamics or electrolyte. Here, we review the spectrum of simulation techniques suitable for electrochemistry, focusing on the key approximations and accuracy considerations for each technique. Additionally, we discuss promising approaches, such as enhanced sampling techniques for atomic configurations and computationally efficient beyond density functional theory (DFT) electronic methods, that will push electrochemical simulations beyond the present frontier.

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

Efficient and universal characterization of atomic structures through a topological graph order parameter

Abstract A graph-based order parameter, based on the topology of the graph itself, is introduced for the characterization of atomistic structures. The order parameter is universal to any material/chemical system and is transferable to all structural geometries. Four sets of data are used to validate both the generalizability and accuracy of the algorithm: (1) liquid lithium configurations spanning up to 300 GPa, (2) condensed phases of carbon along with nanotubes and buckyballs at ambient and high temperature, (3) a diverse set of aluminum configurations including surfaces, compressed and expanded lattices, point defects, grain boundaries, liquids, nanoparticles, all at nonzero temperatures, and (4) eleven niobium oxide crystal phases generated with ab initio molecular dynamics. We compare our proposed method to existing, state-of-the-art methods for the cases of aluminum and niobium oxide. Our order parameter uniquely classifies every configuration and outperforms all studied existing methods, opening the door for its use in a multitude of complex application spaces that can require fine structure-level characterization of atomistic graphs.

36 MATERIALS SCIENCE↗
DOE OSTI2020

In situ high-temperature X-ray diffraction study of Sc-doped titanium oxide nanocrystallites

Titanium dioxide is an inexpensive wide-gap highly ionic semiconductor with striking photocatalytic capabilities in several heterogeneous photoredox reactions. A small crystal size is desirable to maximize the surface area, since photocatalytic reactions occur at the surface of a photocatalyst. Presented in this work are the synthesis and microstructural characterization of 4 at.% Sc-doped TiO 2 (4SDT) prepared by water-based co-precipitation. The crystal structure of 4SDT was examined via in situ high-temperature powder X-ray diffraction experiments from 25 to 1200°C. Rietveld analysis revealed single-phase anatase up to 875°C, while at 900°C the anatase-to-rutile phase transformation occurred and at higher temperatures additional reflections of Sc-rich phases (Sc 2 TiO 5 from 975°C and Ti 3 Sc 4 O 12 or Sc 2 O 3 at 1200°C) were observed. Debye function analysis (DFA) was applied to model the total scattering pattern directly in reciprocal space, allowing the reconstruction of Ti vacancies. Both Rietveld and DFA methods were applied to estimate the nanocrystallite size and shape with consistent growth in crystallite size with temperature: an ellipsoid shape with equatorial ~4.7 nm / axial (001) ~6.9 nm at 25°C to equatorial ~27.9 nm / axial (001) ~39.6 nm at 900°C refined by Rietveld analysis, versus a cylinder shape with D a , b = 4.3 nm and size dispersion σ ab = 1.5 nm, L c = 4.9 nm and σ c = 2.3 nm at 25°C to D a , b = 21.4 nm, σ ab = 8.3 nm, L c = 23.9 and σ c = 10.9 nm at 900°C estimated by DFA. The microstructural changes obtained by Rietveld and DFA methods were supported by high-resolution transmission electron microscopy image analysis, as well as by the less direct nitrogen sorption techniques that provide information on the size of non-agglomerated and dense particles. The Ti site-occupancy factor showed a linear increase from 0.6–0.8 at 25°C to unity at 900°C for anatase, and from ~0.7 at 900°C to unity at 1200°C for rutile, via Rietveld analysis and DFA.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Interaction of extended dislocations with nanovoid clusters

Voids of nanoscale dimensions in irradiated metals can act as obstacles to dislocation motion and cause strengthening. In this work, nanovoid strengthening and the influences of void size, void spacing and material properties, such as stacking fault energies, on dislocation bypass mechanisms are investigated using Phase Field Dislocation Dynamics, a three-dimensional mesoscale model that predicts the minimum energy pathway taken by discrete dislocations. A broad range of face centered cubic metals (copper, nickel, silver, rhodium, and platinum) and nanovoid sizes and spacings are treated, altogether spanning void size–to–dislocation stacking fault width ratios from less than unity to ten. Material γ-surfaces, calculated from ab initio methods, are input directly into the formulation. The analysis reveals that the critical bypass stress scales linearly with the linear void fraction, effective isotropic shear modulus, and ratio of the intrinsic to unstable stacking fault energies. With only a few exceptions, the critical stress is controlled by the stress required for the leading partial to impinge the voids (to move within range of the attractive image stress field of the void). When the void diameter is nearly an order of magnitude greater than the stacking fault width, the mechanism determining critical strength shifts to the stress for the dislocation to breakaway after partially cutting the void. Furthermore, this situation corresponds to that treated by line tension models and is realized here for Pt, with a sub-nanometer stacking fault width.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Analytic reconstruction with massive particles: one-loop amplitudes for $0\to \overline{q} qt\overline{t}H$

We present an analytic reconstruction of one-loop amplitudes for the process $0\to \overline{q} qt\overline{t}H$. Our calculation is a novel use of analytic reconstruction, retaining explicit covariance in the massive spin states through the massive spinor-helicity formalism. The analytic reconstruction relies on embedding the massive five-point kinematics in a fully massless eight-point phase space while still building a minimal ansatz directly in the five-point phase space. In order to obtain compact analytic expressions it is necessary to identify suitable partial fraction decompositions and extract common numerator factors, which we achieve through careful inspection of limits in which pairs of denominators vanish. We find that the resulting amplitudes are more numerically efficient than ones computed using automatic methods but that the gains are not as significant as in the massless case, at least at present. The method opens the door to applications at two-loop order, where numerical efficiency and improvements in the reconstruction methodology are more crucial, especially with regards to the number of free parameters in the ansatz.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Isospin 0 and 2 two-pion scattering at physical pion mass using all-to-all propagators with periodic boundary conditions in lattice QCD

A study of two-pion scattering for the isospin channels, 𝐼 = 0 and 𝐼 = 2, using lattice QCD is presented. Möbius domain-wall fermions, on top of the Iwasaki-DSDR gauge action for gluons with periodic boundary conditions, are used for the lattice computations, which are carried out on two ensembles of gauge field configurations generated by the RBC and UKQCD Collaborations with physical masses, inverse lattice spacings of 1.023 and 1.378 GeV, and spatial extents of 𝐿 = 4.63 and 4.58 fm, respectively. The all-to-all propagator method is employed to compute a matrix of correlation functions of two-pion operators. The generalized eigenvalue problem (GEVP) is solved for a matrix of correlation functions to extract phase shifts with multiple states—two pions with a nonzero relative momentum, as well as two pions at rest. Our results for phase shifts for both the 𝐼 = 0 and 𝐼 = 2 channels are consistent with the Roy equation and chiral perturbation theory, though at this preliminary stage our errors for 𝐼 = 0 are large. An important outcome of this work is that we are successful in extracting two-pion excited states, which are useful for studying 𝐾 → 𝜋⁢𝜋 decay, on physical-mass ensembles using the GEVP.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Toward Accelerating Discovery via Physics-Driven and Interactive Multifidelity Bayesian Optimization

Both computational and experimental material discovery bring forth the challenge of exploring multidimensional and often nondifferentiable parameter spaces, such as phase diagrams of Hamiltonians with multiple interactions, composition spaces of combinatorial libraries, processing spaces, and molecular embedding spaces. Often these systems are expensive or time consuming to evaluate a single instance, and hence classical approaches based on exhaustive grid or random search are too data intensive. This resulted in strong interest toward active learning methods such as Bayesian optimization (BO) where the adaptive exploration occurs based on human learning (discovery) objective. However, classical BO is based on a predefined optimization target, and policies balancing exploration and exploitation are purely data driven. In practical settings, the domain expert can pose prior knowledge of the system in the form of partially known physics laws and exploration policies often vary during the experiment. Here, we propose an interactive workflow building on multifidelity BO (MFBO), starting with classical (data-driven) MFBO, then expand to a proposed structured (physics-driven) structured MFBO (sMFBO), and finally extend it to allow human-in-the-loop interactive interactive MFBO (iMFBO) workflows for adaptive and domain expert aligned exploration. These approaches are demonstrated over highly nonsmooth multifidelity simulation data generated from an Ising model, considering spin–spin interaction as parameter space, lattice sizes as fidelity spaces, and the objective as maximizing heat capacity. Detailed analysis and comparison show the impact of physics knowledge injection and real-time human decisions for improved exploration with increased alignment to ground truth. Here, the associated notebooks allow to reproduce the reported analyses and apply them to other systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Evaluation of Modular and Scalable Downstream Passage Systems for Silver American Eels

Implementation of effective downstream passage for silver-phase American Eel (Anguilla rostrata) at hydropower dams using methods that are economically feasible has been a challenge for many project owners. Due to the potential for high mortality rates during turbine passage, resource agencies have focused on mandatory prescriptions that prevent entrainment by physically excluding eels (i.e., narrow trash rack bar spacing) or eliminating turbines as a passage route (i.e., unit shutdowns at night when eel movement is greatest). Both of these options can significantly impact the economic viability of a project. Consequently, there is a need for alternative technologies that can provide safe and effective downstream passage for silver eels without having a major impact on the ability of a project to generate power during the eel downstream migration period. This need is of particular concern now because about 300 hydropower projects within the range of American Eel freshwater habitats are scheduled to be relicensed by the Federal Energy Regulatory Commission (FERC) in the next 15 years

13 HYDRO ENERGY↗