SEARCH · Engineering Papers

Results for “electrostatic”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 361 records · Page 20

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Ultra-low-level measurements of airborne fission products from the Fukushima Daiichi reactor accident using high volume collection systems at Savannah River National Laboratory

For this study, ultra-low-level measurements of radionuclides in air have been conducted at the Savannah River National Laboratory (SRNL) to determine the atmospheric concentration of fission products released following the Fukushima Daiichi reactor accident on March 11, 2011. Air filter samples were acquired from two high-volume collection systems (a traditional filter-based system and an electrostatic precipitator-based system) to monitor airborne radionuclide concentrations in the period covering from 2 weeks to 3 years after the disaster. The world-wide spread of low-level concentrations of airborne fission products from the Fukushima event provided a unique opportunity to demonstrate SRNL's electrostatic particle collection technology and other improvements in environmental monitoring developed at the Savannah River Site (SRS). Detecting and analyzing the release allowed a comprehensive test of SRS systems for monitoring environmental radioactivity. Gamma-ray-emitting fission products ( 131,132 I, 134,136,137 Cs, and 129,132 Te) and cosmogenic isotopes ( 7 Be and 22 Na) in air were detected and quantified by high-resolution gamma-ray spectroscopy at concentrations as low as 0.07 μBq per standard cubic meter (SCM) (50 mBq total 137 Cs), while plutonium content was quantified by thermal ionization mass spectrometry (TIMS) at concentrations as low as 6.5 × 10 –21 g/SCM (3.0 fg 239+240 Pu). Isotope concentrations measured at SRNL from gamma-ray spectroscopy were compared to independent measurements from Chapel Hill, NC, located approximately 370 km (230 mi) NE of SRNL. Meteorological modeling was also used to predict radionuclide transport from the location of release to both measurement locations.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗

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 OSTI2024

Neutron diffraction from a microgravity-grown crystal reveals the active site hydrogens of the internal aldimine form of tryptophan synthase

Pyridoxal 5′-phosphate (PLP), the biologically active form of vitamin B6, is an essential cofactor in many biosynthetic pathways. The emergence of PLP-dependent enzymes as drug targets and biocatalysts, such as tryptophan synthase (TS), has underlined the demand to understand PLP-dependent catalysis and reaction specificity. The ability of neutron diffraction to resolve the positions of hydrogen atoms makes it an ideal technique to understand how the electrostatic environment and selective protonation of PLP regulates PLP-dependent activities. Facilitated by microgravity crystallization of TS with the Toledo Crystallization Box, we report the 2.1 Å joint X-ray/neutron (XN) structure of TS with PLP in the internal aldimine form. Positions of hydrogens were directly determined in both the α- and β-active sites, including PLP cofactor. The joint XN structure thus provides insight into the selective protonation of the internal aldimine and the electrostatic environment of TS necessary to understand the overall catalytic mechanism.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2021

Large gyro-orbit model of ion velocity distribution in plasma near a wall in a grazing-angle magnetic field

A model is presented for the ion distribution function in a plasma at a solid target with a magnetic field $\boldsymbol {B}$ inclined at a small angle, $\alpha \ll 1$ (in radians), to the target. Adiabatic electrons are assumed, requiring $\alpha \gg \sqrt {Zm_{e}/m_{i}}$ , where $m_{e}$ and $m_{i}$ are the electron and ion mass, respectively, and $Z$ is the charge state of the ion. An electric field $\boldsymbol {E}$ is present to repel electrons, and so the characteristic size of the electrostatic potential $\phi$ is set by the electron temperature $T_{e}$ , $e\phi \sim T_{e}$ , where $e$ is the proton charge. An asymptotic scale separation between the Debye length $\lambda _{D} = \sqrt {\epsilon _0 T_{{e}} / e^{2} n_{{e}} }$ , the ion sound gyro-radius $\rho _{s} = \sqrt { m_{i} ( ZT_{e} + T_{i} ) } / (ZeB)$ and the size of the collisional region $d_{c} = \alpha \lambda _{\textrm {mfp}}$ is assumed, $\lambda _{D} \ll \rho _{s} \ll d_{c}$ . Here $\epsilon _0$ is the permittivity of free space, $n_{e}$ is the electron density, $T_{i}$ is the ion temperature, $B= |\boldsymbol {B}|$ and $\lambda _{\textrm {mfp}}$ is the collisional mean free path of an ion. The form of the ion distribution function is assumed at distances $x$ from the wall such that $\rho _{s} \ll x \ll d_{c}$ , that is, collisions are not treated. A self-consistent solution of the electrostatic potential for $x \sim \rho _{s}$ is required to solve for the quasi-periodic ion trajectories and for the ion distribution function at the target. The large gyro-orbit model presented here allows to bypass the numerical solution of $\phi (x)$ and results in an analytical expression for the ion distribution function at the target. It assumes that $\tau =T_{i}/(ZT_{e})\gg 1$ , and ignores the electric force on the quasi-periodic ion trajectory until close to the target. For $\tau \gtrsim 1$ , the model provides an extremely fast approximation to energy–angle distributions of ions at the target. These can be used to make sputtering predictions.

Physics↗

DOE OSTI2020

Closing the Gap Between Modeling and Experiments in the Self-assembly of Biomolecules at Interfaces and in Solution

Molecular self-assembly is a powerful tool in materials design, wherein non-covalent interactions like electrostatic, hydrophobic, hydrogen bonding, and van der Waals can be exploited to produce supramolecular nanostructures that are functional and highly tunable. Biomolecules are attractive building blocks, as they are biocompatible, biodegradable and adopt a wide array of higher order structures. Moreover, naturally occurring protein systems display a manifold of structures and interactions that can be replicated in synthetic biomolecules. In this perspective, we highlight advances in multiscale simulation techniques across broad spatiotemporal scales that can aid in characterizing self-assembly of hybrid and hierarchical bionanomaterial systems, with an emphasis on physics-based simulation approaches currently employed to study biomolecules at mineral interfaces. The power of these approaches is highlighted across a few recent areas where molecular simulations have advanced our understanding of self-assembly spanning peptides to protein self-assembly. Looking forward, we discuss how in the near future emerging methods in statistical and machine learning will advance this research field in all areas from expanding the capabilities of physics-based simulation methods to enabling new analyses of high throughput experiments. These advances will pave the way for understanding the molecular recognition patterns in systems that are dictated by self-assembly - biomineralizing peptides, hierarchical peptoids, and large protein assemblies, and will aid in the development of a new synthesis science for achieving precise molecular control in materials design

Sampath, Janani↗

DOE OSTI2021

Polymorphism and Optoelectronic Properties in Crystalline Supramolecular Polymers

Supramolecular polymers can emulate some of the physical properties of covalent polymers but offer new opportunities given the possibility of designing monomers that will form highly ordered assemblies with defined shapes. Internally ordered supramolecular polymers formed through nucleation–elongation self-assembly are well-known but highly crystalline examples which exhibit important properties such as light harvesting, charge transport, and ferroelectricity are not common. We report here on a detailed study of supramolecular polymers formed in water by carboxylated naphtho-p-quinodimethane amphiphiles. We found that supramolecular polymerization of these amphiphiles in aqueous media yields crystalline assemblies with morphologies that included ribbons, helically rolled ribbons, and twisted filaments. This polymorphism was found to be controlled exclusively by repulsive electrostatic interactions controlled by the degree of protonation of the carboxylic head groups which also dictates the nature of supramolecular packing. Substoichiometric amounts of base lead to highly crystalline ribbons due to a decreased surface charge density and less electrostatic repulsion. Increasing deprotonation results in helically rolled ribbons with a different polymorph crystal lattice, whereas excessive deprotonation leads to twisted filaments with maximum surface charge density. Ribbons, helical rolled ribbons, and twisted filaments revealed an increasing red shift in their visible absorption maxima. These crystalline assemblies could be potential candidates for solar energy materials and photocatalytic systems.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Covalent Functionalization of Nickel Phosphide Nanocrystals with Aryl-Diazonium Salts

Covalent functionalization of Ni 2 P nanocrystals was demonstrated using aryl-diazonium salts. Spontaneous adsorption of aryl functional groups was observed, with surface coverages ranging from 20 to 96% depending on the native reactivity of the salt as determined by the aryl substitution pattern. Increased coverage was possible for low reactivity species using a sacrificial reductant. Functionalization was confirmed using thermogravimetric analysis, Fourier transform infrared spectroscopy, and X-ray photoelectron spectroscopy. The structure and energetics of this nanocrystal electrocatalyst system, as a function of ligand coverage, were explored with density functional theory calculations. The Hammett parameter of the surface functional group was found to linearly correlate with the change in Ni and P core-electron binding energies and the nanocrystal's experimentally and computationally determined work function. The electrocatalytic activity and stability of the functionalized nanocrystals for hydrogen evolution were also improved when compared to the unfunctionalized material, but a simple trend based on electrostatics was not evident. Density functional theory was used to understand this discrepancy, revealing that H adsorption energies on the covalently functionalized Ni 2 P also do not follow the electrostatic trend and are predictive descriptors of the experimental results.

14 SOLAR ENERGY↗

DOE OSTI2023

Role of Complexation Strength on the Photophysical and Transport Properties of Semiconducting Charged Polymer Complexes

The high polymer fraction in complexes of conjugated and insulating polyelectrolytes offers unique opportunities for the fabrication of conductive thick films and bulk structures. The electrostatic interactions in these systems further provide a handle for controlling their structure and properties. The impact of charge-mediated complexation strength on the photophysical and electronic transport properties in blends of conjugated polyelectrolytes (CPEs) with oppositely charged polymeric ionic liquids (PILs) was examined. Complexes were formed with varying frequency of charged repeat units, from 50 to 100%, on an anionic polythiophene-based CPE and a complimentary cationic PIL. In highly charged complexes, the intimate mixing between the CPE and the PIL reduced the structural disorder along the CPE backbone, enhancing its intrachain conjugation and interchain stacking. In weakly charged complexes (<90%), these chain planarization effects were absent and microphase separation occurred. At all charge fractions examined, the electrical conductivity of an acid-doped complex was higher than that of the unblended constituent CPE. Further, the highest electrical conductivity, near 1 S cm –1 , was found for a charge fraction of 100%. These results demonstrate the potential for designing effective polymeric conductors using electrostatic complexation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

(Fe, Cr)(OH) 3 Coprecipitation in Solution and on Soil: Roles of Surface Functional Groups and Solution pH

The simultaneous precipitation of (Fe, Cr)(OH) 3 nanoparticles in solution (homogeneous) and on soil surfaces (heterogeneous), which controls Cr transport in soil and aquatic systems, was quantified for the first time in the presence of model surfaces, i.e., bare and natural organic matter (NOM)-coated SiO 2 and Al 2 O 3 . Various characterization techniques were combined to explore the surface-ion-precipitate interactions and the controlling mechanisms. (Fe, Cr)(OH) 3 accumulation on negatively charged SiO 2 was mainly governed by electrostatic interactions between hydrolyzed ion species or homogeneous (Fe, Cr)(OH) 3 and surfaces. The elevated pH through protonation of Al 2 O 3 surface hydroxyls resulted in higher Cr/Fe ratios in both homogeneous and heterogeneous coprecipitates. Here, due to ignorable NOM adsorption onto SiO 2 , the amounts of (Fe, Cr)(OH) 3 precipitates on bare/NOM-SiO 2 were similar; contrarily, attributed to favored NOM adsorption onto Al 2 O 3 and consequently carboxyl association with metal ions or (Fe, Cr)(OH) 3 nanoparticles, remarkably more heterogeneous precipitates harvested on NOM-Al 2 O 3 than bare-Al 2 O 3 . With the same solution supersaturation, the total amounts of homogeneous and heterogeneous precipitates were similar irrespective of the substrate type. With lower pH, decreased electrostatic forces between substrates and precipitates shifted (Fe, Cr)(OH) 3 distribution from heterogeneous to homogeneous phases. The quantitative knowledge of (Fe, Cr)(OH) 3 distribution and the controlling mechanisms can assist in better Cr sequestration in natural and engineered settings.

(Fe, Cr)(OH)3 distribution↗

DOE OSTI2022

Harnessing Intermolecular Interactions to Promote Long-Lived Photoinduced Charge Separation from Copper Phenanthroline Chromophores

In this study, facilitating photoinduced electron transfer (PET) while minimizing rapid charge-recombination processes to produce a long-lived charge-separated (CS) state represents a primary challenge associated with achieving efficient solar fuel production. Natural photosynthetic systems employ intermolecular interactions to arrange the electron-transfer relay in reaction centers and promote a directional flow of electrons. This work explores a similar tactic through the synthesis and ground- and excited-state characterization of two Cu(I)bis(phenanthroline) chromophores with homoleptic and heteroleptic coordination geometries and which are functionalized with negatively charged sulfonate groups. The addition of sulfonate groups enables solubility in pure water, and it also induces assembly with the dicationic electron acceptor methyl viologen (MV 2+ ) via bimolecular, dynamic electrostatic interactions. The effect of the sulfonate groups on the ground- and excited-state properties was evaluated by comparison with the unsulfonated analogues in 1:1 acetonitrile/water. The excited-state lifetimes for all sulfonated complexes are similar to what we expect from previous literature, with the exception of the sulfonated heteroleptic complex whose metal-to-ligand charge-transfer (MLCT) lifetime in water has two components that are fit to 10 and 77 ns. For the sulfonated complexes, we detected reduced MV +$\bullet$ in both solvent environments following MLCT excitation, but control measurements in 1:1 acetonitrile/water with the unsulfonated analogues showed no PET to MV 2+ , indicating that electrostatically driven supramolecular assemblies of the sulfonated complexes with MV 2+ facilitate the observed PET. Additionally, the strength of the intermolecular interactions driving the formation of these assemblies changes drastically with the solvent environment. In 1:1 acetonitrile/water, PET occurred from both sulfonated complexes with quantum yields (Φ ET ) of 2-3% but increased to a remarkable 98% for the sulfonated heteroleptic complex with a 3 μs CS-state lifetime in water.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Fast and Accurate Machine Learning Strategy for Calculating Partial Atomic Charges in Metal–Organic Frameworks

Computational high-throughput screening using molecular simulations is a powerful tool for identifying top-performing metal–organic frameworks (MOFs) for gas storage and separation applications. Accurate partial atomic charges are often required to model the electrostatic interactions between the MOF and the adsorbate, especially when the adsorption involves molecules with dipole or quadrupole moments such as water and CO 2 . Although ab initio methods can be used to calculate accurate partial atomic charges, these methods are impractical for screening large material databases because of the high computational cost. We developed a random forest machine learning model to predict the partial atomic charges in MOFs using a small yet meaningful set of features that represent both the elemental properties and the local environment of each atom. The model was trained and tested on a collection of about 320 000 density-derived electrostatic and chemical (DDEC) atomic charges calculated on a subset of the Computation-Ready Experimental Metal–Organic Framework (CoRE MOF-2019) database and separately on charge model 5 (CM5) charges. The model predicts accurate atomic charges for MOFs at a fraction of the computational cost of periodic density functional theory (DFT) and is found to be transferable to other porous molecular crystals and zeolites. In conclusion, a strong correlation is observed between the partial atomic charge and the average electronegativity difference between the central atom and its bonded neighbors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Solvation Thermodynamics of Solutes in Water and Ionic Liquids Using the Multiscale Solvation-Layer Interface Condition Continuum Model

We report molecular assembly processes are generally driven by thermodynamic properties in solutions. Atomistic modeling can be very helpful in designing and understanding complex systems, except that bulk solvent is very inefficient to treat explicitly as discrete molecules. The SLIC/CDC multiscale model combines continuum solvent electrostatics based on the solvent layer interface condition (SLIC) with statistical thermodynamic models for hydrogen bonding and nonpolar modes: cavity formation, dispersion interactions, combinatorial mixing (CDC). The SLIC/CDC model predicts Gibbs energies of solvation for a database of 500 solutes in water with average accuracy better than 1 kcal/mol both for experimental measurements and for explicit-solvent molecular dynamics simulations. The separate SLIC/CDC energy mode values agree quantitatively with those computed from explicit-solvent molecular dynamics. The SLIC/SASA multiscale model combines the SLIC continuum electrostatic model with the solvent accessible surface area (SASA) nonpolar energy mode. The SLIC/SASA model predicts Gibbs energies of solvation with better than 1.4 kcal/mol average accuracy in aqueous systems and better than 1.6 kcal/mol average accuracy in ionic liquids. Both models predict solvation entropies, and are the first implicit-solvation models capable of predicting solvation heat capacities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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 OSTI2023

Theoretical Examination of the Hydroxide Transport in Cobaltocenium-Containing Polyelectrolytes

Polymers incorporating cobaltocenium groups have received attention as promising components of anion-exchange membranes (AEMs), exhibiting a good balance of chemical stability and high ionic conductivity. In this work, we analyze the hydroxide diffusion in the presence of cobaltocenium cations in an aqueous environment based on the molecular dynamics of model systems confined in one dimension to mimic the AEM channels. In order to describe the proton hopping mechanism, the forces are obtained from the electronic structure computed at the density-functional tight-binding level. We find that the hydroxide diffusion depends on the channel size, modulation of the electrostatic interactions by the solvation shell, and its rearrangement ability. Hydroxide diffusion proceeds via both the vehicular and structural diffusion mechanisms with the latter playing a larger role at low diffusion coefficients. The highest diffusion coefficient is observed under moderate water densities (around half the density of liquid water) when there are enough water molecules to form the solvation shell, reducing the electrostatic interaction between ions, yet there is enough space for the water rearrangements during the proton hopping. Furthermore, the effects of cobaltocenium separation, orientation, chemical modifications, and the role of nuclear quantum effects are also discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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 OSTI2023

An Ab Initio -Derived Force Field for Amorphous Silica Interfaces for Use in Molecular Dynamics Simulations

Here, we present a classical interatomic force field, silica-DDEC, to describe the interactions of amorphous and crystalline silica surfaces, parametrized using density functional theory-based charges. Charge schemes for silica surfaces were developed using the density-derived electrostatic and chemical (DDEC) method, which reproduces atomic charges of the periodic models as well as the electrostatic potential away from the atom sites. Lennard–Jones parameters were determined by requiring the correct description of (i) the amorphous silica density, coordination defects, and local coordination geometry, relative to experimental measurements, and (ii) water-silica interatomic distances compared with ab initio results. Deprotonated surface silanol sites are also described within the model based on DDEC charges. The result is a general electronic structure-derived model for describing fully flexible amorphous and crystalline silica surfaces and interactions of liquids with silica surfaces of varying structure and protonation state.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Dynamics of Electric Polarization and Relaxation of Ions at Humid Calcite Surfaces

Mobile ions at mineral surfaces can respond to an applied electric field, adopting a new distribution that effectively represents polarization of the electrical double layer. When the field is released, the ions relax to their equilibrium distribution. In both cases, the dynamics are characteristic of the interface. However, current models of electrokinetic phenomena are not sufficiently robust to accurately predict collective ion dynamics at structurally and chemically specific mineral–water interfaces. Here, in this study, we use electrostatic force microscopy (EFM) to investigate the dynamics of ion relaxation at hydrated calcite (104) surfaces at controlled relative humidity (RH). Electrically biased probes are used to polarize the distributions of calcium and carbonate ions that are intrinsic to this interface across a range of RH values. Polarization kinetics are tracked by monitoring the tip–sample force gradient during charging, and EFM imaging is used to characterize the spatial relaxation dynamics after the applied field is released. Electrostatic finite element modeling of the sample/probe system across length-scales from nanometers to millimeters reproduces the observed stretched exponential charging response. Together, these results allow us to estimate the ion diffusivities at the interface across a wide range of RH values. These diffusivities increase by roughly 5 orders of magnitude as the RH is increased from 5 to 90%, highlighting the critical role of adsorbed water for surface ion solvation that enables ion mobility.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗