Search NASA⌕ Search

SEARCH · Search NASA

Results for “COULOMB POTENTIAL”

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 OSTI2026

On the Trotter Error in Many-body Quantum Dynamics with Coulomb Potentials

Efficient simulation of many-body quantum systems is central to advances in physics, chemistry, and quantum computing, with a key question being whether the simulation cost scales polynomially with the system size. Here, in this work, we analyze many-body quantum systems with Coulomb interactions, which are fundamental to electronic and molecular systems. We prove that Trotterization for such unbounded Hamiltonians achieves a 1/4-order convergence rate, with explicit polynomial dependence on the number of particles. The result holds for all initial wavefunctions in the domain of the Hamiltonian, and the 1/4-order convergence rate is optimal, as previous work has numerically demonstrated that it can be saturated by a specific initial ground state. The main challenges arise from the many-body structure and the singular nature of the Coulomb potential. Our proof strategy differs from prior state-of-the-art Trotter analyses, addressing both difficulties in a unified framework. Our analysis treats the Coulomb potential as an unbounded operator without modification or regularization, and does not rely on spatial discretization, making it compatible with both first- and second-quantized circuit constructions.

Fang, Di [Duke Univ., Durham, NC (United States)]↗
DOE OSTI2025

Proton radiation effects in indium oxide using cascade molecular dynamics simulations

Metal oxide (MO) semiconductors, characterized by their wide band gaps and notable charge transport properties, are promising candidates for electronic applications in extreme environments, including near-Earth space. However, atomistic simulations of radiation–matter interactions in MOs remain challenging due to the limitations of existing interatomic potentials, which often fail to capture both the short-range repulsive forces essential for radiation damage modeling and the long-range electrostatic effects governing defect evolution. In this work, we develop a customized interatomic potential tailored for radiation damage simulations in indium oxide (In 2 O 3 ) as a model system, a representative MO material. Our potential integrates the Ziegler-Biersack-Littmark potential to accurately describe short-range interactions with Buckingham and Coulombic potentials to account for long-range forces. We perform molecular dynamics simulations of low-energy proton irradiation using this custom potential. We employ the primary knock-on atom (PKA) cascade method to study atomic displacements and primary defect formation. Simulations were conducted for 1 keV proton irradiation in four randomly chosen directions, and PKA-driven defect analyses at 5, 10, and 15 keV to examine the effects of direction and energy level on damage generation. Our results provide insight into the impact of irradiation direction and energy level on the cascade evolution and defect formation mechanisms.

Atomistic simulations↗
DOE OSTI2026

Molecular dynamics simulation of hydrodynamic transport coefficients in plasmas

Molecular dynamics simulations are used to calculate transport coefficients in a two-component plasma interacting through a repulsive Coulomb potential. The thermal conductivity, electrical conductivity, electrothermal coefficient, thermoelectric coefficient, and shear viscosity are computed using the Green–Kubo formalism over a broad range of Coulomb coupling strength, 0.01 ≤ Γ ≤ 140. Emphasis is placed on testing standard results of the Chapman–Enskog solution in the weakly coupled regime (Γ ≪ 1) using these first-principles simulations. As expected, the results show good agreement for Γ ≲ 0.1. However, this agreement is only possible if careful attention is paid to the definitions of linear constitutive relations in each of the theoretical models, a point that is often overlooked. For example, the standard Green–Kubo expression for thermal conductivity is a linear combination of thermal conductivity, electrothermal, and thermoelectric coefficients computed in the Chapman–Enskog formalism. Meaningful results for electrical conductivity are obtained over the full range of coupling strengths explored, but it is shown that potential and virial components of the other transport coefficients diverge in the strongly coupled regime (Γ ≫ 1). In this regime, only the kinetic components of the transport coefficients are meaningful for a classical plasma.

Electrical conductivity↗
DOE OSTI2025

Machine learning inversion from small-angle scattering for charged polymers

We develop Monte Carlo simulations for uniformly charged polymers and a machine learning algorithm to interpret the intra-polymer structure factor of the charged polymer system, which can be obtained from small-angle scattering experiments. The polymer is modeled as a chain of fixed-length bonds, where the connected bonds are subject to bending energy, and there is also a screened Coulomb potential for charge interaction between all joints. The bending energy is determined by the intrinsic bending stiffness, and the charge interaction depends on the interaction strength and screening length. All three contribute to the stiffness of the polymer chain and lead to longer and larger polymer conformations. The screening length also introduces a second length scale for the polymer besides the bending persistence length. To obtain the inverse mapping from the structure factor to these polymer conformation and energy-related parameters, we generate a large data set of structure factors by running simulations for a wide range of polymer energy parameters. We use principal component analysis to investigate the intra-polymer structure factors and determine the feasibility of the inversion using the nearest neighbor distance. We employ Gaussian process regression to achieve the inverse mapping and extract the characteristic parameters of polymers from the structure factor with low relative error.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Theory of topological exciton insulators and condensates in flat Chern bands

Excitons are the neutral quasiparticles that form when Coulomb interactions create bound states between electrons and holes. Due to their bosonic nature, excitons are expected to condense and exhibit superfluidity at sufficiently low temperatures. In interacting Chern insulators, excitons may inherit the nontrivial topology and quantum geometry from the underlying electron wavefunctions. We theoretically investigate the excitonic bound states and superfluidity in flat-band insulators pumped with light. We find that the exciton wavefunctions exhibit vortex structures in momentum space, with the total vorticity being equal to the difference of Chern numbers between the conduction and valence bands. Moreover, both the exciton binding energy and the exciton superfluid density are proportional to the Brillouin-zone average of the quantum metric and the Coulomb potential energy per unit cell. Spontaneous emission of circularly polarized light from radiative decay is a detectable signature of the exciton vorticity. We propose that the vorticity can also be experimentally measured via the nonlinear anomalous Hall effect, whereas the exciton superfluidity can be detected by voltage-drop quantization through a combination of quantum geometry and Aharonov–Casher effect. Topological excitons and their superfluid phase could be realized in flat bands of twisted Van der Waals heterostructures.

Science & Technology - Other Topics↗
DOE OSTI2026

Nonadiabatic Strong-Field Photoionization Revisited

We measure strong field ionization of cesium atoms, observing a robust feature near 2⁢𝑈p in the photoelectron spectrum, which we call the intermediate energy structure (IES). Using a Coulomb-corrected strong-field approximation, we show it arises from electrons born with large inward velocities that rapidly undergo forward scattering off the Coulomb potential. The IES is similar to the previously identified low energy structure in that they are both due to forward scattering. However, it is different in that IES requires both nonadiabatic ionization conditions (as defined by the Keldysh parameter) and a more weakly bound initial state. Furthermore, our joint experimental and theoretical study supports the need for an additional parameter in characterization of strong field ionization, which depends both on nonadiabaticity of the process and the principal quantum number of the target. Together with Reiss and Keldysh parameters, this dimensionless parameter expands the description of strong field ionization across a broad range of atomic targets and incorporates ionization from excited states.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2026

Quantifying Outer- and Inner-Coordination Sphere Effects Using Uranium Redox Chemistry in Molten Salt Solutions

Defining the relative influence of intramolecular and intermolecular forces is a fundamental problem in chemistry that is difficult to quantify. To address this challenge, we developed a method to evaluate the relative impact of direct chemical bonding in the inner-coordination sphere vs effects from cations in the outer-coordination sphere by comparative analysis of uranium redox reactivity in various molten salts. We observed that outer-coordination sphere cations (M 1+ ) and inner-coordination sphere anions (X 1– ) both affected uranium redox reactivity, with more polarizing M 1+ and larger X 1– favoring uranium in low oxidation states. Changing M 1+ (Li, Na, K) shifted the U IV + e 1– ⇌ U III (U IV/III ) and U III + 3e 1– → U 0 metal potentials by +330 and +240 mV, respectively. Changing X 1– (Cl, Br, I) caused larger shifts of +440 mV for the U IV/III redox potential and +1060 mV for the U 0 metal deposition potential. Using Coulomb’s Law, we correlated these potentials with electrostatic interactions between UIII and the molten salt. This model provided a facile way of predicting redox chemistry within molten salts.

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

Simulating Crystallization in a Colloidal System Using State Predictive Information Bottleneck Based Enhanced Sampling

Here, we investigate crystal nucleation in supersaturated colloid suspensions using enhanced molecular dynamics simulations augmented with machine learning techniques. The simulations reveal that crystallization in the model colloidal system studied here, with particles interacting through a repulsive screened Coulomb Yukawa potential, proceeds from vapor to dense liquid droplet to crystalline phases across multiple high barriers. Employing a one-dimensional reaction coordinate derived from the State Predictive Information Bottleneck framework, our simulations capture back-and-forth phase transitions across multiple barriers effectively in biased metadynamics simulations. We obtain relative free energy differences between different phases and also quantify the roles of different molecular level features in driving the phase changes.

Chemistry↗
DOE OSTI2024

Fast ground-state-to-ground-state separation of small ion crystals

Rapid separation of linear crystals of trapped ions into different subsets is critical for realizing trapped ion quantum computing architectures where ions are rearranged in trap arrays to achieve all-to-all connectivity between qubits. Here we introduce a general theoretical framework that can be used to describe the separation of same-species and mixed-species crystals into smaller subsets. The framework relies on an efficient description of the evolution of Gaussian motional states under quadratic Hamiltonians that only requires a special solution of the classical equations of motion of the ions to describe their quantum evolution under the influence of a time-dependent applied potential and the ions' mutual Coulomb repulsion. We provide time-dependent applied potentials suitable for separation of a mixed-species three-ion crystal on timescales similar to that of free expansion driven by Coulomb repulsion, with all modes along the crystal axis starting and ending close to their ground states. Three separately confined mixed-species ions can be combined into a crystal held in a single well without energy gain by time-reversal of this separation process.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Mixing effects on spectroscopy and partonic observables of heavy mesons with logarithmic confining potential in a light-front quark model

Using the variational principle, we systematically investigate the mass spectra and wave functions of both 1⁢𝑆 and 2⁢𝑆 state heavy pseudoscalar (𝑃) and vector (𝑉) mesons within the light-front quark model. This approach incorporates a Coulomb plus logarithmic confinement potential to accurately describe the constituent quark and antiquark dynamics. Additionally, spin hyperfine interactions are introduced perturbatively to compute the masses of pseudoscalar and vector mesons. The present analyses of the 1⁢𝑆 and 2⁢𝑆 states require the consideration of mixing between them to account for empirical constraints. These constraints include the mass gap Δ⁢𝑀 𝑃 >Δ⁢𝑀 𝑉 , where Δ⁢𝑀 𝑃⁡(𝑉) =𝑀$^{2⁢𝑆}_{𝑃⁡(𝑉)}$−𝑀$^{1⁢𝑆}_{𝑃⁡(𝑉)}$ and the hierarchy of the decay constants 𝑓 1⁢𝑆 >𝑓 2⁢𝑆 . We find the optimal value of the mixing angle to be 𝜃 =1⁢8°, significantly enhancing the consistency between our spectroscopic predictions and the experimental data compiled by the Particle Data Group. Furthermore, based on the predicted mass, the newly observed resonance 𝐵 𝐽⁡ (5840) could be assigned as a 2 1⁢ 𝑆 0 state in the 𝐵 meson family. The study also reports various pertinent observables, including twist-two distribution amplitudes, electromagnetic form factors, charge radii, 𝜉 moments, and transition form factors that are found to be consistent with both available lattice simulations and experimental data. In addition, our predicted branching ratios for the channels of 𝐵 + →𝜏 + ⁢𝜈 𝜏 as well as rare decays of 𝐵 0 and 𝐵$^0_𝑠$ appear in accordance with experimental data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Ω 3⁢𝑐 ⁢𝑁⁢𝑁 and Ω 3⁢𝑐 ⁢Ω 3⁢𝑐 ⁢𝑁 systems with HAL QCD potentials

Here, this study employs the Faddeev formalism in configuration space to investigate the Ω 3⁢𝑐 ⁢𝑁⁢𝑁 cluster containing a triply charmed Omega baryon (Ω 3⁢𝑐 ). Using the recently reported HAL QCD 𝑆-wave Ω 3⁢𝑐 ⁢𝑁 potentials in the 3 𝑆 1 and 5 𝑆 2 channels, together with the MT-I–III nucleon-nucleon potential and neglecting the Coulomb force, we find no bound state for the Ω 3⁢𝑐 ⁢𝑛⁢𝑝 system. We predict near-threshold resonances in the 𝐽 𝜋 =5/2 + (maximal total spin) and 𝐽 𝜋 =1/2 + (minimal total spin) states, with resonance energies of 1.1 MeV below and 0.0 MeV at the three-body breakup threshold, respectively, at Euclidean time 𝑡/𝑎 =16. A similar analysis of the Ω 3⁢𝑐⁢ Ω 3⁢𝑐⁢ 𝑁 system likewise reveals no bound states, though a possible resonance is indicated.

few-body systems↗
DOE OSTI2024

Particle-Hole Asymmetric Ferromagnetism and Spin Textures in the Triangular Hubbard-Hofstadter Model

In a lattice model subject to a perpendicular magnetic field, when the lattice constant is comparable to the magnetic length, one enters the “Hofstadter regime,” where continuum Landau levels become fractal magnetic Bloch bands. Strong mixing between bands alters the nature of the resulting quantum phases compared to the continuum limit; lattice potential, magnetic field, and Coulomb interaction must be treated on equal footing. Using determinant quantum Monte Carlo and density matrix renormalization group techniques, we study this regime numerically in the context of the Hubbard-Hofstadter model on a triangular lattice. In the field-filling phase diagram, we find a broad wedge-shaped region of ferromagnetic ground states for filling factor ν ≤ 1 , bounded below by filling factor ν = 1 and bounded above by half filling the lowest Hofstadter subband. We observe signatures of SU(2) quantum Hall ferromagnetism at filling factors ν = 1 and ν = 3 . The phases near ν = 1 are particle-hole asymmetric, and we observe a rapid decrease in ground-state spin polarization consistent with the formation of skyrmions only on the electron doped side. At large fields, above the ferromagnetic wedge, we observe a low-spin metallic region with spin correlations peaked at small momenta. We argue that the phenomenology of this region likely results from exchange interaction mixing fractal Hofstadter subbands. The phase diagram derived beyond the continuum limit points to a rich landscape to explore interaction effects in magnetic Bloch bands. Published by the American Physical Society 2024

Ding, Jixun K.↗
DOE OSTI2025

Automated workflow for non-empirical Wannier-localized optimal tuning of range-separated hybrid functionals

Here, we introduce an automated workflow for generating non-empirical Wannier-localized optimally-tuned screened range-separated hybrid (WOT-SRSH) functionals. WOT-SRSH functionals have been shown to yield highly accurate fundamental band gaps, band structures, and optical spectra for bulk and 2D semiconductors and insulators. Our workflow automatically and efficiently determines the WOT-SRSH functional parameters for a given crystal structure and composition, approximately enforcing the correct screened long-range Coulomb interaction and an ionization potential ansatz. In contrast to previous manual tuning approaches, our tuning procedure relies on a new search algorithm that only requires a few hybrid functional calculations with minimal user input. We demonstrate our workflow on 23 previously studied semiconductors and insulators, reporting the same high level of accuracy. By automating the tuning process and improving its computational efficiency, the approach outlined here enables applications of the WOT-SRSH functional to compute spectroscopic and optoelectronic properties for a wide range of materials.

Gant, Stephen E. [University of California, Berkel↗
DOE OSTI2025

Visualizing the Three-Dimensional Arrangement of Hydrogen Atoms in Organic Molecules by Coulomb Explosion Imaging

Structure-sensitive methods based on femtosecond light or electron pulses are now making it possible to measure how molecular structures change during light-induced processes. Despite significant progress, high-fidelity imaging of nuclear positions remains a challenge even for relatively small molecular systems and, notably, regarding the positions of hydrogen atoms. As demonstrated in recent work, X-ray-induced Coulomb explosion imaging (CEI) may overcome this obstacle, as its sensitivity does not depend on the mass of the imaged atoms. The photoinduced ring opening of the heterocyclic molecule 2(5 H )-thiophenone has attracted recent interest. Here, in this work, we show that CEI offers a powerful route to imaging the peripheral H atoms in this molecule and thus, more generally, to tracking detailed nuclear motions (e.g., isomerizations) in organic molecules on ultrafast time scales. Specifically, we record momentum-space Coulomb explosion images that report on the three-dimensional positioning of all nuclei within the molecule, for instance, distinguishing H atoms in C–H bonds that lie within or are directed out of the plane defined by the heavy atoms. The prospect of imaging peripheral H atoms to probe photochemical dynamics is explored by coupling ab initio molecular dynamics with classical Coulomb explosion simulations, thereby differentiating potential photoproduct isomers, including those whose structures primarily differ in the position of the hydrogens.

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

Alloying Interlayers for Anode‐Less Lithium‐Metal Polymer Batteries

Anode‐less lithium‐metal batteries potentially offer further increased energy densities. However, the Coulombic efficiency (CE) of lithium plating and stripping, as a classical measure of Li inventory reversibility, is commonly still insufficient for achieving long‐lasting rechargeable batteries. Herein, the potential benefits of employing thin (20 nm) metal interlayers of Ag, Pt, and Au on Cu to alloy with Li in Cu||Li half‐cells and induce homogeneous Li plating with poly(ethylene oxide)‐based electrolytes are investigated. Interestingly, not all alloying interlayers enable a higher CE compared to neat Cu foil with 84%—specifically Ag@Cu with only 81%, while the best performing one, Au@Cu, provides a substantially increased CE of 91%. While generally the formation of “dead lithium” is found to be the major source of CE, this appears to be less pronounced in the case of Au@Cu, indeed. Further improvement can be achieved by carefully adjusting the cell voltage to a region in which the continuous de‐/alloying is suppressed, yielding a further enhanced CE of 94%, thus highlighting the need for a comprehensive approach to design suitable electrode chemistries and designs beyond a “simple” material improvement.

anode-less battery↗
DOE OSTI2025

The Chemistry, Formation and Passivation of Solid-Electrolyte Interphase (SEI) from FSI- Breakdown at Li-metal Potential

LiFSI-based liquid electrolytes are promising for realizing high Coulombic efficiency (CE) and long cycle life of next-generation high-energy Li metal battery. Nevertheless, the fundamentals of both the chemistry and formation of solid-electrolyte interphase (SEI) in these electrolytes have remained elusive. Herein, we extensively combine electrochemistry and X-ray photoelectron spectroscopy (XPS) to illustrate the reaction pathways of electrolyte decomposition, especially FSI- breakdown, as well as the dissolution trend of various inorganic SEI components. The SEI-forming reactions at Cu/electrolyte interface manifest a potential (E) dependence, driven by direct electron-transfer (ET) and Li underpotential-deposition (Li UPD) at high and low E regions, respectively. By closely correlating the XPS data obtained with and without solvent washing, we demonstrate that not all the products derived from SEI reactions are incorporated into forming a passivating surface layer, with a notable portion dissolving into liquid electrolyte instead. Importantly, high-CE electrolyte dictated by the chemical identity of solvent, is found to exhibit minimized interfacial reactivity of SEI formation enabling more effective interphasial passivation, via preferential integration of passivating ingredients such as LiF. Overall, this work formulates a systematic understating of bridging the SEI chemistry, formation and passivation while identifying key characteristics of high-performance electrolyte for guiding future design.

Bao, Zhenan [SLAC National Accelerator Laboratory ↗
DOE OSTI2025

16 O + 16 O fusion at sub-coulomb energies

A simple optical model (OPM) method using non-monotonic (NM) potentials characterized by a repulsive core from the microscopic theory of the Pauli-led energy-density functional (EDF) has been developed to investigate the 16 O + 16 O fusion at sub-Coulomb energies relevant to the oxygen burning. Here, the study involves the analysis of the experimental fusion cross-section (FCS) data in energy range 6.92 ≤ E cm ≤ 13.83 MeV, which includes the Coulomb barrier region E cm = 10.0–11.92 MeV. Apart from an excellent description of the existing FCS data in the studied energy range, the associated hindrance, characterized by the S-factor reaching a maximum and then gradually falling off at lower energies, so far suggested empirically for the system, is reproduced down to 4 MeV for the first time in the simple OPM. An S-factor maximum of S 0 = 3.15 × 10 25 MeV.b at E 0 = 7.14 MeV is observed with T ≃ 2.6 GK, which conforms to the values reported for quiescent and explosive burning. Our reaction rate, deduced from the NM potential, compares well with the Caughlan and Fowler data. Dominant partial waves implicit in the observed maximum S-factor in the studied Gamow range are also explored. Our present findings, with the success of NM potentials, suggest that the nucleus–nucleus potential is non-monotonic.

16O+16O↗
DOE OSTI2024

Coulomb confinement in the Hamiltonian limit

The Gribov-Zwanziger scenario attributes the phenomenon of confinement to the instantaneous interaction term in the QCD Hamiltonian in the Coulomb gauge. For a static quark-antiquark pair, it leads to a potential energy that increases linearly with the distance between them. Lattice studies of the SU(2) Yang-Mills theory determined the corresponding (Coulomb) string tension for sources in the fundamental representation, 𝜎 𝐶 , to be about 3 times larger than the Wilson loop string tension, 𝜎 𝐹 . It is far above the Zwanziger variational bound, 𝜎 𝐶 ≥ 𝜎 𝐹 . We argue that the value often reported in the literature is artificially inflated. We examine the lattice definition of the instantaneous potential, find the source of the string tension’s enhancement, and perform its improved determination in SU(2) lattice gauge theory. We report our conservative estimate for the value of the Coulomb string tension as 𝜎 𝐶 /𝜎 𝐹 = 2.0 ± 0.4 and discuss its phenomenological implications.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗