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 OSTI2024

Quantum Drude oscillators coupled with Coulomb potential as an efficient model for bonded and non-covalent interactions in atomic dimers

The quantum Drude oscillator (QDO) model has been widely used as an efficient surrogate to describe the electric response properties of matter as well as long-range interactions in molecules and materials. Most commonly, QDOs are coupled within the dipole approximation so that the Hamiltonian can be exactly diagonalized, which forms the basis for the many-body dispersion method [Phys. Rev. Lett. 108, 236402 (2012)]. The dipole coupling is efficient and allows us to study non-covalent many-body effects in systems with thousands of atoms. However, there are two limitations: (i) the need to regularize the interaction at short distances with empirical damping functions and (ii) the lack of multipolar effects in the coupling potential. In this work, we convincingly address both limitations of the dipole-coupled QDO model by presenting a numerically exact solution of the Coulomb-coupled QDO model by means of quantum Monte Carlo methods. We calculate the potential-energy surfaces of homogeneous QDO dimers, analyzing their properties as a function of the three tunable parameters: frequency, reduced mass, and charge. We study the coupled-QDO model behavior at short distances and show how to parameterize this model to enable an effective description of chemical bonds, such as the covalent bond in the H2 molecule.

Chemistry↗
DOE OSTI2023

Charged-Particle Bound States in Periodic Boxes

We consider the binding energy of a two-body system with a repulsive Coulomb interaction in a finite periodic volume. We define the finite-volume Coulomb potential as the usual Coulomb potential, except that the distance is defined the shortest separation between the two bodies in the periodic volume. We investigate this problem in one and three-dimensional periodic boxes and derive the asymptotic behavior of the volume dependence for bound states with zero angular momentum in terms of Whittaker functions. We benchmark our results against numerical calculations and show how the method can be used to extract asymptotic normalization coefficients for charged-particle bound states. Furthermore, the results we derive here have immediate applications for calculations of atomic nuclei in finite periodic volumes for the case where the leading finite-volume correction is associated with two charged clusters.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

A numerical Poisson solver with improved radial solutions for a self-consistent locally scaled self-interaction correction method

Abstract The universal applicability of density functional approximations is limited by self-interaction error made by these functionals. Recently, a novel one-electron self-interaction-correction (SIC) method that uses an iso-orbital indicator to apply the SIC at each point in space by scaling the exchange-correlation and Coulomb energy densities was proposed. The locally scaled SIC (LSIC) method is exact for the one-electron densities, and unlike the well-known Perdew–Zunger SIC (PZSIC) method recovers the uniform electron gas limit of the uncorrected density functional approximation, and reduces to PZSIC method as a special case when isoorbital indicator is set to the unity. Here, we present a numerical scheme that we have adopted to evaluate the Coulomb potential of the electron density scaled by the iso-orbital indicator required for the self-consistent LSIC calculations. After analyzing the behavior of the finite difference method (FDM) and the green function solution to the radial part of the Poisson equation, we adopt a hybrid approach that uses the FDM for the Coulomb potential due to the monopole and the GF for all higher-order terms. The performance of the resultant hybrid method is assessed using a variety of systems. The results show improved accuracy than earlier numerical schemes. We also find that, even with a generic set of radial grid parameters, accurate energy differences can be obtained using a numerical Coulomb solver in standard density functional studies.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Anomalous ellipticity dependence of the generation of near-threshold harmonics in noble gases

We measured the ellipticity dependence of the harmonic yield of He, Ne, Ar, Kr, and Xe. It is found that the harmonic yield in general decreases monotonically with increasing laser ellipticity but anomalous dependence does occur for each atom for some harmonic orders. Additionally, we found that the degree of anomaly is weakest in helium but is stronger for the heavier atoms. The anomaly also depends on laser intensity. To explain these features, we developed a quantum trajectory Monte Carlo (QTMC) model based on the quantum path integral theory to study high-order harmonic generation (HHG). The model includes the effect of Coulomb potential from the target ion, thus it is capable of calculating near-threshold harmonics quantitatively. This model reveals that the presence of Coulomb potential would generate orbiting trajectories which are responsible for the anomalous ellipticity dependence near the threshold.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS1983

Modification of the photoionization of hydrogen by a low-frequency laser

The modification of the photoionization cross section due to the presence of an intense low-frequency laser is obtained for the hydrogen-atom target. Previous calculations of this effect have neglected the effects of the Coulomb potential in the final state. A sum rule was obtained which indicates that the photoelectric current is unaffected by the presence of the low-frequency laser to lowest order in the laser frequency. It is found that the Coulomb potential changes this result and that the photoelectric current is changed by a factor which depends upon the laser parameters and the electron momentum. A simple experiment is suggested to confirm this prediction.

Fiordilino, E.↗
DOE OSTI2021

Excited Rydberg states in MoSe 2 /WSe 2 heterostructures

The functional form of Coulomb interactions in the transition metal dichalcogenides (TDMs) and other van der Waals solids is critical to many of their unique properties, e.g. strongly-correlated electron states, superconductivity and emergent ferromagnetism. This paper presents measurements of key excitonic energy levels in MoSe 2 /WSe 2 heterostructures. These measurements are obtained from resonance Raman experiments on specific Raman peaks only observed at excited states of the excitons. This data is used to validate a model of the Coulomb potential in these structures which predicts the exciton energies to within ~5 meV. This model is used to determine the effect of heterostructure formation on the single-particle band gaps of the layers and will have a wide applicability in designing the next generation of more complex TDM structures.

Materials Science↗
DOE OSTI2020

Tunable discrete scale invariance in transition-metal pentatelluride flakes

Log-periodic quantum oscillations discovered in transition-metal pentatelluride give a clear demonstration of discrete scale invariance (DSI) in solid-state materials. The peculiar phenomenon is convincingly interpreted as the presence of two-body quasi-bound states in a Coulomb potential. However, the modifications of the Coulomb interactions in many-body systems having a Dirac-like spectrum are not fully understood. Here, we report the observation of tunable log-periodic oscillations and DSI in ZrTe 5 and HfTe 5 flakes. By reducing the flakes thickness, the characteristic scale factor is tuned to a much smaller value due to the reduction of the vacuum polarization effect. The decreasing of the scale factor demonstrates the many-body effect on the DSI, which has rarely been discussed hitherto. Furthermore, the cut-offs of oscillations are quantitatively explained by considering the Thomas-Fermi screening effect. Finally, our work clarifies the many-body effect on DSI and paves a way to tune the DSI in quantum materials.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A method to predict texture effect on ion beam channeling analysis of polycrystals and the application to study the mosaic spreading effect in highly oriented pyrolytic graphite

We propose a method to convert the channeling Rutherford backscattering spectrum yield map of a single crystal to a polycrystal through a matrix rotation technique. The rotation matrix is determined by the deviation of the crystal axial direction from the original z axis. The final yield map is created after averaging the rotated yields using the texture function as the weight factor. For highly oriented pyrolytic graphite (HOPG) exhibiting mosaic spread, the method leads to a Gaussian kernel averaging of the map obtained from a single crystal. The yield map of a single crystal is obtained by a simulation of ion trajectories in a potential field described by Moliere screened Coulomb potentials. Yield maps are calculated under various σ values (standard deviations of mosaic spread). The simulated results are compared with experimental results obtained using 1.2 MeV alpha particle. σ is extracted through the best fitting, demonstrating that the method can be used to obtain texture details. The effects of mosaic spread on minimum yield χ min and the half-width at half maximum of angular scans ψ ½ are systematically modeled and compared with previous theoretical equations. The study also shows that previous theoretical equations are valid only at small σ values. The proposed method can be applied to any type of polycrystal and is not limited to HOPG. It provides near-surface mosaic spread and crystallography information with a longitudinal depth resolution of tens of nanometers and is not influenced by grain shapes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1990

The effects of vacuum polarization on thermonuclear reaction rates

Added to the pure Coulomb potential, the contribution from vacuum polarization increases the barrier, reducing the wave function (u) for reacting nuclei within the range of nuclear forces. The cross section and reaction rate are then reduced accordingly by a factor proportional to u squared. The effect is treated by evaluating the vacuum polarization potential as a small correction to the Coulomb term, then computing u in a WKB formulation. The calculation is done analytically employing the small r power-series expansion for the Uehling potential to express the final result in terms of convenient parameters. At a temperature of 1.4 x 10 to the 7th K the (negative) correction is 1.3 percent for the fundamental fusion process p + p yields d + e(+) + nu.

Gould, Robert J.↗