Search NASA⌕ Search

SEARCH · Search NASA

Results for “electronic wave function”

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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Extreme Heating of Minor Ions in Imbalanced Solar-wind Turbulence

Minor ions in the solar corona are heated to extreme temperatures, far in excess of those of the electrons and protons that comprise the bulk of the plasma. These highly nonthermal distributions make minor ions sensitive probes of the collisionless processes that heat the corona and power the solar wind. The recent discovery of the "helicity barrier" offers a mechanism in which imbalanced Alfvénic turbulence in low-β plasmas preferentially heats protons over electrons, generating high-frequency, proton-cyclotron-resonant fluctuations. We use the hybrid-kinetic particle-in-cell code Pegasus++ to drive imbalanced Alfvénic turbulence in a 3D low-β plasma with additional passive ion species, He 2+ and O 5+ . A helicity barrier naturally develops, followed by clear phase-space signatures of oblique proton-cyclotron-wave heating and Landau-resonant heating from the imbalanced Alfvénic fluctuations. The former results in characteristically arced ion velocity distribution functions, whose non-bi-Maxwellian features are shown by linear ALPS calculations to be critical to the heating process. Additional features include a steep transition-range electromagnetic spectrum, proton-cyclotron waves propagating in the direction of the imbalance, significantly enhanced proton-to-electron heating ratios, ion temperatures that are considerably more perpendicular with respect to magnetic field, and extreme heating of heavier species in a manner consistent with mass scalings inferred from spacecraft measurements. None of these features are realized in an otherwise equivalent simulation of balanced turbulence. If seen simultaneously in the fast solar wind, these signatures of the helicity barrier would testify to the necessity of incorporating turbulence imbalance in a complete theory for the evolution of the solar wind.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Hybrid Kohn-Sham+Thomas-Fermi scheme for high-temperature density functional theory

We report density functional theory has revolutionised our understanding of dense, strongly coupled systems in the condensed-matter limit. At elevated temperatures, such as those found in warm dense matter and dense plasma, the computational cost of calculating enough states to accurately describe a thermal ensemble of electrons is a major obstacle to practical Kohn-Sham calculations. This problem has recently been tackled from a number of approaches, including a plane-wave approximation for electrons at high energy. In this work, we demonstrate how the Thomas-Fermi density of states, equivalent to a local plane-wave approximation, can yield improved results for thermodynamic quantities, electron density profiles, and interatomic forces while requiring fewer states than existing methods.

36 MATERIALS SCIENCE↗
NASA NTRS1989

Generation models of electron conics

Electron distribution functions (EDFs) with a peak oblique to the magnetic field, adjacent to but distinct from loss-cone features, have been observed by the DE 1, Viking, and S3-3 spacecraft in passes through the nightside auroral zone, polar cap, dayside cusp, and extended dayside auroral oval. Using particle simulations, two types of wave excitation and particle acceleration mechanisms which may contribute to producing these electron conic distributions are investigated. The first involves excitation of upper-hybrid waves by the electron loss cone, and the subsequent perpendicular heating of the background and thermal electrons. The second involves excitation of downward propagating parallel modes by an auroral electron beam which frequently accompanies the upflowing electron conics. These modes provide parallel acceleration, which modifies the GDF. Those electrons which are not lost to the atmosphere and mirror back up the magnetic field line give rise to enhancements in the GDF at the edge of the loss cone.

Roth, I.↗
DOE OSTI2021

Ab initio molecular dynamics on quantum computers

Ab initio molecular dynamics (AIMD) is a valuable technique for studying molecules and materials at finite temperatures where the nuclei evolve on potential energy surfaces obtained from accurate electronic structure calculations. In this work, we present an approach to running AIMD simulations on noisy intermediate-scale quantum (NISQ)-era quantum computers. The electronic energies are calculated on a quantum computer using the variational quantum eigensolver (VQE) method. Algorithms for computation of analytical gradients entirely on a quantum computer require quantum fault-tolerant hardware, which is beyond NISQ-era. Therefore, we compute the energy gradients numerically using finite differences, the Hellmann-Feynman theorem, and a correlated sampling technique. This method only requires additional classical calculations of electron integrals for each degree of freedom without any additional computations on a quantum computer beyond the initial VQE run. As a proof of concept, AIMD simulations are demonstrated for the H-2 molecule on IBM quantum devices. In addition, we demonstrate the validity of the method for larger molecules using full configuration interaction wave functions. As quantum hardware and noise mitigation techniques continue to improve, the method can be utilized for studying larger molecular systems.

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

Ionization state and dielectric constant in cold rarefied hydrocarbon plasmas of inertial confinement fusion

A combined approach to study cold rarefied matter is introduced that includes a semianalytical method based on the free-energy minimization and ab initio calculations based on the finite-temperature density-functional theory. The approach is used to calculate the ionization state of hydrocarbon (CH) under the shock-release conditions in inertial confinement fusion. Furthermore, the dielectric constant of CH is calculated using the Kubo-Greenwood formulation and contribution from atomic polarizabilities is found to be as important as the free-electron contribution. Using the ionization state and dielectric constant, the electron density profile in the rarefaction wave of the shock-release plasma is obtained.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Analysis of diagonal G and subspace W approximations within fully self-consistent GW calculations for bulk semiconducting systems

Fully self-consistent GW (sc-GW) methods are now available to evaluate quasiparticle and spectral properties of various molecular and bulk systems. However, such techniques based on the full matrix of G and W are computationally demanding. Additionally, the routinely used single-shot GW approximation (G 0 W 0 ) has an undesirable dependency on the choice of initial exchange-correlation functional. In the literature, many so-called self-consistent GW methods are based on diagonal approximation of G and low-ranking approximation of W. It is thus worth checking how good such approximations are in comparison with the full matrix method. In this work, we consider AlAs, AlP, GaP, and ZnS as the prototype systems to perform sc-GW calculations by expressing the full G matrix using a plane-wave basis set. We compared our sc-GW results with the diagonal G and subspace W approximated sc-GW results (sc-GW-diagG and sc-GW-subW methods). In the sc-GW-diagG method, interacting G is expanded in the eigenvectors of noninteracting G such that only diagonal elements are retained, whereas the number of eigenmodes is truncated in sc-GW-subW calculations. A systematic analysis of the results obtained from the above techniques is presented. The differences in the quasiparticle band gap between the approximated and the full matrix sc-GW approaches are mostly less than 1.7%, which validates such widely adopted approximations, and also shows how such low-ranking approximation can be used to include higher-order terms such as the vertex correction.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Self-organization of photoionized plasmas via kinetic instabilities

Abstract Self-organization in an unmagnetized collisionless plasma (in this paper) refers to formation of transient coherent structures such as collective oscillations (electrostatic waves) or magnetic fields resulting from so-called kinetic effects in the plasma. This topical review provides a comprehensive analysis of the self-organization of strong-field photoionized, non-equilibrium plasmas through kinetic instabilities. The authors propose and demonstrate a novel experimental platform that enables the formation of dense plasmas with known highly anisotropic and non-thermal electron velocity distribution functions on a timescale on the order of an inverse electron plasma frequency. We then show that such plasmas are highly susceptible to a hierarchy of kinetic instabilities, including two-stream, current filamentation and Weibel, that convert a fraction of the electron kinetic energy into electric and/or magnetic energy stored in self-organized structures. The electrostatic waves so produced are measured using a collective light (Thomson) scattering technique with femtosecond resolution as the kinetic instabilities aided by collisions eventually thermalize the plasma electrons. In addition, we describe a novel experimental technique that has made it possible to map the temporal evolution of the wavenumber spectrum of the thermal Weibel instability with picosecond resolution, which leads to the formation of quasi-static coherent magnetic fields with different topologies in photoionized plasmas. Finally, the paper summarizes the important results and discusses future directions on this topic.

Zhang, Chaojie (ORCID:0000000180353014)↗
NASA NTRS1988

Growth of electron plasma waves above and below f(p) in the electron foreshock

This paper investigates the conditions required for electron beams to drive wave growth significantly above and below the electron plasma frequency, f(p), by numerically solving the linear dispersion equation. It is shown that kinetic growth well below f(p) may occur over a broad range of frequencies due to the beam instability, when the electron beam is slow, dilute, and relatively cold. Alternatively, a cold or sharp feature at low parallel velocities in the distribution function may drive kinetic growth significantly below f(p). Kinetic broadband growth significantly above f(p) is explained in terms of faster warmer beams. A unified qualitative theory for the narrow-band and broad-band waves is proposed.

Cairns, Iver H.↗
DOE OSTI2021

Effects of multiply charged ions on microturbulence-driven electron transport in partially magnetized plasmas

Herein, nonlinear interaction between kinetic instabilities in partially magnetized plasmas in the presence of multiply-charged ion streams is investigated using kinetic simulations. It was observed by Hara and Tsikata [Phys. Rev. E 102, 023202 (2020)] that the axial ion-ion two-stream instability (IITSI) due to singly and doubly charged ion streams, coupled with the azimuthal electron cyclotron drift instability (ECDI), enhances cross-field electron transport. In the present study, it is observed that the addition of triply charged ions (as a third ion species) contributes to damping of the excited modes, leading to a reduction in the cross-field electron transport. The net instability-driven electron transport is shown to be a function not only of the azimuthal modes, such as the ECDI, but of the multiple ion species that dictate the development of additional plasma waves. It is found that trapping of the higher ion charge states within the plasma waves results in a broadening of the ion velocity distribution functions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Systematic comparison and cross-validation of fixed-node diffusion Monte Carlo and phaseless auxiliary-field quantum Monte Carlo in solids

Quantum Monte Carlo (QMC) methods are some of the most accurate methods for simulating correlated electronic systems. Here, we investigate the compatibility, strengths, and weaknesses of two such methods, namely, diffusion Monte Carlo (DMC) and auxiliary-field quantum Monte Carlo (AFQMC). The multideterminant trial wave functions employed in both approaches are generated using the configuration interaction using a perturbative selection made iteratively (CIPSI) technique. Complete basis-set full configuration interaction energies estimated with CIPSI are used as a reference in this comparative study between DMC and AFQMC. By focusing on a set of canonical finite–size solid-state systems, we show that both QMC methods can be made to systematically converge towards the same energy once basis-set effects and systematic biases have been removed. AFQMC shows a much smaller dependence on the trial wave function than DMC while simultaneously exhibiting a much larger basis-set dependence. We outline some of the remaining challenges and opportunities for improving these approaches.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Confinement-Induced Chiral Edge Channel Interaction in Quantum Anomalous Hall Insulators

In quantum anomalous Hall (QAH) insulators, the interior is insulating but electrons can travel with zero resistance along one-dimensional (1D) conducting paths known as chiral edge channels (CECs). These CECs have been predicted to be confined to the 1D edges and exponentially decay in the two-dimensional (2D) bulk. Here, in this Letter, we present the results of a systematic study of QAH devices fashioned in a Hall bar geometry of different widths under gate voltages. At the charge neutral point, the QAH effect persists in a Hall bar device with a width of only ~7⁢2 nm, implying the intrinsic decaying length of CECs is less than ~3⁢6 nm. In the electron-doped regime, we find that the Hall resistance deviates quickly from the quantized value when the sample width is less than 1 μ⁢m. Our theoretical calculations suggest that the wave function of CEC first decays exponentially and then shows a long tail due to disorder-induced bulk states. Therefore, the deviation from the quantized Hall resistance in narrow QAH samples originates from the interaction between two opposite CECs mediated by disorder-induced bulk states in QAH insulators, consistent with our experimental observations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Spin-orbit-proximitized ferromagnetic metal by monolayer transition metal dichalcogenide: Atlas of spectral functions, spin textures, and spin-orbit torques in Co/MoSe 2 , Co/WSe 2 and Co/TaSe 2 heterostructures

In this work, the heterostructures composed of an ultrathin ferromagnetic metal (FM) and a material hosting strong spin-orbit (SO) coupling are the principal resource for SO torque and spin-to-charge conversion nonequilibrium effects in spintronics. The key quantity in theoretical description of these effects is nonequilibrium spin density, which can appear on any monolayer of the heterostructure through which the current is flowing as long as the monolayer bands are affected by the native or proximity induced SO coupling. Here we demonstrate how hybridization of wave functions of Co layer and a monolayer of transition metal dichalcogenides (TMDs)—such as semiconducting MoSe 2 and WSe 2 or metallic TaSe 2 —can lead to dramatic transmutation of electronic and spin structure of Co within some distance away from its interface with TMD, when compared to the bulk of Co or its surface in contact with a vacuum. This is due to proximity induced SO splitting of Co bands encoded in the spectral functions and spin textures on its monolayers, which we obtain using noncollinear density functional theory (ncDFT) combined with equilibrium Green function (GF) calculations. In fact, SO splitting is present due to structural inversion asymmetry of the bilayer—i.e., just the presence of the interface—even if SO coupling within TMD monolayer is artificially switched off in ncDFT calculations. However, switching it on makes the effects associated with proximity SO coupling within Co layer about five times larger. Injecting spin-unpolarized charge current through SO-proximitized monolayers of Co generates nonequilibrium spin density over them, so that its cross product with the magnetization of Co determines SO torque. The SO torque computed via first-principles quantum transport methodology, which combines ncDFT with nonequilibrium GF calculations, can be used as the screening parameter to identify optimal combination of materials and their interfaces for applications in spintronics. In particular, we identify heterostructure two-monolayer-Co/monolayer-WSe 2 as the most optimal one, at least in the clean limit.

2-dimensional systems↗
NASA NTRS1972

Theory on the precipitation of magnetospheric electrons at the time of a sudden commencement.

At the time of a sudden commencement or even of a sudden impulse, very-low-frequency waves can be excited. In turn, the wave turbulence drives the pitch-angle diffusion of energetic electrons trapped in the magnetosphere, causing them to be precipitated into the lower ionosphere. The distribution of this precipitation as a function of geomagnetic latitude and local time is shown to be directly related to the characteristics of the magnetic disturbances.

Perona, G. E.↗
NASA NTRS1984

Path-integrated growth of auroral kilometric radiation

Using Poeverlein's graphical method, three dimensional ray path calculations are performed to evaluate the path-integrated growth of auroral kilometric radiation (AKR). The ray tracing results indicate that waves whose initial wave vector lie in the local meridian plane continue to propagate in that plane and that among these waves, those with frequencies near the cutoff frequency (f sub R = 0) refract substantially, where as those with frequencies well above the cutoff frequency suffer little refraction. It is also shown that waves whose initial wave vector lie outside of the local meridian plane propagate in the longitudinal as well as the radial and the latitudinal directions. The refraction of these waves is also highly dependent upon the wave frequency, i.e., waves with frequencies near f sub R = 0 refract substantially, whereas waves with frequencies much above f sub R = 0 undergo little refraction. In order to test the electron cyclotron maser mechanisms as a method for generation of AKR, a typical electron distribution function measured in the auroral zone by the S3-3 satellite, is used to calculate path-integrated growths of representative rays. The results of this study indicate that electron distribution functions like those measured by the S3-3 satellite are not capable of amplifying cosmic noise background to the observed intensities of auroral kilometric radiation, and that much steeper slopes at the edges of the loss cone are required. The presence of such distribution functions in the auroral zone is plausible if one assumes that backscattered electrons in this region have energies less than a few hundred eV.

Omidi, N.↗
NASA NTRS1986

The synchrotron-maser theory of type II solar radio emission processes - The physical model and generation mechanism

A theory is proposed to explain the generation mechanism of type II solar radio bursts. It is suggested that the shock wave formed at the leading edge of a coronal transient can accelerate electrons. Because of the nature of the acceleration process, the energized electrons can possess a 'hollow-beam' type distribution function. When the electron beam propagates along the ambient magnetic field to lower altitudes and attains larger pitch angles, a synchrotron-maser instability can set in. This instability leads to the amplification of unpolarized or weakly polarized radiation. The present discussion incorporates a model which describes the ambient magnetic field and background plasma by means of MHD simulation. The potential emission regions may be located approximately, according to the time-dependent MHD simulation. Since the average local plasma frequency in the source region can be evaluated from the MHD model, the frequent drift associated with the radiation may be estimated. The result seems to be in good agreement with that derived from observations.

Wu, C. S.↗
NASA NTRS1997

Effective Collision Strengths for Electron Impact Excitation of Inelastic Transitions in S III

We have calculated electron collisional excitation strengths for all electric dipole forbidden, semi-forbidden, and allowed transitions among the lowest 17 LS states 3s(exp 2)3p(exp 2) P-3, D-1, S-1, 3s3p(exp 3)S-5(exp 0), D-3(exp 0), P-3(exp 0), P-1(exp 0), S-3(exp 0), D-1(exp 0), 3S(exp 2)3p3d D-1(exp 0), F-3(exp 0), P-3(exp 0), D-3(exp 0), F-3(exp 0), P-1(exp 0), and 3S(exp 2)3p4S P-3(exp 0), P-l(exp 0) of S III using the R-matrix method. These S m states are represented by fairly extensive configuration-interaction wave functions that yield excited state energies in close agreement with recent laboratory measurements. Rydberg series of resonances converging to the excited state thresholds are explicitly included in the scattering calculation. The effective collision strengths are determined assuming Maxwellian distribution of electron energies. These are listed over a wide temperature range ([0.5-10] x 10(exp 4) K) and compared, where possible, with other available calculations. Subject headings: atomic data - atomic processes

Tayal, S. S.↗
DOE OSTI2023

Trends in electronic structures and s ± -wave pairing for the rare-earth series in bilayer nickelate superconductor R 3 Ni 2 O 7

The recent discovery of pressure-induced superconductivity in the bilayer La 3 Ni 2 O 7 (LNO) has opened a new platform for the study of unconventional superconductors. In this publication, we investigate theoretically the whole family of bilayer 327-type nickelates R 3 Ni 2 O 7 (R = rare-earth elements) under pressure. From La to Lu, the lattice constants and volume decrease, leading to enhanced in-plane and out-of-plane hoppings, resulting in an effectively reduced electronic correlation U / W. Furthermore, the Ni's t 2g states shift away from the e g states, while the crystal-field splitting between d 3z 2 -r 2 and d x 2 -y 2 is almost unchanged. In addition, six candidates were found to become stable in the Fmmm phase, with increasing values of critical pressure as the atomic number increases. Similar to the case of LNO, the s ± -wave pairing tendency dominates in all candidates, due to the nesting between the M = (π , π) and the X = (π , 0) and Y = (0 , π) points in the Brillouin zone. Then, T c is expected to decrease as the radius of rare-earth (RE) ions decreases. In conclusion, our results suggest that LNO is already the “optimal” candidate, with Ce a close competitor, among the whole of the RE bilayer nickelates, and to increase T c we suggest growing on special substrates with larger in-plane lattice spacings.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗