Search NASA⌕ Search

SEARCH · Search NASA

Results for “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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Improving the efficiency of G 0 W 0 calculations with approximate spectral decompositions of dielectric matrices

Recently, it was shown that the calculation of quasiparticle energies using the G 0 W 0 approximation can be performed without computing explicitly any virtual electronic states, by expanding the Green function and screened Coulomb interaction in terms of the eigenstates of the static dielectric matrix. Avoiding the evaluation of virtual electronic states leads to improved efficiency and ease of convergence of G 0 W 0 calculations. Here, we propose a further improvement of the efficiency of these calculations, based on an approximation of density-density response functions of molecules and solids. The approximation relies on the calculation of a subset of eigenvectors of the dielectric matrix using the kinetic operator instead of the full Hamiltonian, and it does not lead to any substantial loss of accuracy for the quasiparticle energies. The computational savings introduced by this approximation depend on the system, and they become more substantial as the number of electrons increases.

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

Nanoshell quantum dots: Quantum confinement beyond the exciton Bohr radius

Nanoshell quantum dots (QDs) represent a novel class of colloidal semiconductor nanocrystals (NCs), which supports tunable optoelectronic properties over the extended range of particle sizes. Traditionally, the ability to control the bandgap of colloidal semiconductor NCs is limited to small-size nanostructures, where photoinduced charges are confined by Coulomb interactions. A notorious drawback of such a restricted size range concerns the fact that assemblies of smaller nanoparticles tend to exhibit a greater density of interfacial and surface defects. This presents a potential problem for device applications of semiconductor NCs where the charge transport across nanoparticle films is important, as in the case of solar cells, field-effect transistors, and photoelectrochemical devices. Furthermore, the morphology of nanoshell QDs addresses this issue by enabling the quantum-confinement in the shell layer, where two-dimensional excitons can exist, regardless of the total particle size. Such a geometry exhibits one of the lowest surface-to-volume ratios among existing QD architectures and, therefore, could potentially lead to improved charge-transport and multi-exciton characteristics. The expected benefits of the nanoshell architecture were recently demonstrated by a number of reports on the CdS bulk /CdSe nanoshell model system, showing an improved photoconductivity of solids and increased lifetime of multi-exciton populations. Along these lines, this perspective will summarize the recent work on CdS bulk /CdSe nanoshell colloids and discuss the possibility of employing other nanoshell semiconductor combinations in light-harvesting and lasing applications.

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

Extended Lagrangian Born–Oppenheimer molecular dynamics for orbital-free density-functional theory and polarizable charge equilibration models

We report extended Lagrangian Born–Oppenheimer molecular dynamics (XL-BOMD) is formulated for orbital-free Hohenberg–Kohn density-functional theory and for charge equilibration and polarizable force-field models that can be derived from the same orbital-free framework. The purpose is to introduce the most recent features of orbital-based XL-BOMD to molecular dynamics simulations based on charge equilibration and polarizable force-field models. These features include a metric tensor generalization of the extended harmonic potential, preconditioners, and the ability to use only a single Coulomb summation to determine the fully equilibrated charges and the interatomic forces in each time step for the shadow Born–Oppenheimer potential energy surface. The orbital-free formulation has a charge-dependent, short-range energy term that is separate from long-range Coulomb interactions. This enables local parameterizations of the short-range energy term, while the long-range electrostatic interactions can be treated separately. The theory is illustrated for molecular dynamics simulations of an atomistic system described by a charge equilibration model with periodic boundary conditions. The system of linear equations that determines the equilibrated charges and the forces is diagonal, and only a single Ewald summation is needed in each time step. The simulations exhibit the same features in accuracy, convergence, and stability as are expected from orbital-based XL-BOMD.

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

Understanding supercooling mechanism in sodium sulfate decahydrate phase-change material

Salt hydrate-based phase-change materials are considered promising for future heat storage applications in residential heating/cooling systems. Smooth phase transition from the liquid to solid phase and vice versa is essential for effective heat exchanger; however, supercooling in salt hydrates delays the onset of liquid–solid phase transition. We investigate the molecular level mechanism of supercooling in sodium sulfate decahydrate (SSD). SSD is a complex salt hydrate whose properties are governed by electrostatic forces that include pure Coulombic interactions as well as hydrogen bonds. Experimentally, we examine the importance of a nucleator in reducing supercooling temperatures. We investigated the effect of various mass concentrations of a borax nucleator on a decrease of supercooling temperatures. Molecular dynamics simulation techniques are used to obtain a basic understanding of supercooling in SSD. We observe that by introducing borax as a nucleator, there is a decrease in the supercooling temperature before nucleation. Our molecular dynamics simulations show that long-range electrostatics between sodium and sulfate ion pairs and that with polar water molecules is responsible for delayed nucleation in SSD that results in supercooling, and also, dynamics of charged molecules slows down. The lack of crystallization leads to amorphous structures in supercooled SSD.

36 MATERIALS SCIENCE↗
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 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 OSTI2024

Ultrafast laser triggered electron emission from ultrananocrystalline diamond pyramid tip cathode

Nitrogen-incorporated ultrananocrystalline diamond [(N)UNCD] pyramid tip cathode has been considered as a next-generation high peak current electron source for dielectric laser accelerators as well as other high peak current particle accelerator applications. In this work, we study non-linear photoemission from an (N)UNCD pyramid tip cathode using an ultrafast laser with the pulse length of 150 fs with the central wavelength of 800 nm in the peak intensity range of 10 9 –10 10 W/cm 2 . We demonstrated that as the incident laser intensity increases, the current emitted from the nano-tip first increases as a power function with an exponent of about 5 and then starts to roll over to an exponent of 3. This roll over is attributed to the Coulomb interaction between electrons emitted from the tip also known as the space charge. We also measured the photoemission electron energy spectra that show electrons with energies as high as ~10 eV. Based on the shape of the electron energy spectra, we conclude that the high-energy electrons are thermally emitted electrons due to ultrafast laser heating at the tip of the (N)UNCD pyramid tip cathode.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

GW with hybrid functionals for large molecular systems

A low-cost approach for stochastically sampling static exchange during time-dependent Hartree–Fock-type propagation is presented. This enables the use of an excellent hybrid density functional theory (DFT) starting point for stochastic GW quasiparticle energy calculations. Generalized Kohn–Sham molecular orbitals and energies, rather than those of a local-DFT calculation, are used for building the Green function and effective Coulomb interaction. The use of an optimally tuned hybrid diminishes the starting point dependency in one-shot stochastic GW, effectively avoiding the need for self-consistent GW iterations.

Chemistry↗
DOE OSTI2024

Lanthanide binding peptide surfactants at air–aqueous interfaces for interfacial separation of rare earth elements

Rare earth elements (REEs) are critical materials to modern technologies. They are obtained by selective separation from mining feedstocks consisting of mixtures of their trivalent cation. We are developing an all-aqueous, bioinspired, interfacial separation using peptides as amphiphilic molecular extractants. Lanthanide binding tags (LBTs) are amphiphilic peptide sequences based on the EF-hand metal binding loops of calcium-binding proteins which complex selectively REEs. We study LBTs optimized for coordination to Tb 3+ using luminescence spectroscopy, surface tensiometry, X-ray reflectivity, and X-ray fluorescence near total reflection, and find that these LBTs capture Tb 3+ in bulk and adsorb the complex to the interface. Molecular dynamics show that the binding pocket remains intact upon adsorption. We find that, if the net negative charge on the peptide results in a negatively charged complex, excess cations are recruited to the interface by nonselective Coulombic interactions that compromise selective REE capture. If, however, the net negative charge on the peptide is −3, resulting in a neutral complex, a 1:1 surface ratio of cation to peptide is achieved. Surface adsorption of the neutral peptide complexes from an equimolar mixture of Tb 3+ and La 3+ demonstrates a switchable platform dictated by bulk and interfacial effects. The adsorption layer becomes enriched in the favored Tb 3+ when the bulk peptide is saturated, but selective to La 3+ for undersaturation due to a higher surface activity of the La 3+ complex.

Ortuno Macias, Luis E. (ORCID:0000000284342192)↗
DOE OSTI2022

Sub Coulomb barrier d+ 208 Pb scattering in the time-dependent basis function approach

We employ the non-perturbative time-dependent basis function (tBF) approach to study the scattering of the deuteron on 208 Pb below the Coulomb barrier. We obtain the bound and discretized scattering states of the projectile, which form the basis representation of the tBF approach, by diagonalizing a realistic Hamiltonian in a large harmonic oscillator basis. We find that the higher-order inelastic scattering effects are noticeable for sub barrier scatterings with the tBF method. We have successfully reproduced experimental sub Coulomb barrier elastic cross section ratios with the tBF approach by considering only the electric dipole (E1) component of the Coulomb interaction between the projectile and the target during scatterings. Here, we find that the correction of the polarization potential to the Rutherford trajectory is dominant in reproducing the data at very low bombarding energies, whereas the role of internal transitions of the deuteron projectile induced by the E1 interaction during the scattering becomes increasingly significant at higher bombarding energies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Development of a new quantum trajectory molecular dynamics framework

An extension to the wave packet description of quantum plasmas is presented, where the wave packet can be elongated in arbitrary directions. A generalized Ewald summation is constructed for the wave packet models accounting for long-range Coulomb interactions and fermionic effects are approximated by purpose-built Pauli potentials, self-consistent with the wave packets used. We demonstrate its numerical implementation with good parallel support and close to linear scaling in particle number, used for comparisons with the more common wave packet employing isotropic states. Ground state and thermal properties are compared between the models with differences occurring primarily in the electronic subsystem. Especially, the electrical conductivity of dense hydrogen is investigated where a 15% increase in DC conductivity can be seen in our wave packet model compared with other models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Significance of the structural configuration of B2 disorder in Co and Ti based Heusler alloys

We investigate here structural (at local and global levels) and transport properties for 𝑋 2 ⁢MnAl (𝑋= Co and Ti). Additionally, the magnetic properties were also studied for Ti 2 ⁢MnAl. Our x-ray diffraction results show that both the compounds stabilize in B2 disordered phase with cubic structure of 𝑃⁢𝑚⁢$\overline{3}$⁢𝑚 space group. Further, the structural configuration of the above disordered phase for both the compounds was identified using combined studies of x-ray absorption spectroscopy and multiple scattering calculations at the transition metal 𝐾 edges. Upon such identification, in the case of Co 2 ⁢Mn 1−𝑦⁢ Cr 𝑦 ⁢Al (𝑦= 0, 0.05, 0.1, 0.2) with change in 𝑦, we are able to establish a better connection quantitatively between the inverse of Mn-Co bonds and peak in the temperature-dependent resistivity. This highlights the crucial importance of a detailed understanding of the nature of B2 disorder. In the case of 𝑦=0, in the temperature range of study, the resistivity is driven by the functional form associated with (a) three-dimensional enhanced electron-electron Coulomb interaction scattering mechanism and (b) an unconventional one-magnon process. For Ti 2 ⁢MnAl, the transport shows metallic glasslike behavior at high temperature, while at low temperature it follows both the Cote-Meisel's model and quantum correction model. In this compound, the magnetic studies suggest the formation of superparamagnetic clusters in the paramagnetic matrix at low temperatures. Our density functional theory results are in line with the transport and magnetic properties. In literature, the spin polarization percentage (𝑃) for 𝑋= Co in B2 disordered phase is 76%. However, the present results emphasize the fact that in B2 disordered phase, the value of 𝑃 can range from 90% to 71% depending on the structural configuration introduced by swapping of the atomic positions of Mn and Al. For the compound under study, the value of percentage spin polarization obtained ranges between 82% to 85%. In addition, we also identify the origin of the difference in the shape of the Mn 3⁢𝑑 density of states for both the alloys. In conclusion, our results for 𝑋= Co alloy highlights the importance of identifying the specific structural configuration associated with a particular disorder category especially in the estimation of 𝑇 𝑐 and 𝑃 and for 𝑋= Ti, the physical properties can be tuned by varying the position of 𝐸 𝐹 and thereby its utilization in device applications.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Comparative analysis of plasmon modes in layered Lindhard metals and strange metals

The enigmatic strange metal remains one of the central unsolved problems of 21st century science. Understanding this phase of matter requires knowledge of the momentum- and energy-resolved dynamic charge susceptibility 𝜒⁡(𝑞,𝜔), especially at finite momentum. Inelastic electron scattering (EELS), performed in either transmission or reflection geometry, is a powerful probe of 𝜒⁡(𝑞,𝜔). For the prototypical strange metal Bi 2 ⁢Sr 2 ⁢CaCu 2 ⁢O 8+𝑥 , transmission- and reflection EELS, and infrared (IR) spectroscopy agree at 𝑞∼0, all revealing a highly damped plasmon near 1 eV. At larger 𝑞, however, EELS results show unresolved discrepancies. Since IR data are highly reproducible, it is advantageous to use IR data to calculate what the expected EELS response should be at modest 𝑞. Building on prior momentum-resolved reflection geometry M-EELS work [Chen et al., Phys. Rev. B 109, 045108 (2024)], we extend this approach to transmission EELS for finite stacks of metallic layers, comparing a “textbook” Lindhard metal to a strange metal. In the Lindhard case, the low-𝑞 response is dominated by long-lived, standing wave plasmon modes arising from interlayer Coulomb coupling, with in-plane dispersions that resemble the well-known Fetter modes of layered metals. This behavior depends only on the geometry and the long-range nature of the Coulomb interaction and is largely insensitive to layer details. At larger 𝑞, the response reflects the microscopic properties of individual layers. For the strange metal, calculations based on IR data predict a highly damped plasmon with weak dispersion and no distinct surface mode. While our results match IR and M-EELS at low 𝑞, they do not reproduce any published EELS spectra at large 𝑞, highlighting unresolved discrepancies that demand further experimental investigation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

New scaling and nuclear structure aspects in heavy-ion fusion reactions

Three new behaviors have been found in comparisons of fusion cross sections for different collision systems. root (1) Replacing the energy E with a scaling one, E scal = (E-V g )/($\sqrt{2}$W g ), is successful for washing out the Coulomb interaction in the spectra of fusion cross sections, where V g and W g are barrier height and width of the single-Gaussian barrier distribution model. (2) In a representation of σE vs the scaling energy, E scal , all data sets display in parallel. Here, the ratio for sigma E from any two fusion systems over the whole range is a constant value. That behavior is also studied in another representation, in which the data sets display as parallel horizontal lines for any heavy-ion fusion system. (3) The constant ratio value is the ratio of parameter products, $R^2_gW_g$, of the two systems; where R g is the barrier radius obtained in the single-Gaussian barrier distribution model. Moreover, when comparing neighboring collision systems at the same E scal , the ratio of sigma is near a constant value within a few percent over the whole range. Thus a quantitative comparison for the fusion enhancement for neighboring systems is developed. The present finding could be beneficial for predicting unmeasured fusion cross sections.

Jiang, C. L. [Argonne National Laboratory (ANL), A↗
DOE OSTI2026

Interlayer Exciton Condensates between Second Landau Level Orbitals in Double Bilayer Graphene

Here, we present Coulomb-drag measurements on a heterostructure comprising two Bernal-stacked bilayer graphene (BLG) sheets separated by a 2.5 nm hexagonal boron nitride (hBN) spacer in the quantum Hall (QH) regime. Using top and bottom gate control, together with an interlayer bias, we independently tune the two BLG layers into either the lowest (N = 0) or second (N = 1) Landau level (LL) orbital and probe their interlayer QH states. When both layers occupy the N = 0 orbital, we observe both interlayer exciton condensates (ECs) at integer total filling and interlayer fractional QH states, echoing the results in double monolayer graphene. In contrast to previous studies, however, when both BLG layers occupy the N = 1 orbital, we also observe quantized drag signals, signifying an interlayer exciton condensate formed between the second LLs. By tuning the layer degree of freedom, we find that this N = 1 EC state arises only when the N = 1 wave function in each BLG is polarized toward the hBN interface to maximize the interlayer Coulomb interaction.

Electrical conductivity↗
DOE OSTI2020

Dual-cone variational calculation of the two-electron reduced density matrix

The computation of strongly correlated quantum systems is challenging because of its potentially exponential scaling in the number of electron configurations. Variational calculation of the two-electron reduced density matrix (2-RDM) without the many-electron wave function exploits the pairwise nature of the electronic Coulomb interaction to compute a lower bound on the ground-state energy with polynomial computational scaling. Recently, a dual-cone formulation of the variational 2-RDM calculation was shown to generate the ground-state energy, albeit not the 2-RDM, at a substantially reduced computational cost, especially for higher $\textit{N}$-representability conditions such as the T2 constraint. Here we generalize the dual-cone variational 2-RDM method to compute not only the ground-state energy but also the 2-RDM. The central result is that we can compute the 2-RDM from a generalization of the Hellmann-Feynman theorem. Specifically, we prove that in the Lagrangian formulation of the dual-cone optimization the 2-RDM is the Lagrange multiplier. Further, we apply the method to computing the energies and properties of strongly correlated electrons—including atomic charges, electron densities, dipole moments, and orbital occupations—in an illustrative hydrogen chain and the nitrogen-fixation catalyst FeMoco. The dual variational computation of the 2-RDM with T2 or higher N -representability conditions provides a polynomially scaling approach to strongly correlated molecules and materials with significant applications in atomic and molecular and condensed-matter chemistry and physics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗