SEARCH · Engineering Papers

Results for “Electrostatic Potential”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 145 records · Page 8

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Quantum physical reality of polar-nonpolar oxide heterostructures

Conducting interfaces between polar and nonpolar insulating oxides, e.g., LaAlO 3 /SrTiO 3 , have generated interest for both fundamental physics and oxide-electronics applications. Current understanding is based on an amalgamation of a classical electrostatic model (polar catastrophe model) that was originally derived for semi-infinite solids and quantum density-functional-theory (DFT) results on ultrathin films. Here we report comprehensive DFT calculations that unveil a very different purely quantum physical reality. We show that, for ultrathin polar films, the interfacial dipole does not control the electrostatic potential in the polar film—the surface and interface play equal roles, and the absence or presence of centrosymmetry in the physical LAO film results in different, purely quantum mechanisms for the generation of a conducting interface, neither involving physical-charged transfer. Predictions are made that can be tested and can guide technology development.

36 MATERIALS SCIENCE↗

DOE OSTI2021

First-principles modeling of conductivity at the (001), (110), and (111) SrTiO 3 /LaAlO 3 heterointerfaces

Here, the complex polar oxide heterojunction of SrTiO 3 /LaAlO 3 (STO/LAO) is of great interest due to the emergent physical phenomena observed at the interface. STO and LAO separately are wide band-gap insulators. However, upon joining them at the 001, 110, and 111 crystallographic planes, the interface undergoes a transition to a conductive state. Although first-principles modeling of the 001 plane interface has been widely studied, there is a lack of reports regarding the 110 and 111. This paper expands the theoretical model of the STO/LAO heterointerface to the three crystallographic planes (001, 110, and 111) where the conductivity has been experimentally reported. The calculations showed that whereas at the 001 interface the conductivity appears at a critical thickness of 4 monolayers of LAO, the 110 and 111 planes have no clear critical thickness; these two interfaces were always conductive. Nevertheless, the number of conductive electrons per unit cell increases with the thickness of the LAO layer in the 110 and 111 interfaces. This is related to the energy levels downshifting due to the electrostatic potential buildup (which was in the opposite direction respect to the 001 interface), increasing the number of conductive sates below the Fermi level. Given the absence of a critical thickness and the fact that chemical intermixing and oxygen vacancies at the interface were not considered, the main mechanism responsible for the conductivity in the 110 and 111 planes was attributed to the large structure reconstruction that locally changes the energy levels at the interface causing charge transfer and accumulation at the layers close to the interface.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2022

Unified ab initio description of Fröhlich electron-phonon interactions in two-dimensional and three-dimensional materials

Ab initio calculations of electron-phonon interactions including the polar Fröhlich coupling have advanced considerably in recent years. The Fröhlich electron-phonon matrix element is by now well understood in the case of bulk three-dimensional (3D) materials. In the case of two-dimensional (2D) materials, the standard procedure to include Fröhlich coupling is to employ Coulomb truncation, so as to eliminate artificial interactions between periodic images of the 2D layer. While these techniques are well established, the transition of the Fröhlich coupling from three to two dimensions has not been investigated. Furthermore, it remains unclear what error one makes when describing 2D systems using the standard bulk formalism in a periodic supercell geometry. In this work, we generalize previous work on the ab initio Fröhlich electron-phonon matrix element in bulk materials by investigating the electrostatic potential of atomic dipoles in a periodic supercell consisting of a 2D material and a continuum dielectric slab. We obtain a unified expression for the matrix element, which reduces to the existing formulas for three-dimensional and 2D systems when the interlayer separation tends to zero or infinity, respectively. This expression enables an accurate description of the Fröhlich matrix element in 2D systems without resorting to Coulomb truncation. We validate our approach by direct ab initio density-functional perturbation theory calculations for monolayer BN and MoS 2 , and we provide a simple expression for the 2D Fröhlich matrix element that can be used in model Hamiltonian approaches. The formalism outlined in this work may find applications in calculations of polarons, quasiparticle renormalization, transport coefficients, and superconductivity, in 2D and quasi-2D materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2020

Lower-Hybrid Drift Waves Driving Electron Nongyrotropic Heating and Vortical Flows in a Magnetic Reconnection Layer

Here, we report measurements of lower-hybrid drift waves driving electron heating and vortical flows in an electron-scale reconnection layer under a guide field. Electrons accelerated by the electrostatic potential of the waves exhibit perpendicular and nongyrotropic heating. The vortical flows generate magnetic field perturbations comparable to the guide field magnitude. The measurements reveal a new regime of electron-wave interaction and how this interaction modifies the electron dynamics in the reconnection layer.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Ablating Ion Velocity Distributions in Short-Pulse-Heated Solids via X-Ray Doppler Shifts

Solids ablate under laser irradiation, but experiments have not previously characterized the initiation of this process at ultrarelativistic laser intensities. Here, we present first measurements of bulk ion velocity distributions as ablation begins, captured as a function of depth via Doppler-shifted x-ray line emission from two viewing angles. Bayesian analysis indicates that bulk ions are either nearly stationary or flowing outward at the plasma sound speed. The measurements quantitatively constrain the laser-plasma ablation mechanism, suggesting that a steplike electrostatic potential structure drives solid disassembly.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Interplay between charge ordering and geometric ferroelectricity in LuFe 2 O 4 /LuFeO 3 superlattices

Oxide superlattices have drawn great attentions owing to the intriguing coupling among elastic, electrical, and magnetic orderings at the interfaces and the emerging of improper ferroelectricity. Here, superlattices composed of hexagonal LuFeO 3 (h-LuFeO 3 ) and LuFe 2 O 4 are investigated via density functional theory calculations. h-LuFeO 3 is a well-known multiferroic material that is stable only in thin-film or doped bulk state, while LuFe 2 O 4 is a charge ordered material where the existence of ferroelectricity is still in controversy. We have found that the charge ordering (CO) induced polarizations in LuFe 2 O 4 layers coexist with the geometric polarizations in h-LuFeO 3 layers in the (LuFe 2 O 4 ) m /(LuFeO 3 ) n superlattices with different periodicities, and the ferroelectric states are generally preferred over the antiferroelectric states for LuFe 2 O 4 in superlattices. The out-of-plane polarizations in h-LuFeO 3 and LuFe 2 O 4 layers tend to be aligned in parallel, and the overall polarization increases with the ratio of h-LuFeO 3 . The influence of layered polarizations on the local electrostatic potential is not significant except the detected small trend caused by the CO-induced polarization within a FeO bilayer. Besides, the local electronic structures show that the Fermi level position in a certain layer can be tuned by the valences of Fe in this layer and the polarization distributions in neighboring layers. LuFe 2 O 4 layers sandwiched between thick h-LuFeO 3 layers are more susceptible. The calculated configurations of the superlattices are supported by atomic-resolution transmission electron microscopy experiments. Our results pave a way for tunable ferroelectricity in superlattice systems and create a new playground for manipulating the coupling between various degrees of freedom.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2026

pyDiSCaMB : enabling the use of multipolar scattering factors in Phenix

Multipolar scattering models, such as the transferable aspherical atom model, account for atomic chemical interactions and provide a more accurate representation of experimental data. However, the simpler independent atom model (IAM), which assumes non-interacting atoms, is the only model available in the most widely used macromolecular refinement programs. This is primarily because IAM offers a hard-to-beat combination of computational efficiency and modelling power at typical macromolecular resolutions. By contrast, more accurate multipolar modelling has historically been limited due to its computational cost and the absence of an interface between software capable of calculating structure factors and gradients based on multipolar models and software designed for macromolecular refinement. This work introduces pyDiSCaMB , a Python software package designed to integrate between the computational crystallography toolbox ( cctbx ) and the quantum crystallography library DiSCaMB ( Densities in Structural Chemistry and Molecular Biology ), thus enabling multipolar scattering models in Phenix 's toolkit. The implementation, features and capabilities of pyDiSCaMB are presented, the runtimes for the calculation of structure factor and target gradients with respect to atomic parameters are explored, and Fourier images of electrostatic potential, electron density and deformation maps are computed as illustrative examples. The pyDiSCaMB library will make multipolar modelling widely available to the structural biology community, potentially transforming refinement and model-building for both crystallography and cryogenic electron microscopy (cryoEM).

MATTS data bank↗

DOE OSTI2020

Operando Microscopy Characterization of Perovskite Solar Cells

In this work, we developed operando Kelvin probe force microscopy (KPFM) to study the electrostatic potential distribution across perovskite cells under light and forward bias to gain a deeper understanding of device operation physics. As a case study, we selected perovskite cells with a SnO 2 -based electron-selective layer (ESL), which showed great potential for fabricating high-efficiency, hysteresis-free devices due to the deeper conduction band and higher electron mobility of SnO 2 . The as-made device showed a main junction at the perovskite/spiro interface. After light soaking and applying forward bias, the junction quality improved, possibly explained by filling trap states at the interfaces and by the perovskite absorber perhaps having a self-poling effect; the main junction is observed at the ESL/perovskite interface. The results are consistent with current-voltage measurements, device performance improves mainly with fill factor enhancement. The operando KPFM results should more closely reflect the real case during current density-voltage measurements or solar cell operation. The operando KPFM technique that we have developed can be a powerful tool to provide a deeper understanding of the device operation mechanism and to further optimize the device.

14 SOLAR ENERGY↗

DOE OSTI2020

Advanced modeling of field enhanced thermionic emission

Shaped emitters are of interest to a broad range of applications in vacuum electronic devices. In particular, thermionic energy converters (TECs) take advantage of shaped emitters to increase the local surface field, thereby extracting more current for a given cathode temperature and applied voltage. However, modeling these devices is challenging; Warp [J.-L. Vay, D. P. Grote, R. H. Cohen, and A. Friedman, Comput. Sci. Discov. 5, 014019 (2012)] is a fully 3D particle-in-cell code capable of handling a wide range of physics problems and is well suited to modeling TECs. Additionally, recent improvements to Warp have enabled the accurate modeling of emitters with arbitrary curved surfaces. Specifically, the inclusion of subgrid resolution for computing the electrostatic potential and the ability to apply mesh refinement for specific areas of interest allow for a more accurate solution to the fields on these surfaces. These improvements coupled with Warp’s ability to handle variable particle weights make it an ideal candidate for simulating these complex devices. In this paper, the authors study the applicability of different subgrid configurations for simulating shaped emission surfaces and field convergence for different mesh-refinement techniques. They then implement a custom weighting algorithm that allows for uniform sampling of emission surfaces with a large variation in the surface electric field. They then use this algorithm to study emission for curved emitters in both the field-enhancement regime and the space-charge regime.

Edelen, Jonathan P. (ORCID:0000000215180652)↗

DOE OSTI2020

Strain-dependent luminescence and piezoelectricity in monolayer transition metal dichalcogenides

The modification of optical and electronic properties of transition metal dichalcogenides via mechanical deformation has been widely studied. Their ability to withstand large deformations before rupture has enabled large tunability of the bandgap, and further, the spatially varying strain has been shown to control the spatial distribution of the bandgap and lead to effects such as carrier funneling. Monolayer transition metal dichalcogenides exhibit a significant piezoelectric effect that could couple to a spatially inhomogeneous strain distribution to influence electronic and optical behavior. We investigate both experimentally and theoretically an example case of photoluminescence in structures with a strain distribution similar to that employed in single-photon emitters but generated here via nanoindentation. Using a mechanical model for strain induced by nanoindentation, we show that piezoelectricity can result in charge densities reaching 10 12 e/cm 2 and can generate electrostatic potential variations on the order of ±0.1 V across the suspended monolayer. We analyze the implications of these results for luminescence and exciton transport in monolayer transition metal dichalcogenides with spatially varying strain.

42 ENGINEERING↗

DOE OSTI2021

Determining valence band offsets in heterojunctions using a single core-level x-ray photoelectron spectrum

The behavior and functionality of semiconductor heterojunctions depend critically on the alignments of the valence and conduction bands at the various interfaces. Traditionally band alignment has been measured by x-ray photoelectron spectroscopy using pairs of distinct core levels, one from each side of the interface, to track the valence band maxima in the two materials. Here, we demonstrate that it is possible to determine band alignment using a single core-level photoelectron spectrum for an element that is common to both materials. The energy splitting between the photoemission features originating in the two materials is shown to be dominated by the difference in electrostatic potential across the interface, leading to a reliable determination of the band alignment.

Wang, Le↗

DOE OSTI2022

Gyrokinetic simulations of momentum flux parasitic to free-energy transfer

Ion Landau damping interacts with a portion of the E×B drift to cause a non-diffusive outward flux of co-current toroidal angular momentum. Quantitative evaluation of this momentum flux requires nonlinear simulations to determine fL, the fraction of fluctuation free energy that passes through ion Landau damping, in fully developed turbulence. Nonlinear gyrokinetic simulations with the GKW code confirm the presence of the systematic symmetry-breaking momentum flux. For simulations with adiabatic electrons, fL scales inversely with the ion temperature gradient, because only the ion curvature drift can transfer free energy to the electrostatic potential. Although kinetic electrons should in principle relax this restriction, the ion Landau damping measured in collisionless kinetic-electron simulations remained at low levels comparable with ion-curvature-drift transfer, except when magnetic shear was strong. A set of simulations scanning the electron pitch-angle scattering rate showed only a weak variation of fL with the electron collisionality. However, collisional-electron simulations with electron temperature greater than ion temperature unambiguously showed electron-curvature-drift transfer supporting ion Landau damping, leading to a corresponding enhancement of the symmetry-breaking momentum flux.

intrinsic rotation↗

DOE OSTI2024

Kinetic Plasma Simulation in the MOOSE Framework: Verification of Electrostatic Particle In Cell Capabilities

In magnetic confinement nuclear fusion reactors, the interaction between the plasma edge and plasma facing components is extremely important. At the plasma edge, a kinetic representation such as particle-in-cell (rather than a fluid representation) is required to accurately capture the plasma behavior. General purpose particle-in-cell plasma simulation capabilities have been developed in the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework. This new capability is a part of the development of a new MOOSE-based framework for modeling plasma facing components, the Fusion ENergy Integrated multiphys-X (FENIX) framework. In this work, the verification of foundational particle-in-cell capabilities in FENIX is presented. This new plasma simulation capability has three main components: moving particles in discrete steps on the finite element mesh, mapping charge density from the particle's location to the finite element mesh, and solving for the electrostatic potential based on the charge density mapped from particles to the mesh. In this paper, simple verification problems demonstrating each of these new capabilities are presented, and future work includes electromagnetic capabilities and Monte Carlo collisions with neutral gas particles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Close-Out Project Description for Koepke's Dept of Energy grant DE-SC0021405

Objectives: To establish, for lab & space conditions, EM-IEDDI’s (electromagnetic shear-driven instability's) dispersion relation, unstable range, instability threshold, and mode characteristics, we need LAPD’s Alfven-wave-favorable electromagnetic-style conditions, including higher "beta" (0.001 < beta ≤ 0.3) and low-collisionality. Also, we attempted to intentionally launch or spontaneously destabilize compressional and shear Alfven waves in the strong, localized, perpendicular-velocity-shear region at the interface between coaxial plasmas (one plasma cylinder inside an outer, otherwise hollow, tube, each having a different, controllable, value of plasma electrostatic potential, i.e., “space” potential (not to be confused with the temperature-dependent “floating” potential of an object immersed in the plasma). Nonlinear wave-wave interactions between same-family (EM-IEDD or Alfven) and cross-family (EM-IEDD-with-Alfven) fluctuations were targeted for documentation over a range of spectral overlap. Although laboratory experiments were conducted, the following theoretical work was left unfinished: Analytical non-modal prediction Computational non-modal prediction Check to see if Mikhailenko’s theory formulation leads to his published graphs

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Acid-Base Equilibrium and Dielectric Environment Regulate Charge in Supramolecular Nanofibers

Peptide amphiphiles are a class of molecules that can self-assemble into a variety of supramolecular structures, including high-aspect-ratio nanofibers. It is challenging to model and predict the charges in these supramolecular nanofibers because the ionization state of the peptides are not fixed but liable to change due to the acid-base equilibrium that is coupled to the structural organization of the peptide amphiphile molecules. Here, we have developed a theoretical model to describe and predict the amount of charge found on self-assembled peptide amphiphiles as a function of pH and ion concentration. In particular, we computed the amount of charge of peptide amphiphiles nanofibers with the sequence C 16 − V 2 A 2 E 2 . In our theoretical formulation, we consider charge regulation of the carboxylic acid groups, which involves the acid-base chemical equilibrium of the glutamic acid residues and the possibility of ion condensation. The charge regulation is coupled with the local dielectric environment by allowing for a varying dielectric constant that also includes a position-dependent electrostatic solvation energy for the charged species. We find that the charges on the glutamic acid residues of the peptide amphiphile nanofiber are much lower than the same functional group in aqueous solution. There is a strong coupling between the charging via the acid-base equilibrium and the local dielectric environment. Our model predicts a much lower degree of deprotonation for a position-dependent relative dielectric constant compared to a constant dielectric background. Furthermore, the shape and size of the electrostatic potential as well as the counterion distribution are quantitatively and qualitatively different. These results indicate that an accurate model of peptide amphiphile self-assembly must take into account charge regulation of acidic groups through acid–base equilibria and ion condensation, as well as coupling to the local dielectric environment.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Visualization techniques for the gyrokinetic tokamak simulation code

Gyrokinetic simulations of plasma microturbulence in tokamaks are challenging to visualize because the compute grid follows the magnetic field lines that spiral around the torus. We have overcome this challenge by developing three new approaches that improve visualization of gyrokinetics. Our techniques work directly with the topology of magnetic flux surfaces where the simulation stores variables in concentric rings on poloidal planes (vertical cross sections of the torus). Our visualization preview step triangulates each consecutive pair of rings to display the data on a poloidal plane. The second visualization technique follows spiral field lines around the torus and constructs polygons to visualize a flux surface. Third, the poloidal triangles are connected between planes to form prisms that compose a 3-D model of the entire torus. The visualization workflow produces detailed geometry that matches the high resolution, irregular compute grid for every time step. The surface and solid models are displayed in scientific visualization programs to effectively explore and communicate the results, including fluctuation of electron density, ion temperature, and electrostatic potential. Highly detailed renderings verify plasma behavior along magnetic field lines over time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Structural and Mechanistic Basis for the Inactivation of Human Ornithine Aminotransferase by (3 S ,4 S )-3-Amino-4-fluorocyclopentenecarboxylic Acid

Ornithine aminotransferase (OAT) is overexpressed in hepatocellular carcinoma (HCC), and we previously showed that inactivation of OAT inhibits the growth of HCC. Recently, we found that (3S,4S)-3-amino-4-fluorocyclopentenecarboxylic acid (5) was a potent inactivator of γ-aminobutyric acid aminotransferase (GABA-AT), proceeding by an enamine mechanism. Here we describe our investigations into the activity and mechanism of 5 as an inactivator of human OAT. We have found that 5 exhibits 10-fold less inactivation efficiency (k inact /K I ) against hOAT than GABA-AT. A comprehensive mechanistic study was carried out to understand its inactivation mechanism with hOAT. pK a and electrostatic potential calculations were performed to further support the notion that the α,β-unsaturated alkene of 5 is critical for enhancing acidity and nucleophilicity of the corresponding intermediates and ultimately responsible for the improved inactivation efficiency of 5 over the corresponding saturated analogue (4). Intact protein mass spectrometry and the crystal structure complex with hOAT provide evidence to conclude that 5 mainly inactivates hOAT through noncovalent interactions, and that, unlike with GABA-AT, covalent binding with hOAT is a minor component of the total inhibition which is unique relative to other monofluoro-substituted derivatives. Furthermore, based on the results of transient-state measurements and free energy calculations, it is suggested that the α,β-unsaturated carboxylate group of PLP-bound 5 may be directly involved in the inactivation cascade by forming an enolate intermediate. Overall, compound 5 exhibits unusual structural conversions which are catalyzed by specific residues within hOAT, ultimately leading to an enamine mechanism-based inactivation of hOAT through noncovalent interactions and covalent modification.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2024

Image-Based Digital Twin for Assessing the Coupled Electro-Chemo-Mechanical Behavior of Li-Ion Batteries

In energy storage materials, strong electrochemical-mechanical coupling and highly anisotropic material properties contribute to the formation and propagation of micro-cracking during charge/discharge cycling, resulting in reduced performance and service life. In this work, a digital twin is created to investigate the performance of a Li-ion battery cathode and simulate degradation accumulation. Pixel-based model construction is used to represent the complex material geometries from microstructural images supplied by the National Renewable Energy Laboratory (NREL). Because of the expected large deformation and crack opening, the reproducing kernel particle method (RKPM), a meshfree method with discretization at the image pixels, is used to approximate the field variables: electrostatic potential, concentration, and displacement. An interface modified reproducing kernel (IM-RK) is constructed by scaling a smooth kernel function with an interface-distance function to achieve strategic discontinuity types (i.e. weak discontinuities for strain discontinuities and strong discontinuities for cracks) and alleviate Gibbs oscillations near these transition zones. The discrete nature of the images' pixel points is employed throughout the model and approximation construction. The mechanical model is verified using an image-based microstructure under tensile loading. A transient electro-chemo-mechanical coupled simulation is performed to evaluate the potential micro-cracking induced degradation of the battery cathode material subjected to charge/discharge cycling.

battery degradation↗