Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Coulomb”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,189 records · Page 66
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

An optimally tuned range-separated hybrid starting point for ab initio GW plus Bethe–Salpeter equation calculations of molecules

The ab initio GW plus Bethe–Salpeter equation (GW-BSE, where G is the one particle Green's function and W is the screened Coulomb interaction) approach has emerged as a leading method for predicting excitations in both solids and molecules with a predictive power contingent upon several factors. Among these factors are the (1) generalized Kohn–Sham eigensystem used to construct the GW self-energy and to solve the BSE and (2) the efficacy and suitability of the Tamm–Dancoff approximation. Here, we present a detailed benchmark study of low-lying singlet excitations from a generalized Kohn–Sham (gKS) starting point based on an optimally tuned range-separated hybrid (OTRSH) functional. We show that the use of this gKS starting point with one-shot G0W0 and G0W0-BSE leads to the lowest mean absolute errors (MAEs) and mean signed errors (MSEs), with respect to high-accuracy reference values, demonstrated in the literature thus far for the ionization potentials of the GW100 benchmark set and for low-lying neutral excitations of Thiel’s set molecules in the gas phase, without the need for self-consistency. The MSEs and MAEs of one-shot G0W0-BSE@OTRSH excitation energies are comparable to or lower than those obtained with other functional starting points after self-consistency. Additionally, we compare these results with linear-response time-dependent density functional theory (TDDFT) calculations and find GW-BSE to be superior to TDDFT when calculations are based on the same exchange-correlation functional. This work demonstrates tuned range-separated hybrids used in combination with GW and GW-BSE can greatly suppress starting point dependence for molecules, leading to accuracy similar to that for higher-order wavefunction-based theories for molecules without the need for costlier iterations to self-consistency.

McKeon, Caroline A. (ORCID:0000000217373503)↗
DOE OSTI2022

Perspective on solid-phase epitaxy as a method for searching novel topological phases in pyrochlore iridate thin films

Pyrochlore iridates exhibit various novel topological phenomena due to their topology and electron correlation. Notably, pyrochlore iridate was the first proposed material system shown to host the time-reversal symmetry broken Weyl semimetal phase. After this profound theoretical proposal, extensive experimental attempts have been made to synthesize high-quality samples in both bulk and film forms. In particular, adjusting thin film geometry represents a practical way to tune the U/ t ( U: Coulomb interaction and t: for hopping parameter) and identify the various topological phases in the regime U ~ λ ( λ: spin–orbit interaction). However, the instability of pyrochlore iridates at high temperature and low oxygen pressure has long been a barrier to growing pyrochlore iridate thin films by conventional film growth methods. To overcome this, pyrochlore iridate films have often been grown by the solid-phase epitaxy method, which uses a metastable amorphous layer grown at low temperature. During a high temperature post-annealing process, the layer in contact with a single crystalline substrate will crystallize epitaxially in the solid state by rearranging atoms at the interface. Here, we present a perspective on the solid-phase epitaxy as a method to synthesize epitaxial pyrochlore iridate thin films and a way to search for novel correlated phenomena.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Quasilinear theory of general electromagnetic fluctuations including discrete particle effects for magnetized plasmas: General analysis

The general quasilinear Fokker–Planck kinetic equation for the gyrophase-averaged plasma particle distribution functions in magnetized plasmas is derived, making no restrictions on the energy of the particles and on the frequency of the electromagnetic fluctuations and avoiding the often made Coulomb approximation of the electromagnetic interactions. Here, the inclusion of discrete particle effects breaks the dichotomy of nonlinear kinetic plasma theory divided into the test particle and the test fluctuation approximation because it provides expression of both the non-collective and collective electromagnetic fluctuation spectra in terms of the plasma particle distribution functions. Within the validity of the quasilinear approach, the resulting full quasilinear transport equation can be regarded as a determining nonlinear equation for the time evolution of the plasma particle distribution functions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Global gyrokinetic study of shaping effects on electromagnetic modes at NSTX aspect ratio with ad hoc parallel magnetic perturbation effects

Plasma shaping may have a stronger effect on global turbulence in tight-aspect-ratio tokamaks than in conventional-aspect-ratio tokamaks due to the higher toroidicity and more acute poloidal asymmetry in the magnetic field. In addition, previous local gyrokinetic studies have shown that it is necessary to include parallel magnetic field perturbations in order to accurately compute growth rates of electromagnetic modes in tight-aspect-ratio tokamaks. In this work, the effects of elongation and triangularity on global, ion-scale, linear electromagnetic modes are studied at National Spherical Torus Experiment (NSTX) aspect ratio and high plasma β using the global gyrokinetic particle-in-cell code XGC. The effects of compressional magnetic perturbations are approximated via a well-known modification to the particle drifts that was developed for flux-tube simulations [Joiner et al., Phys. Plasmas 17, 072104 (2010)], without proof of its validity in a global simulation, with the gyrokinetic codes GENE and GEM being used for local verification and global cross-verification. Magnetic equilibria are re-constructed for each distinct plasma profile that is used. Coulomb collision effects are not considered. Within the limitations imposed by the present study, it is found that linear growth rates of electromagnetic modes (collisionless microtearing modes and kinetic ballooning modes) are significantly reduced in a high-elongation and high-triangularity NSTX-like geometry compared to a circular NSTX-like geometry. For example, growth rates of kinetic ballooning modes at high- β are reduced to the level of that of collisionless trapped electron modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Self-consistent implementation of locally scaled self-interaction-correction method

Recently proposed local self-interaction correction (LSIC) method is a one-electron self-interaction-correction (SIC) method that uses an iso-orbital indicator to apply the SIC at each point in space by scaling the exchange–correlation and Coulomb energy densities. The LSIC method is exact for the one-electron densities, also recovers the uniform electron gas limit of the uncorrected density functional approximation, and reduces to the well-known Perdew–Zunger SIC (PZSIC) method as a special case. This article presents the self-consistent implementation of the LSIC method using the ratio of Weizsäcker and Kohn–Sham kinetic energy densities as an iso-orbital indicator. The atomic forces as well as the forces on the Fermi-Löwdin orbitals are also implemented for the LSIC energy functional. Results show that LSIC with the simplest local spin density functional predicts atomization energies of the AE6 dataset better than some of the most widely used generalized-gradient-approximation (GGA) functional and barrier heights of the BH6 database better than some of the most widely used hybrid functionals. The LSIC method [a mean absolute error (MAE) of 0.008 Å] predicts bond lengths of a small set of molecules better than the PZSIC-LSDA (MAE 0.042 Å) and LSDA (0.011 Å). Furthermore, this work shows that accurate results can be obtained from the simplest density functional by removing the self-interaction-errors using an appropriately designed SIC method.

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

Resolving the mystery of electron perpendicular temperature spike in the plasma sheath

A large family of plasmas has collisional mean-free-path much longer than the non-neutral sheath width, which scales with the plasma Debye length. The plasmas, particularly the electrons, assume strong temperature anisotropy in the sheath. The temperature in the sheath flow direction (Te∥) is lower and drops toward the wall as a result of the decompressional cooling by the accelerating sheath flow. The electron temperature in the transverse direction of the flow field (Te⊥) not only is higher but also spikes up in the sheath. This abnormal behavior of Te⊥ spike is found to be the result of a negative gradient of the parallel heat flux of transverse degrees of freedom (qes) in the sheath. The non-zero heat flux qes is induced by pitch-angle scattering of electrons via either their interaction with self-excited electromagnetic waves in a nearly collisionless plasma or Coulomb collision in a collisional plasma, or both in the intermediate regime of plasma collisionality.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Atomistic insight into the effects of electrostatic fields on hydrocarbon reaction kinetics

Reactive Molecular Dynamics (MD) and Density Functional Theory (DFT) computations are performed to provide insight into the effects of external electrostatic fields on hydrocarbon reaction kinetics. By comparing the results from MD and DFT, the suitability of the MD method in modeling electrodynamics is first assessed. Results show that the electric field-induced polarization predicted by the MD charge equilibration method is in good agreement with various DFT charge partitioning schemes. Then, the effects of oriented external electric fields on the transition pathways of non-redox reactions are investigated. Results on the minimum energy path suggest that electric fields can cause catalysis or inhibition of oxidation reactions, whereas pyrolysis reactions are not affected due to the weaker electronegativity of the hydrogen and carbon atoms. MD simulations of isolated reactions show that the reaction kinetics is also affected by applied external Lorentz forces and interatomic Coulomb forces since they can increase or decrease the energy of collision depending on the molecular conformation. In addition, electric fields can affect the kinetics of polar species and force them to align in the direction of field lines. These effects are attributed to energy transfer via intermolecular collisions and stabilization under the external Lorentz force. The kinetics of apolar species is not significantly affected mainly due to the weak induced dipole moment even under strong electric fields. The dynamics and reaction rates of species are studied by means of large-scale combustion simulations of n-dodecane and oxygen mixtures. Results show that under strong electric fields, the fuel, oxidizer, and most product molecules experience translational and rotational acceleration mainly due to close charge transfer along with a reduction in their vibrational energy due to stabilization. This study will serve as a basis to improve the current methods used in MD and to develop novel methodologies for the modeling of macroscale reacting flows under external electrostatic fields.

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

Distributed memory, GPU accelerated Fock construction for hybrid, Gaussian basis density functional theory

With the growing reliance of modern supercomputers on accelerator-based architecture such a graphics processing units (GPUs), the development and optimization of electronic structure methods to exploit these massively parallel resources has become a recent priority. While significant strides have been made in the development GPU accelerated, distributed memory algorithms for many modern electronic structure methods, the primary focus of GPU development for Gaussian basis atomic orbital methods has been for shared memory systems with only a handful of examples pursing massive parallelism. Here in this work, we present a set of distributed memory algorithms for the evaluation of the Coulomb and exact exchange matrices for hybrid Kohn–Sham DFT with Gaussian basis sets via direct density-fitted (DF-J-Engine) and seminumerical (sn-K) methods, respectively. The absolute performance and strong scalability of the developed methods are demonstrated on systems ranging from a few hundred to over one thousand atoms using up to 128 NVIDIA A100 GPUs on the Perlmutter supercomputer.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Faraday effect in collisional magnetized plasmas

Faraday rotation is a valuable diagnostic tool for investigating laboratory and astrophysical plasmas, but it has almost exclusively been treated in the collisionless limit (despite the non-negligible role of Coulomb collisions in many laboratory plasmas). Here, we show that Faraday rotation can occur in collisional plasmas and that the usual effect is curiously altered by collisions. Namely, an initially linearly polarized light wave—propagating parallel to a uniform magnetic field—becomes elliptically polarized in the collisional plasma, and collisional absorption of the wave also occurs. Moreover, the ellipticity and rotation angle are quantifiably sensitive to the fidelity of the collisional transport coefficients. As we will demonstrate with particle-in-cell EPOCH simulations, these effects offer a stringent diagnostic tool for benchmarking multi-scale plasma simulation codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Impact of solvation on the GW quasiparticle spectra of molecules

First-principles calculations for electrochemistry require accurate treatment of both electronic structure and solvation. The perturbative GW approximation starting from density functional theory (DFT) calculations accurately models materials systems with varying dimensionality. Continuum solvation models enable efficient treatment of solvation effects in DFT calculations, but their applications with beyond-DFT electronic structure methods such as GW have been limited. Here, we introduce the frequency-dependent liquid polarizability from a nonlocal continuum solvation model in the screened Coulomb interaction of full-frequency GW calculations with a solvated DFT starting point. In this study, we show that the liquid screening contributions substantially reduce the HOMO–LUMO gap of molecules by 3–5 eV, while solvent effects on the DFT starting point negligibly impact the GW gap. The resulting framework facilitates the simultaneous electronic and solvation accuracy needed for first-principles electrochemistry.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Atomistic models of In and Ga diffusion in Cu(In,Ga)Se 2

To improve the performance of Cu(In,Ga)Se 2 thin-film photovoltaic devices, a robust understanding of the dominant diffusion pathways of the alloy species In and Ga is needed. Here, the most probable defect complexes and mechanisms for In and Ga diffusion are identified with the aid of density functional theory. The binding energies and migration barriers for these complexes are calculated in bulk CuInSe 2 and CuGaSe 2 . Analytic models and kinetic lattice Monte Carlo simulations are employed to predict the diffusivity of In and Ga under variations in composition and temperature. Here, we find that a model based on coulombic interactions between group III antisites and vacancies on the Cu-sublattice produces results that match well with experiment.

14 SOLAR ENERGY↗
DOE OSTI2023

Influence of trap-assisted and intrinsic Auger–Meitner recombination on efficiency droop in green InGaN/GaN LEDs

Here, we study the impact of deep-level defects on trap-assisted Auger–Meitner recombination in c-plane InGaN/GaN LEDs using a small-signal electroluminescence (SSEL) method and deep-level optical spectroscopy (DLOS). Carrier dynamics information, including carrier lifetime, recombination rate, and carrier density, is obtained from SSEL, while DLOS is used to obtain the deep-level defect density. Through fitting the nonradiative recombination rates of wafers with different deep-level defect densities, we obtain the Shockley–Read–Hall (SRH) and trap-assisted Auger–Meitner recombination (TAAR) coefficients. We show that defect-related nonradiative recombination, including both SRH and TAAR, accounts for a relatively small fraction of the total nonradiative recombination, which is dominated by intrinsic Auger–Meitner recombination. The interplay between carrier localization and Coulomb enhancement has a different impact on radiative and intrinsic Auger–Meitner recombination. Evidence is presented that the imbalance between the change of radiative and intrinsic Auger–Meitner recombination is the primary cause of the efficiency droop at high carrier densities in the samples studied.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Development of spinning-disk solid sample delivery system for high-repetition rate x-ray free electron laser experiments

X-ray free-electron lasers (XFELs) deliver intense x-ray pulses that destroy the sample in a single shot by a Coulomb explosion. Experiments using XFEL pulse trains or the new generation of high-repetition rate XFELs require rapid sample replacement beyond those provided by the systems now used at low repletion-rate XFELs. We describe the development and characterization of a system based on a spinning disk to continuously deliver a solid sample into an XFEL interaction point at very high speeds. We tested our system at the Linac Coherent Light Source and European XFEL hard x-ray nano-focus instruments, employing it to deliver a 25 μm copper foil sample, which can be used as a gain medium for stimulated x-ray emission for the proposed x-ray laser oscillator.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

The difference between molecules and materials: Reassessing the role of exact conditions in density functional theory

Exact conditions have long been used to guide the construction of density functional approximations. However, hundreds of empirical-based approximations tailored for chemistry are in use, of which many neglect these conditions in their design. We analyze well-known conditions and revive several obscure ones. Two crucial distinctions are drawn: that between necessary and sufficient conditions and that between all electronic densities and the subset of realistic Coulombic ground states. Simple search algorithms find that many empirical approximations satisfy many exact conditions for realistic densities and non-empirical approximations satisfy even more conditions than those enforced in their construction. The role of exact conditions in developing approximations is revisited.

Chemistry↗
DOE OSTI2023

Use of FLOSIC for understanding anion-solvent interactions

An Achille’s heel of lower-rung density-functional approximations is that the highest-occupied-molecular-orbital energy levels of anions, known to be stable or metastable in nature, are often found to be positive in the worst case or above the lowest-unoccupied-molecular-orbital levels on neighboring complexes that are not expected to accept charge. A trianionic example, [Cr(C2O4)3]3−, is of interest for constraining models linking Cr isotope ratios in rock samples to oxygen levels in Earth’s atmosphere over geological timescales. Here we describe how crowd sourcing can be used to carry out self-consistent Fermi–Löwdin–Orbital-Self-Interaction corrected calculations (FLOSIC) on this trianion in solution. The calculations give a physically correct description of the electronic structure of the trianion and water. In contrast, uncorrected local density approximation (LDA) calculations result in approximately half of the anion charge being transferred to the water bath due to the effects of self-interaction error. Use of group-theory and the intrinsic sparsity of the theory enables calculations roughly 125 times faster than our initial implementation in the large N limit reached here. By integrating charge density densities and Coulomb potentials over regions of space and analyzing core-level shifts of the Cr and O atoms as a function of position and functional, we unambiguously show that FLOSIC, relative to LDA, reverses incorrect solute-solvent charge transfer in the trianion-water complex. In comparison to other functionals investigated herein, including Hartree–Fock and the local density approximation, the FLOSIC Cr 1s eigenvalues provide the best agreement with experimental core ionization energies.

Chemistry↗
DOE OSTI2023

Shivalik Plasma Device-I, a glow discharge device to study the collective dynamics of dusty plasma

Using a uniquely configured glow discharge-based Shivalik Plasma Device-I, we demonstrate a variety of collective phenomena in dusty plasma away from the glow discharge region. The cylindrical glass device produces plasma using parallel disc-shaped electrodes with a smaller anode size than the cathode. The dust microparticles are initially sprinkled over the grounded cathode. These particles acquire a significant negative charge upon plasma formation, resulting in their levitation due to the balance between the Coulomb force and gravity. The new device supports the levitation of a big-sized (10 × 8 × 5 cm 3 ) three-dimensional dust cloud over the glass surface. It contrasts the dusty plasma formations in-between electrodes reported earlier. As the discharge voltage varies from high to low, the dust cloud travels from over the glass surface to between the electrodes. A complex interplay of dust void over the cathode, a sharp density gradient, and gravity lead to self-excitation of collective dust phenomena. It includes dust density waves (phase velocity, v ph ~4 cm/s), dust cloud oscillation (frequency, f = 5 Hz), sheared flow (flow velocity, v f ~ 1 cm/s), and multiple-sized dust vortices. These dust vortices provided an excellent platform for studying turbulent mixing phenomena. The power spectrum analysis agreed with two-dimensional Kolmogorov power-law scaling. This is an ideal dusty plasma apparatus where we can create or move the dust cloud to a location of choice from the glass surface to in-between the electrodes and excite one among many collective dust dynamics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Implementation of a drive cylinder for low collisional experiments on magnetic reconnection

For most laboratory plasma experiments, Coulomb collisions between the particle species are sufficiently frequent that the particle distribution functions are relaxed to a near-Maxwellian form. This hampers the applicability of such experiments to phenomena observed in tenuous and near-collisionless space plasma. The Terrestrial Reconnection EXperiment (TREX) at the Wisconsin Plasma Physics Laboratory aims to study collisionless reconnection for parameters relevant to the Earth’s magnetosphere. To reduce the role of collisional effects, a reconnection Drive Cylinder has been developed, which increases both the effective system size of the TREX configuration and the rate at which reconnection can be driven. These two effects now permit TREX to reach a kinetic reconnection regime where collisional effects are minimized. The Drive Cylinder is comprised of 12 single loop drive-coils connected in parallel to a 10 kV capacitor bank. Insulated sheets of aluminum are applied to smooth the magnetic fields and enhance the drive efficiency. Following is a description of the technical details and performance of the Drive Cylinder.

Instruments & Instrumentation↗
DOE OSTI2024

Numerical thermalization in 2D PIC simulations: Practical estimates for low-temperature plasma simulations

The process of numerical thermalization in particle-in-cell (PIC) simulations has been studied extensively. It is analogous to Coulomb collisions in real plasmas, causing particle velocity distributions (VDFs) to evolve toward a Maxwellian as macroparticles experience polarization drag and resonantly interact with the fluctuation spectrum. This paper presents a practical tutorial on the effects of numerical thermalization in 2D PIC applications. Scenarios of interest include simulations, which must be run for many thousands of plasma periods and contain a population of cold electrons that leave the simulation space very slowly. This is particularly relevant to many low-temperature plasma discharges and materials processing applications. We present numerical drag and diffusion coefficients and their associated timescales for a variety of grid resolutions, discussing the circumstances under which the electron VDF is modified by numerical thermalization. Though the effects described here have been known for many decades, direct comparison of analytically derived, velocity-dependent numerical relaxation timescales to those of other relevant processes has not often been applied in practice due to complications that arise in calculating thermalization rates in 1D simulations. Using these comparisons, we estimate the impact of numerical thermalization in several examples of low-temperature plasma applications including capacitively coupled plasma discharges, inductively coupled plasma discharges, beam plasmas, and hollow cathode discharges. Finally, we discuss possible strategies for mitigating numerical relaxation effects in 2D PIC simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗