Search NASA⌕ Search

SEARCH · Search NASA

Results for “first principles simulations”

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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Surface Transfer Doping in MoO 3– x /Hydrogenated Diamond Heterostructure

Surface transfer doping is proposed to be a potential solution for doping diamond, which is hard to dope for applications in high-power electronics. While MoO 3 is found to be an effective surface electron acceptor for hydrogen-terminated diamond with a negative electron affinity, the effects of commonly existing oxygen vacancies remain elusive. We have performed reactive molecular dynamics simulations to study the deposition of MoO 3–x on a hydrogenated diamond (111) surface and used first-principles calculations based on density functional theory to investigate the electronic structures and charge transfer mechanisms. We find that MoO 3–x is an effective surface electron acceptor and the spatial extent of doped holes in hydrogenated diamond is extended, promoting excellent transport properties. Charge transfer is found to monotonically decrease with the level of oxygen vacancy, providing guidance for engineering of the surface transfer doping process.

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

Understanding (La,Sr)(Co,Fe)O$_{3-δ}$ Phase Instability within SOECs Using a Combined Experimental and Atomistic Modeling Approach

Understanding the onset of degradation in the air electrode within solid oxide electrolysis cells (SOECs), and the subsequent impact on cell performance, is a critical step in mitigating the performance losses and stability issues of SOECs. In an effort to identify early onset degradation phenomena, SOECs were characterized as fabricated and after testing potentiostatically at 1.3 V for 1000 h at 750 °C. SOEC air electrodes composed of a 1:1 composite of La 0.6 Sr 0.4 Co 0.2 Fe 0.8 O 3–δ (6428-LSCF) and Gd 0.1 Ce 0.9 O 1.95 (GDC) were studied using synchrotron X-ray diffraction (XRD), scanning transmission electron microscopy coupled with energy dispersive X-ray spectroscopy (STEM-EDS), and X-ray absorption near-edge spectroscopy (XANES) to evaluate the changes in the air electrode structurally and chemically. These techniques show the migration of Sr species from the air electrode through pores in the GDC barrier layer, progressing to the electrolyte boundary, where it accumulates and reacts with (Zr 0.84 Y 0.16 )O 2–δ (YSZ) to form SrZrO 3 . Microscopy results are paired with atomistic simulations to better understand the relationship between the thermodynamic instability of 6428-LSCF and cell fabrication/testing conditions. First-principles calculations reveal that LSCF-6428 is not stable during cell manufacturing and testing conditions, which supports the experimental identification of secondary phases in both as-fabricated and tested cells. Together, these results demonstrate that the challenging environments encountered by SOECs during cell manufacturing and operation lead to instabilities of the target 6428-LSCF anode material and underscore the need for more durable, high-performing SOEC components.

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

Nuclear Theory from First Principles to Forefront Experiments (Final Report for DE-SC0018638)

The Lee research group is a part of the Nuclear Lattice Effective Field Theory Collaboration, which has developed and performed ab initio lattice simulations of nuclear structure, scattering, reactions, and many-body systems. The other senior members of the collaboration include Ulf-G. Meißner at Bonn/Julich, Evgeny Epelbaum and Hermann Krebs at Bochum, Timo Lahde and Thomas Luu at Julich, and Gautam Rupak at Mississippi State. Our letter “Ab initio alpha-alpha scattering” was featured in a Nature News and Views article. Another letter “Nuclear binding near a quantum phase transition” was highlighted in a Viewpoint article in the online APS journal Physics as well as a news article in the IOP publication Physics World (September 21, 2016). Our letter “Eigenvector continuation with subspace learning” was also highlighted a Synopsis article in Physics. There are many promising ab initio approaches being used to calculate the properties of few-and many-nucleon systems. This includes the no-core shell model, symmetry-adapted no-core shell model quantum Monte Carlo, auxiliary-field diffusion Monte Carlo, self-consistent Green’s functions, many-body perturbation theory, in-medium similarity renormalization group, and coupled cluster methods. Lattice effective field theory is another ab initio approach which combines the framework of effective field theory with lattice Monte Carlo methods to allow favorable scaling from few- to many-body systems. Perhaps the most important aspect of lattice effective field theory is that its strengths and weaknesses are orthogonal to that of other ab initio methods. For example, lattice effective field theory has little difficulty in probing cluster structures and collectivity, while such features are much more difficult using other methods. Furthermore it can be used to compute superfluid condensate fractions as well as the phase diagram of strongly matter and the density and temperature dependence of clustering. Lattice effective field theory was first used in simulations of infinite nuclear matter and infinite neutron matter with pions and without pions. In addition to the efforts by our collaboration, there have been recent efforts by other groups as well.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Effects of cluster expansion on the locations of phase transition boundary as a first step to quantify uncertainty in first principles statistical mechanics framework

Predicting phase diagrams from first principle calculations eliminates the need of tedious experimental trials and errors. Fully automating first principle phase diagram calculations without any sort of human intervention has been a long daunting task and troubling scientific communities for decades. This grand problem remains not fully resolved, largely due to the vastly high-dimensional parameter space associated with density functional theory, cluster expansion, lattice Monte Carlo, and the substantial uncertainty propagating through a set of complex simulations. As a first step to tackle this grand problem, we reported a first demonstration of how sensitive phase boundary locations can be to various cluster expansion fittings and input DFT training data. To the best knowledge of the authors, this study reported the first ever attempt to quantify uncertainty in first principles statistical mechanics framework. In addition, a semi-automated phase transition detection algorithm has been devised in this paper to deal with the associated statistical errors and uncertainties from Monte Carlo method and its predecessor DFT calculations and cluster expansions. This algorithm has been applied in a classical cluster expansion Mg-Cd binary alloy system to detect phase transitions at various locations of phase diagram using different cluster expansions to demonstrate its predictive power and quantify uncertainties. The results suggested that using Chebyshev basis function shifted the transition locations toward the dilute solid solution phase and expanded the phase stability of concentrated ordered phases to wider composition range. The addition of perturbed defect configurations into training data set lowered the order-disorder transition temperature to be away from true transition temperature, suggesting the transition nature is indeed configurational disorder dominated rather than defect assisted. Finally, we found that the weighting has negligible effect on the transition locations, except for the case of using Chebyshev basis function to fit all non-weighted configurations that can be susceptible to Monte Carlo sampling hysteresis.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Change of Electrical and Transport Properties of Nickel Oxide by Carrier Concentration and Temperature through First-Principle Calculations

Abstract Nickel is typically used as one of the main components in electrical contact devices or connectors. Nickel oxide (NiO) is usually formed on the surfaces of electrodes and can negatively impact system performance by introducing electrical contact resistance. The thermal, electrical, and transport properties of NiO, as a Mott insulator or a p-type semiconductor, can be altered by operating and environmental conditions such as temperature and stress/strain by contact. In this study, we investigate the fundamental material properties of NiO through the first-principle calculations. First, we obtain and compare the lattice parameter, magnetic moment, and electronic structure for NiO via the WIEN2K simulations with four different potentials (i.e., GGA, GGA + U, LSDA, and LSDA + U). Then, using the WIEN2K simulation results with LSDA + U potential that produces a highly accurate bandgap for NiO, we calculate the electrical conductivity and electrical part of the thermal conductivity of nickel and NiO as a function of temperature and carrier concentration through the BoltzTraP simulations. Systematic simulation results revealed that the electrical conductivity relative to the relaxation time for NiO increases with the carrier concentration, while it shows a slightly decreasing trend with temperature under a fixed carrier concentration. By contrast, the electrical part of the thermal conductivity shows an increasing trend considering carrier concentration and temperature.

Shah, Nayem Md. Reza↗
DOE OSTI2026

Characterizing Defect Dynamics in Silicon Carbide Using Symmetry-Adapted Collective Variables and Machine Learning Interatomic Potentials

Silicon carbide (SiC) divacancies are attractive candidates for spin-defect qubits possessing long coherence times and optical addressability. The high activation barriers associated with SiC defect formation and motion pose challenges for their study by first-principles molecular dynamics. In this work, we develop and deploy machine learning interatomic potentials (MLIPs) to accelerate defect dynamics simulations while retaining ab initio accuracy. We employ an active learning strategy comprising symmetry-adapted collective variable discovery and enhanced sampling to compile configurationally diverse training data, calculation of energies and forces using density functional theory (DFT), and training of an E(3)-equivariant MLIP based on the Allegro model. Here, the trained MLIP reproduces DFT-level accuracy in defect transition activation free energy barriers, enables the efficient and stable simulation of multidefect 216-atom supercells, and permits an analysis of the temperature dependence of defect thermodynamic stability and formation/annihilation kinetics to propose an optimal annealing temperature to maximally stabilize VV divacancies.

Computer simulations↗
DOE OSTI2026

Ab Initio Polariton Spectra of ZnTPP Molecules Collectively Coupled Inside an Optical Cavity

Exciton-polaritons are quasi-particles formed by the quantum mechanical hybridization of electronic and photonic excitations. Despite extensive investigations, a fundamental understanding of molecular polariton spectra and the polariton delocalization from an ab initio theoretical perspective remains elusive. We simulate experimentally measured linear transmission spectroscopy of many Zinc(II) tetraphenylporphyrin (ZnTPP) molecules collectively coupled to a cavity from first principles. Our theoretical approach incorporates many low-lying electronic excitations in ZnTPP molecules, as well as collective light-matter couplings between ZnTPP and the quantized radiation modes, both of which are shown to be the key to accurately recovering the experimental spectra. We further analyzed to what extent the polariton and dark states are delocalized over many molecules, for the first time, using fully ab initio descriptions of the molecules. We finally investigate the line width as a function of detuning, providing new theoretical insights into the experimentally observed motional narrowing behavior. Our work presents first-ofits- kind theoretical studies on molecular polariton spectra, offering a new perspective on molecular polariton formation in realistic ab initio molecular systems whose rich, many-state nature provides spectral features enabled by the high density of electronic states beyond simple quantum optics models.

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

Toward exascale whole-device modeling of fusion devices: Porting the GENE gyrokinetic microturbulence code to GPU

GENE solves the five-dimensional gyrokinetic equations to simulate the development and evolution of plasma microturbulence in magnetic fusion devices. The plasma model used is close to first principles and computationally very expensive to solve in the relevant physical regimes. In order to use the emerging computational capabilities to gain new physics insights, several new numerical and computational developments are required. Here, we focus on the fact that it is crucial to efficiently utilize GPUs (graphics processing units) that provide the vast majority of the computational power on such systems. In this paper, we describe the various porting approaches considered and given the constraints of the GENE code and its development model, justify the decisions made, and describe the path taken in porting GENE to GPUs. We introduce a novel library called gtensor that was developed along the way to support the process. Performance results are presented for the ported code, which in a single node of the Summit supercomputer achieves a speed-up of almost 15× compared to running on central processing unit (CPU) only. Typical GPU kernels are memory-bound, achieving about 90% of peak. Our analysis shows that there is still room for improvement if we can refactor/fuse kernels to achieve higher arithmetic intensity. We also performed a weak parallel scalability study, which shows that the code runs well on a massively parallel system, but communication costs start becoming a significant bottleneck.

Germaschewski, K. (ORCID:0000000284956354)↗
DOE OSTI2024

High-pressure characterization of Ag 3 AuTe 2 : Implications for strain-induced band tuning

Recent band structure calculations have suggested the potential for band tuning in the chiral semiconductor Ag 3 AuTe 2 to zero upon application of negative strain. In this study, we report on the synthesis of polycrystalline Ag 3 AuTe 2 and investigate its transport and optical properties and mechanical compressibility. Transport measurements reveal the semiconducting behavior of Ag 3 AuTe 2 with high resistivity and an activation energy E a of 0.2 eV. The optical bandgap determined by diffuse reflectance measurements is about three times wider than the experimental E a ⁠. Despite the difference, both experimental gaps fall within the range of predicted bandgaps by our first-principles density functional theory (DFT) calculations employing the Perdew–Burke–Ernzerhof and modified Becke–Johnson methods. Furthermore, our DFT simulations predict a progressive narrowing of the bandgap under compressive strain, with a full closure expected at a strain of –4% relative to the lattice parameter. To evaluate the feasibility of gap tunability at such substantial strain, the high-pressure behavior of Ag 3 AuTe 2 was investigated by in situ high-pressure x-ray diffraction up to 47 GPa. Mechanical compression beyond 4% resulted in a pressure-induced structural transformation, indicating the possibility of substantial gap modulation under extreme compression conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Measurement of the small-scale 3D Lyman- α forest power spectrum

Small-scale correlations measured in the Lyman-α (Lyα) forest encode information about the intergalactic medium and the primordial matter power spectrum. In this article, we present and implement a simple method to measure the 3-dimensional power spectrum, P 3D , of the Lyα forest at wavenumbers k corresponding to small, ~ Mpc scales. In order to estimate P 3D from sparsely and unevenly distributed data samples, we rely on averaging 1-dimensional Fourier Transforms, as previously carried out to estimate the 1-dimensional power spectrum of the Lyα forest, P 1D . Further, this methodology exhibits a very low computational cost. We confirm the validity of this approach through its application to Nyx cosmological hydrodynamical simulations. Subsequently, we apply our method to the eBOSS DR16 Lyα forest sample, providing as a proof of principle, a first P 3D measurement averaged over two redshift bins z = 2.2 and z = 2.4. This work highlights the potential for forthcoming P 3D measurements, from upcoming large spectroscopic surveys, to untangle degeneracies in the cosmological interpretation of P 1D .

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Thermomagnetic properties of Bi 2 Te 3 single crystal in the temperature range from 55 K to 380 K

Magnetothermoelectric transport provides an understanding of coupled electron-hole-phonon current in topological materials and has applications in energy conversion and cooling. In this paper, we study the Nernst coefficient, the magneto-Seebeck coefficient, and the magnetoresistance of single-crystalline Bi 2 Te 3 under external magnetic field in the range of –3T to 3T and in the temperature range of 55 K to 380 K. Moreau's relation is employed to justify both the overall trend of the Nernst coefficient and the temperature at which this coefficient changes sign. We observe a nonlinear relationship between the Nernst coefficient and the applied magnetic field in the temperature range of 55 K to 255 K. An increase in both the Nernst coefficient and the magneto-Seebeck coefficient is observed as the temperature is reduced which can be attributed to the increased mobility of the carriers at lower temperatures. First-principles density functional theory calculations were carried out to physically model the experimental data including electronic and transport properties. Simulation findings agreed with the experiments and provide a theoretical insight to justify the measurements.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Modeling and Simulation of Electrostatics of Ge$_{\text{1-x}}$Sn$_{\text{x}}$ Layers Grown on Ge Substrates

This work introduces a comprehensive simulation tool that provides a robust 1D Schrödinger – Poisson solver for modeling the electrostatics of heterostructures with an arbitrary number of layers, and non-uniform doping profiles along with the treatment of partial ionization of dopants at low temperatures. The effective masses are derived from the first-principles calculations. The solver is used to characterize three Ge 1-x Sn x /Ge heterostructures with non-uniform doping profiles and determine the subband structure at various temperatures. Here, the simulation results of the sheet carrier densities show excellent agreement with the experimentally extracted data, thus demonstrating the capabilities of the solver.

42 ENGINEERING↗
DOE OSTI2022

Fully Kinetic Shearing-box Simulations of Magnetorotational Turbulence in 2D and 3D. I. Pair Plasmas

The magnetorotational instability (MRI) is a fundamental mechanism determining the macroscopic dynamics of astrophysical accretion disks. In collisionless accretion flows around supermassive black holes, MRI-driven plasma turbulence cascading to microscopic (i.e., kinetic) scales can result in enhanced angular-momentum transport and redistribution, nonthermal particle acceleration, and a two-temperature state where electrons and ions are heated unequally. However, this microscopic physics cannot be captured with standard magnetohydrodynamic (MHD) approaches typically employed to study the MRI. In this work, we explore the nonlinear development of MRI turbulence in a pair plasma, employing fully kinetic particle-in-cell (PIC) simulations in two and three dimensions. First, we thoroughly study the axisymmetric MRI with 2D simulations, explaining how and why the 2D geometry produces results that differ substantially from 3D MHD expectations. We then perform the largest (to date) 3D simulations, for which we employ a novel shearing-box approach, demonstrating that 3D PIC models can reproduce the mesoscale (i.e., MHD) MRI dynamics in sufficiently large runs. With our fully kinetic simulations, we are able to describe the nonthermal particle acceleration and angular-momentum transport driven by the collisionless MRI. Since these microscopic processes ultimately lead to the emission of potentially measurable radiation in accreting plasmas, our work is of prime importance to understand current and future observations from first principles, beyond the limitations imposed by fluid (MHD) models. While in this first study we focus on pair plasmas for simplicity, our results represent an essential step toward designing more realistic electron–ion simulations, on which we will focus in future work.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Autonomous Microgrid Restoration Using Grid-Forming Inverters and Smart Circuit Breakers

The proliferation of distributed inverter-based resources (IBRs) raises the questions if these IBRs can be used to blackstart microgrids and distribution feeders after major outages. In this paper, we propose and evaluate an autonomous microgrid restoration concept using grid-forming (GFM) IBRs and smart circuit creakers (SCBs). The concept is first explored in simulation platform and then a hardware testbed containing actual GFM inverters is developed to demonstrate these functionalities. A combination of dispatchable virtual oscillator control (dVOC) and droop-based control schemes have been designed for GFM inverter controls in the software simulations and hardware testbed. Subsequently, operation of SCBs using two distinct principles have been demonstrated that can restore or connect portions of the network.

black start↗
DOE OSTI2020

Shear Softening in a Metallic Glass: First-Principles Local-Stress Analysis

Metallic glasses deform elastically under stress. However, the atomic-level origin of elastic properties of metallic glasses remain unclear. In this Letter using ab initio molecular dynamics simulations of the Cu 5 0Zr 50 metallic glass under shear strain, we show that the heterogeneous stress relaxation results in the increased charge transfer from Zr to Cu atoms, enhancing the softening of the shear modulus. Changes in compositional short-range order and atomic position shifts due to the nonaffine deformation are discussed. It is shown that the Zr subsystem exhibits a stiff behavior, whereas the displacements of Cu atoms from their initial positions, induced by the strain, provide the stress drop and softening.

36 MATERIALS SCIENCE↗
DOE OSTI2020

A Flexible and Secure Evaluation Platform for Overvoltage Protection in Power Electronics Systems

This paper summarizes the fault mechanisms specifically for the overvoltage protection schemes in power electronic converters and the resulted power electronics enabled systems. A flexible and secure evaluation platform for overvoltage protection, has also been presented based on this conceptual mechanism. The proposed platform aims at testing the protection performance of distinctive overvoltage protection schemes. It would be applied to verify the effectiveness, identify potential strengths and weaknesses, and obtain the most critical protection response time to provide sufficient engineering insights for the designer. The operation principle and detailed design guideline have been presented first. The comprehensive component stress has been analyzed. An intensive set of LTSpice simulations have been conducted to show the validity of the proposed assessment approach.

Zheng, Sheng↗
DOE OSTI2022

Parallel Simulation of Beam Dynamics in Particle Accelerators [Slides]

LW3D was used to explore the expected CSR from an undulator-based diagnostic using a first-principles approach. It was also used for comparison with an analytic Green-function approach. LW3D was extended to model transport systems involving dipoles and quadrupoles. It was used to explore the radiation produced in a 4-dipole bunch compressor. Radiation ahead of the electron bunch at the beamline exit was shown to be associated with wavefronts produced in upstream dipoles overtaking the bunch.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Code interoperability extends the scope of quantum simulations

The functionality of many materials is critically dependent on the integration of dissimilar components and on the interfaces that arise between them. The description of such heterogeneous components requires the development and deployment of first principles methods, coupled to appropriate dynamical descriptions of matter and advanced sampling techniques, in order to capture all the relevant length and time scales of importance to the materials’ performance. It is thus essential to build simple, streamlined computational schemes for the prediction and design of multiple properties of broad classes of materials, by developing interoperable codes which can be efficiently coupled to each other to perform complex tasks. We discuss the use of interoperable codes to simulate the structural and spectroscopic characterization of materials, including chemical reactions for catalysis, the description of defects for quantum information science, and heat and charge transport.

36 MATERIALS SCIENCE↗