Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ELECTROSTATIC CHARGING”

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

Water Dielectric Function at Finite Wavelength and the Convergence of Its Large Wavelength Limit (Static Dielectric) by Fluctuation Formulas

Here we provide a general analysis of the longitudinal dielectric function ε l (k) via fluctuation formulas at finite wavelength k, including the static dielectric constant ε w = ε l (k → 0), and analyze the different sources of errors, deriving explicit formulas. Simulations with the SPC/E water model show that the convergence of fluctuation formulas to compute the static dielectric constant is slow and requires long simulations. The analysis of ε l (k) allows us to identify the long- and short-wavelength limits and provide a precise determination of the location of the poles as well as the zeros in the complex plane, which determine the short-distance form of the electrostatic (Coulomb) interaction among charged particles.

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

Adsorption and Catalytic Activity of Gold Nanoparticles in Mesoporous Silica: Effect of Pore Size and Dispersion Salinity

The assembled state of nanoparticles (NPs) within porous matrices plays a governing role in directing their biological, electronic, and catalytic properties. However, the effects of the spatial confinement and environmental factors, such as salinity, on the NP assemblies within the pores are poorly understood. In this study, we use adsorption isotherms, spectrophotometry, and small-angle neutron scattering to develop a better understanding of the effect of spatial confinement on the assembled state and catalytic performance of gold (Au) NPs in propylamine-functionalized SBA-15 and MCM-41 mesoporous silica materials (mSiO 2 ). We carry out a detailed investigation of the effect of pore diameter and ionic strength on the packing and spatial distribution of AuNPs within mSiO 2 to get a comprehensive insight into the structure, functioning, and activity of these NPs. We demonstrate the ability of the adsorbed AuNPs to withstand aggregation under high salinity conditions. We attribute the observed preservation of the adsorbed state of AuNPs to the strong electrostatic attraction between oppositely charged pore walls and AuNPs. The preservation of the structure allows the AuNPs to retain their catalytic activity for a model reaction in high salinity aqueous solution, here, the reduction of p-nitrophenol to p-aminophenol, which otherwise is significantly diminished due to bulk aggregation of the AuNPs. This fundamental study demonstrates the critical role of confinement and dispersion salinity on the adsorption and catalytic performance of NPs.

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

Effect of Flexibility in Molecular Simulations of Carbon Dioxide Adsorption and Diffusion in a Cuprous Triazolate Framework

Using fixed atom force fields to model gas adsorption in flexible metal organic frameworks (MOFs) is known to pose difficulties in accurately reproducing and predicting experimental results. This paper studies the difference in accuracy between flexible and fixed atom force fields in reproducing CO 2 adsorption measurements in MAF-2 ([Cu(etz)]∞ (MAF-2, Hetz) 3,5-diethyl-1,2,4-triazole), an NbO-type triazolate scaffold with a bcu cavity system and attached ethyl groups. The flexible force field used to run the hybrid molecular dynamics and grand canonical Monte Carlo calculations were generated using the QuickFF software incorporating van der Waals parameters from the Universal Force Field (UFF) and density derived electrostatic and chemical (DDEC) charges. The fixed atom force field used was composed of UFF van der Waals parameters together with DDEC charges. The calculations were run at 298 K and at pressures of 0.1, 0.3, 0.5, 0.8, and 1 bar. It was observed that for this MOF the rigid force field overpredicted gas adsorption, whereas the flexible force field values closely matched experimental data. In the flexible structure, the freely moving ethyl groups of MAF-2 hindered adsorption, reducing the interaction energy between CO 2 and the N atoms of the triazolate framework as well as reducing the size of the largest cavity diameter. Here, the combination of these factors led to improved prediction of adsorption values with the flexible forcefield as compared to the rigid forcefield, demonstrating the need for modeling MOFs in a way more indicative of their behavior.

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

Synthesis and Characterization of Uranium Complexes Supported by Substituted Aryldimethylsilylanilide Ligands

Here, we report the synthesis and characterization of substituted aryldimethylsilyldiisopropylanilide ligands and their respective bisamido complexes of U(III), (3,5-R 2 -PhMe 2 SiNDipp) 2 UI(dioxane) x (1, R = H, x= 0; 2, R = Me, x = 0; 3, R = t Bu, x = 1). We found that the steric bulk of the 3,5-R 2 -Ph ring affects the hapticity of the U–arene interaction. In the solid-state, 1 is a U–(η 6 -arene) complex, while 2 is a bis(U–(η 1 -arene)) complex. Theoretically calculated bond orders at PBE0 and PBE0-D3 levels of theory support these hapticity assignments. The 3,5- t Bu 2 -Ph rings of 3 are too bulky to interact with U and solid-state metrical parameters initially suggested a U–(η 1 -arene) interaction with one of the Dipp rings. However, bond order calculations show that this interaction is even weaker than in the previously reported ((PhMe 2 Si) 2 N) 3 U complex, leading to the conclusion that 3 is best described as a U–(η o -arene) complex. Molecular orbital analyses in conjunction with electron localization methods reveal that the U–(η 6 -arene) bonding in 1 is primarily electrostatic in nature. Some charge transfer takes place from the arene π orbitals to the U 6d/5f hybrid orbitals in addition to subtle δ-back-bonding. In 2 and 3, both π and δ interactions are substantially weaker, in agreement with the differences in the U–arene coordination modes. Surprisingly, attempts to generate less sterically bulky (3,5-R 2 -PhMe 2 SiNPh) 2 UI complexes results in disproportionation to homoleptic tetraamido (3,5-R 2 -PhMe 2 SiNPh) 4 U(IV) (4, R = H; 5, R = Me) complexes.

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

Photochemical Upconversion in Water Using Cu(I) MLCT Excited States: Role of Energy Shuttling at the Micellar/Water Interface

Photochemical upconversion (UC) through triplet–triplet annihilation (TTA), which employs a visible absorbing triplet photosensitizer and an annihilator, is a process that generates a high energy photon from two lower energy photons. TTA-UC has been largely developed in pure organic solvents and solid-state polymeric constructs while featuring near exclusive use of rare and expensive metals within the photosensitizer. In this current investigation, we demonstrate that TTA-UC from the long lifetime earth-abundant photosensitizer [Cu(dsbtmp) 2 ](PF) 6 (dsbtmp = 2,9-di(sec-butyl)-3,4,7,8-tetramethyl-1,10-phenanthroline), abbreviated as Cu-PS, functions in water through encapsulation within a cationic-based assembly. Cetyltrimethylammonium bromide was the surfactant of choice as it electrostatically binds the negatively charged water-soluble 10-phenylanthracene-9-carboxylate (PAC) acceptor/annihilator and ultimately facilitates energy transfer across the interface. Efficient triplet–triplet energy transfer (TTET) from Cu-PS to the PAC acceptor was achieved in this aqueous assembly. Unfortunately, the hindered mobility of the PAC moieties ultimately hampered the annihilation process, and this was reflected in attenuated TTA rates and efficiencies. The combined experimental data illustrated that the water-soluble PAC acceptor was able to vectorially deliver the excited-state energy stored in Cu-PS across the interface into the bulk aqueous solution by engaging in excited-state electron transfer with methyl viologen acceptors. Lastly, these results are important for remotely operating photoredox reactions in water while rendering a photosensitizer spatially isolated in the hydrophobic core of a micelle.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Construction of an Anion-Pillared MOF Database and the Screening of MOFs Suitable for Xe/Kr Separation

The separation of xenon/krypton (Xe/Kr) mixtures is a challenging process. Many porous materials allow the adsorption of both Xe and Kr but only with low selectivity. Anion-pillared metal–organic frameworks (MOFs), featuring the anion groups as structural pillars, show potential in gas separations, but only a limited number of them have been synthesized. Here, we describe a collection of 936 anion-pillared MOFs based on 22 experimentally available structures. We performed density functional theory (DFT) optimization and then assigned density-derived electrostatic and chemical (DDEC) charges for each MOF to make them well suited to molecular simulations. The structural properties of the MOFs vary more strongly with the choice of the organic ligand than with other aspects like fluorine groups and metal centers. We then screened the entire collection of MOFs in the context of Xe/Kr separation at room temperature. Compared with previously reported MOFs, the interpenetrated MOF SIFSIX-6-Cd-i is predicted to perform better for Xe/Kr separations, with a good balance between working capacity (1.62 mmol/g) and separation selectivity (16.4) at 298 K and 100 kPa. Here, we also found that the heterogeneity of fluorine groups within a MOF can help to enhance Xe working capacity without reducing the Xe/Kr selectivity, suggesting that synthesis of anion-pillared MOFs with mixed fluorine groups may lead to improved Xe/Kr separation performance.

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

Ethanol-Induced Condensation and Decondensation in DNA-Linked Nanoparticles: A Nucleosome-like Model for the Condensed State

Inspired by the conventional use of ethanol to induce DNA precipitation, ethanol condensation has been applied as a routine method to dynamically tune “bond” lengths (i.e., the surface-to-surface distances between adjacent nanoparticles that are linked by DNA) and thermal stabilities of colloidal crystals involving DNA-linked nanoparticles. However, the underlying mechanism of how the DNA bond that links gold nanoparticles changes in this class of colloidal crystals in response to ethanol remains unclear. Here, we conducted a series of all-atom molecular dynamic (MD) simulations to explore the free energy landscape for DNA condensation and decondensation. Our simulations confirm that DNA condensation is energetically much more favorable under 80% ethanol conditions than in pure water, as a result of ethanol’s role in enhancing electrostatic interactions between oppositely charged species. Moreover, the condensed DNA adopts B-form in pure water and A-form in 80% ethanol, which indicates that the higher-order transition does not affect DNA’s conformational preferences. We further propose a nucleosome-like supercoiled model for the DNA condensed state, and we show that the DNA end-to-end distance derived from this model matches the experimentally measured DNA bond length of about 3 nm in the fully condensed state for DNA where the measured length is 16 nm in water. Altogether, this study provides an atomistic understanding of the mechanism underlying ethanol-induced condensation and water-induced decondensation, while our proposed nucleosome-like model allows the design of new strategies for interpreting experimental studies of DNA condensation.

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

Alkali-Metal Interlocking of 2D V 4 O 10 Sheets Defines Discretized Interlayer Shear Relationships

Low-dimensional materials manifest structural anisotropy, quantum confinement, and tightly bound excitonic states, which make them attractive building blocks that can be assembled within three-dimensional laterally stitched heterostructures, stacked van der Waals solids, and complex moiré superlattices. Ion intercalation in the galleries between layered materials provides a means of modifying interlayer separation and coupling, but it is also known to drive the shearing of the layers. In this article, we explore the distinct ligand coordination environments afforded by vanadyl oxygens of singular [V 4 O 10 ] sheets and examine how the size, polarizability, and stoichiometry of Group I cations sandwiched between such layers determine the interlocking of the sheets in stacked structures. Based on the topochemical insertion of alkali-metal ions into the layered λ-V 2 O 5 , we identify seven types of guest ion coordination sites discretized into four distinct regimes of interlayer shear in units of half octahedral widths. The coordination preferences of intercalated cations govern how they interlock 2D [V 4 O 10 ] sheets and engender specific shear conformations. We present evidence that static and dynamic disorder in guest ion arrangement modulate the magnetic structure of the intercalated compounds based on electrostatic polarization, localization of charge and spin density, and lattice distortion. The results illustrate the use of topochemical ion insertion to modulate stacking relationships and magnetic transition characteristics.

36 MATERIALS SCIENCE↗
DOE OSTI2022

“Did you eat a MOLEcule today?” An Interactive Demonstration of Molecular Weight with Real-Time Breath Analysis Using Mass Spectrometry for All Ages

Although mass spectrometry is a widely used analytical tool, age-appropriate, interactive outreach activities for laboratory visitors, especially children, are lacking. Furthermore, the presented interactive demonstration, “Did you eat a MOLEcule today?”, introduces all ages to molecular weight concepts and mass spectrometry in a research laboratory, while connecting the concepts to real-world applications. Through real-time breath analysis, participants explore the concepts of molecular weight, electrostatic field manipulation of charged molecules, and analyte identification by mass analysis. This module is rapid and highly adaptable for outreach activities but also includes age- or classroom-appropriate variations to decrease or increase difficulty levels. The presented interactive demonstration has repeatedly been implemented, with over 2300 participants during six annual “Take Our Daughters & Sons to Work Day” and two corporate “Family Day” outreach activities, successfully engaging, exciting, and educating both kids and parents.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Unraveling the core of fuel cell performance: engineering the ionomer/catalyst interface

The biggest obstacle to the widespread implementation of polymer electrolyte membrane fuel cells (PEMFCs) is their cost, primarily due to the use of platinum catalysts. The high intrinsic catalyst activity exhibited on a rotating disk electrode (RDE) is rarely realized in a membrane electrode assembly (MEA), which is a long-standing challenge for PEMFCs and a cause of low catalyst utilization. To translate the high RDE performance of a catalyst into a MEA, the design of an ideal ionomer/catalyst interface is proposed: a thin, conformal ionomer film covers the maximum surface of a Pt nanoparticle and thus simultaneously maximizes catalyst utilization, (i.e., high mass activity and electrochemically active surface area) and O 2 diffusion rate (i.e., high current density performance) without compromising proton conduction. Building such an interface is a long-standing challenge due to the lack of interaction between the ionomer and catalyst particles, resulting in large ionomer agglomerates and inhomogeneous ionomer coverage over the catalyst nanoparticle, with consequent poor fuel cell performance. In this work, this ionomer/catalyst interface has been engineered, utilizing the electrostatic attraction between positively charged catalyst and negatively charged ionomer particles in a catalyst ink and preserved in a solid catalyst layer. As a result, this interface leads to previously unachieved proton exchange membrane fuel cell performance in terms of both catalyst utilization (75% vs. 45%) and peak/rated power density (i.e., 1.430/0.930 W cm -2 , H 2 /air, cathode Pt loading: 0.1 mgPt cm -2 ) for pure Pt catalysts, even better than those of Pt alloy catalysts. This work demonstrates the formation of an interface in the liquid phase (using ultra-small-angle X-ray scattering in combination with cryo-TEM, isothermal–titration–calorimetry) and the preservation of the interface in the solid catalyst layer (using TEM) and estimates the effective coverage and thickness of the ionomer film (using limiting current density, RDE and fuel cell performance).

25 ENERGY STORAGE↗
DOE OSTI2022

Highly stretchable ionically crosslinked acrylate elastomers inspired by polyelectrolyte complexes

Dynamic bonds are a powerful approach to tailor the mechanical properties of elastomers and introduce shape-memory, self-healing, and recyclability. Among the library of dynamic crosslinks, electrostatic interactions among oppositely charged ions have been shown to enable tough and resilient elastomers and hydrogels. In this work, we investigate the mechanical properties of ionically crosslinked ethyl acrylate-based elastomers assembled from oppositely charged copolymers. Using both infrared and Raman spectroscopy, we confirm that ionic interactions are established among polymer chains. We find that the glass transition temperature of the complex is in between the two individual copolymers, while the complex demonstrates higher stiffness and more recovery, indicating that ionic bonds can strengthen and enhance recovery of these elastomers. We compare cycles to increasing strain levels at different strain rates, and hypothesize that at fast strain rates ionic bonds dynamically break and reform while entanglements do not have time to slip, and at slow strain rates ionic interactions are disrupted and these entanglements slip significantly. Further, we show that a higher ionic to neutral monomer ratio can increase the stiffness, but its effect on recovery is minimal. Lastly, taking advantage of the versatility of acrylates, ethyl acrylate is replaced with the more hydrophilic 2-hydroxyethyl acrylate, and the latter is shown to exhibit better recovery and self-healing at a cost of stiffness and strength. The design principles uncovered for these easy-to-manufacture polyelectrolyte complex-inspired bulk materials can be broadly applied to tailor elastomer stiffness, strength, inelastic recovery, and self-healing for various applications.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Control of core–shell nanoparticles properties through plasma synthesis: a computational study

The improved properties of core–shell nanoparticles (CSNPs) over homogeneous nanoparticles (NPs) have expanded and diversified the applications of these nanomaterials. However, controlling the properties of CSNPs can be a challenging task. Low temperature plasmas have proven to be an effective method of producing NPs with uniform size and morphology, and high yield. That said, NP transport and growth dynamics are sensitive to LTP properties. We report on a computational investigation of the evolution of Ge–Si CSNP properties as a function of operating conditions through the modeling of a flowing, two-zone inductively coupled plasma (ICP) reactor. Ar/GeH 4 and Ar/SiH 4 gas mixtures were supplied to separate plasma zones at a pressure of 1 Torr to promote growth of Ge cores and Si shells. The negatively charged CSNPs are trapped electrostatically in the vicinity of the antennas where the plasma is generated and where the majority of particle growth occurs. Particles that grow to a critical size are then de-trapped by fluid drag due to neutral gas flow. A two-dimensional hybrid plasma model coupled with a three-dimensional kinetic NP transport model were utilized to resolve plasma chemistry and NP growth processes that take place on distinct timescales. The trends in CSNP properties and trapping mechanisms associated with flow rate, applied ICP power and inlet precursor fraction are discussed. While the spatial distribution of plasma produced radical species can have significant impact on the NP growth process, the NP transport dynamics are what ultimately dictates the growth environment that is unique to each particle and so determines their final dimension and composition. The key to optimizing reactor conditions involves controlling the spatial density of growth species and plasma profile as a means to tailor particle trapping dynamics suitable to produce CSNPs for a specific application.

36 MATERIALS SCIENCE↗
DOE OSTI2025

High-resolution X-ray structure of Gln143Asn manganese superoxide dismutase captures multiple hydrogen peroxide-binding sites

Human mitochondrial manganese superoxide dismutase (MnSOD) converts superoxide into hydrogen peroxide and molecular oxygen, serving as a key defence against oxidative damage. Despite extensive studies, the full structural characterization of H 2 O 2 -binding sites in MnSOD remains largely unexplored. Previous H 2 O 2 -soaked MnSOD structures have identified two distinct H 2 O 2 -binding sites: one directly ligated to the catalytic manganese (LIG position) and another at the active-site gateway (PEO position) between the second-shell residues Tyr34 and His30. In this study, a kinetically impaired Gln143Asn MnSOD variant is used to trap and explore additional H 2 O 2 -binding sites beyond the second-shell solvent gate. In the wild-type enzyme, Gln143 mediates proton transfers with the manganese-bound solvent (WAT1) to drive redox cycling of the metal, which is necessary for effective superoxide dismutation. Substitution with Asn stalls catalysis because the increased distance from WAT1 disrupts critical proton-coupled electron-transfer (PCET) events, and the redox cycling of the active-site metal is impaired. This, in turn, stalls the electrostatic cycling of positive charge on the enzyme surface and enhances the likelihood of trapping transient H 2 O 2 -bound states in this variant. The results reveal several H 2 O 2 molecules leading up to the active site, in addition to the canonical LIG and PEO positions.

Dasgupta, Medhanjali [University of Nebraska Medic↗
DOE OSTI2020

Chemically Cross-Linked Cellulose Nanocrystal Aerogels for Effective Removal of Cation Dye

In this study, porous aerogels were prepared by directional freeze-drying via cross-linking cellulose nanocrystals (CNCs) with poly(methyl vinyl ether-co-maleic acid) (PMVEMA) and poly(ethylene glycol) (PEG). The thermal properties and physical adsorption performance toward cation methylene blue dye of the obtained CNC aerogels were investigated. The maximum degradation temperature was increased from 324°C of CNCs to 355°C of cross-linked CNC aerogels. The dye adsorption isotherm results showed that the maximum methylene blue adsorption capacity of CNC aerogels was 116.2 mg g −1 , according to the Langmuir model, which was mainly due to the electrostatic attractions between negatively charged carboxyl groups or sulfonate groups on the CNC aerogles and cation MB molecules. The reusability test showed that the CNC aerogels contained the same dye adsorption performance in five adsorption/desorption cycles. Overall, this study described an ideal alternative for water purification with high dye adsorption capacity and enhanced physical performance.

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

Intramolecular Electrostatic Effects on O 2 , CO 2 , and Acetate Binding to a Cationic Iron Porphyrin

Noncovalent electrostatic interactions are important in many biological and chemical reactions, especially those that involve charged intermediates. There has been a growing interest in using electrostatic ligand designs—placing charges in the second coordination sphere—to improve molecular reactivity, catalysis, and electrocatalysis. For instance, an iron porphyrin bearing four cationic ortho -trimethylanilinium groups, Fe( o -TMA), has been reported to be an exceptional electrocatalyst for both the carbon dioxide reduction reaction (CO 2 RR) and the oxygen reduction reaction (ORR). These reactions involve many different steps, and it is not evident which steps are affected by the four positive charges, or why. By comparing Fe( o -TMA) with the related iron-tetraphenylporphyrin, this work examines how covalently positioned charged groups affect substrate binding and other key pre-equilibria of both the ORR and CO 2 RR, specifically acetate, dioxygen, and carbon dioxide binding. This study is among the first to directly measure the effects of electrostatics on ligand-binding. The results show that adding electrostatic groups to a catalyst design often results in a complex interplay of multiple effects, including changes in pre-equilibria prior to substrate binding, combinations of through-space and inductive contributions, and effects of ionic strength and solution dielectric. The inverse half-order dependence of binding constant on ionic strength is proposed as a clear marker for an electrostatic effect. Overall, the conclusions provide guidance for the increasingly popular electrostatic ligand designs in catalysis and other reactivity.

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

Electrostatic Asymmetry of Wurtzite Nanocrystals and Resulting Photocatalytic Properties

Efficient light absorption and high energy of charge carriers of zinc cadmium sulfide (ZCS) make this semiconductor attractive for many photocatalytic reactions. Despite marked successes in shape-controlled synthesis of ZCS central to their photocatalytic performance, recombination of charge carriers as they migrate through the nanoscale particles results in losses of excitation energy, markedly reducing the photocatalytic activity of ZCS and other heterogeneous photocatalysts. Here, in this study, we show that the electrostatic asymmetry of nanostructures, previously discovered for nearly spherical nanoparticles, also manifests in wurtzite ZCS with planar geometry. The electrostatic asymmetry assists charge separation and substantially increases the yield of photocatalytic reactions in monocrystalline ZCS. The synthesized ZCS nanorods and nanoplates with identical chemical composition were found to have markedly different photocatalytic activity for evolution of hydrogen in water. Despite much smaller specific surface areas, the ~500 nm wide nanoplates displayed a hydrogen evolution rate 12 times higher than the ~35 nm long nanorods, also outperforming other ZCS photocatalysts. Experimental and computational data indicate that the homo- and heterojunction-free ZCS nanoplates with continuous wurtzite lattice behave essentially as nanoscale dipoles with intrinsic dipole moment as high as 48.39 D per unit cell. Electric-field-directed migration of charge carriers stimulates their localization on opposite parts of the nanoplates. Direct imaging of the intraparticle electrical field using off-axis electron holography confirmed their electrostatic asymmetry. Polarization-enhanced charge separation provides a new pathway to efficient and stable photocatalysts for sustainable energy technologies.

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