Search NASA⌕ Search

SEARCH · Search NASA

Results for “GO”

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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

DFTB+, a software package for efficient approximate density functional theory based atomistic simulations

DFTB+ is a versatile community developed open source software package offering fast and efficient methods for carrying out atomistic quantum mechanical simulations. By implementing various methods approximating density functional theory (DFT), such as the density functional based tight binding (DFTB) and the extended tight binding method, it enables simulations of large systems and long timescales with reasonable accuracy while being considerably faster for typical simulations than the respective ab initio methods. Based on the DFTB framework, it additionally offers approximated versions of various DFT extensions including hybrid functionals, time dependent formalism for treating excited systems, electron transport using non-equilibrium Green’s functions, and many more. DFTB+ can be used as a user-friendly standalone application in addition to being embedded into other software packages as a library or acting as a calculation-server accessed by socket communication. We give an overview of the recently developed capabilities of the DFTB+ code, demonstrating with a few use case examples, discuss the strengths and weaknesses of the various features, and also discuss on-going developments and possible future perspectives.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Two-color scattering for the measurement of neutrals at the edge of fusion devices

Laser two-color scattering (TCS) is proposed to detect the neutral species in the edge of fusion devices, namely, tokamaks. Here, TCS uses two wavelengths to probe both the laser Rayleigh scattering and Thomson scattering of the neutral-electron bath, with emphasis on neutral density measurements such as that of hydrogen and deuterium. Modeling of the Rayleigh scattering of tokamak neutral species under various plasma conditions (electron density and temperature) shows that, with an appropriate filtering of the Thomson signal and by going to ultraviolet-region wavelengths, identification of the Rayleigh signal can be achieved. Photon count and signal fractions are calculated in two test cases, one in the midplane region of the National Spherical Torus Experiment and one in the divertor region of DIII-D. An uncertainty analysis and discussion of the feasibility of the TCS diagnostic is also presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Parameter space for magnetization effects in high-energy-density plasmas

Magnetic fields are well known to affect the evolution of fluids via the J × B force, where J is the current density and B is the magnetic field. This force leads to the influence of magnetic fields on hydrodynamics (magnetohydrodynamics). Magnetic fields are often neglected in modeling of high-energy-density plasmas, since J × B is very small compared with the plasma pressure gradients. However, many experiments lie in a separate part of parameter space where the plasma is indirectly affected via magnetization of the heat flux and charged particle transport. This is true even for initially unmagnetized plasmas, since misaligned density and temperature gradients can self-generate magnetic fields. By comparing terms in the induction equation, we go on to estimate the regions of parameter space where these self-generated fields are strong enough to affect the hydrodynamics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Challenges of relaxed n -type GaP on Si and strategies to enable low threading dislocation density

We directly show that doping type strongly affects the threading dislocation density (TDD) of relaxed GaP on Si, with n-type GaP having a TDD of ~3.1 × 10 7 cm -2 , nearly 30× higher than both p-type and unintentionally doped GaP at ~1.1 × 10 6 cm -2 . Such a high TDD is undesirable since n-GaP on Si serves as the starting point for the growth of epitaxial III-V/Si multi-junction solar cells. After highlighting additional challenges for highly n-doped GaP on Si including increased surface roughness, anisotropic strain relaxation, and inhomogeneous TDD distributions from blocking of the dislocation glide, we go on to show that the TDD of n-GaP on Si rises by 10× as the doping concentration increases from ~5 × 10 16 to ~2 × 10 18 cm -3 . Next, we investigate the effects of additional dopant choices on the TDD, determining that electronic effects dominate over solute effects on the dislocation velocity at these concentrations. Finally, we demonstrate the respective roles of compressively strained superlattices, low-temperature initiation, and lowered n-type doping concentration on reducing the TDD for n-GaP on Si. By combining all three, we attain relaxed n-GaP on Si with a TDD of 1.54(±0.20) × 10 6 cm -2 , approaching parity with p-GaP on Si. Such high-quality n-GaP on Si will play an important role in boosting the efficiency of epitaxial III-V/Si multi-junction solar cells.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Charge breeders: Development of diagnostic tools to probe the underlying physics

Charge breeders were developed more than 20 years ago in the context of radioactive ion beam (RIB) production. The main goal is to boost the charge state of a singly charged RIB to a higher value matching the A/Q of a post-accelerator. In that way, the RIB produced at some tens of keV can efficiently be accelerated to energies in a range of few MeV/u up to several tens of MeV/u, which is of interest for nuclear structure and nuclear astrophysics experiments. Two families of charge breeders have been developed: one based on Electron Cyclotron Resonance Ion Sources (ECRISs) suitable for Continuous Wave (CW) operation and another built around Electron Beam Ion Sources (EBISs) used for pulsed post-accelerators. For many years, experimental studies have been carried out to enhance the charge breeding process and improve the extracted beam quality, i.e., purity, intensity, emittance, and time structure. The evolution of the charge breeders is also closely related to the emergence of new facilities delivering even more exotic beams. Diagnostic tools, such as Faraday cups, beam profile monitors, and emittance scanners, are mandatory to tune the machines and evaluate their performances. However, to go beyond that, a better understanding of plasma physics (in ECRIS) and electron-ion interactions (in EBIS), as well as the 1+ beam characteristics, for instance, energy spread transverse emittance, is needed. Therefore, simulation codes modeling the physics phenomena in both type of charge breeders were developed to support the understanding of experimental outputs and to gain insights into non-observable parameters. Nevertheless, more sophisticated diagnostic tools are required to improve the charge breeding performance and to extend the potential applications toward new RIB facilities. This Review will present the two kinds of charge breeders and their technical evolutions. It will review the diagnostic tools and simulation codes employed for operation and evaluation of the internal physics processes.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Electron-only reconnection and associated electron heating and acceleration in PHASMA

In this work, using incoherent Thomson scattering, electron heating and acceleration at the electron velocity distribution function (EVDF) level are investigated during electron-only reconnection in the PHAse Space MApping (PHASMA) facility. Reconnection arises during the merger of two kink-free flux ropes. Both push and pull type reconnection occur in a single discharge. Electron heating is localized around the separatrix, and the electron temperature increases continuously along the separatrix with distance from the X-line. The local measured gain in enthalpy flux is up to 70% of the incoming Poynting flux. Notably, non-Maxwellian EVDFs comprised of a warm bulk population and a cold beam are directly measured during the electron-only reconnection. The electron beam velocity is comparable to, and scales with, electron Alfvén speed, revealing the signature of electron acceleration caused by electron-only reconnection. The observation of oppositely directed electron beams on either side of the X-point provides “smoking-gun” evidence of the occurrence of electron-only reconnection in PHASMA. 2D particle-in-cell simulations agree well with the laboratory measurements. The measured conversion of Poynting flux into electron enthalpy is consistent with recent observations of electron-only reconnection in the magnetosheath [Phan et al., Nature 557, 202 (2018)] at similar dimensionless parameters as in the experiments. The laboratory measurements go beyond the magnetosheath observations by directly resolving the electron temperature gain.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Stochastic evaluation of four-component relativistic second-order many-body perturbation energies: A potentially quadratic-scaling correlation method

A second-order many-body perturbation correction to the relativistic Dirac-Hartree-Fock energy is evaluated stochastically by integrating 13-dimensional products of four-component spinors and Coulomb potentials. The integration in the real space of electron coordinates is carried out by the Monte Carlo (MC) method with the Metropolis sampling, whereas the MC integration in the imaginary-time domain is performed by the inverse-CDF (cumulative distribution function) method. The computational cost to reach a given relative statistical error for spatially compact but heavy molecules is observed to be no worse than cubic and possibly quadratic with the number of electrons or basis functions. This is a vast improvement over the quintic scaling of the conventional, deterministic second-order many-body perturbation method. The algorithm is also easily and efficiently parallelized with demonstrated 92% strong scalability going from 64 to 4096 processors for a fixed job size.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

A review on recent machine learning applications for imaging mass spectrometry studies

Imaging mass spectrometry (IMS) is a powerful analytical technique widely used in biology, chemistry, and materials science fields that continue to expand. IMS provides a qualitative compositional analysis and spatial mapping with high chemical specificity. The spatial mapping information can be 2D or 3D depending on the analysis technique employed. Due to the combination of complex mass spectra coupled with spatial information, large high-dimensional datasets (hyperspectral) are often produced. Therefore, the use of automated computational methods for an exploratory analysis is highly beneficial. The fast-paced development of artificial intelligence (AI) and machine learning (ML) tools has received significant attention in recent years. These tools, in principle, can enable the unification of data collection and analysis into a single pipeline to make sampling and analysis decisions on the go. There are various ML approaches that have been applied to IMS data over the last decade. Here, in this review, we discuss recent examples of the common unsupervised (principal component analysis, non-negative matrix factorization, k-means clustering, uniform manifold approximation and projection), supervised (random forest, logistic regression, XGboost, support vector machine), and other methods applied to various IMS datasets in the past five years. The information from this review will be useful for specialists from both IMS and ML fields since it summarizes current and representative studies of computational ML-based exploratory methods for IMS.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Quantum water desalination: Water generation through separate pathways for protons and hydroxide ions in membranes

Much of the water desalination strategies has focused on designing pores and membranes that transport water and reject ions and other molecules at a high rate. In this paper, we discuss an approach where protons (H+) and hydroxide (OH−) ions are transported via different mechanisms through a porous membrane, and subsequently, once they have been transported through the membrane, they recombine to generate water. 2D materials such as graphene and MoS2 have generated significant interest for applications such as desalination. Here, we explore the applicability of one such 2D material—a cubic Ti2C MXene membrane—in desalination by creating a OH− ion selective pore, which significantly suppresses protons but allows OH− ions and water to go through. The catalytic properties of MXenes enable the dissociation of water on the surface, and the dissociated protons translocate through the membrane via quantum-dominated phenomena such as hopping from interstitial-to-interstitial. OH− ions translocate through a positively charged pore and recombine with protons on the other side of the membrane to form water. Our results indicate that water molecules generated via quantum processes can significantly enhance the overall transport of water across the membrane.

Physics↗
DOE OSTI2023

Correlated Dirac–Coulomb–Breit multiconfigurational self-consistent-field methods

The fully correlated frequency-independent Dirac–Coulomb–Breit Hamiltonian provides the most accurate description of electron–electron interaction before going to a genuine relativistic quantum electrodynamics theory of many-electron systems. In this work, we introduce a correlated Dirac–Coulomb–Breit multiconfigurational self-consistent-field method within the frameworks of complete active space and density matrix renormalization group. In this approach, the Dirac–Coulomb–Breit Hamiltonian is included variationally in both the mean-field and correlated electron treatment. Here, we also analyze the importance of the Breit operator in electron correlation and the rotation between the positive- and negative-orbital space in the no-virtual-pair approximation. Atomic fine-structure splittings and lanthanide contraction in diatomic fluorides are used as benchmark studies to understand the contribution from the Breit correlation.

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

Polarization consistent dielectric screening in polarizable continuum model calculations of solvation energies

A polarization consistent framework, where dielectric screening is affected consistently in polarizable continuum model (PCM) calculations, is employed for the study of solvation energies. The computational framework combines a screened range-separated-hybrid functional (SRSH) with PCM calculations, SRSH-PCM, where dielectric screening is imposed in both PCM self-consistent reaction field (SCRF) iterations and the electronic structure Hamiltonian. We begin by demonstrating the impact of modifying the Hamiltonian to include such dielectric screening in SCRF iterations by considering the solutions of electrostatically embedded Hartree–Fock (HF) exact exchange equations. Long-range screened HF-PCM calculations are shown to capture properly the linear dependence of gap energy of frontier orbitals on the inverse of the dielectric constant, whereas unscreened HF-PCM orbital energies are fallaciously semi-constant with respect to the dielectric constant and, therefore, inconsistent with the ionization energy gaps. Similar trends affect density functional theory (DFT) calculations that aim to achieve predictive quality. Importantly, the dielectric screened calculations are shown to significantly affect DFT- and HF PCM-based solvation energies, where screened solvation energies are smaller compared to the unscreened values. Importantly, SRSH-PCM, therefore, appears to reduce the tendency of DFT-PCM to overestimate solvation energies, where we find the effect to increase with the dielectric constant and the polarity of the molecular solute, trends that enhance the quality of DFT-PCM calculations of solvation energy. Understanding the relationship of dielectric screening in the Hamiltonian and DFT-PCM calculations can ultimately benefit on-going efforts for the design of predictive and parameter free descriptions of solvation energies.

Chemistry↗
DOE OSTI2023

Time dependent model for non-inductive ECH X-I current ramp-up for SHPD tokamak facility

Non-inductive (NI) plasma current start-up and ramp-up is an important research topic for spherical tokamak (ST) based reactors and fusion pilot plant (FPP). For a compact FPP, the OH flux availability is highly restricted due to its compact geometry. Efficient fundamental extraordinary mode (X-I) electron cyclotron heating (ECH) current start-up and ramp-up regime was identified for a reactor-like high toroidal magnetic field range which has more than a hundred times higher current drive efficiency compared to more conventional ECH methods for the relevant start-up temperature range. High current drive efficiency is possible due to the strong X-I fundamental ECH interaction only with unidirectional passing electrons constrained by the wave accessibility conditions. Here, we extend the X-I electron cyclotron current drive (ECCD) investigation to a time dependent model to simulate the non-inductive current ramp-up to 10 MA for the Sustained High-Power Density (SHPD) facility. As the X-I ECCD driven current I EC rises, due to the back EMF driven negative current, the net plasma current Ip rises more slowly with the current resistive time scale. For tokamak confinement time (both L-mode and H-Mode) which tends to rise with I p , a positive feedback results and even with constant applied ECH power, T e0 , I EC , and I p can continue to rise to very high values. However, in a realistic situation, the T e0 rise should saturate due to a number of factors such as enhanced core radiation and increased power loss at high temperature. To simulate this effect, we adopt a maximum T e0 model which would limit the temperature rise to certain T e0 . With this model, we investigated the current ramp-up for various ECH power levels and the maximum T e0 of 15, 20, and 25 keV. We find that while the power required is reduced with increasing T e0 limit due to increased current drive efficiency, the time to reach 10 MA tends to go up due to the reduced plasma resistivity for the higher T e0 limit. We also find for a given Te0 limit, the time to reach 10 MA tends to be reduced by increasing the applied ECH power by over driving the current ramp-up where I EC is driven at significantly higher level than 10 MA. While the current ramp-up time may not be an issue for the steady-state reactor systems, if it is desirable to minimize the current ramp-up time, it is prudent to have a sufficient ECH power for current over-drive and have some T e0 limiting tools such as impurity seeding for enhanced radiation. In conclusion, a well-controlled NI ECH start-up also has a potential of improving the tokamak start-up reliability and avoid run-away electrons while the NI off-axis current drive could enhance MHD stability and plasma performance improvements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Stability of Kuramoto networks subject to large and small fluctuations from heterogeneous and spatially correlated noise

Oscillatory networks subjected to noise are broadly used to model physical and technological systems. Due to their nonlinear coupling, such networks typically have multiple stable and unstable states that a network might visit due to noise. In this article, we focus on the assessment of fluctuations resulting from heterogeneous and spatially correlated noise inputs on Kuramoto model networks. We evaluate the typical, small fluctuations near synchronized states and connect the network variance to the overlap between stable modes of synchronization and the input noise covariance. Going beyond small to large fluctuations, we introduce the indicator mode approximation that projects the dynamics onto a single amplitude dimension. Such an approximation allows for estimating rates of fluctuations to saddle instabilities, resulting in phase slips between connected oscillators. Statistics for both regimes are quantified in terms of effective noise amplitudes that are compared and contrasted for several noise models. In conclusion, bridging the gap between small and large fluctuations, we show that a larger network variance does not necessarily lead to higher rates of large fluctuations.

42 ENGINEERING↗
DOE OSTI2024

Learning the stable and metastable phase diagram to accelerate the discovery of metastable phases of boron

Boron, an element of captivating chemical intricacy, has been surrounded by controversies ever since its discovery in 1808. The complexities of boron stem from its unique position between metals and insulators in the Periodic Table. Recent computational studies have shed light on some of the stable boron allotropes. However, the demand for multifunctionality necessitates the need to go beyond the stable phases into the realm of metastability and explore the potentially vast but elusive metastable phases of boron. Traditional search for stable phases of materials has focused on identifying materials with the lowest enthalpy. Here, we introduce a workflow that uses reinforcement learning coupled with decision trees, such as Monte Carlo tree search, to search for stable and metastable boron phases, with enthalpy as the objective. We discover new boron metastable phases and construct a phase diagram that locates their phase space (T, P) at different levels of metastability (ΔG) from the ground state and provides useful information on the domains of relative stability of the various stable and metastable boron phases.

42 ENGINEERING↗
DOE OSTI2024

A method for examining ensemble averaging forms during the transition to turbulence in HED systems for application to RANS models

This paper discusses a strategy to initialize a two-dimensional (2D) Reynolds-averaged Navier–Stokes model [LANL's Besnard–Harlow–Rauenzahn (BHR) model] in order to describe an unsteady transitional Richtmyer–Meshkov (RM)-induced flow observed in on-going high-energy-density ensemble experiments performed on the OMEGA-EP facility. The experiments consist of a nominal single-mode perturbation (initial amplitude a 0 ≈ 10 and wavelength $λ$ = 100μm) with target-to-target variations in the surface roughness subjected to the RM instability with delayed Rayleigh–Taylor in a heavy-to-light configuration. Our strategy leverages high-resolution three-dimensional (3D) implicit large eddy simulations (ILES) simulations to initialize BHR-relevant parameters and subsequently validate the 2D BHR results against the 3D ILES simulations. A suite of five 3D ILES simulations corresponding to five experimental target profiles is undertaken to generate an ensemble dataset. Using ensemble averages from the 3D simulations to initialize the turbulent kinetic energy in the BHR model ( K 0 ) demonstrates the ability of the model to predict the time evolution of the interface as well as the density-specific-volume covariance, b . To quantify the sensitivity of the BHR results to the choice of K 0 and the initial turbulent length scale, S 0 , we execute a parameter sweep spanning four orders of magnitude for both S 0 and K 0 , generating a parameter space consisting of 26 simulations. The Pearson's correlation coefficient is used as a measure of discrepancy between the 2D BHR and 3D ILES simulations and reveals that the ranges 8≲S 0 ≲20 μm and 10 9 ≲K 0 ≲10 10 cm 2 /s 2 produce predictions that agree best with the 3D ILES results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Symmetry breaking and self-interaction correction in the chromium atom and dimer

Density functional approximations to the exchange–correlation energy can often identify strongly correlated systems and estimate their energetics through energy-minimizing symmetry-breaking. In particular, the binding energy curve of the strongly correlated chromium dimer is described qualitatively by the local spin density approximation (LSDA) and almost quantitatively by the Perdew–Burke–Ernzerhof generalized gradient approximation (PBE-GGA), where the symmetry breaking is antiferromagnetic for both. Here, we show that a full Perdew–Zunger self-interaction-correction (SIC) to LSDA seems to go too far by creating an unphysical symmetry-broken state, with effectively zero magnetic moment but non-zero spin density on each atom, which lies ~4 eV below the antiferromagnetic solution. A similar symmetry-breaking, observed in the atom, better corresponds to the 3d↑↑4s↑3d↓↓4s↓ configuration than to the standard 3d↑↑↑↑↑4s↑. For this new solution, the total energy of the dimer at its observed bond length is higher than that of the separated atoms. In conclusion, these results can be regarded as qualitative evidence that the SIC needs to be scaled down in many-electron regions.

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

Experimental evidence of enhanced radial transport in small ELM regimes at DIII-D

Small/type-II edge-localized-modes (ELMs), carrying 1% of the plasma stored energy, are found to deposit only 45 ± 5% of the ELM power near the strike point, and the remaining 55 ± 5% to the far scrape-off-layer (SOL). Small ELMs spread their power over a larger area compared to type-I ELMs, where such a ratio is about 60% and 40% for near- and far-SOL regions, respectively. The larger spread is reflected in the heat flux width (λq) in the SOL for the intra-small ELMs profile of 6.0 mm, almost a factor 2 larger than that of type-I ELMs of 3.15 mm, for similar plasma conditions and magnetic configuration. At the ELM peak, the small ELMs λq is found to be up to 4 times larger than for the type-I ELMs, going from 2 to 7.9 mm, indicating enhanced radial transport in the neon-seeded small ELM scenario. Inter-ELM λqs have been also calculated at the secondary outer divertor in quasi-double-null (QDN) discharges. It is found that, on average, λq is 2.2 times larger in the high-separatrix-density small ELM regime, compared to a reference type-I ELM one. These findings are supported by small ELMs radial velocity profiles, measured at the outer midplane with a fast reciprocating probe, showing a decay length (λvr) in the SOL of 12.8 cm, which is 3.3 times larger than that for the type-I ELMs of 3.9 cm. This analysis shows that small ELMs, although attractive for future machines due to low peak heat flux and large λq, might be of concern for the larger flux to the outer wall.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

A static quantum embedding scheme based on coupled cluster theory

Here, we develop a static quantum embedding scheme that utilizes different levels of approximations to coupled cluster (CC) theory for an active fragment region and its environment. To reduce the computational cost, we solve the local fragment problem using a high-level CC method and address the environment problem with a lower-level Møller–Plesset (MP) perturbative method. This embedding approach inherits many conceptual developments from the hybrid second-order Møller–Plesset (MP2) and CC works by Nooijen [J. Chem. Phys. 111, 10815 (1999)] and Bochevarov and Sherrill [J. Chem. Phys. 122, 234110 (2005)]. We go beyond those works here by primarily targeting a specific localized fragment of a molecule and also introducing an alternative mechanism to relax the environment within this framework. We will call this approach MP-CC. We demonstrate the effectiveness of MP-CC on several potential energy curves and a set of thermochemical reaction energies, using CC with singles and doubles as the fragment solver, and MP2-like treatments of the environment. The results are substantially improved by the inclusion of orbital relaxation in the environment. Using localized bonds as the active fragment, we also report results for N=N bond breaking in azomethane and for the central C–C bond torsion in butadiene. We find that when the fragment Hilbert space size remains fixed (e.g., when determined by an intrinsic atomic orbital approach), the method achieves comparable accuracy with both a small and a large basis set. Additionally, our results indicate that increasing the fragment Hilbert space size systematically enhances the accuracy of observables, approaching the precision of the full CC solver.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗