Search NASA⌕ Search

SEARCH · Search NASA

Results for “Screened coulomb interaction”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

A Filon-like integration strategy for calculating exact exchange in periodic boundary conditions: a plane-wave DFT implementation

An efficient and accurate approach for calculating exact exchange and other two-electron integrals has been developed for periodic electronic structure methods. Traditional approaches used for integrating over the Brillouin zone in band structure calculations, e.g. trapezoidal or Monkhorst-Pack, are not accurate enough for two-electron integrals. This is because their integrands contain multiple singularities over the double integration of the Brillouin zone, which with simple integration methods lead to very inaccurate results. A common approach to this problem has been to replace the Coulomb interaction with a screened Coulomb interaction that removes singularities from the integrands in the two-electron integrals, albeit at the inelegance of having to introduce a screening factor which must precomputed or guessed. Instead of introducing screened Coulomb interactions in an ad hoc way, the method developed in this work derives an effective screened potential using a Filon-like integration approach that is based only on the lattice parameters. This approach overcomes the limitations of traditionally defined screened Coulomb interactions for calculating two-electron integrals, and makes chemistry many-body calculations tractable in periodic boundary conditions. This method has been applied to several systems for which conventional DFT methods do not work well, including the reaction pathways for the addition of H 2 to phenol and Au$^{–}_{20}$ nanoparticle, and the electron transfer of a charge trapped state in the Fe(II) containing mica, annite.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Downfolding from ab initio to interacting model Hamiltonians: comprehensive analysis and benchmarking of the DFT+cRPA approach

Abstract Model Hamiltonians are regularly derived from first principles to describe correlated matter. However, the standard methods for this contain a number of largely unexplored approximations. For a strongly correlated impurity model system, here we carefully compare a standard downfolding technique with the best possible ground-truth estimates for charge-neutral excited-state energies and wave functions using state-of-the-art first-principles many-body wave function approaches. To this end, we use the vanadocene molecule and analyze all downfolding aspects, including the Hamiltonian form, target basis, double-counting correction, and Coulomb interaction screening models. We find that the choice of target-space basis functions emerges as a key factor for the quality of the downfolded results, while orbital-dependent double-counting corrections diminish the quality. Background screening of the Coulomb interaction matrix elements primarily affects crystal-field excitations. Our benchmark uncovers the relative importance of each downfolding step and offers insights into the potential accuracy of minimal downfolded model Hamiltonians.

Chemistry↗
DOE OSTI2026

Quantitative SANS and multi-model analysis of spacer-dependent micellization of urea-based gemini surfactants

The micellization behavior of urea-based cationic gemini surfactants was investigated using small-angle neutron scattering (SANS) with multi-model form factor analysis. A homologous series of surfactants with urea group included in the hydrophobic tail and polymethylene spacers consisting of two to ten methylene units was analyzed using three form factor models: a core–shell ellipsoid and two variants of homogeneous ellipsoids. The results from all models show a consistent trend of the micelle structures, confirming that the spacer length critically influences micellar geometry, aggregation number, and hydration. The surfactant with four CH 2 groups in the spacer formed the largest micelles with the highest aggregation number, while longer spacers led to progressively smaller, more compact aggregates. The shell hydration—quantified as the volume fraction of heavy water within the hydrophilic region—decreased systematically with increasing spacer length due to enhanced hydrophobicity of the headgroup-spacer region. Intermicellar interactions, modeled as screened Coulomb interaction using the rescaled mean spherical approximation (RMSA), revealed the strongest electrostatic repulsion for the case of four methylene groups in the spacer, corresponding to the highest micellar charge and largest interparticle spacing. The observed spacer-dependent trends were robust across all modeling approaches, demonstrating that the spacer length serves as a key structural determinant of self-assembly in this type of urea-based gemini systems. These findings provide insight into the design of gemini surfactants with tailored aggregation behavior for applications in drug delivery, nanostructure templating, and solubilization technologies.

Core–shell ellipsoid model↗
DOE OSTI2023

Improved quasiparticle self-consistent electronic band structure and excitons in β – LiGaO 2

The band structure of β–LiGaO 2 is calculated using the quasiparticle self-consistent QSGWˆ method where the screened Coulomb interaction Wˆ is evaluated including electron-hole interaction ladder diagrams and G is the one-electron Green's function. Improved convergence compared to previous calculations leads to a significantly larger band gap of about 7.0 eV. However, exciton binding energies are found to be large and lead to an exciton gap of about 6.0 eV if also a zero-point-motion correction of about –0.4 eV is included. Furthermore, these results are in excellent agreement with recent experimental results on the onset of absorption. Besides the excitons observed thus far, the calculations indicate the existence of a Rydberg-like series of exciton excited states, which is however modified from the classical Wannier exciton model by the anisotropies of the material and the more complex mixing of Bloch states in the excitons resulting from the Bethe-Salpeter equation. The exciton fine structure and the exciton wave functions are visualized and analyzed in various ways.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Screening of Coulombic Interactions To Achieve a Higher Power Factor in Conjugated Polymers

Thermoelectric properties of conducting polymers typically suffer from molecular chain disordering, as charge transport is predominantly controlled by morphology. This is especially more problematic when counterions are introduced to tune the carrier concentration for optimal thermoelectric performance, which disturbs the morphology further. Here, in this work, we introduce a new avenue for enhancing thermoelectric properties without needing to regulate the morphology, namely, by controlling the coulombic interaction between polarons and counterions. We perform in situ de-doping thermoelectric experiments over 3 orders of magnitude change in electrical conductivity of three distinct thermoelectric polymers, namely, poly(3-hexylthiophene-2,5-diyl) (P3HT), poly[2,5-bis(3-dodecylthiophen-2-yl)thieno[3,2-b]thiophene] (PBTTT-C 12 ), and poly[2,5-(2-octyldodecyl)-3,6-diketopyrrolopyrrole-alt-5,5-(2,5-di(thien-2-yl)thieno[3,2-b]thiophene)] (OD-PDPP2T-TT) conjugated polymers, followed by grazing-incidence wide-angle X-ray scattering (GIWAXS) to study their respective morphologies. We demonstrate a 9-fold enhancement in the thermoelectric power factor in OD-PDPP2T-TT compared to PBTTT-C 12 and link it to the coulombic screening of charge carriers, including in the optimally doped regime. We support this hypothesis using Boltzmann transport equations and show that, in both P3HT and PBTTT-C 12 , as the polymer is doped, impurity scattering remains the dominant scattering mechanism, while in OD-PDPP2T-TT, the scattering mechanism changes from impurity to acoustic phonon limited, resulting in more effective screening of ionized counterions. Our results provide an additional knob to enhance the fundamental understanding of thermoelectric physics of conducting polymers and provide a pathway to achieve higher performance in the field of organic thermoelectrics.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Real-space representation of the quasiparticle self-consistent GW self-energy and its application to defect calculations

The quasiparticle self-consistent (QS) GW (G for Green's function, W for screened Coulomb interaction) approach incorporates the corrections of the quasiparticle energies from their Kohn-Sham density functional theory (DFT) eigenvalues by means of an energy-independent and Hermitian self-energy matrix usually given in the basis set of the DFT eigenstates. By expanding these into an atom-centered basis set (specifically here the linearized muffin-tin orbitals) a real space representation of the self-energy corrections becomes possible. In this work, We show that this representation is relatively short-ranged. This offers opportunities to construct the self-energy of a complex system from parts of the system by a cut-and-paste method. Specifically for a point defect, represented in a large supercell, the self-energy can be constructed from those of the host and a smaller defect-containing cell. The self-energy of the periodic host can be constructed simply from a GW calculation for the primitive cell. We show for the case of the As Ga in GaAs that the defect part can already be well represented by a minimal eight-atom cell and allows us to construct the self-energy for a 64-atom cell in good agreement with direct QSGW calculations for the large cell. Using this approach to an even larger 216-atom cell shows the defect band approaches an isolated defect level. The calculations also allow us to identify a second defect band which appears as a resonance near the conduction band minimum. The results on the extracted defect levels agree well with Green's function calculations for an isolated defect and with experimental data.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Dynamic response functions of two-dimensional Dirac fermions with screened Coulomb and short-range interactions

Here, we consider a screened Coulomb interaction between electrons in graphene and determine their dynamic response functions, such as a longitudinal and a transverse electric conductivity and a polarization function and compare them to the corresponding quantities in the short-range interaction model. The calculations are performed to all orders for short-range interaction by taking into account the self-energy renormalization of the electron velocity and using a ladder approximation to account for the vertex corrections, ensuring that the Ward identity (charge conservation law) is satisfied. Our findings predict a resonant response of interacting electron-hole pairs at a particular frequency below the threshold $\textit{qv = ω}$ and further predict an instability for sufficiently strong interactions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1993

Plasma-screening effects on the electron-impact excitation of hydrogenic ions in dense plasmas

Plasma-screening effects are investigated on electron-impact excitation of hydrogenic ions in dense plasmas. Scaled cross sections Z(exp 4) sigma for 1s yields 2s and 1s yields 2p are obtained for a Debye-Hueckel model of the screened Coulomb interaction. Ground and excited bound wave functions are modified in the screened Coulomb potential (Debye-Hueckel model) using the Ritz variation method. The resulting atomic wave functions and their eigenenergies agree well with the numerical and high-order perturbation theory calculations for the interesting domain of the Debye length not less than 10. The Born approximation is used to describe the continuum states of the projectile electron. Plasma screening effects on the atomic electrons cannot be neglected in the high-density cases. Including these effects, the cross sections are appreciably increased for 1s yields 2s transitions and decreased for 1s yields 2p transitions.

Jung, Young-Dae↗
DOE OSTI2024

Liquid state theory of the structure of model polymerized ionic liquids

We employ polymer integral equation theory to study a simplified model of semiflexible polymerized ionic liquids (PolyILs) that interact via hard core repulsions and short range screened Coulomb interactions. The multi-scale structure in real and Fourier space of PolyILs (ions chosen to mimic Li, Na, K, Br, PF 6 , and TFSI) are determined as a function of melt density, Coulomb interaction strength, and ion size. Comparisons with a homopolymer melt, a neutral polymer–solvent-like athermal mixture, and an atomic ionic liquid are carried out to elucidate the distinct manner that ions mediate changes of polymer packing, the role of excluded volume effects, and the influence of chain connectivity, respectively. The effect of Coulomb strength depends in a rich manner on ion size and density, reflecting the interplay of steric packing, ion adsorption, and charge layering. Ion-mediated bridging of monomers is found, which intensifies for larger ions. Intermediate range charge layering correlations are characterized by a many-body screening length that grows with PolyIL density, cooling, and Coulomb strength, in disagreement with Debye–Hückel theory, but in accord with experiments. Qualitative differences in the collective structure, including an ion-size-dependent bifurcation of the polymer structure factor peak and pair correlation function, are predicted. The monomer cage order parameter increases significantly, but its collective ion counterpart decreases, as ions become smaller. Such behaviors allow one to categorize PolyILs into two broad classes of small and large ions. Furthermore, dynamical implications of the predicted structural results are qualitatively discussed.

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

Robustness of electronic screening effects in electron spectroscopies: Example of V 2 O 5

In bulk and low-dimensional extended systems, the screening of excitations by the electron cloud is a key feature governing spectroscopic properties. Widely used computational approaches, especially in the framework of many-body perturbation theory, such as the GW approximation and the resulting approximate Bethe-Salpeter equation, are explicitly formulated in terms of the screened Coulomb interaction. In the present work, we explore the effect of screening in absorption and electron energy loss spectroscopy, concentrating on the effect of local distortions on the screening and elucidating the resulting changes in the various spectra. Using the layered bulk oxide V 2 O 5 as prototype material, we show in which way local distortions affect the screening, and in which way changes in the screening impact electron energy loss and absorption spectra including excitons. We highlight cancellations that make many-body effects in the spectra very robust with respect to structural modifications, while the band structure undergoes significant changes and the nature of the excitations may also be affected. Here, this yields insight concerning the structure-properties relations that are crucial for the use of V 2 O 5 as energy storage material, and more generally, that may be used to optimize the analysis and the calculation of electronic spectra in complex materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

GPU-Accelerated Solution of the Bethe–Salpeter Equation for Large and Heterogeneous Systems

We present a massively parallel GPU-accelerated implementation of the Bethe–Salpeter equation (BSE) for the calculation of the vertical excitation energies (VEEs) and optical absorption spectra of condensed and molecular systems, starting from single-particle eigenvalues and eigenvectors obtained with density functional theory. The algorithms adopted here circumvent the slowly converging sums over empty and occupied states and the inversion of large dielectric matrices through a density matrix perturbation theory approach and a low-rank decomposition of the screened Coulomb interaction, respectively. Further computational savings are achieved by exploiting the nearsightedness of the density matrix of semiconductors and insulators to reduce the number of screened Coulomb integrals. We scale our calculations to thousands of GPUs with a hierarchical loop and data distribution strategy. The efficacy of our method is demonstrated by computing the VEEs of several spin defects in wide-band-gap materials, showing that supercells with up to 1000 atoms are necessary to obtain converged results. We discuss the validity of the common approximation that solves the BSE with truncated sums over empty and occupied states. In conclusion, we then apply our GW-BSE implementation to a diamond lattice with 1727 atoms to study the symmetry breaking of triplet states caused by the interaction of a point defect with an extended line defect.

Absorption spectra↗
DOE OSTI2024

Comparison of interband related optical transitions and excitons in ZnGeN 2 and GaN

The optical dielectric function of ZnGeN 2 is calculated from the interband transitions using the energy bands calculated in the quasiparticle self-consistent (QS)GW^ method using two different levels of approximation: the independent particle approximation (IPA) and the Bethe-Salpeter Equation (BSE) approach. The first allows us to relate peaks in ε 2 to specific bands and k-points but does not include electron-hole interaction effects. The second includes electron-hole interaction or excitonic effects. The corresponding changes in the shape of ε 2 (ω) are found to be similar to those in GaN. The screened Coulomb interaction W^ is here calculated including electron-hole interactions in the polarization function and gives a band structure already going beyond the random phase approximation. The static dielectric constants including only electronic screening, commonly referred to as ε∞ , were calculated separately by extrapolating the wave vector dependent macroscopic dielectric function in the long-wavelength limit. Below the quasiparticle gap, we find three bound excitons optically active for different polarization. The convergence of these bound excitons with respect to the density of the k-mesh used in the BSE is studied and found to require a fine mesh. It is also found that these bound excitons originate from only the lowest conduction band and the top three valence bands. To incorporate the lattice screening, we include a scaling factor (ε∞/ε 0 ) 2 , which allows us to obtain exciton binding energies of the correct order of magnitude similar to those in GaN. Here, the excitons are related to each of the three fold split valence bands and the splittings of the latter are also studied as function of strain. Finally, a relation between the anisotropic effective masses and the valence band splitting is pointed out and explained.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Superconductivity in a doped valley coherent insulator in magic angle graphene: Goldstone-mediated pairing and Kohn-Luttinger mechanism

We consider magic angle graphene in the doping regime around charge neutrality and study the connection between a recently proposed inter-valley coherent insulator at zero doping and the neighboring superconducting domes. The breaking of the valley U(1) symmetry leads to massless Goldstone modes, which couple to the doped charge carriers. We derive the effective interaction between these Goldstone modes and the conduction electrons and study its role in mediating superconductivity. Combining it with the screened Coulomb potential, we find weak-coupling superconducting instabilities in the two-component p- and d-wave channels. The competition between the two channels is set by the distance between the bilayer graphene device and the metallic gates. We find that the p-wave instability originates from the attraction mediated by the Goldstone modes, while the d-wave pairing is caused purely by the screened Coulomb interaction, similarly to the Kohn-Luttinger mechanism.

Kozii, Vladyslav↗
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 OSTI2021

Optical response and band structure of LiCoO 2 including electron-hole interaction effects

The optical response functions and band structures of LiCoO 2 are studied at different levels of approximation, from density functional theory (DFT) in the generalized gradient approximation (GGA) to quasiparticle self-consistent QSGW (with G for Green's function and W for screened Coulomb interaction) without and with ladder diagrams ($QSG\hat{W}$) and the Bethe Salpeter Equation (BSE) approach. The QSGW method is found to strongly overestimate the band gap and electron-hole or excitonic effects are found to be important. They lower the quasiparticle gap by only about 11~\% but the lowest energy peaks in absorption are found to be excitonic in nature. The contributions from different band to band transitions and the relation of excitons to band-to-band transitions are analyzed. Finally, the excitons are found to be strongly localized. Finally, a comparison to experimental data is presented.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Ten-valley excitonic complexes in charge-tunable monolayer WSe 2

Excitons dominate the optical response of two-dimensional (2D) semiconductors. Strong interactions produce peculiar excitonic complexes, which provide a testing ground for exciton and quantum many-body theories. Here, we report a hitherto unobserved many-body exciton that emerges upon filling both the K and Q valleys of WSe 2 . We optically probe the exciton landscape using charge-tunable devices with unusually thin dielectrics that facilitate doping up to several 10 13 cm −2 . We observe the emergence of the thermodynamically stable complex when 10 valleys are electrostatically filled. We gain insight into its physics using magneto-optical measurements. Our results are well-described by a model where the number of distinguishable Fermi seas interacting with the photoexcited electron-hole pair defines the complex’s behavior. In addition to expanding the repertoire of excitons in 2D semiconductors, this complex could probe the limit of exciton models and answer open questions about screened Coulomb interactions in 2D semiconductors.

M exciton↗
DOE OSTI2022

Full radiator-perturber interaction in computer simulations of hydrogenic spectral line broadening by plasmas

Spectral line broadening by plasmas can be computed by solving the equation of motion for the dipole of the radiating system perturbed by a fluctuating potential obtained from computer simulations. Such calculations have relied on the multipole expansion of the radiator-plasma interaction often keeping only the dipole term. With increasing density, however, higher multipoles as well as plasma perturbers overlapping the bound electron wave functions are expected to become important. For hydrogenic systems, the atomic matrix elements of the full Coulomb and screened Coulomb interactions are given by analytical formulas. Using these results, a computer simulation approach that accounts for the full radiator-plasma interaction is developed. One benefit is the removal of inherent strong collision divergences in the multipole expansion approximation. Furthermore, it yields the plasma polarization shift produced by perturbers penetrating the wave function of the radiator bound electrons. Here, the model was applied to hydrogenic argon Ly-α , Ly- β , and Ly-$\gamma$ spectral lines in a dense argon plasma at free electron densities of 10 24 or 10 25 cm ⁻3 and temperature of $800$ $eV$ relevant to plasma diagnostic techniques for inertial confinement fusion implosions.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Quantum Embedding Theory for Strongly Correlated States in Materials

Quantum embedding theories are promising approaches to investigate strongly correlated electronic states of active regions of large-scale molecular or condensed systems. Notable examples are spin defects in semiconductors and insulators. We present a detailed derivation of a quantum embedding theory recently introduced, which is based on the definition of effective Hamiltonians. The effect of the environment on a chosen active space is accounted for through screened Coulomb interactions evaluated using density functional theory. Importantly, the random phase approximation is not required, and the evaluation of virtual electronic orbitals is circumvented with algorithms previously developed in the context of calculations based on many-body perturbation theory. In addition, we generalize the quantum embedding theory to active spaces composed of orbitals that are not eigenstates of Kohn–Sham Hamiltonians. Finally, we report results for spin defects in semiconductors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗