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 865 records · Page 48
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Coupled-cluster based approach for core-level states in condensed phase: Theory and application to different protonated forms of aqueous glycine

A theoretical approach for calculating core-level states in condensed phase is presented. The approach is based on equation-of-motion coupled-cluster theory (EOMCC) and effective fragment potential (EFP) method. By introducing an approximate treatment of double excitations in the EOM-CCSD (EOM-CC with single and double substitutions) ansatz, we address poor convergence issues that are encountered for the core-level states and significantly reduce computational costs. While the approximations introduce relatively large errors in the absolute values of transition energies, the errors are systematic. Consequently, chemical shifts, changes in ionization energies relative to reference systems, are reproduced reasonably well. By using different protonation forms of solvated glycine as a benchmark system, we show that our protocol is capable of reproducing the experimental chemical shifts with a quantitative accuracy. The results demonstrate that chemical shifts are very sensitive to the solvent interactions and that explicit treatment of solvent, such as EFP, is essential for achieving quantitative accuracy.

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

Accelerated steady-state electrostatic particle-in-cell simulation of Langmuir probes

First-principles particle-in-cell (PIC) simulation is a powerful tool for understanding plasma behavior, but this power often comes at great computational expense. Artificially reducing the ion/electron mass ratio is a time-honored practice to reduce simulation costs. Usually, this is a severe approximation. However, for steady-state collisionless, electrostatic (Vlasov–Poisson) systems, the solution with reduced mass ratio can be scaled to the solution for the real mass ratio, with no approximation. This “scaled mass” method, which works with already-existing PIC codes, can reduce the computation time for a large class of electrostatic PIC simulations by the square root of the mass ratio. The particle distributions of the resulting steady state must be trivially rescaled to yield the true distributions, but the self-consistent electrostatic field is independent of the mass ratio. This method is equivalent to “numerical timestepping,” an approach that evolves electron and ion populations with different time steps. Numerical timestepping can be viewed as a special case of the speed-limited PIC (SLPIC) method, which is not restricted to steady-state phenomena. Although the scaled-mass approach is simplest, numerical timestepping and SLPIC more easily generalize to include other effects, such as collisions. The equivalence of these new approaches is demonstrated by applying them to simulate a cylindrical Langmuir probe in electron–argon plasma, speeding up simulation by two orders of magnitude. Methods such as SLPIC can therefore play an invaluable role in interpreting probe measurements by including geometric effects, collisions, secondary emission, and non-Maxwellian distributions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Origin of the injection-dependent emission blueshift and linewidth broadening of III-nitride light-emitting diodes

III-nitride light-emitting diodes (LEDs) exhibit an injection-dependent emission blueshift and linewidth broadening that is severely detrimental to their color purity. By using first-principles multi-scale modeling that accurately captures the competition between polarization-charge screening, phase-space filling, and many-body plasma renormalization, we explain the current-dependent spectral characteristics of polar III-nitride LEDs fabricated with state-of-the-art quantum wells. Our analysis uncovers a fundamental connection between carrier dynamics and the injection-dependent spectral characteristics of light-emitting materials. For example, polar III-nitride LEDs offer poor control over their injection-dependent color purity due to their poor hole transport and slow carrier-recombination dynamics, which forces them to operate at or near degenerate carrier densities. Designs that accelerate carrier recombination and transport and reduce the carrier density required to operate LEDs at a given current density lessen their injection-dependent wavelength shift and linewidth broadening.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Sources of error and methods to improve accuracy in interface state density analysis using quasi-static capacitance–voltage measurements in wide bandgap semiconductors

Characterizing interface trap states in commercial wide bandgap devices using frequency-based measurements requires unconventionally high probing frequencies to account for both fast and slow traps associated with wide bandgap materials. The C - ψs technique has been suggested as a viable quasi-static method for determining the interface trap state densities in wide bandgap systems, but the results are shown to be susceptible to errors in the analysis procedure. Here this work explores the primary sources of errors present in the C - ψs technique using an analytical model that describes the apparent response for wide bandgap MOS capacitor devices. Measurement noise is shown to greatly impact the linear fitting routine of the 1/C$^{*2}_s$ vs Ψ s plot to calibrate the additive constant in the surface potential/gate voltage relationship, and an inexact knowledge of the oxide capacitance is also shown to impede interface trap state analysis near the band edge. In addition, a slight nonlinearity that is typically present throughout the 1/C$^{*2}_s$ vs Ψ s plot hinders the accurate estimation of interface trap densities, which is demonstrated for a fabricated n-SiC MOS capacitor device. Methods are suggested to improve quasi-static analysis, including a novel method to determine an approximate integration constant without relying on a linear fitting routine.

42 ENGINEERING↗
DOE OSTI2025

On finite-dimensional smoothed-particle Hamiltonian reductions of the Vlasov equation

The inclusion of spatial smoothing in finite-dimensional particle-based Hamiltonian reductions of the Vlasov equation and related models is considered. Here, this work investigates the underlying Hamiltonian structure of such smoothed particle-based methods for Hamiltonian systems and the small-scale regularization such methods implicitly make in approximating the continuum theory. In the context of the Vlasov–Poisson equation and other mean-field Lie–Poisson systems, of which Vlasov–Poisson is a special case, smoothing amounts to a convolutive regularization of the Hamiltonian. This regularization may be interpreted as a change of the inner product structure used to identify the dual space in the Lie–Poisson Hamiltonian formulation. In particular, the shape function used for spatial smoothing may be identified as the kernel function of a reproducing kernel Hilbert space whose inner product is used to define the Lie–Poisson Hamiltonian structure. It is likewise possible to introduce smoothing in the Vlasov–Maxwell system, but in this case the Poisson bracket must be modified rather than the Hamiltonian. The smoothing applied to the Vlasov–Maxwell system is incorporated by inserting smoothing in the map from canonical to kinematic coordinates. In the filtered system, the Lorentz force law and the current, the two terms coupling the Vlasov equation with Maxwell’s equations, are spatially smoothed.

Hamiltonian mechanics↗
DOE OSTI2015

Some new function spaces of variable smoothness

A new Besov space of variable smoothness is introduced on which the norm is defined in terms of difference relations. This space is shown to be the trace of a weighted Sobolev space with a weight in the corresponding Muckenhoupt class. Methods of nonlinear spline approximation are applied to derive an atomic decomposition theorem for functions in a Besov space of variable smoothness. A complete description of traces on the hyperplane of a Besov space of variable smoothness and of a weighted Besov space with a weight in the corresponding Muckenhoupt class is given. Bibliography: 27 titles.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

DNS of turbulent flow with hemispherical wall roughness

The present study of the effect of roughness density on the mean flow turbulence parameters is motivated by the need for new generation of boundary conditions for multiphase computational multiphase fluid dynamics (CMFD) models applied to boiling flows. Effect of roughness element density on the turbulent flow in a channel is quantified through direct numerical simulations (DNSs). The Navier--Stokes equations are solved using finite element method and bubbles are approximated as rigid near-hemispherical obstacles at the wall. Six different cases were analysed including channel flow with smooth wall and channel flow with rough wall for five different bubble nucleation site densities. Friction factor and the law of the wall was calculated and compared with the previously published results. Existing correlations for nucleating bubble site density dependency on a wall heat flux were used to obtain a relation between the heat flux and the friction factor, leading to the law of the wall dependency on the heat flux. This separate effect study provides new guidelines on how the heat flux in subcooled boiling regime affects the turbulence behaviour near the wall and guides the computational fluid dynamics model development for boiling two-phase flows.

Mechanics↗
DOE OSTI2017

Radiation damage buildup by athermal defect reactions in nickel and concentrated nickel alloys

We develop a new method using binary collision approximation simulating the Rutherford backscattering spectrometry in channeling conditions (RBS/C) from molecular dynamics atom coordinates of irradiated cells. The approach allows comparing experimental and simulated RBS/C signals as a function of depth without fitting parameters. The simulated RBS/C spectra of irradiated Ni and concentrated solid solution alloys (CSAs, NiFe and NiCoCr) show a good agreement with the experimental results. The good agreement indicates the damage evolution under damage overlap conditions in Ni and CSAs at room temperature is dominated by defect recombination and migration induced by irradiation rather than activated thermally.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Atomic oxygen density determination in the effluent of the COST reference source using in situ effective lifetime measurements in the presence of a liquid interface

Abstract Spatially resolved, absolute densities of atomic oxygen are measured for several helium-based admixtures in the effluent of the COST Reference Microplasma Jet using two-photon absorption laser induced fluorescence (TALIF). Admixtures investigated include a helium-only admixture, four helium/oxygen admixtures (0.1 % , 0.5 % , 0.6 % , and 1.0 % oxygen), and a helium/water admixture (2500 ppm water), chosen to coincide with previously published characterizations of plasma-treated liquid. The atomic oxygen TALIF signal is calibrated for density using the noble gas xenon, which possesses a very similar two-photon excitation and fluorescence scheme. Measurements are conducted for the jet operating in ambient air with both an open effluent and a liquid surface present, allowing for comparison with liquid phase measurements conducted under similar conditions. The presence of a water surface does not appear to alter the background chemistry in the effluent but reduces O densities close to the liquid interface when compared to a similar distance from the nozzle in an open effluent case. This may be the result of a reduction in flow velocity caused by the liquid obstructing the gas flow. Additionally, measurements near the liquid surface revealed a region of atomic oxygen well outside of where the core of the effluent impinges on the liquid. This is likely relevant for applications as it considerably expands the surface area subject to O absorption. Critically, in situ measurements of the effective lifetimes of the laser-excited 3p 3 P J state of atomic oxygen were recorded in the effluent by employing a picosecond (ps) laser and a nanosecond (ns) ICCD. By experimentally determining the contribution from collisional quenching via the in situ effective lifetime measurements, significant improvements in the accuracy of the atomic oxygen density calibration were made, with differences of approximately 30 % from existing methods of estimating quenching rates at atmospheric pressure. Finally, spatially resolved atomic oxygen densities allow for an investigation of O formation and extinction pathways in the effluent and a comparison between admixtures.

laser diagnostics↗
DOE OSTI2020

Stability and response of trapped solitary wave solutions of coupled nonlinear Schrödinger equations in an external, $\mathcal{PT}$- and supersymmetric potential

In this work, we present trapped solitary wave solutions of a coupled nonlinear Schrödinger (NLS) system in 1 + 1 dimensions in the presence of an external, supersymmetric and complex $\mathcal{PT}$-symmetric potential. The Schrödinger system this work focuses on possesses exact solutions whose existence, stability, and spatio-temporal dynamics are investigated by means of analytical and numerical methods. Two different variational approximations are considered where the stability and dynamics of the solitary waves are explored in terms of eight and twelve time-dependent collective coordinates (CCs). We find regions of stability for specific potential choices as well as analytic expressions for the small oscillation frequencies in the CC approximation. Our findings are further supported by performing systematic numerical simulations of the NLS system.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Deblending galaxies with variational autoencoders: A joint multiband, multi-instrument approach

ABSTRACT Blending of galaxies has a major contribution in the systematic error budget of weak-lensing studies, affecting photometric and shape measurements, particularly for ground-based, deep, photometric galaxy surveys, such as the Rubin Observatory Legacy Survey of Space and Time (LSST). Existing deblenders mostly rely on analytic modelling of galaxy profiles and suffer from the lack of flexible yet accurate models. We propose to use generative models based on deep neural networks, namely variational autoencoders (VAE), to learn probabilistic models directly from data. We train a VAE on images of centred, isolated galaxies, which we reuse, as a prior, in a second VAE-like neural network in charge of deblending galaxies. We train our networks on simulated images including six LSST bandpass filters and the visible and near-infrared bands of the Euclid satellite, as our method naturally generalizes to multiple bands and can incorporate data from multiple instruments. We obtain median reconstruction errors on ellipticities and r-band magnitude between ±0.01 and ±0.05, respectively, in most cases, and ellipticity multiplicative bias of 1.6 per cent for blended objects in the optimal configuration. We also study the impact of decentring and prove the method to be robust. This method only requires the approximate centre of each target galaxy, but no assumptions about the number of surrounding objects, pointing to an iterative detection/deblending procedure we leave for future work. Finally, we discuss future challenges about training on real data and obtain encouraging results when applying transfer learning.

Arcelin, Bastien↗
DOE OSTI2022

The luminosity function of ringed galaxies

ABSTRACT We perform an analysis of the luminosity functions (LFs) of two types of ringed galaxies – polar-ring galaxies and collisional ring galaxies – using data from the Sloan Digital Sky Survey (SDSS). Both classes of galaxies were formed as a result of interaction with their environment and they are very rare objects. We constructed LFs of galaxies by different methods and found their approximations by the Schechter function. The luminosity functions of both types of galaxies show a systematic fall-off at low luminosities. The polar structures around bright (Mr ≤ −20m) and red (g − r > +0.8) galaxies are about twice as common as around blue ones. The LF of collisional rings is shifted towards brighter luminosities compared to polar-ring galaxies. We analysed the published data on the ringed galaxies in several deep fields and confirmed the increase in their volume density with redshift: up to z ∼ 1 their density grows as (1 + z)m, where m ≳ 5.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Reconstructing the gravitational field of the local Universe

Tests of gravity at the galaxy scale are in their infancy. As a first step to systematically uncovering the gravitational significance of galaxies, we map three fundamental gravitational variables – the Newtonian potential, acceleration and curvature – over the galaxy environments of the local Universe to a distance of approximately 200 Mpc. Our method combines the contributions from galaxies in an all-sky redshift survey, haloes from an N-body simulation hosting low-luminosity objects, and linear and quasi-linear modes of the density field. We use the ranges of these variables to determine the extent to which galaxies expand the scope of generic tests of gravity and are capable of constraining specific classes of model for which they have special significance. In conclusion, we investigate the improvements afforded by upcoming galaxy surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Ab initio calculations of monopole sum rules: From finite nuclei to infinite nuclear matter

We compute moments of the isoscalar monopole response of 𝑁 = 𝑍 closed-shell nuclei based on chiral nucleon-nucleon plus three-nucleon interactions. We employ the random-phase approximation (RPA) and two ab initio many-body approaches, the in-medium similarity renormalization group (IMSRG) and coupled-cluster theory (CC). In the IMSRG framework, the moments are obtained as ground-state expectation values, whereas in the CC approach, they are evaluated through excited-state calculations. We find good agreement between the IMSRG and CC results across all nuclei studied. RPA provides a reasonable approximation to the correlated methods if the interaction is soft. From the calculated moments, we extract average energies of the monopole response, compute finite-nucleus incompressibilities, and estimate the incompressibility of symmetric nuclear matter by a fit to a leptodermous expansion. Our extrapolated values are lower than those obtained in nuclear-matter calculations with the same interactions, but the values are consistent with phenomenological ranges.

Bonaiti, Francesca [Michigan State Univ., East Lan↗
DOE OSTI2020

Selective acoustic control of photon-mediated qubit-qubit interactions

Quantum technologies such as quantum sensing, quantum imaging, quantum communications, and quantum computing rely on the ability to actively manipulate the quantum state of light and matter. Quantum emitters, such as color centers trapped in solids, are a useful platform for the realization of elementary building blocks (qubits) of quantum information systems. In particular, the modular nature of such solid-state devices opens up the possibility to connect them into quantum networks and create nonclassical states of light shared among many qubits. The function of a quantum network relies on efficient and controllable interactions among individual qubits. In this context, we present a scheme where optically active qubits of differing excitation energies are mutually coupled via a dispersive interaction with a shared mode of an optical cavity. This generally off-resonant interaction is prohibitive of direct exchange of information among the qubits. However, we propose a scheme in which by acoustically modulating the qubit excitation energies it is, in fact, possible to tune to resonance a preselected pair of qubits and thus open a communication channel between them. This method potentially enables fast (approximately nanoseconds) and parallelizable on-demand control of a large number of physical qubits. In this work, we develop an analytical and a numerical theoretical model demonstrating this principle and suggest feasible experimental scenarios to test the theoretical predictions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Capturing many-body correlation effects with quantum and classical computing

Theoretical descriptions of excited states of molecular systems in high-energy regimes are crucial for supporting and driving many experimental efforts at light source facilities. However, capturing their complicated correlation effects requires formalisms that provide a hierarchical infrastructure of approximations. These approximations lead to an increased overhead in classical computing methods and, therefore, decisions regarding the ranking of approximations and the quality of results must be made on purely numerical grounds. The emergence of quantum computing methods has the potential to change this situation. Here, in this study, we demonstrate the efficiency of the quantum phase estimator (QPE) in identifying core-level states relevant to x-ray photoelectron spectroscopy. We compare and validate the QPE predictions with exact diagonalization and real-time equation-of-motion coupled-cluster formulations, which are some of the most accurate methods for states dominated by collective correlation effects.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2016

Strong-field ionization and gauge dependence of nonlocal potentials

Nonlocal potential models have been used in place of the Coulomb potential in the Schrödinger equation as an efficient means of exploring high-field laser-atom interaction in previous works. Although these models have found use in modeling phenomena including photoionization and ejected electron momentum spectra, they are known to break electromagnetic gauge invariance. This paper examines if there is a preferred gauge for the linear field response and photoionization characteristics of nonlocal atomic binding potentials in the length and velocity gauges. Overall, it is found that the length gauge is preferable for a wide range of parameters.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2017

Kinetic-energy release of fragments from electron-impact dissociation of the molecular hydrogen ion and its isotopologues

We calculate here the kinetic-energy release distributions of fragments produced for electron-impact dissociation of the vibrationally excited molecular hydrogen ion ${{\mathrm{H}}_{2}}^{+}$ and its isotopologues ${{\mathrm{D}}_{2}}^{+}$ and ${{\mathrm{T}}_{2}}^{+}$. Here we apply the adiabatic-nuclei convergent close-coupling method and compare results with several different methods, including the ${\delta}$ approximation. Results are presented for a number of dissociative excitation transitions and dissociative ionization as a function of the initial vibrational state of the molecule. We confirm that the square root approximation is a good approximation for the adiabatic-nuclei kinetic-energy release cross sections of ${{\mathrm{H}}_{2}}^{+}$. Agreement with experiment, where available, is good.

74 ATOMIC AND MOLECULAR PHYSICS↗