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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Metastable Bound States of the Two-Dimensional Bimagnetoexcitons in the Lowest Landau Levels Approximation

The possible existence of the bound states of the interacting two-dimensional (2D) magnetoexcitons in the lowest Landau levels (LLLs) approximation was investigated using the Landau gauge description. The magnetoexcitons taking part in the formation of the bound state with resultant wave vector k=0 have opposite in-plane wave vectors k and −k and look as two electric dipoles with the arms oriented in-plane perpendicularly to the corresponding wave vectors. The bound state of two antiparallel dipoles moving with equal probability in any direction of the plane with equal but antiparallel wave vectors is characterized by the variational wave function of the relative motion φ{sub n}(k) depending on the modulus |k|. The spins of two electrons and the effective spins of two holes forming the bound states were combined separately in the symmetric or in the antisymmetric forms (↑↓ +η ↓↑) with the same parameter η=±1 for electrons and holes. In the case of the variational wave function φ{sub 2}(k)=(8α{sup 3}){sup 1/2}k{sup 2}l{sub 0}{sup 2}exp [−αk{sup 2}l{sub 0}{sup 2}] the maximum density of the magnetoexcitons in the momentum space representation is concentrated on the in-plane ring with the radius k{sub r}=1/(√αl{sub 0}). The stable bound states of the bimagnetoexciton molecule do not exist for both spin orientations. Instead of them, a deep metastable bound state with the activation barrier comparable with the ionization potential I{sub l} of the magnetoexciton with k=0 was revealed in the case η=1 and α=0.5. In the case η=−1 and α=3.4 only a shallow metastable bound state can appear.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

TurboRVB: A many-body toolkit for ab initio electronic simulations by quantum Monte Carlo

TurboRVB is a computational package for ab initio Quantum Monte Carlo (QMC) simulations of both molecular and bulk electronic systems. The code implements two types of well established QMC algorithms: Variational Monte Carlo (VMC) and diffusion Monte Carlo in its robust and efficient lattice regularized variant. A key feature of the code is the possibility of using strongly correlated many-body wave functions (WFs), capable of describing several materials with very high accuracy, even when standard mean-field approaches [e.g., density functional theory (DFT)] fail. The electronic WF is obtained by applying a Jastrow factor, which takes into account dynamical correlations, to the most general mean-field ground state, written either as an antisymmetrized geminal power with spin-singlet pairing or as a Pfaffian, including both singlet and triplet correlations. This WF can be viewed as an efficient implementation of the so-called resonating valence bond (RVB) Ansatz, first proposed by Pauling and Anderson in quantum chemistry [L. Pauling, The Nature of the Chemical Bond (Cornell University Press, 1960)] and condensed matter physics [P.W. Anderson, Mat. Res. Bull 8, 153 (1973)], respectively. The RVB Ansatz implemented in TurboRVB has a large variational freedom, including the Jastrow correlated Slater determinant as its simplest, but nontrivial case. Moreover, it has the remarkable advantage of remaining with an affordable computational cost, proportional to the one spent for the evaluation of a single Slater determinant. Therefore, its application to large systems is computationally feasible. The WF is expanded in a localized basis set. Several basis set functions are implemented, such as Gaussian, Slater, and mixed types, with no restriction on the choice of their contraction. The code implements the adjoint algorithmic differentiation that enables a very efficient evaluation of energy derivatives, comprising the ionic forces. Thus, one can perform structural optimizations and molecular dynamics in the canonical NVT ensemble at the VMC level. For the electronic part, a full WF optimization (Jastrow and antisymmetric parts together) is made possible, thanks to state-of-the-art stochastic algorithms for energy minimization. In the optimization procedure, the first guess can be obtained at the mean-field level by a built-in DFT driver. The code was efficiently parallelized by using a hybrid MPI-OpenMP protocol, which is also an ideal environment for exploiting the computational power of modern Graphics Processing Unit accelerators.

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

SchrödingerNet: A Universal Neural Network Solver for the Schrödinger Equation

Recent advances in machine learning have facilitated numerically accurate solution of the electronic Schrödinger equation (SE) by integrating various neural network (NN)-based wave function ansatzes with variational Monte Carlo methods. Nevertheless, such NN-based methods are all based on the Born–Oppenheimer approximation (BOA) and require computationally expensive training for each nuclear configuration. In this work, we propose a novel NN architecture, SchrödingerNet, to solve the full electronic-nuclear SE by defining a loss function designed to equalize local energies across the system. This approach is based on a translationally, rotationally and permutationally symmetry-adapted total wave function ansatz that includes both nuclear and electronic coordinates. Furthermore, this strategy not only allows for an efficient and accurate generation of a continuous potential energy surface at any geometry within the well-sampled nuclear configuration space, but also incorporates non-BOA corrections, through a single training process. Comparison with benchmarks of atomic and small molecular systems demonstrates its accuracy and efficiency.

Chemical calculations↗
DOE OSTI2020

Ground state wave functions for single-band Hubbard models from the Gutzwiller conjugate gradient minimisation theory

The Gutzwiller conjugate gradient minimisation (GCGM) theory is an ab initio quantum many-body theory for computing the ground-state properties of infinite systems. Previous applications of GCGM provides satisfying accuracy of ground-state energy of Hubbard models. In the current work, we address the problem of whether the obtained wave function is a good approximation for the true ground state by comparing the correlation functions with the benchmark data. Additionally, our results confirms the accuracy of the reproduced ground state of the regular Hubbard model, but with some exception for the frustrated Hubbard model.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2018

Locally coupled open subsystems: A formalism for affordable electronic structure calculations featuring fractional charges and size consistency

This manuscript introduces a methodology (within the Born-Oppenheimer picture) to compute electronic ground-state properties of molecules and solids/surfaces with fractionally occupied components. Given a user-defined division of the molecule into subsystems, our theory uses an auxiliary global Hamiltonian that is defined as the sum of subsystem Hamiltonians, plus the spatial integral of a second-quantized local operator that allows the electrons to be transferred between subsystems. This electron transfer operator depends on a local potential that can be determined using density functional approximations and/or other techniques such as machine learning. The present framework employs superpositions of tensor-product wave functions, which can satisfy size consistency and avoid spurious fractional charges at large bond distances. The electronic population of each subsystem is in general a positive real number and is obtained from wave-function amplitudes, which are calculated by means of ground-state matrix diagonalization (or matrix propagation in the time-dependent case). Furthermore, our method can provide pathways to explore charge-transfer effects in environments where dividing the molecule into subsystems is convenient and to develop computationally affordable electronic structure algorithms.

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

Application of the Screened Hydrogenic Model to Light Atoms

The purpose of this work was to develop a Screened Hydrogenic Model (SHM) to accurately calculate the electron energies for light atoms and ions with up to ten electrons for atomic numbers up to 18. The total energy of an atom or ion was calculated with effective nuclear charges and screening parameters for each electron type (1s, 2s and 2p) within a specific electron configuration. Multiple energy states, centered at the total energy, were calculated for elec-tron configurations that have Russell Saunders coupling. The energy of each electron included its relativistic energy, EREL, but close overall agreement between the calculated and experimental energies of multi-electron configura-tions required that the one-electron expression for EREL be modified in a simple manner. In the present work, 98% of the 587 calculated energies for light atoms/ions have a relative error within ±0.1% of the corresponding experimental energies. The effective nuclear charges described in this work allow hydrogen-like wave functions to be defined for the electrons within a multi-electron configuration. The SHM, described in this work, is available for future calculations involving light atoms and ions.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Shakeup and shakeoff spectra in the electron capture decay of atomic 7 Be

The most stringent laboratory-based experimental limits on the existence of submegaelectronvolt sterile neutrinos are currently set by decay spectroscopy of radioactive 7 Be embedded into superconducting sensors. The systematic uncertainties are dominated by the modeling of the electron shakeup and shakeoff spectra that are not based on state-of-the-art atomic theory and do not include electron correlations or relativistic effects. We have used the multiconfiguration Dirac-Fock formalism to obtain correlated wave functions ab initio and compute all single and double shake processes in the electron capture decay of atomic 7 Be . The simulations can explain some but not all of the observed spectral features, likely because the wave functions are modified by the Ta sensor material that the 7 Be is embedded into. The new models also show that the 𝐿/𝐾 electron capture ratio of 7 Be in Ta has previously been slightly underestimated revising the previous value of 0.070(7) to a new value of 0.0756(20).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

A parallel orbital-updating based plane-wave basis method for electronic structure calculations

Highlights: • Propose three parallel orbital-updating based plane-wave basis methods for electronic structure calculations. • These new methods can avoid the generating of large scale eigenvalue problems and then reduce the computational cost. • These new methods allow for two-level parallelization which is particularly interesting for large scale parallelization. • Numerical experiments show that these new methods are reliable and efficient for large scale calculations on modern supercomputers. - Abstract: Motivated by the recently proposed parallel orbital-updating approach in real space method , we propose a parallel orbital-updating based plane-wave basis method for electronic structure calculations, for solving the corresponding eigenvalue problems. In addition, we propose two new modified parallel orbital-updating methods. Compared to the traditional plane-wave methods, our methods allow for two-level parallelization, which is particularly interesting for large scale parallelization. Numerical experiments show that these new methods are more reliable and efficient for large scale calculations on modern supercomputers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Whistler wave generation by electron temperature anisotropy during magnetic reconnection at the magnetopause

Two magnetopause reconnection events of the Magnetospheric Multiscale mission with whistler wave activity are presented. The whistler mode around half of the electron cyclotron frequency is excited near the magnetospheric separatrix. In both events, there are positive correlations between the whistler wave and the lower hybrid drift instability (LHDI) activities, suggesting a possible role of LHDI in the whistler wave generation. A sudden change in the electron pitch angle distribution (PAD) function of energetic electrons is observed right after intense LHDI activity. This change in the PAD leads to temperature anisotropy in energetic electrons which is responsible for the whistler wave excitation. The measured dispersion relation demonstrates that the whistler wave propagates toward the X line nearly parallel to the magnetic field line. Furthermore, a linear analysis with the measured distribution function verifies that the whistler mode is excited by the temperature anisotropy in energetic electrons.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

An efficient random-sampling method for calculating double occupancy of Gutzwiller wave function in single-band 1D and 2D lattices

In this paper, we report a random sampling method for computing the expectation value of physical quantities based on the Gutzwiler variational wave function. As the first application, we calculated the double occupancy, which is a critical quantity for under- standing the correlation effects in many-body systems, for single-band 1D and 2D lattices. We demonstrated that the random sampling scheme is more efficient than an existing Metropolis Monte-Carlo algorithm. For the 1D Hubbard model with only nearest-neighbor hopping, our results are almost identical to the exact analytic solution. We have also studied systems to which analytic solutions are not available, including the 1D lattices with next-nearest-neighbor hopping and 2D lattices. In addition, constraints on real-space con gurations can be easily implemented in the current scheme to further improve the Gutzwiller wave function. As an example, we calculated the double occupancy for 1D Hubbard model by applying the constraint that all double-occupied sites are paired with an empty site. With enhanced correlation between double-occupied and empty sites, the constraint results in much improved ground-state energy for 1D Hubbard model with strong on-site repulsion.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2019

Spin–Orbit Coupling Constants in Atoms and Ions of Transition Elements: Comparison of Effective Core Potentials, Model Core Potentials, and All-Electron Methods

The spin–orbit coupling constants (SOCC) in atoms and ions of the first- through third-row transition elements were calculated for the low-lying atomic states whose main electron configuration is [ n d] q ( q = 1–4 and 6–9, n = the principal quantum number), using four different approaches: (1) a nonrelativistic Hamiltonian used to construct multiconfiguration self-consistent field (MCSCF) wave functions utilizing effective core potentials and their associated basis sets within the framework of second-order configuration interaction (SOCI) to calculate spin–orbit couplings (SOC) using one-electron Breit–Pauli Hamiltonian (BPH), (2) a nonrelativistic Hamiltonian used to construct MCSCF wave functions utilizing model core potentials and their associated basis sets within the framework of SOCI to calculate SOC using the full BPH, (3) nonrelativistic and spin-independent relativistic Hamiltonians used to construct MCSCF wave functions utilizing all-electron (AE) basis sets within the framework of SOCI to calculate SOC using the full BPH, and (4) a relativistic Hamiltonian given by the exact two-component (X2C) transformation for construction of Kramers-restricted relativistic configuration interaction wave functions. In this investigation, these four approaches are referred to as ECP, MCP, AE, and X2C methods, respectively. The ECP, MCP, and AE methods are so-called two-step approaches (TSA), while the X2C method is a one-step approach (OSA). In the AE method, three different calculations—relativistic elimination of small components (RESC), third-order Douglas–Kroll–Hess (DKH3), and infinite-order two-component (IOTC) relativistic correction—were performed for the estimation of the scalar relativistic components in addition to those of the nonscalar relativistic (NSR) contributions. The calculated SOCC are compared to the available experimental data via the Landé interval rule. Although there are several exceptions, including states whose main configuration is [ n d] 5 , the average differences between the ECP and AE (IOTC) SOCC and between the ECP and the X2C SOCC are mostly less than 20%. The differences between the ECP and the experimental SOCC are even smaller. No serious discrepancy was found between the TSA and OSA predictions of SOCC for the first- and second-row transition elements in comparison to experiment. For atoms and ions of the third-row transition elements, the SOCC calculated through the Landé interval rule are not reliable. The low-energy spin-mixed (SM) states originating from a [5d]q configuration ( q = 2–4) have a larger energy lowering due to the SOC effects, in comparison with those for atoms and ions of the first- and second-row transition elements. For the spin-mixed (SM) states originating from a [5d]q configuration ( q = 6–8), the energy lowering of all 4 F 7/2 , 5 D 1 , and 5 D 3 states due to the SOC effects is smaller than those of the other SM states. Furthermore this difficulty, which also arises for the MCP, AE, and X2C (OSA) approaches, suggests that the LS-coupling scheme is inappropriate.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Electron transfer calculations between edge sharing octahedra in hematite, goethite, and annite

A key reaction underlying the charge transport in iron containing oxides, clays, micas is the Fe$^{2+}$-Fe$^{3+}$ exchange reaction between edge-sharing iron octahedra. These reactions facilitate conduction in these minerals by the thermally-activated hopping of small polarons across the lattice. Depending on the mineral and local charge state the small polaron can either encase an electron or hole. The probability for conduction of small polarons depends strongly on the height and adiabicity of the reaction barrier, with larger and more diabatic barriers yielding slow conduction associated with either weak coupling or a large prerequisite rearrangement of the lattice during charge transport. To model these reactions, a first principle electron transfer (ET) method was developed to model the small polaron hopping between the edge-sharing octahedra sites in hematite ($e^{-}$ polaron), goethite ($e^{-}$ polaron), and annite ($h^{+}$ polaron) bulk structures. The ET method is based on electronic structure methods (i.e., plane-wave Density Functional Theory) capable of performing calculations with periodic cells and large size systems efficiently while at the same time being accurate enough to be used in the estimation of the electron-transfer coupling matrix element, $V_{AB}$, and the electron transfer transmission factor, $\kappa_{el}$. Additionally, the calculations confirmed the existence of small polarons in all three minerals, and the reactions were predicted to be strongly adiabatic. It was found that transfer of a hole in the octahedral layer of annite had an adiabatic barrier of $0.311$ eV, and the transfer of an extra electron in hematite and goethite had adiabatic barriers of $0.242$ eV and $0.232$ eV respectively. The electronic coupling parameters,$V_{AB}$, were found to be $0.188$ eV, $0.196$ eV, and $0.102$ eV respectively for hematite, goethite, and annite. While similar bonding topologies pertain, the underlying basis for the differences is the subtle differences in local structures.

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

The role of high-order electron correlation effects in a model system for non-valence correlation-bound anions

The diffusion Monte Carlo (DMC), auxiliary field quantum Monte Carlo (AFQMC), and equation-of-motion coupled cluster (EOM-CC) methods are used to calculate the electron binding energy (EBE) of the non-valence anion state of a model (H2O)4 cluster. Two geometries are considered, one at which the anion is unbound and the other at which it is bound in the Hartree–Fock (HF) approximation. It is demonstrated that DMC calculations can recover from the use of a HF trial wave function that has collapsed onto a discretized continuum solution, although larger EBEs are obtained when using a trial wave function for the anion that provides a more realistic description of the charge distribution and, hence, of the nodal surface. For the geometry at which the cluster has a non-valence correlation-bound anion, both the inclusion of triples in the EOM-CC method and the inclusion of supplemental diffuse d functions in the basis set are important. DMC calculations with suitable trial wave functions give EBE values in good agreement with our best estimate EOM-CC result. AFQMC using a trial wave function for the anion with a realistic electron density gives a value of the EBE nearly identical to the EOM-CC result when using the same basis set. For the geometry at which the anion is bound in the HF approximation, the inclusion of triple excitations in the EOM-CC calculations is much less important. The best estimate EOM-CC EBE value is in good agreement with the results of DMC calculations with appropriate trial wave functions.

Upadhyay, Shiv (ORCID:0000000285010501)↗
DOE OSTI2024

Unveiling the Redox Noninnocence of Metallocorroles: Exploring K-Edge X-ray Absorption Near-Edge Spectroscopy with a Multiconfigurational Wave Function Approach

X-ray absorption near-edge spectroscopy (XANES) is an advanced technique for probing the local electronic structure of catalysts, effectively identifying the noninnocent nature of ligands in transition-metal complexes. Metallocorroles with noninnocent corrole rings exhibit unusual electronic structures that challenge traditional density functional theory (DFT) methods, necessitating more rigorous approaches to describe electron correlation accurately. We explored K-edge XANES spectra of Fe, Mn, and Co metallocorroles using TDDFT and wave function-based methods. This is the first investigation employing multireference methods, specifically RASSCF, RASPT2, and MC-PDFT, to analyze the redox noninnocent nature of metallocorroles reflected in their XANES spectra. We quantified the noninnocent character of the corrole and the oxidation states of the metals, capturing more than singly excited excitations responsible for the pre-edge peak. Our findings demonstrate the importance of these advanced computational techniques for accurately predicting XANES spectra, providing a reliable understanding of the electronic properties of such complexes. In conclusion, this study offers a new strategy for investigating ligand redox noninnocence via integrated experimental and computational XANES.

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

New Dimensions in the Theory of Excited States and X-ray Spectra (Final Report)

This Final Technical Report briefly summarizes the achievements during the lifetime of our DOE BES grant DE-FG02-97ER45623. The long-term goal of this project has been the development of quantitative theories of the interaction between radiation and matter, with a focus on x-ray spectroscopies. X-ray spectra have long been among the most important probes of atomic-scale structure and properties of matter, ranging from atoms and molecular systems to condensed matter and exotic states. These spectroscopies are widely used in investigations at the major DOE synchrotron x-ray facilities and related centers world-wide. In addition to fundamental theory, a major goal of our project has been the development of computational software that implements the theory for calculations of x-ray spectra of various materials throughout the periodic table. Due to the complex nature of x-ray spectra, quantitative theory is essential for its interpretation. The theory is challenging since it involves excited state electronic structure and many-body correlation effects that go beyond independent particle approximations like DFT or Hartree-Fock. Moreover, the experimental investigations typically involve a broad range of energy, time, and temperature scales, from the UV-Vis to hard x-ray energies of order 10 4 eV, and temperatures T from ambient up to the warm-dense-matter regime where the Fermi energy kBTF is of order a few eV, i.e., temperatures of order 105 Kelvin. This broad range of experimental conditions has fostered many novel theoretical approaches and computational techniques, many of which we have developed systematically over the duration of the grant. In contrast to the traditional wave-function approach of quantum theory and electronic structure methods, our theoretical approach is based on modern Green's function techniques. This approach is better suited for aperiodic structures, excited states, and broad spectral ranges, since it avoids the computational bottlenecks of sum-over-states approaches, as in the Fermi golden rule. This theoretical framework has been incorporated into efficient, user-friendly x-ray spectroscopy software which is now used routinely worldwide to simulate and analyze spectra. These theoretical tools provide an essential complement to synchrotron and next-generation light sources, which are used to investigate complex materials with ever increasing precision. Moreover, the synergism between theory, computation and experiment contributed by our research enhances scientific understanding and creates opportunities for innovations in materials and energy science and in many fields. As documented in this Report, this research grant has been remarkably successful in achieving these goals. In particular, this grant has supported the development of the x-ray spectroscopy software suite known as FEFF (named for an effective scattering amplitude f eff in the theory). The FEFF codes have become one of the premier tools for quantitative simulations of x-ray spectra as documented by many thousands of citations in the Web of Science and Google-Scholar.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Taming the fixed-node error in diffusion Monte Carlo via range separation

By combining density-functional theory (DFT) and wave function theory via the range separation (RS) of the interelectronic Coulomb operator, we obtain accurate fixed-node diffusion Monte Carlo (FN-DMC) energies with compact multi-determinant trial wave functions. In particular, we combine here short-range exchange-correlation functionals with a flavor of selected configuration interaction known as configuration interaction using a perturbative selection made iteratively (CIPSI), a scheme that we label RS-DFT-CIPSI. One of the take-home messages of the present study is that RS-DFT-CIPSI trial wave functions yield lower fixed-node energies with more compact multi-determinant expansions than CIPSI, especially for small basis sets. Indeed, as the CIPSI component of RS-DFT-CIPSI is relieved from describing the short-range part of the correlation hole around the electron-electron coalescence points, the number of determinants in the trial wave function required to reach a given accuracy is significantly reduced as compared to a conventional CIPSI calculation. Importantly, by performing various numerical experiments, we evidence that the RS-DFT scheme essentially plays the role of a simple Jastrow factor by mimicking short-range correlation effects, hence avoiding the burden of performing a stochastic optimization. Considering the 55 atomization energies of the Gaussian-1 benchmark set of molecules, we show that using a fixed value of mu = 0.5 bohr(-1) provides effective error cancellations as well as compact trial wave functions, making the present method a good candidate for the accurate description of large chemical systems.

Scemama, Anthony↗
DOE OSTI2019

Computational investigation of ion cyclotron heating on Proto-MPEX

Ion cyclotron heating (ICH) on the Prototype Material Plasma Exposure eXperiment (Proto-MPEX) is to be accomplished using the “beach-heating” technique. Beach heating has not been previously demonstrated to efficiently heat core ions at the high electron density values present in Proto-MPEX. This work numerically investigates the wave propagation characteristics of the ICH region on Proto-MPEX to explore avenues for efficient core ion heating. The analysis reveals that finite electron temperature effects are required to predict core ion heating. Cold plasma dispersion analysis and full-wave simulations show that the inertial Alfvén wave (IAW) is restricted from coupling power into the core plasma because (1) the group velocity is too shallow for its energy to penetrate into the core before damping in the periphery and (2) when operating in a magnetic field where ω/ω ci ≳ 0.7, the IAW is cut off from the core plasma by the Alfvén resonance. However, including kinetic effects shows that the kinetic Alfvén wave (KAW) can propagate in the electron temperature regime in Proto-MPEX. Full-wave simulations show that when the electron temperature is increased to T e > 2 eV and the edge electron density is sufficiently high n e edge > 1 × 10 17 m -3 , ion power absorption in the core increases substantially (≈25% of total power). Here, the increase in ion power absorption in the core is attributed to the propagation of the KAW. Calculations of electron and ion power absorption show that the electron heating is localized around the Alfvén resonance, while the ion heating is localized at the fundamental ion cyclotron resonance.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

The Nature of Electron Mobility in Hybrid Perovskite CH 3 NH 3 PbI 3

The peroskovite CH 3 NH 3 PbI 3 one of the most promising candidates for cheap and high-efficiency solar cells. One of its unique features is the long carrier diffusion length (>100 μm), but its carrier mobility is rather modest. The nature of the mobility is unclear. Here, using nonadiabatic wave function dynamics simulations, we show that the random rotations of the CH 3 NH 3 cations play an important role in the carrier mobility. Our previous work showed that the electron and hole wave functions were localized and spatially separated due to the random orientations of the CH 3 NH 3 cations in the tetragonal phase. We find that the localized carriers are able to conduct random walks due to the electrostatic potential fluctuation caused by the CH 3 NH 3 random rotations. The calculated electron mobilities are in the experimentally measured range. We thus conclude that the carrier mobility of CH 3 NH 3 PbI 3 is likely driven by the dynamic disorder that causes the fluctuation of the electrostatic potential.

14 SOLAR ENERGY↗