Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Coulomb”

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 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 OSTI2024

Combination of Counterion Size and Doping Concentration Determines the Electronic and Thermoelectric Properties of Semiconducting Polymers

In both chemical and electrochemical doping of organic semiconductors (OSCs), a counterion, either from the electrolyte or ionized dopant, balances the charge introduced to the OSC. Despite the large influence of this counterion on OSC optical and electronic response, there remains substantial debate on how a fundamental parameter, ion size, impacts these properties. Here, this work resolves much of this debate by accounting for two doping regimes. In the low-doping regime, the Coulomb binding energies between charge carriers on the OSC and the counterions are significant, and larger counterions lead to decreased Coulomb interactions, more delocalized charge carriers, and higher electrical conductivities. In the high-doping regime, the Coulomb binding energies become negligible due to the increased dielectric constant of the films and a smoothing of the energy landscape; thereby, the electrical conductivities depend primarily on the extent of morphological disorder in the OSC. Moreover, in regioregular poly(3-hexylthiophene), rr-P3HT, smaller counterions lead to greater bipolaron concentrations in the low-doping regime due to the increased Coulomb interactions. Emphasizing the impact of the counterion size, it is shown that larger counterions can lead to increased thermoelectric power factors for rr-P3HT.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Critical review of the Mogi failure criterion based on true-triaxial laboratory data analysis and theoretical considerations

The influence of the intermediate principal stress, σ 2 , is disregarded when rock strength is analyzed using the Mohr-Coulomb or Hoek-Brown criteria. Several polyaxial strength criteria have been proposed to account for the influence of σ 2 . One of the methods that received considerable attention is the Mogi criterion. The criterion modifies the Nadai or Drucker-Prager approaches by exchanging the normal octahedral stress with what is introduced as the mean effective stress, σ m,2 . This change is based on experimental evidence that the influence of the intermediate principal stress on strength is less pronounced than the influence of the minimum principal stress. However, the Mogi criterion was derived somewhat empirically and consequently it is essential to analyze its potential theoretical deficiencies. In the current note, we analyzed the Mogi failure criterion, describing its characteristics in both τ oct –σ m,2 and σ 1 -σ 2 spaces. We investigated the relation between Mohr-Coulomb and Mogi criteria, which marks the basis for the commonly used Mogi-Coulomb approach. Here, we compared the performance of the Mogi criterion against available polyaxial strength data using three fitting approaches: Mogi-Coulomb based on triaxial compression results only, linear Mogi based on all polyaxial strength results, and power law Mogi also based on all data. The analysis reveals several serious deficiencies of the Mogi approach which all suggest that the Mogi criterion cannot properly represent the σ 2 dependence of rock strength in many cases.

58 GEOSCIENCES↗
DOE OSTI2022

Understanding the Origin of the Nonpassivating Behavior of Si-Based Anodes during the Initial Cycles

In this contribution, we combined electrochemical cycling and X-ray photoelectron spectroscopy (XPS) to understand the nonpassivating behavior of the solid electrolyte interphase (SEI) on Si anodes during the first cycles. Based on galvanostatic measurements, we show that the irreversible capacity loss is reduced after the first cycle, and it stays almost constant from the second cycle onwards. XPS was used to determine the root causes of the Coulombic inefficiency, showing that the rate of decomposition of the organic solvents strongly decreased after the first cycle, whereas the rate of salt decomposition is almost unchanged between cycles. We determine that the inhibition of the decomposition reaction of the organic solvent is responsible for the lower Coulombic loss during the second electrochemical cycle in comparison to the first, whereas the nonpassivating behavior toward the salt decomposition is one of the main causes of capacity loss upon cycling. Here, we further revisit the role of cracking in contributing to capacity loss. Whereas high volumetric expansion remains an issue plaguing the performance of Si anodes, our chronoamperometry studies reveal that the SEI formed on Si anodes does not passivate even when the electrode is fully expanded, and no additional surface is exposed. Overall, our work establishes the need to address the chemical and electrochemical instability of the SEI on the Si anode in addition to the more notorious issue of cracking.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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 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 OSTI2021

Distribution function of nuclei from $\mathcal{e}^±$ scattering in the presence of a strong primordial magnetic field

The amplitude of the primordial magnetic field (PMF) is constrained from observational limits on primordial nuclear abundances. Within this constraint, it is possible that nuclear motion is regulated by Coulomb scattering with electrons and positrons (e ± s), while e ± s are affected by a PMF rather than collisions. For example, at a temperature of 10 9 K, thermal nuclei typically experience ~ 10 21 scatterings per second that are dominated by very small angle scattering leading to minuscule changes in the nuclear kinetic energy of order $\mathcal{O}$(1) eV. In this paper the upper limit on the effects of a possible discretization of the e ± momenta by the PMF on the nuclear momentum distribution is estimated under the extreme assumptions that the momentum of the e ± is relaxed before and after Coulomb scattering to Landau levels, and that during Coulomb scattering the PMF is neglected. This assumption explicitly breaks the time reversal invariance of Coulomb scattering, and the Maxwell-Boltzmann distribution is not a trivial steady state solution of the Boltzmann equation under these assumptions. Here we numerically evaluate the collision terms in the Boltzmann equation, and show that the introduction of a special direction in the e ± distribution by the PMF generates no directional dependence of the collisional destruction term of nuclei. Large anisotropies in the nuclear distribution function are then constrained from big bang nucleosynthesis. Ultimately, we conclude that a PMF does not significantly affect the isotropy or big bang nucleosynthesis.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Viscosity of the magnetized strongly coupled one-component plasma

The viscosity tensor of the magnetized one-component plasma, consisting of five independent shear viscosity coefficients, a bulk viscosity coefficient, and a cross coefficient, is computed using equilibrium molecular dynamics simulations and the Green-Kubo relations. A broad range of Coulomb coupling and magnetization strength conditions are studied. Magnetization is found to strongly influence the shear viscosity coefficients when the gyrofrequency exceeds the Coulomb collision frequency. Three regimes are identified as the Coulomb coupling strength and magnetization strength are varied. Here, the Green-Kubo relations are used to separate kinetic and potential energy contributions to each viscosity coefficient, showing how each contribution depends upon the magnetization strength. The shear viscosity coefficient associated with the component of the pressure tensor parallel to the magnetic field, and the two coefficients associated with the component perpendicular to the magnetic field, are all found to merge to a common value at strong Coulomb coupling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Charging and ion ejection dynamics of large helium nanodroplets exposed to intense femtosecond soft xray pulses

Ion ejection from charged helium nanodroplets exposed to intense femtosecond soft X-ray pulses is studied by single-pulse ion time-of-flight (TOF) spectroscopy in coincidence with small-angle X-ray scattering. Scattering images encode the droplet size and absolute photon flux incident on each droplet, while ion TOF spectra are used to determine the maximum ion kinetic energy, Ekin, of He+ j fragments (j = 1-4). Measurements span HeN droplet sizes between N similar to 10(7) and similar to 10(10) (radii R-0 = 78-578 nm), and droplet charges between similar to 9x10(-5) and similar to 3x10(-3) e/atom. Conditions encompass a wide range of ionization and expansion regimes, from departure of all photoelectrons from the droplet, leading to pure Coulomb explosion, to substantial electron trapping by the electrostatic potential of the charged droplet, indicating the onset of hydrodynamic expansion. The unique combination of absolute Xray intensities, droplet sizes, and ion Ekin on an event-by-event basis reveals a detailed picture of the correlations between the ionization conditions and the ejection dynamics of the ionic fragments. The maximum E-kin of He+ is found to be governed by Coulomb repulsion from unscreened cations across all expansion regimes. The impact of ion-atom interactions resulting from the relatively low charge densities is increasingly relevant with less electron trapping. The findings are consistent with the emergence of a charged spherical shell around a quasineutral plasma core as the degree of ionization increases. The results demonstrate a complex relationship between measured ion E-kin and droplet ionization conditions that can only be disentangled through the use of coincident single-pulse TOF and scattering data.

Helium↗
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 OSTI2023

Limits on the Range and Rate of Change in Power Take-Off Load in Ocean Wave Energy Conversion: A Study Using Model Predictive Control

Previous work comparing power take-off (PTO) architectures for ocean wave-powered reverse osmosis suggests that variable displacement in the wave energy converter (WEC)-driven pump does not offer a significant performance advantage. A limitation of that study is that the WEC was subject to a constant load within a given sea state (“Coulomb damping”) and did not account for controlled, moment-to-moment variation of the PTO load enabled by a variable displacement pump. This study explores the potential performance advantage of a variable PTO load over Coulomb damping. Model predictive control is used to provide optimal load control with constraints on the PTO load. The constraints include minimum and maximum loads and a limit on the rate of load adjustment. Parameter studies on these constraints enable conclusions about PTO design requirements in addition to providing an estimated performance advantage over Coulomb damping. Numerical simulation of the Oyster 1 WEC is carried out with performance weighted by historical sea state data from Humboldt Bay, CA. The results show a performance advantage of up to 20% higher yearly-average power absorption over Coulomb damping. Additionally, the parameter studies suggest that the PTO load should be adjustable down to at least 25% of the maximum load and should be adjustable between the minimum and maximum loads within a few seconds.

16 TIDAL AND WAVE POWER↗
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 OSTI2021

A three-step model of high harmonic generation using complex classical trajectories

Highlights: • Rigorous trajectory formulation of strong field physics and high harmonic generation. • The Coulomb ground state is represented by an ensemble of complex trajectories. • Trajectories transition seamlessly from the ground state to ionization or recollision. • The ionized wavepacket is reproduced well from the complex trajectories. • Classically allowed and forbidden processes are treated in a single framework. We present a new trajectory formulation of high harmonic generation that treats classically allowed and classically forbidden processes within a single dynamical framework. Complex trajectories orbit the nucleus, producing the stationary Coulomb ground state. When the field is turned on, these complex trajectories continue their motion in the field-dressed Coulomb potential and therefore tunnel ionization, unbound evolution and recollision are described within a single, seamless framework. The new formulation can bring mechanistic understanding to a broad range of strong field physics effects.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Deuteron-Deuteron Elastic and Three and Four-Body Breakup Scattering Using the Faddeev-Yakubovskii Equations

The deuteron-deuteron elastic and three and four-body breakup scattering cross section have been calculated using the Faddeev-Yakubovskii (FY) chain-of-partition momentum-space equations. In this calculation the initial two-cluster potential is split into separable and non-separable components, and the effective potential is reduced to elastic two-body and three and four-body breakup open-channels and a closed-channel many-body contribution. The closed-channel contribution is determined by minimizing a variational bound. The Coulomb interaction was included by expanding the initial and final Coulomb states in a Coulomb-Sturmian basis. The three sets of chain-of-partition integral equations were solved for the elastic and three and four-body breakup scattering amplitudes. The calculations were performed for the S = 2 spin/L = 0 angular-momentum state. The elastic and double-breakup calculations were performed for energies up to E = 5.48 MeV, while the single-breakup calculations were performed for energies up to E = 4.17 MeV. In the case of elastic scattering, the calculated scattering length of 5 a dd = 7.8 ± 0.3 fm is in good agreement with a FY cluster reduction calculation. The calculated phase shift is smaller than that predicted by the resonating group model and this difference is believed to be due to the differences in the potential and calculation methods. The breakup cross sections were calculated as a function of initial deuteron momentum and fragmented-deuteron momentum. Here, the d+d → d+n+p cross sections were compared with neutron yield measurements and, while the measurements also included the L > 0 components, the general features were consistent. Estimates of the calculational uncertainties/bias are provided.

42 ENGINEERING↗
DOE OSTI2022

Exploring the strong-coupling region of SU( N ) Seiberg-Witten theory

We consider the Seiberg-Witten solution of pure N = 2 gauge theory in four dimensions, with gauge group SU(N). A simple exact series expansion for the dependence of the 2(N – 1) Seiberg-Witten periods a I (u), a DI (u) on the N – 1 Coulomb-branch moduli un is obtained around the Z 2N -symmetric point of the Coulomb branch, where all u n vanish. This generalizes earlier results for N = 2 in terms of hypergeometric functions, and for N = 3 in terms of Appell functions. Using these and other analytical results, combined with numerical computations, we explore the global structure of the Kähler potential K = 1/2Σ I Im(a¯ I a DI ), which is single valued on the Coulomb branch. Evidence is presented that K is a convex function, with a unique minimum at the Z 2N -symmetric point. Finally, we explore candidate walls of marginal stability in the vicinity of this point, and their relation to the surface of vanishing Kähler potential.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Fluorinated co-solvent promises Li-S batteries under lean-electrolyte conditions

Practical lithium-sulfur batteries stipulate the use of a lean-electrolyte and a high Coulombic efficiency of the lithium metal anode. In this work, we employ 1,1,2,2-tetrafluoroethyl 2,2,2-trifluoroethyl ether as a co-solvent in the electrolyte of Li-S batteries to meet the demands. The co-solvent fosters the formation of a LiF-rich solid electrolyte interphase on lithium metal anode, as revealed by Ab initio molecular dynamics and Femtosecond Stimulated Raman spectroscopy. The co-solvent results in an average Coulombic efficiency of 99.4% for lithium plating/stripping cycles. Full cells with a capacity ratio of 2 between lithium anode and the S@PAN cathode (3 gs cm –2 ) exhibited a stable cycle life over 100 cycles at an electrolyte/sulfur ratio of 2 µL mg –1 , validating the high Coulombic efficiency of the lithium anode and demonstrating the compatibility of the electrolyte with both electrodes. To enhance the energy density, we prepared a hybrid cathode composed of 45 wt.% VS 2 mixed with ZnS-coated Li 2 S@graphene as the cathode. Based on the mass of both electrodes and the electrolyte, the full cell delivers an energy density of 483 Wh kg –1 , which demonstrates viable Li-S batteries with the lean-electrolyte.

25 ENERGY STORAGE↗
DOE OSTI2022

Phosphoric acid pre-treatment to tailor polybenzimidazole membranes for vanadium redox flow batteries

We repot vanadium redox flow batteries (VRFBs) use ion-selective membranes for transporting ionic species while separating the positive and negative electrolytes. In this paper, we report phosphoric acid doped polybenzimidazole (PBI) membranes that yield high ionic selectivity and conductivity in VRFBs. The phosphoric acid pre-treatment swells the PBI matrix irreversibly and increases its sulfuric acid doping in VRFB electrolytes. The pre-treated membranes show comparable area resistance with Nafion-212 with better selectivity towards vanadium ions. A low resistance was obtained with a reduced membrane thickness, which can reduce the overall cost of materials. The VRFB using an optimized phosphoric acid pre-treated PBI membrane demonstrates coulombic, voltage, and energy efficiencies of 99.7, 90.3, and 90.0%, respectively, at 40 mA cm -2 while achieving ~40% higher discharge capacity compared to Nafion-212. Furthermore, the VRFB cell shows excellent cycling stability, i.e., only 6.8% or 0.068% per cycle capacity decay for 100 cycles while maintaining coulombic efficiency at >99.9% (98.5% coulombic efficiency for Nafion-212 with a capacity decay of 0.42% per cycle) proving the effectiveness of pre-treatment of phosphoric acid on PBI membranes.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Systematic Evaluation of Carbon Hosts for High-Energy Rechargeable Lithium-Metal Batteries

Rechargeable lithium batteries (RLBs) have attracted wide attention for achieving high energy and low cost, but the poor stability of the Li anode has remained as a large challenge. In the last several years, great efforts have been made to develop carbon anode structures for RLB applications. A wide range of carbon structures as well as many chemical modification strategies have been investigated to improve the Coulombic efficiency and extend the cycle life of the anodes or the cells. However, currently there is no clear conclusion on key factors that determine the performance of such the types of carbon structures due to the wide range of experimental conditions used in such studies. In this study, we systematically evaluated and compared the electrochemical properties of representative carbon anode materials as well as chemically functionalized carbon materials in different electrolytes. We have found that chemical functionalization, metallic loading, and the types of micro or nano-structures all have an effect on Li deposition/stripping process, but by far the most critical factor is the nature of the electrolytes. All carbon materials produced poor Coulombic efficiency and cycle life in conventional carbonate or ether-based electrolytes. However, when a Li metal compatible is used, the average Coulombic efficiency of Li metal deposition on all carbon structures is increased to 99% or higher. This study provides critical lessons on the directions for carbon anodes for RLBs. Based on these findings, we constructed Li metal coin cells using prelithiated hard carbon, LiNi0.6Mn0.2Co0.2O2 (NMC622) cathode, and compatible electrolyte and tested the cells under conditions of lean electrolyte, lean Li metal and high cathode loading required for more than 300 Wh kg-1 specific energy. More than 200 stable cycles were demonstrated with over 80% capacity retention. Furthermore, the existing challenges and potential approaches have been discussed to further extend the cyclability of lithium-metal batteries under practical conditions.

Liu, Yao↗