Search NASA⌕ Search

SEARCH · Search NASA

Results for “Approximation 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 1,135 records · Page 63
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

First-Principles Approach to the Conductance of Covalently Bound Molecular Junctions

In this study, we develop a first-principles approach for accurate conductance calculations of covalently bound molecular junctions. Our approach extends the DFT+$Σ$ method, an approximate GW -based self-energy correction scheme acting on a tractable molecular subspace (based on a gas-phase reference of the same dimension) that corrects level alignment in the junction relative to density functional theory (DFT). We introduce a new extended gas-phase reference system, consisting of the molecule and several lead atoms, whose frontier orbitals maximally project onto the conducting orbitals of the junction. With this choice of reference, our self-energy correction to the Kohn–Sham Hamiltonian takes into account mixing of the gas-phase reference orbitals upon the formation of the junction. We apply our generalized DFT+$Σ$ approach to a series of alkane–chain junctions in which the molecules are covalently bound to the leads via carboxyl terminal groups. Our results lead to conductance values in quantitative agreement with experiment. We also revisit the well-studied Au–bipyridine–Au junction and show that we recover the original DFT+$Σ$ approach for relatively weak donor–acceptor molecule–lead binding.

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

Post-Fermentation Recovery of Biobased Carboxylic Acids

Carboxylic acids are common products produced from the bioconversion of renewable feedstocks. In these processes the separation of the acid product from fermentation broth is the most energy and cost intensive unit operation. Thus, the development of robust, scalable separation approaches that can be applied to a variety of carboxylates is of critical importance to the development of processes that utilize carboxylic acids as platform chemicals. Here we report a batch separation method that includes cell and particulate removal, cation exchange, activated carbon treatment, dewatering with a polymer resin, and product recovery. This method is demonstrated on two unique fermentation broths both derived from corn stover hydrolysate to separate neat succinic and propionic acid. For succinic acid, a crystallization yield of 91% with a product purity of 99.93% was achieved. To our knowledge this is the highest reported crystallization yield and purity for the recovery of succinic acid. Additionally, the method requires approximately 50% less energy compared to standard evaporative crystallization approaches. For propionic acid, neat liquid product was obtained with a distillation yield of 80% and purity of 98%. Finally, these excellent results achieved in terms of yield and purity for succinic and propionic acid, two acids with widely different physical properties, from chemically complex hydrolysate broth demonstrates the effective and robust nature of this approach.

09 BIOMASS FUELS↗
DOE OSTI2022

Accurate temperature diagnostics for matter under extreme conditions

The experimental investigation of matter under extreme densities and temperatures, as in astrophysical objects and nuclear fusion applications, constitutes one of the most active frontiers at the interface of material science, plasma physics, and engineering. The central obstacle is given by the rigorous interpretation of the experimental results, as even the diagnosis of basic parameters like the temperature T is rendered difficult at these extreme conditions. Here, we present a simple, approximation-free method to extract the temperature of arbitrarily complex materials in thermal equilibrium from X-ray Thomson scattering experiments, without the need for any simulations or an explicit deconvolution. Our paradigm can be readily implemented at modern facilities and corresponding experiments will have a profound impact on our understanding of warm dense matter and beyond, and open up a variety of appealing possibilities in the context of thermonuclear fusion, laboratory astrophysics, and related disciplines.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Mapping brucellosis risk in Kenya and its implications for control strategies in sub-Saharan Africa

In Sub-Saharan Africa (SSA), effective brucellosis control is limited, in part, by the lack of long-term commitments by governments to control the disease and the absence of reliable national human and livestock population-based data to inform policies. Therefore, we conducted a study to establish the national prevalence and develop a risk map for Brucella spp. in cattle to contribute to plans to eliminate the disease in Kenya by the year 2040. We randomly generated 268 geolocations and distributed them across Kenya, proportionate to the area of each of the five agroecological zones and the associated cattle population. Cattle herds closest to each selected geolocation were identified for sampling. Up to 25 cattle were sampled per geolocation and a semi-structured questionnaire was administered to their owners. We tested 6,593 cattle samples for Brucella immunoglobulin G (IgG) antibodies using an Enzyme-linked immunosorbent assay (ELISA). We assessed potential risk factors and performed spatial analyses and prevalence mapping using approximate Bayesian inference implemented via the integrated nested Laplace approximation (INLA) method. The national Brucella spp. prevalence was 6.8% (95% CI: 6.2–7.4%). Exposure levels varied significantly between agro-ecological zones, with a high of 8.5% in the very arid zone with the lowest agricultural potential relative to a low of 0.0% in the agro-alpine zone with the highest agricultural potential. Additionally, seroprevalence increased with herd size, and the odds of seropositivity were significantly higher for females and adult animals than for males or calves. Similarly, animals with a history of abortion, or with multiple reproductive syndromes had higher seropositivity than those without. At the herd level, the risk of Brucella spp. transmission was higher in larger herds, and herds with a history of reproductive problems such as abortion, giving birth to weak calves, or having swollen testes. Geographic localities with high Brucella seroprevalence occurred in northern, eastern, and southern regions of Kenya all primarily characterized by semi-arid or arid agro-ecological zones dominated by livestock pastoralism interspersed with vast areas with mixed livestock-wildlife systems. The large spatial extent of our survey provides compelling evidence for the widespread geographical distribution of brucellosis risk across Kenya in a manner easily understandable for policymakers. Our findings can provide a basis for risk-stratified pilot studies aiming to investigate the cost-effectiveness and efficacy of singular and combined preventive intervention strategies that seek to inform Kenya’s Brucellosis Control Policy.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Martensite decomposition during rapid heating of Ti-6Al-4V studied via in situ synchrotron X-ray diffraction

Martensite, α‘, commonly appears in Ti-6Al-4V upon rapid cooling from above the β-transus temperature. It is known that α‘ decomposes into α and β at high temperatures but well below the β-transus temperature. Here, we study the decomposition of martensitic Ti-6Al-4V under rapid laser heating, employing in situ synchrotron X-ray diffraction. A comparison is made with post-annealed Ti-6Al-4V under heating to elucidate changes without martensite decomposition. The fast acquisition of X-ray diffraction data at 250 Hz temporally resolves the decomposition process initiated by annihilating dislocations in α‘. The recovery process is accompanied by structural changes in martensite, followed by the phase transformation to β. Thermal profiles estimated from the lattice parameter data reveal the influence of heating rates and dislocation densities on the decomposition process. Throughout the analysis of the diffraction profiles with respect to estimated temperature, we propose a straightforward method for approximating the initiation temperature of martensite decomposition.

42 ENGINEERING↗
DOE OSTI2025

Electric dipole excitations near the neutron separation energies in 96 Mo

Electric dipole strength near the neutron separation energy significantly impacts nuclear structure properties and astrophysical scenarios. These excitations are complex in nature and may involve the so-called pygmy dipole resonance (PDR). Transition densities play a crucial role in understanding the nature of nuclear excited states, including collective excitations, as well as in constructing transition potentials in DWBA or coupled-channels equations. In this work, we focus on electric dipole excitations in spherical molybdenum isotopes, particularly 96 Mo, employing fully consistent Hartree-Fock-Bogoliubov (HFB) and Quasiparticle Random Phase Approximation (QRPA) methods. We analyze the dipole strength near the neutron separation energy, which represents the threshold for neutron capture processes, and examine the isospin characteristics of PDR states through transition density calculations. Examination of proton and neutron transition densities reveals distinctive features of each dipole state, indicating their isoscalar and isovector nature. We observe that the primary component in the enhanced low-energy region exhibits isovector character. The PDR displays a mixture of isoscalar and isovector nature, distinguishing it from the isovector giant dipole resonance (IVGDR). These findings lay the groundwork for future investigations into the role of transition densities in reaction models and for their application to inelastic scattering calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2015

Accurate thermoelastic tensor and acoustic velocities of NaCl

Despite the importance of thermoelastic properties of minerals in geology and geophysics, their measurement at high pressures and temperatures are still challenging. Thus, ab initio calculations are an essential tool for predicting these properties at extreme conditions. Owing to the approximate description of the exchange-correlation energy, approximations used in calculations of vibrational effects, and numerical/methodological approximations, these methods produce systematic deviations. Hybrid schemes combining experimental data and theoretical results have emerged as a way to reconcile available information and offer more reliable predictions at experimentally inaccessible thermodynamics conditions. Here we introduce a method to improve the calculated thermoelastic tensor by using highly accurate thermal equation of state (EoS). The corrective scheme is general, applicable to crystalline solids with any symmetry, and can produce accurate results at conditions where experimental data may not exist. We apply it to rock-salt-type NaCl, a material whose structural properties have been challenging to describe accurately by standard ab initio methods and whose acoustic/seismic properties are important for the gas and oil industry.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Highly parallelisable simulations of time-dependent viscoplastic fluid flow with structured adaptive mesh refinement

We present the extension of an efficient and highly parallelisable framework for incompressible fluid flow simulations to viscoplastic fluids. The system is governed by incompressible conservation of mass, the Cauchy momentum equation, and a generalised Newtonian constitutive law. In order to simulate a wide range of viscoplastic fluids, we employ the Herschel-Bulkley model for yield-stress fluids with nonlinear stress-strain dependency above the yield limit. We utilise Papanastasiou regularisation in our algorithm to deal with the singularity in apparent viscosity. The resulting system of partial differential equations is solved using the IAMR (Incompressible Adaptive Mesh Refinement) code, which uses second-order Godunov methodology for the advective terms and semi-implicit diffusion in the context of an approximate projection method to solve adaptively refined meshes. By augmenting the IAMR code with the ability to simulate regularised Herschel-Bulkley fluids, we obtain efficient numerical software for time-dependent viscoplastic flow in three dimensions, which can be used to investigate systems not considered previously due to computational expense. We validate results from simulations using this new capability against previously published data for Bingham plastics and power-law fluids in the two-dimensional lid-driven cavity. In doing so, we expand the range of Bingham and Reynolds numbers which have been considered in the benchmark tests. Moreover, extensions to time-dependent flow of Herschel-Bulkley fluids and three spatial dimensions offer new insights into the flow of viscoplastic fluids in this test case, and we provide missing benchmark results for these extensions.

42 ENGINEERING↗
DOE OSTI2019

Simulating conical intersection dynamics in the condensed phase with hybrid quantum master equations

We introduce a framework for simulating relaxation dynamics through a conical intersection of an open quantum system that combines methods to approximate the motion of degrees of freedom with disparate time and energy scales. In the vicinity of a conical intersection, a few degrees of freedom render the nuclear dynamics nonadiabatic with respect to the electronic degrees of freedom. We treat these strongly coupled modes by evolving their wavepacket dynamics in the absence of additional coupling exactly. The remaining weakly coupled nuclear degrees of freedom are partitioned into modes that are fast relative to the nonadiabatic coupling and those that are slow. The fast degrees of freedom can be traced out and treated with second-order perturbation theory in the form of the time-convolutionless master equation. The slow degrees of freedom are assumed to be frozen over the ultrafast relaxation and treated as sources of static disorder. Similarly, we adopt the recently developed frozen-mode extension to second-order quantum master equations. We benchmark this approach to numerically exact results in models of pyrazine internal conversion and rhodopsin photoisomerization. We use this framework to study the dependence of the quantum yield on the reorganization energy and the characteristic time scale of the bath in a two-mode model of photoisomerization. We believe that the yield is monotonically increasing with reorganization energy for a Markovian bath but monotonically decreasing with reorganization energy for a non-Markovian bath. This reflects the subtle interplay between dissipation and decoherence in conical intersection dynamics in the condensed phase.

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

Twisted bilayer graphene as topological heavy fermion: II. Analytical approximations of the model parameters

The recently-introduced topological heavy fermion (THF) model [1] of twisted bilayer graphene (TBG) aims to reconcile the quantum-dot-like electronic structure of the latter observed by scanning tunneling microscopy, with its electron delocalization seen in transport measurements. The THF model achieves this by coupling localized (heavy) fermions with anomalous conduction electrons. Originally, the parameters of the THF model were obtained numerically from the Bistritzer–Macdonald (BM) model of TBG [1]. In this work, we derive analytical expressions for the THF model parameters as a function of the twist angle, the ratio between the tunneling amplitudes at the AA and AB regions (w 0 /w 1 ), and the screening length of the interaction potential. By numerically computing the THF model parameters across an extensive experimentally-relevant parameter space, we show that the resulting approximations are remarkably good, i.e., within the 30% relative error for almost the entire parameter space. At the single-particle level, the THF model accurately captures the energy spectrum of the BM model over a large phase space of angles and tunneling amplitude ratios. When interactions are included, we also show that the THF description of TBG is good around the magic angle for realistic values of the tunneling amplitude ratios (0.6 ≤ w 0 /w 1 ≤ 1.0), for which the hybridization between the localized and conduction fermions γ is smaller than the onsite repulsion of the heavy fermions U 1 (i.e., |γ| < U 1 ).

Approximation methods↗
DOE OSTI2022

Electronic states and transitions of PrO and PrO + probed by threshold ionization spectroscopy and spin–orbit multiconfiguration perturbation theory

The precise ionization energy of praseodymium oxide (PrO) seeded in supersonic molecular beams is measured with mass-analyzed threshold ionization (MATI) spectroscopy. Here, a total of 33 spin–orbit (SO) states of PrO and 23 SO states of PrO + are predicted by second-order multiconfigurational quasi-degenerate perturbation (MCQDPT2) theory. Electronic transitions from four low-energy SO levels of the neutral molecule to the ground state of the singly charged cation are identified by combining the MATI spectroscopic measurements with the MCQDPT2 calculations. The precise ionization energy is used to reassess the ionization energies and the reaction enthalpies of the Pr + O → PrO + + e – chemi-ionization reaction reported in the literature. An empirical formula that uses atomic electronic parameters is proposed to predict the ionization energies of lanthanide monoxides, and the empirical calculations match well with available precise experimental measurements.

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

Nonadiabatic dynamics with classical trajectories: The problem of an initial coherent superposition of electronic states

Advances in coherent light sources and development of pump–probe techniques in recent decades have opened the way to study electronic motion in its natural time scale. When an ultrashort laser pulse interacts with a molecular target, a coherent superposition of electronic states is created and the triggered electron dynamics is coupled to the nuclear motion. A natural and computationally efficient choice to simulate this correlated dynamics is a trajectory-based method where the quantum-mechanical electronic evolution is coupled to a classical-like nuclear dynamics. These methods must approximate the initial correlated electron–nuclear state by associating an initial electronic wavefunction to each classical trajectory in the ensemble. Different possibilities exist that reproduce the initial populations of the exact molecular wavefunction when represented in a basis. We show that different choices yield different dynamics and explore the effect of this choice in Ehrenfest, surface hopping, and exact-factorization-based coupled-trajectory schemes in a one-dimensional two-electronic-state model system that can be solved numerically exactly. Furthermore, this work aims to clarify the problems that standard trajectory-based techniques might have when a coherent superposition of electronic states is created to initialize the dynamics, to discuss what properties and observables are affected by different choices of electronic initial conditions and to point out the importance of quantum-momentum-induced electronic transitions in coupled-trajectory schemes.

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

Introduction to spectral line shape theory

Spectral line-shape models are an important part of understanding high-energy-density (HED) plasmas. Models are needed for calculating opacity of materials and can serve as diagnostics for astrophysical and laboratory plasmas. However, much of the literature on line shapes is directed toward specialists. This perspective makes it difficult for non-specialists to enter the field. We have two broad goals with this topical review. First, we aim to give information so that others in HED physics may better understand the current field. Additionally, this first goal may help guide future experiments to test different aspects of the theory. Second, we provide an introduction for those who might be interested in line-shape theory, and enough materials to be able to navigate the field and the literature. We give a high-level overview of line broadening process, as well as dive into the formalism, available methods, and approximations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Real-time observation of toroidal current redistributions induced by three-dimensional MHD phenomena triggering vertical displacement events in tokamak plasmas

Three-dimensional MHD instabilities, including edge-localized modes (ELMs) and internal reconnection events (IREs), have been observed to precipitate loss of vertical stability in tokamak plasmas, resulting in vertical displacement events (VDEs). This vertical destabilization can occur due to toroidal current redistributions and/or shape changes resulting from these phenomena. Using a recently introduced method for rapidly reconstructing the two-dimensional toroidal plasma current density profile in real-time, results are presented that demonstrate the specific current distribution changes that occur during ELMs (on KSTAR) and IREs (on MAST-U) that lead to loss of vertical control. The method most efficiently reconstructs the toroidal current density profile by doing so on a basis of principal components of historical profiles. These principal components isolate dominant current profile dynamics, improving interpretability, increasing speed, and reducing dimensionality of the profile computation. On KSTAR, this computation is executed in the real-time plasma control system at a rate of 10 kHz (limited by available CPU cycle times), allowing the current profile evolution to be assessed at several times over the course of each ELM event. Further, by incorporating the reconstructions into a novel vertical stability metric, the contribution of specific current profile dynamics to the loss of vertical stability can be assessed in real-time for VDE avoidance and improved understanding of the causal relationship between three-dimensional MHD phenomena and VDEs. The success of this method in approximating toroidal current density profiles from kinetic equilibrium reconstructions is also presented ($R^2=0.990$), along with its capability to produce other equilibrium quantities of interest in real-time at high time resolution.

edge-localized modes↗
DOE OSTI2022

Anisotropic angular scattering models of elastic electron-neutral collisions for Monte Carlo plasma simulations

Many laboratory and industrial plasma applications require accurate modeling techniques to understand the interplay between microscopic and macroscopic processes. A prime example of this interplay is how particle and Monte Carlo (MC) simulation codes describe angular scattering of electrons following elastic scattering events. The forward peaked nature of high energy electron elastic scattering is relatively trivial to accurately describe in plasma simulations. However, for lower energy collisions, which produce near isotropic or backward peaked differential cross sections, there is not a strong consensus among the plasma modeling community on how to best describe these angular scattering trends. Here, in this study, we propose a systematic method to approximate the aforementioned non-trivial angular scattering behavior with a formula that can be readily implemented in particle-in-cell (PIC) and/or MC plasma simulation codes. The present approach is specifically applied to fusion relevant atomic hydrogen and helium, as well as for molecular hydrogen, and results are also applicable to the atomic isotopes and homonuclear molecular isotopologues of these species. Comparisons between the present angular distribution function and benchmark scattering data were used to validate the proposed models. In addition, two-term Boltzmann calculations and PIC direct simulation MC simulations revealed that the proposed angular distribution function is accurate, agreeing very well with benchmark convergent close-coupling scattering calculations, and electron transport measurements. These studies confirmed that the present angular distribution function model can be utilized without the need of renormalization to the momentum transfer cross section (as opposed to using the elastic scattering integrated cross section), which has been suggested by several studies in order to correct for deficient angular scattering models, and to agree with transport measurements. Hence, the present anisotropic angular scattering model can be utilized to accurately model the momentum transfer as well as the electron trajectories of elastic collisions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Phoebe: a high-performance framework for solving phonon and electron Boltzmann transport equations

Understanding the electrical and thermal transport properties of materials is critical to the design of electronics, sensors, and energy conversion devices. Computational modeling can accurately predict material properties but, in order to be reliable, requires accurate descriptions of electron and phonon states and their interactions. While first-principles methods are capable of describing the energy spectrum of each carrier, using them to compute transport properties is still a formidable task, both computationally demanding and memory intensive, requiring integration of fine microscopic scattering details for estimation of macroscopic transport properties. To address this challenge, we present Phoebe—a newly developed software package that includes the effects of electron–phonon, phonon–phonon, boundary, and isotope scattering in computations of electrical and thermal transport properties of materials with a variety of available methods and approximations. This open source C++ code combines MPI-OpenMP hybrid parallelization with GPU acceleration and distributed memory structures to manage computational cost, allowing Phoebe to effectively take advantage of contemporary computing infrastructures. We demonstrate that Phoebe accurately and efficiently predicts a wide range of transport properties, opening avenues for accelerated computational analysis of complex crystals.

36 MATERIALS SCIENCE↗