Search NASA⌕ Search

SEARCH · Search NASA

Results for “ELECTROSTATICS”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Functionalized Porous Polymer Networks as High-Performance PFAS Adsorbents

Toxic per- and polyfluoroalkyl substances (PFAS) are now found in nearly every water source on the planet. Exposure to these molecules can have negative health consequences, but the low concentration of PFAS relative to other solutes in water makes their removal challenging. Adsorbents offer a promising treatment route, but often exhibit low selectivities and removal capacities, as well as slow kinetics. The performance in these metrics can be improved by chemically optimizing PFAS binding sites and maximizing PFAS-adsorbent interactions. To explore how to achieve this, a porous polymer network solid (PPN-6, also known as PAF-1) was postsynthetically modified with various chemical moieties capable of leveraging unique combinations of electrostatic, hydrogen-bonding, hydrophobic, and fluorophilic interactions with PFAS molecules. Batch adsorption experiments and computational studies revealed that electrostatic and hydrogen-bonding interactions drive short-chain PFAS adsorption, while hydrophobic and fluorophilic interactions improve long-chain PFAS adsorption. In complex water matrices, a combination of electrostatic and fluorophilic interactions led to the greatest total PFAS removal. The best-performing material, functionalized with a fluorinated alkylammonium (PPN-6-FNDMB), selectively adsorbs PFAS with high capacity (up to 4.0 mmol/g) and rapid kinetics (equilibrium reached in <30 s). Furthermore, PPN-6-FNDMB outperforms several commercial adsorbents, achieving near-complete removal of 21 different PFAS from a groundwater sample collected at a US Air Force base. The PFAS could subsequently be desorbed from PPN-6-FNDMB, concentrating them by a factor of over 50 times. The recycled PPN-6-FNDMB could then be reused with minimal losses in long-chain PFAS adsorption capacity over four cycles.

Pezoulas, Ethan R↗
DOE OSTI2026

How Does Water Dissociation Work in Bipolar Membranes?

Bipolar membranes (BPMs) create counteracting spatial gradients of pH and electrostatic potential in electrochemical systems, enabling applications in pH regulation, electrocatalysis, and separations. At the polarized junction of a BPM the water dissociation (WD, 2H2O ⇌ H3O+ + OH-) reaction can be driven, but it remains poorly understood. In this Perspective, we integrate molecular insights from bulk-water autoionization and the associated field effects with continuum descriptions of BPM electrostatics and experimental WD kinetic analyses to describe possible mechanisms of voltage-driven WD. Pristine BPM junctions highlight both the limits of primarily electric-field-driven WD and the practical challenges of junction stability at extreme reverse bias. Introducing heterogeneous catalyst layers, commonly metal oxides and graphene oxides, accelerates WD by orders of magnitude through hypothesized coupled effects in which surface acid-base functionality and high-density hydroxyl sites mediate proton-transfer steps, and catalyst mobile electronic/ionic charges redistribute the junction electric potential drop to shape the local electric fields and reactive microenvironments. Kinetic analyses suggest two regimes of heterogeneous WD mechanism, including field-driven ordering of interfacial water and a Second-Wien-Effect dissociation-barrier lowering. We conclude by defining the key unknown variables (local pH, electrostatic potential, catalyst charge state and relationships among mechanisms) and outlining experimental and multiscale modeling strategies needed for predictive WD catalysis and for controlling related ion-transfer reactions.

Wu, Yifan↗
DOE OSTI2025

Why the electron chooses one of two symmetry-related paths in the type II bacterial photosynthetic reaction centers

A series of electron and proton transfer reactions are carried out in photosynthetic reaction centers that convert the energy of light into high energy reduced products and add to the proton gradient. All photosynthetic reaction centers (photosystems I, II, and different bacterial reaction centers (bRCs)) have their charge separating cofactors arranged with two, c 2 -symmetric paths from a primary donor to a terminal acceptor. In Type II reaction centers (PSII and bRCs), electrons utilize only one active branch. In bRCs the electron transfer occurs in the A branch from the excited state of the bacteriochlorophyll (BChl) dimer P 860 through BChl A , bacteriopheophytin A (BPh A ), to ubiquinone Q A , and then to the B-branch Q B . B-branch BChl B and BPh B do not participate. A similar path is found in PSII. These proteins challenge our understanding of how proteins can turn on or turn off the activity of bound ligands. The shift of the electrochemical midpoint potential (E m ) of each cofactor by the protein environment is calculated using continuum electrostatics within the Multi-Conformation Continuum Electrostatic (MCCE) program. In bRCs from Rhodobacter sphaeroides, the electrostatic potentials favor reduction of the active branch cofactors. The residues that contribute to a more positive potential on the A branch are identified. These are often distant from the cofactors. The large difference of the potential at the competing BChls may direct the electron transfer to the A branch. Residues TyrM210 and PheL181 are at symmetrical positions on the A and B branches, near P 860 , between BChl A and BChl B . Mutated bRCs, that swap the Tyr and Phe (YF mutant), have been shown to support some electron transfer to the B-branch cofactors. Here, the calculated E m s with the YF swap show more positive B-branch potentials, which help explain why there is now electron transfer along the B branch.

Bacterial Reaction Center↗
DOE OSTI2026

Search for active and inactive ion insertion sites in organic crystalline materials

The position of mobile active and inactive ions, specifically ion insertion sites, within organic crystals, significantly affects the properties of organic materials used for energy storage and ionic transport. Identifying the positions of these atomic (and ionic) sites in an organic crystal is challenging, especially when the element has low X-ray scattering power, such as lithium (Li) and hydrogen, which are difficult to detect by powder X-ray diffraction. First-principles calculations, exemplified by density functional theory (DFT), are very practical for confirming the relative stability of ion positions in materials. However, the lack of effective strategies to identify ion sites in these organic crystalline frameworks renders this task extremely challenging. This work presents two algorithms: the (i) efficient location of ion insertion sites from extrema in electrostatic local potential and charge density (ELIISE), and the (ii) ElectRostatic InsertioN (ERIN), which leverage charge density and electrostatic potential fields accessed from first-principles calculations, combined with the Simultaneous Ion Insertion and Evaluation (SIIE) workflow. SIIE inserts all ions simultaneously—to determine ion positions in organic crystals. We demonstrate that these methods accurately reproduce known ion positions in 16 organic materials and identify previously overlooked low-energy sites in tetralithium 2,6-naphthalenedicarboxylate (Li4NDC), an organic electrode material, highlighting the importance of inserting all ions simultaneously, as in the SIIE workflow. These algorithms are also integrated with off-the-shelf machine learning potentials, yielding promising results comparable to first-principles findings.

Gopidi, Harshan Reddy [Argonne National Laboratory↗
DOE OSTI2025

Gate control of electron correlations towards Mott field effect transistors (MottFET)

The modulation of channel conductance in field-effect transistors (FETs) via metal-insulator-semiconductor structures has revolutionized information processing and storage. However, the limitations of silicon-based FETs in electrical switching have driven the search for new materials capable of overcoming these constraints. Electrostatic gating of competing electronic phases in a Mott material near its metal-to-insulator transition (MIT) offers the prospects of substantial modulation of the conducting electrons and electrical resistivity through small changes in band filling. While electrostatic control of the MIT has been previously reported, the advancement of Mott materials toward novel Mott-based transistors, MottFETs requires the realization of their unique physical properties in a solid-state device. In this study, we present gate control of electron correlation using a solid-state device utilizing the oxide Mott system La 1-x Sr x VO 3 as a correlated FET channel. We report on a gate resistance response that cannot be explained in a purely electrostatic framework. This behavior suggests an enhancement in effective mass and a reduction in effective carrier density as the system approaches its insulating state, consistent with theoretical predictions for Mott systems, suggesting at least 100× charge gain originating from the correlated behavior. These preliminary results pave the way toward the development of highly efficient, low-power electronic devices that could surpass the performance bottlenecks of conventional FETs by leveraging the electronic phase transitions of correlated electron systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Shadow molecular dynamics for flexible multipole models

Shadow molecular dynamics provide an efficient and stable atomistic simulation framework for flexible charge models with long-range electrostatic interactions. Shadow molecular dynamics simulations are driven by approximate “shadow” Born–Oppenheimer potentials for which the exact charges and forces are directly accessible without relying on costly (and approximate) iterative solvers. While previous implementations have been limited to atomic monopole charge distributions, we extend this approach to flexible multipole models. We derive detailed expressions for the shadow energy functions, potentials, and force terms, explicitly incorporating monopole–monopole, dipole–monopole, and dipole–dipole interactions. In our formulation, both atomic monopoles and atomic dipoles are treated as extended dynamical variables alongside the propagation of the nuclear degrees of freedom. We demonstrate that introducing the additional dipole degrees of freedom preserves the stability and accuracy previously seen in monopole-only shadow molecular dynamics simulations. In addition, we present a shadow molecular dynamics scheme where the monopole charges are held fixed while the dipoles remain flexible. Our extended shadow dynamics provide a framework for stable, computationally efficient, and versatile molecular dynamics simulations involving long-range interactions between flexible multipoles. This is of particular current interest in combination with machine-learned interatomic potentials, including long-range electrostatic interactions.

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

Probing the Nanoscale Excitonic Landscape and Confinement of Excitons in Gated Monolayer Semiconductors

Engineering and probing excitonic properties at the nanoscale remains a central challenge in quantum photonics and optoelectronics. While exciton confinement via electrical control and strain engineering has been demonstrated in 2D semiconductors, substantial nanoscale heterogeneity limits the scalability of 2D quantum photonic device architectures. Here, in this work, we use cathodoluminescence spectroscopy to probe the excitonic landscape of monolayer WS 2 under electrostatic gating. Exploiting the high spatial resolution of the converged electron beam, we resolve a homojunction arising between gated and ungated regions. Moreover, we reveal an exciton confinement channel arising from an unconventional doping mechanism driven by the interplay between the electron beam and the applied gate fields. These findings offer new insights into the optoelectronic behavior of monolayer semiconductors under the combined influence of electron-beam excitation and electrostatic gating. Our approach provides a pathway for exciton manipulation at the nanoscale and opens opportunities for controlling quantum-confined exciton transport in 2D materials.

2D materials↗
DOE OSTI2024

Long‐Range Non‐Coulombic Coupling at the LaAlO 3 /SrTiO 3 Interface

The LaAlO 3 /SrTiO 3 interface hosts a plethora of gate‐tunable electronic phases. Gating of LaAlO 3 /SrTiO 3 interfaces is usually assumed to occur electrostatically. However, increasing evidence suggests that non‐local interactions can influence and, in some cases, dominate the coupling between applied gate voltages and electronic properties. Here, quasi‐1D ballistic electron waveguides are sketched at the LaAlO 3 /SrTiO 3 interface as a probe to understand how gate tunability varies as a function of spatial separation. Gate tunability measurements reveal the scaling law to be at odds with the pure electrostatic coupling observed in traditional semiconductor systems. The non‐Coulombic gating at the interface is attributed to a long‐range nanoelectromechanical coupling between the gate and electron waveguide, possibly mediated by the ferroelastic domains in SrTiO 3 . The long‐range interactions at the LaAlO 3 /SrTiO 3 interface add unexpected richness and complexity to this correlated electron system.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Oversaturating Liquid Interfaces with Nanoparticle‐Surfactants

Abstract Assemblies of nanoparticles at liquid interfaces hold promise as dynamic “active” systems when there are convenient methods to drive the system out of equilibrium via crowding. To this end, we show that oversaturated assemblies of charged nanoparticles can be realized and held in that state with an external electric field. Upon removal of the field, strong interparticle repulsive forces cause a high in‐plane electrostatic pressure that is released in an explosive emulsification. We quantify the packing of the assembly as it is driven into the oversaturated state under an applied electric field. Physiochemical conditions substantially affect the intensity of the induced explosive emulsification, underscoring the crucial role of interparticle electrostatic repulsion.

Wu, Xuefei↗
DOE OSTI2025

GPU-accelerated kinetic Hall thruster simulations in WarpX

Abstract Two-dimensional (axial-azimuthal) simulations of a Hall thruster are performed using the open-source particle-in-cell code WarpX. The simulation conditions are chosen to match those of the axial-azimuthal benchmark first reported by Charoy et al. in 2019. A range of numerical and solver parameters is investigated in order to find those which yield the best performance. It is found that WarpX completes the benchmark case in 3.8 days on an Nvidia V100 GPU, and in as low as 1.5 days on a more recent Nvidia H100 GPU. Of the numerical parameters investigated, it is determined that the field-solve tolerance and particle resampling thresholds have the largest effect on the simulation wall time and that particle resampling may artificially widen electron velocity distribution functions, leading to unphysical heating. A semi-implicit scheme for the electrostatic field solve is tested and is found to produce results consistent to within 10% of the benchmark in less than twelve hours. The scaling properties of the electrostatic solver to multiple GPUs are also assessed on a uniform plasma test problem. The results of this work are discussed in the context of advancements in GPU hardware and the suitability of kinetic Hall thruster simulations for engineering applications.

Marks, Thomas A.↗
DOE OSTI2024

High-throughput measurements of CO 2 permeance and solubility in ionic liquid reveal a synergistic role of ionic interactions and void fractions

The factors that govern CO 2 solubility in ionic liquids (ILs) are of great interest for the development of new materials for CO 2 capture and utilization. The cationic functional group (i.e., imidazolium, pyrrolidinium, pyridinium, etc.), alkyl chain length of cation, degree of fluorination of anion, anion size, and the void fraction in IL are known to influence CO 2 solubility. However, a comprehensive explanation of how these factors collectively affect CO 2 solubility has not been developed yet. This knowledge gap is largely attributed to the lack of CO 2 solubility data for IL structures other than imidazolium based ILs. We report here an automated high-throughput (HT) setup for the measurement of CO 2 solubility in room-temperature ILs (RTILs) combining six different anions and nine different cations for a total of 19 different specific ranges of RTILs. The HT setup first dispenses up to 200 µL of RTILs in a 96-well microtiter plate and then utilizes a robotic arm to measure cyclic voltammogram (CV) in each well using maneuverable Ag electrodes. The Cottrell analysis of the CO 2 reduction CV peak provides a direct measurement of CO 2 permeance in RTILs, which yields Henry’s constant from the estimated diffusion coefficient of CO 2 . Henry’s constants thus obtained are in very good agreement with those reported earlier. The measured CO 2 permeance and Henry’s constant of all RTILs seem to follow a first-order dependence on void fraction and a second-order dependence on electrostatic interaction between anion and cation of IL, with some synergistic dependence on the product of a void fraction and electrostatic interaction, making them two important descriptors for the design of novel ILs.

CO2 Solubility↗
DOE OSTI2026

Theoretical model of dielectrophoretic water vapor condensation during multiphase flow of electrically charged water droplets in air

Water droplet growth in air, influenced by vapor diffusion and electrical charging, has significant implications for indoor air quality, cloud formation, and aerosol dynamics. Here, this study presents a comprehensive theoretical model that describes the growth dynamics of electrically charged water droplets under external electric fields, incorporating diffusion, electrostatic, and dielectrophoretic forces to enhance understanding and practical applications in air dehumidification and environmental control. A theoretical model was established to analyze the growth dynamics of the charged water droplets. The approach extended the thermodynamic framework of the Modified Kelvin-Thomson equation by integrating electrostatic and dielectrophoretic effects, enabling a comprehensive evaluation of vapor concentration near the surface of charged droplets. A sensitivity analysis was performed to assess the roles of droplet size, electric field intensity, surface charge, and vapor concentration on condensation rates. The study progressed beyond the classical diffusion-limited growth of neutral droplets, addressing the behavior of charged droplets in electrohydrodynamic environments. The model highlighted the significance of dielectrophoretic drift flux for dipolar vapor molecules, especially in micron-sized droplets ranging from 10 to 100 µm. Findings demonstrated that, for smaller droplets subjected to elevated electrical charge, average growth rates varied from 0.3 to 1.6 µm/s if the air relative humidity increased from 50 to 100 %. Analytical results are corroborated by numerical simulations, which exhibited strong concordance, thereby validating the proposed theoretical approach. This research introduced an innovative and scalable analytical methodology for electric field-enhanced droplet growth, offering valuable insights into the advancement of energy-efficient air dehumidification, water desalination, and atmospheric water harvesting technologies.

42 ENGINEERING↗
DOE OSTI2026

Cosolvent-tuned interactions in ionic liquids: A vibrational and quantum-chemical study of ethylene glycol ratio effects

Ionic liquids (ILs) are attractive media for CO 2 capture but remain limited by viscosity and cost. Blending ILs with ethylene glycol (EG) is a practical route to mitigate these constraints, yet the molecular origins of cosolvent effects and their dependence on composition are not well resolved. We combine Fourier-transform infrared (FT-IR) spectroscopy with quantum-chemical (DFT) analysis to elucidate how the IL:EG molar ratio modulates intermolecular interactions and electronic structure. Computed vibrational frequencies enable mode assignment and deconvolution of overlapping bands, revealing systematic, ratio-dependent shifts and broadenings in (i) EG O–H stretching, (ii) cation and EG C–H stretchings (imidazolium C2–H, C4–H, C5–H, methyl and ethyl groups, -CH2 of EG), (iii) anion signature modes (e.g., CN motifs), and (iv) EG C–O and C–C stretchings, consistent with the redistribution of hydrogen-bonding networks. Molecular electrostatic potential (MESP) maps quantify attenuation of extreme potential regions with increasing EG, indicating progressive screening of cation–anion electrostatic interactions. Quantum Theory of Atoms in Molecules (QTAIM) identifies emergent bond critical points between EG and the IL ions, while Reduced Density Gradient–Noncovalent Interaction (RDG–NCI) analysis differentiates strong directional hydrogen bonds from dispersive contacts across compositions. Together, these results show that EG fraction controls a switch from predominantly ion–ion to mixed ion–EG coordination, altering local polarity and polarizability that underlie the observed FT-IR trends. The framework provides composition–structure–spectrum relationships that can guide rational selection of IL:EG ratios to balance favorable molecular interactions with practical performance targets in scalable CO 2 capture systems.

DAC↗
DOE OSTI2025

Liquid Phase Modeling in Porous Media: Adsorption of Methanol and Ethanol in H-MFI in Condensed Water

Zeolites are used in the chemical and separation industries for their exceptional selectivity, adsorption capacity, regenerability, and stability in gas and liquid phase processing. Here, we developed an explicit solvation method for predicting solvent/condensed phase effects on adsorption free energies in microporous media such as zeolites based on the hybrid quantum mechanical/molecular mechanical free energy perturbation (QM/MM-FEP) technique. Our explicit solvation method for zeolite systems, called eSZS, aims to capture site-specific interactions during the adsorption process at the Brønsted acid sites of H-MFI zeolite while still considering the diverse configuration space of the solvent molecules. This strategy is ideal for chemical reactions or adsorbates that interact with the microporous medium in few distinct adsorbate/transition state configurations, i.e., the harmonic or similar approximations are acceptable for the adsorbate/transition state while such approximations break down for the solvent molecules that require extensive configuration space sampling. In this way, our approach effectively overcomes the limitations of implicit solvation models and classical force field methods for describing solvation effects on chemical reactions within porous materials such as zeolites. Specifically, in this study, we investigated various aspects of our hybrid QM/MM approach, including QM cluster size dependencies in a periodic electrostatically embedded cluster model (PEECM), rules for link atoms at the QM/MM boundary, and functional and basis set considerations for converged and reasonably accurate gas and aqueous phase methanol and ethanol adsorption free energy predictions in H-MFI. For gas phase adsorption of methanol and ethanol in H-MFI at a Brønsted acid site in T12 position, we compute adsorption free energies at 298 K of −0.61 and −0.75 eV, respectively, using a PEECM containing 50 Si and 1 Al atom with ωB97x-D/def2-TZVP level of theory. For solvent effect calculations, we sample the aqueous phase using grand canonical Monte Carlo (GCMC) simulations to (1) obtain a mean field of electrostatic interactions in the reaction system and (2) perform a rigorous free energy perturbation calculation. Similar to the experimentally and computationally observed endergonic solvation effects observed for hydrocarbon adsorption on metal surfaces, we also observe that a condensed aqueous environment destabilizes methanol and ethanol at these acid sites in H-MFI at 298 K. Specifically, the computed solvation free energies of adsorption (ΔΔG solv ) for methanol and ethanol are +0.44 and +0.54 eV, respectively. From this study, it is evident that adsorbates (methanol and ethanol) are competing with water for adsorption space inside the H-MFI zeolite, leading to an endergonic solvation effect. Here, we expect that the endergonic, aqueous solvent effect during adsorption in microporous zeolites is highly tunable by changing the pore size and hydrophobicity of the microporous material as this will affect the water density inside the pore structure.

Adsorption↗
DOE OSTI2025

STEPs-SOL, a Peptoid Force Field Parameterization to Include Solvent Effects

As peptoids (N-substituted glycines) continue to gain popularity as a class of biomimetic polymers, the importance and demand for accurate force fields in molecular simulations also grow. Building on the vacuum-optimized Systematic and Extensible Force Field for Peptoids (STEPs) force field, here we present STEPs-SOL, a novel peptoid force field parametrization that effectively incorporates solvent effects to enhance the accuracy of peptoid simulations. The development of STEPs-SOL is based on the need for precise electrostatic modeling achieved through solvent-specific partial charge optimization. Here, our systematic approach significantly improves agreement with experimental measurements, reducing the mean absolute error in cis/trans ratio predictions (ΔG c/t ) by an average of 38% across multiple peptoid residues and solvent environments. This improved parametrization addresses computational challenges associated with nonbonded energies while maintaining a workflow that relies on high-level quantum mechanical data rather than depending solely on limited experimental equilibrium properties. By evaluating the effects of conformational bias in restrained electrostatic potential (RESP) charge generation and examining their impact on peptoid conformations in various solvents, we enhance our understanding of peptoid structural dynamics while providing a more accurate modeling framework.

force field↗
DOE OSTI2025

Investigating Electron Conductivity Regimes in the Bacterial Cytochrome Wire OmcS

The anaerobic bacterium Geobacter sulfurreducens produces extracellular, electronically conductive cytochrome polymer wires that are conductive over micron length scales. Structure models from cryo-electron microscopy data show OmcS wires form a linear chain of hemes along the protein wire axis, which is proposed as the structural basis supporting their electronic properties. However, the mechanism by which this heme arrangement supports long-range electronic conduction remains unknown. Structure models from cryo-electron microscopy data show these wires form a linear chain of hemes along the protein wire axis, which is proposed as the structural basis supporting their electronic properties. Existing computational models using static heme redox potentials and coupling energies fail to explain experimental observations, predicting conductances 10,000 to 100,000 times lower than measured values. Here, we investigate how dynamic disorder affects site energies, interheme coupling, and long-range electronic conductivity within these cytochrome wires. We introduce an approach to extract charge carrier site information directly from Kohn–Sham density functional theory, without employing projector schemes, and show that site and coupling energies are highly sensitive to changes in interheme geometry and the surrounding electrostatic environment. Unlike models that incorporate dynamic disorder as a thermally averaged quantity, our quantum charge carrier model incorporates proxies for dynamic disorder through decoherence corrections, yielding predicted diffusion coefficient closer to what is expected from experiment and comparable with other organic-based electronic materials. Based on these simulations, we propose that the instantaneous fluctuations of the local electrostatic environment can transiently lift energy degeneracies and delocalize charge carriers. Furthermore, these studies reveal how incorporating dynamic fluctuations associated with the environment resolves the discrepancy between theory and experiment in microbial cytochrome wires and highlight design principles for bioinspired, heme-based conductive materials.

Bioinorganic chemistry↗
DOE OSTI2024

Voltage-Dependent First-Principles Barriers to Li Transport within Li-Ion Battery Solid Electrolyte Interphases

Charging a Li-ion battery requires Li-ion transport between the cathode and the anode. This Li-ion transport is dependent on (among other factors) the electrostatic environment that the ion encounters within the solid electrolyte interphase (SEI), which separates the anode from the surrounding electrolyte. A previous first-principles work has illuminated the reaction barriers through likely atomistic SEI environments but has had difficulty accurately reflecting the larger electrostatic potential landscape that an ion encounters moving through the SEI. In this work, we apply the recently developed quantum continuum approximation (QCA) technique to provide an equilibrium electronic potentiostat for first-principles interface calculations. Using QCA, we calculate the potential barrier for Li-ion transport through LiF, Li 2 O, and Li 2 CO 3 SEIs along with LiF–LiF and LiF–Li 2 O grain boundaries, all paired with Li metal anodes. Here, we demonstrate that the SEI potential barrier is dependent on the electrochemical potentials of the anode in each system. Finally, we use these techniques to estimate the change in the diffusion barrier for a Li ion moving in a LiF SEI as a function of the anode potential. We find that properly accounting for interface and electronic voltage effects significantly lowers reaction barriers compared with previous literature results.

25 ENERGY STORAGE↗
DOE OSTI2024

Predicting Partial Atomic Charges in Metal–Organic Frameworks: An Extension to Ionic MOFs

Molecular simulation is an invaluable tool to predict and understand the usage of metal–organic frameworks (MOFs) for gas storage and separation applications. Accurate partial atomic charges, commonly obtained from density functional theory (DFT) calculations, are often required to model the electrostatic interactions between the MOF and adsorbates, especially when the adsorbates have dipole or quadrupole moments, such as water and CO 2 . Machine learning (ML) models have been previously employed to predict partial charges and avoid the computational cost associated with DFT calculations. However, previous ML models suffer from small training data sets, which limit their scope of application. In this work, we introduce two novel machine learning models, PACMOF2-neutral and PACMOF2-ionic, aimed at predicting the density-derived electrostatic and chemical (DDEC6) partial atomic charges for both neutral and ionic MOFs. These models not only yield DFT-level accuracy at a fraction of the computational cost but also demonstrate a remarkable improvement in prediction of adsorption, as validated with grand canonical Monte Carlo simulations. Furthermore, the robustness and fast computational time of the PACMOF2 models, along with their transferability to other porous materials such as covalent organic frameworks and zeolites, underscores their potential in high-throughput screening of MOFs for diverse applications.

36 MATERIALS SCIENCE↗