SEARCH · Engineering Papers

Results for “interface diffusion”

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 1,423 records · Page 79

All sources NASA NTRS DOE OSTI

NASA NTRS2020

NOAA-20 Visible Infrared Imaging Radiometer Suite day–night band on-orbit calibration and performance

The NOAA-20 Visible Infrared Imaging Radiometer Suite (VIIRS) instrument has been successfully operating on orbit since November 28, 2017. The day–night band (DNB) onboard NOAA-20 VIIRS is a panchromatic channel covering wavelengths from 0.5 to 0.9 μm, capable of observing the Earth scene in visible/near-infrared spectral range at a spatial resolution of 750 m. The DNB operates at low-, mid-, or high-radiometric gain stages, and it uses an onboard solar diffuser (SD) panel for low-gain stage calibration. The SD observations also provide a means to compute gain ratios between low-to-mid and mid-to-high-gain stages. With their large dynamic range and high sensitivity, the DNB detectors can make observations during both daytime and nighttime. We provide an assessment of the DNB on-orbit performance and behavior in the first two-year mission period and beyond. The calibration methodology used by the VIIRS Characterization Support Team in support of the NASA Earth science community has been described. The trending of on-board calibrators dark-offsets, SD gains, and gain ratios, and signal-to-noise ratio at minimum radiance have been analyzed, especially during key events such as the nadir and cryocooler doors opening. Furthermore, we performed intercomparison studies between Suomi National Polar-orbiting Partnership and NOAA-20 instruments and evaluated DNB radiometric calibration and characterization, including the SD degradation, detector gains, gain ratios, and straylight correction, as well as the calibration comparison between the NOAA Interface Data Processing Segment look-up-tables and our delivery results.

VIIRS↗

DOE OSTI2024

Optical and microstructural characterization of Er 3+ doped epitaxial cerium oxide on silicon

Rare-earth ion dopants in solid-state hosts are ideal candidates for quantum communication technologies, such as quantum memories, due to the intrinsic spin–photon interface of the rare-earth ion combined with the integration methods available in the solid state. Erbium-doped cerium oxide (Er:CeO 2 ) is a particularly promising host material platform for such a quantum memory, as it combines the telecom-wavelength (~ 1.5 μm) 4f–4f transition of erbium, a predicted long electron spin coherence time when embedded in CeO 2 , and a small lattice mismatch with silicon. In this work, we report on the epitaxial growth of Er:CeO 2 thin films on silicon using molecular beam epitaxy, with controlled erbium concentration between 2 and 130 parts per million (ppm). We carry out a detailed microstructural study to verify the CeO 2 host structure and characterize the spin and optical properties of the embedded Er 3+ ions as a function of doping density. In as-grown Er:CeO 2 in the 2–3 ppm regime, we identify an EPR linewidth of 245(1) MHz, an optical inhomogeneous linewidth of 9.5(2) GHz, an optical excited state lifetime of 3.5(1) ms, and a spectral diffusion-limited homogeneous linewidth as narrow as 4.8(3) MHz. We test the annealing of Er:CeO 2 films up to 900 °C, which yields narrowing of the inhomogeneous linewidth by 20% and extension of the excited state lifetime by 40%.

36 MATERIALS SCIENCE↗

DOE OSTI2020

In Operando Calorimetric Measurements for Activated Carbon Electrodes in Ionic Liquid Electrolytes under Large Potential Windows

This study aims to investigate the effect of the potential window on heat generation in carbon-based electrical double layer capacitors (EDLCs) with ionic-liquid (IL)-based electrolytes using in operando calorimetry. The EDLCs consisted of two identical activated-carbon electrodes with either neat 1-butyl-1-methylpyrrolidinium bis(trifluoromethane-sulfonyl)imide ([Pyr 14 ][TFSI]) electrolyte or 1.0 m [Pyr 14 ][TFSI] in propylene carbonate (PC) as electrolyte. The instantaneous heat generation rate at each electrode was measured under galvanostatic cycling for different potential windows ranging from 1 to 4 V. First, the heat generation rates at the positive and negative electrodes differed significantly in neat IL owing to the differences in the ion sizes and diffusion coefficients. However, these differences were minimized when the IL was diluted in PC. Second, for EDLC in neat [Pyr 14 ][TFSI] at high potential window (4 V), a pronounced endothermic peak was observed at the beginning of the charging step at the positive electrode owing to TFSI - intercalation in the activated carbon. On the other hand, for EDLC in 1.0 m [Pyr 14 ][TFSI] in PC at potential window above 3 V, an endothermic peak was observed only at the negative electrode owing to the decomposition of PC. Third, for both neat and diluted [Pyr 14 ][TFSI] electrolytes, the irreversible heat generation rate increased with increasing potential window and exceeded Joule heating. This was attributed to the effect of potential-dependent charge redistribution resistance. Additionally, a further increase in the irreversible heat generation rate was observed for the largest potential windows owing to the degradation of the PC solvent. Finally, for both types of electrolyte, the reversible heat generation rate increased with increasing potential window because of the increase in the amount of ion adsorbed/desorbed at the electrode/electrolyte interface.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

From suppressed void growth to significant void swelling in NiCoFeCr complex concentrated solid-solution alloy

Void swelling can result in dimensional instability and undermine the safe operation of nuclear reactors. Current strategies to inhibit void swelling mainly focus on enhancing defect absorption and recombination by introducing high-density defect sinks. Complex concentrated solid-solution alloys (CSAs), including high-entropy alloys, can withstand severe radiation damage due to their inherent chemical complexity without interfaces. However, the underlying mechanisms for void suppression in CSAs are far from clear. In this research, we studied the void evolution with respect to irradiation depths, doses, and temperatures in equiatomic NiCoFeCr under 3 MeV Ni ion irradiations. At relatively low doses (16 and 54 displacements per atom, dpa), voids form mainly outside of the ion-damaged region, and void formation in the peak damage region is suppressed, leading to negligible swelling. However, with further increase of dose (86 up to 250 dpa), significant void growth occurs in the peak damage region and extended dislocation lines dominate instead of short dislocation lines and loops formed at lower doses. From 500 to 700 °C, the dislocation density decreases while dislocations grow. Although the overall void swelling increases dramatically from 500 to 580 °C at 54 dpa, void growth in the peak damage region is still suppressed. The transition from suppressed void growth to significant void swelling is attributed to dislocation evolution and local chemical inhomogeneity (enrichment of Fe/Cr in the matrix) at higher doses. Our study shows that controlling element diffusion and defect evolution through tuning chemical complexity can further enhance the swelling resistance of CSAs.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Effects of temperature and dose rate on ion-irradiated γ-LiAlO 2 pellets

Defect accumulation and microstructural evolution during ion irradiation at elevated temperatures are governed by competing processes of defect production, driven by the dose rate, and defect recovery, controlled by diffusion, interaction, and annihilation. Here, this study investigates the effects of irradiation temperature and the dose rate on microstructural evolution, deuterium retention, and lithium volatilization in γ-LiAlO 2 pellets subjected to sequential He + and D + ion irradiation. Experiments were performed to a total fluence of 3 × 10 17 (He + + D + )/cm 2 at 623, 673, 723, and 773 K with an average He + dose rate of 7.7 × 10 −4 dpa/s, and to 2 × 10 17 (He + + D + )/cm 2 at 773 K with dose rates of 6.8 × 10 −5 , 2.9 × 10 −4 , and 7.3 × 10 −4 dpa/s. At 623 K, the microstructure was dominated by cavities and fractures with no observable precipitate formation, while small precipitates emerged at 673 K. Increasing the irradiation temperature to 723–773 K promoted the formation of larger, faceted LiAl 5 O 8 precipitates, and surface amorphization, accompanied by pronounced lithium depletion and H–D isotopic exchange. At 773 K, medium and high dose rates produced an amorphized surface layer over a crystalline subsurface containing LiAl 5 O 8 precipitates and blisters at the crystalline–amorphous interface, whereas low-dose-rate irradiation preserved surface crystallinity with cavities distributed in the matrix, around precipitates, and along grain boundaries. Precipitate morphology was anisotropic with limited size dependence on the dose rate. These results elucidate the coupled effects of temperature and the dose rate and demonstrate that sequential He + and D 2 + irradiation at 773 K reproduces key microstructural features and H isotope behavior observed in neutron-irradiated γ-LiAlO 2 at 573 K.

dose rate effects↗

NASA NTRSDate not supplied

Detached Melt and Vapor Growth of InI in SUBSA hardware

Indium iodide (InI) single crystals, grown by the Bridgman process, have shown a significant promise as room temperature detector material for X-rays and γ-rays due to the large, close to ideal, bandgap of Eg = 2 eV, and high electrical resistivity. Yet the transport properties of the charge and the detector performance remain inadequate. The Bridgman growth of semiconductor crystals in microgravity have typically resulted in partial dewetting and Marangoni convection, which disturbed the diffusion-controlled growth and produced nonuniform dopant distribution. However, unconfined or partially confined solidification drastically reduced the number of stress-induced defects compared to the confined growth in crucibles on Earth.The goal of the present investigation was to explore melt growth of a heavy metal halide in microgravity, focusing on dewetting, in order to improve the crystalline perfection, i.e., to reduce the number of stress-related defects. InI was chosen because it is non-toxic and has a low melting point of only 365 C (~150 C below the melting point of InSb, used in the SUBSA experiments in 2002). Furthermore, InI melts and evaporates congruently, making it suitable for both growth by directional solidification and growth from the vapor phase. In preparation for the flight experiments, ground-based experiments were conducted in the SUBSA ground unit, consisting of 4 melt growth and 2 vapor growth experiments. The crystals were grown by a gradient freeze method by lowering the setpoint of the furnace. The transparent section of the SUBSA furnace did not allow observations of the solid-liquid interface. Contrary to our expectations, the possible occurrence of dewetting during growth could not be observed. We will present i) the ampoule design and preparation, ii) the experiments conducted in the SUBSA ground unit, iii) the images obtained during the microgravity experiments and iv) the results of the characterization of the grown crystals and detector fabrication.

Crystal Growth↗

DOE OSTI2021

Zinc binding alters the conformational dynamics and drives the transport cycle of the cation diffusion facilitator YiiP

YiiP is a secondary transporter that couples Zn 2+ transport to the proton motive force. Structural studies of YiiP from prokaryotes and Znt8 from humans have revealed three different Zn 2+ sites and a conserved homodimeric architecture. These structures define the inward-facing and outward-facing states that characterize the archetypal alternating access mechanism of transport. To study the effects of Zn 2+ binding on the conformational transition, we use cryo-EM together with molecular dynamics simulation to compare structures of YiiP from Shewanella oneidensis in the presence and absence of Zn 2+ . To enable single-particle cryo-EM, we used a phage-display library to develop a Fab antibody fragment with high affinity for YiiP, thus producing a YiiP/Fab complex. To perform MD simulations, we developed a nonbonded dummy model for Zn 2+ and validated its performance with known Zn 2+ -binding proteins. Using these tools, we find that, in the presence of Zn 2+ , YiiP adopts an inward-facing conformation consistent with that previously seen in tubular crystals. After removal of Zn 2+ with high-affinity chelators, YiiP exhibits enhanced flexibility and adopts a novel conformation that appears to be intermediate between inward-facing and outward-facing states. This conformation involves closure of a hydrophobic gate that has been postulated to control access to the primary transport site. Comparison of several independent cryo-EM maps suggests that the transition from the inward-facing state is controlled by occupancy of a secondary Zn 2+ site at the cytoplasmic membrane interface. This work enhances our understanding of individual Zn 2+ binding sites and their role in the conformational dynamics that govern the transport cycle.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2022

Development and Applications of an eReaxFF Force Field for Graphitic Anodes of Lithium-Ion Batteries

Graphene is one of the most promising materials for lithium-ion battery anodes due to its superior electronic conductivity, high surface area for lithium intercalation, fast ionic diffusivity and enhanced specific capacity. A reliable description of many battery processes requires an explicit description of electrochemical interactions involving electrons. A detailed atomistic modeling of electronic conduction and non-zero voltage simulations of graphitic materials require the inclusion of an explicit electronic degree of freedom. To enable large length- and time-scale simulations of electron conduction in graphitic anodes, we developed an eReaxFF force field concept describing graphitic materials with an explicit electron. The newly developed force field, verified against quantum chemistry-based data describing, amongst others, electron affinities and equation of states, reproduces the qualitative behavior of electron conductivity in pristine and imperfect graphitic materials at different applied temperatures and voltages. In addition, excess electron localization near a defect site estimated from eReaxFF simulations agree quite well with the corresponding density functional theory calculations. Here, our eReaxFF simulations show the initiation of lithium-metal-plating driven by electron transfer from the graphene surface to the exposed lithium ions demonstrating the method’s potential for studying lithium-graphene interactions with explicit electrons and explain many unresolved electrode and electrode-electrolyte interface processes.

25 ENERGY STORAGE↗

DOE OSTI2021

Assessment of Cladding-Cladding Delamination Failures in U-Mo Fuel Plates

Monolithic fuel is a candidate fuel form for the conversion of high-performance research reactors. This plate-type fuel consists of a high-density, U-Mo alloy-based fuel in a monolithic form that is sandwiched between zirconium diffusion barriers and encapsulated in an aluminum cladding. To qualify this fuel system, the program must demonstrate that the fuel meets the safety standards, has a predictable behavior, and performs well in a reactor. Although many plates have shown satisfactory thermo-mechanical performance in reactor, a small number of experimental plates exhibited failures, rendering performance of these plates unacceptable. Post-irradiation examinations of failed plates have revealed that various failures modes are present, including pillowing, blistering, buckling, warping, delamination, fracture, rupture, oxide spallation, etc. Previous studies have implied that such various operational and loading conditions could generate a complex stress-strain state, facilitating distinct failure modes that could impact the structural integrity. To understand such various failure modes, a comprehensive set of studies focusing on failure modes is being performed. As a part of these studies, this work studies the cladding-cladding separation failures that were observed in some of previous irradiation tests. To understand the root cause of cladding separation failures, selected plates from previous irradiation experiments were simulated under comparable operational conditions. Stress-strain state of the plates at the interface were examined. The results were comparatively evaluated to understand possible effects of stress state of the bond-line on delamination failures. This paper presents these findings, discusses such failure modes, and examines the influence of various factors on previously observed cladding delamination failures.

42 ENGINEERING↗

DOE OSTI2023

Modeling laser interactions with aluminum and tantalum targets using a hybrid atomistic-continuum model

A hybrid atomistic-continuum method can model the microstructure evolution of metals subjected to laser irradiation. This method combines classical molecular dynamics (MD) simulations with the two-temperature model (TTM) to account for the laser energy absorption and heat diffusion behavior. Accurate prediction of the temperature evolution in the combined MD-TTM method requires reliable accuracy in electron heat capacity, electron thermal conductivity, and electron–phonon coupling factor across the temperatures generated. This study uses the electronic density of states (DOS) obtained from first-principle calculations. The calculated electron temperature-dependent parameters are used in MD-TTM simulations to study the laser metal interactions in FCC and BCC metals and the phenomenon of laser shock loading and melting. This study uses FCC Al and BCC Ta as model systems to demonstrate this capability. When subjected to short pulsed laser shocks, the dynamic failure behavior predicted using temperature-dependent parameters is compared with the experimentally reported single-crystal and nanocrystalline Al and Ta systems. The MD-TTM simulations also investigate laser ablation and melting behavior of Ta to compare with the ablation threshold reported experimentally. In conclusion, this manuscript demonstrates that integrating the temperature-dependent parameters into MD-TTM simulations leads to the accurate modeling of the laser–metal interaction and allows the prediction of the kinetics of the solid–liquid interface.

36 MATERIALS SCIENCE↗

NASA NTRS1991

Creep and fracture of dispersion-strengthened materials

The creep and fracture of dispersion strengthened materials is reviewed. A compilation of creep data on several alloys showed that the reported values of the stress exponent for creep varied between 3.5 and 100. The activation energy for creep exceeded that for lattice self diffusion in the matrix in the case of some materials and a threshold stress behavior was generally reported in these instances. The threshold stress is shown to be dependent on the interparticle spacing and it is significantly affected by the initial microstructure. The effect of particle size and the nature of the dispersoid on the threshold stress is not well understood at the present time. In general, most studies indicate that the microstructure after creep is similar to that before testing and very few dislocations are usually observed. It is shown that the stress acting on a dispersoid due to a rapidly moving dislocation can exceed the particle yield strength of the G sub p/1000, where G sub p is the shear modulus of the dispersoid. The case when the particle deforms is examined and it is suggested that the dislocation creep threshold stress of the alloy is equal to the yield strength of the dispersoid under these conditions. These results indicate that the possibility that the dispersoid creep threshold stress is determined by either the particle yield strength or the stress required to detach a dislocation from the dispersoid matrix interface. The conditions under which the threshold stress is influenced by one or the other mechanism are discussed and it is shown that the particle yield strength is important until the extent of dislocation core relaxation at the dispersoid matrix interface exceeds about 25 pct. depending on the nature of the particle matrix combination. Finally, the effect of grain boundaries and grain morphology on the creep and fracture behavior of dispersoid strengthened alloys is examined.

Raj, Sai V.↗

DOE OSTI2024

Toward a high-fidelity tritium transport modeling for retention and permeation experiments

We report Tritium Migration Analysis Program version 8 (TMAP8), the latest version of TMAP, was developed within the framework of the Multiphysics Object-Oriented Simulation Environment (MOOSE). Created at Idaho National Laboratory (INL), MOOSE is an open-source, dimension-agnostic, fully coupled, and fully implicit multiphysics platform featuring massively parallel computation capabilities. Using TMAP8, tritium transport in a divertor monoblock was analyzed to elucidate the effects of pulsed operation (up to fifty 1,600 s plasma discharge and cool-down cycles) on the tritium in-vessel inventory source term and ex-vessel release term (i.e., tritium retention and permeation) for safety analysis. With its built-in Message Passing Interface capability, TMAP8 can, in under 2 h, simulate tritium transport in three different layered materials (i.e., tungsten, copper, and copper-chromium-zirconium alloy) in 2D geometry, using a single device/computer with 10 cores. The MOOSE-based TMAP8 code can leverage other MOOSE tools developed under the Nuclear Energy Advanced Modeling and Simulation program to perform tritium and thermal transport in complex geometries and multiphysics environments. And via its massively parallel computation, MOOSE will enable the fusion pilot plant designers to conduct high-fidelity multiphysics modeling for the design of the divertor and blanket systems as well as for the safety analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS2003

Particle Engulfment and Pushing By Solidifying Interfaces

The study of particle behavior at solid/liquid interfaces (SLI s) is at the center of the Particle Engulfment and Pushing (PEP) research program. Interactions of particles with SLI s have been of interest since the 1960 s, starting with geological observations, i.e., frost heaving. Ever since, this field of research has become significant to such diverse areas as metal matrix composite materials, fabrication of superconductors, and inclusion control in steels. The PEP research effort is geared towards understanding the fundamental physics of the interaction between particles and a planar SLI. Experimental work including 1-g and mu-g experiments accompany the development of analytical and numerical models. The experimental work comprised of substantial groundwork with aluminum (Al) and zinc (Zn) matrices containing spherical zirconia particles, mu-g experiments with metallic Al matrices and the use of transparent organic metal-analogue materials. The modeling efforts have grown from the initial steady-state analytical model to dynamic models, accounting for the initial acceleration of a particle at rest by an advancing SLI. To gain a more comprehensive understanding, numerical models were developed to account for the influence of the thermal and solutal field. Current efforts are geared towards coupling the diffusive 2-D front tracking model with a fluid flow model to account for differences in the physics of interaction between 1-g and -g environments. A significant amount of this theoretical investigation has been and is being performed by co-investigators at NASA MSFC.

Source record↗

DOE OSTI2021

Evaluating the stability and activity of dilute Cu-based alloys for electrochemical CO 2 reduction

Cu-based catalysts currently offer the most promising route to actively and selectively produce value-added chemicals via electrochemical reduction of CO 2 (eCO 2 R); yet further improvements are required for their wide-scale deployment in carbon mitigation efforts. Here, we systematically investigate a family of dilute Cu-based alloys to explore their viability as active and selective catalysts for eCO 2 R through a combined theoretical–experimental approach. Using a quantum–classical modeling approach that accounts for dynamic solvation effects, we assess the stability and activity of model single-atom catalysts under eCO 2 R conditions. Our calculations identify that the presence of eCO2R intermediates, such as CO*, H*, and OH*, may dynamically influence the local catalyst surface composition. Additionally, we identify through binding energy descriptors of the CO*, CHO*, and OCCO* dimer intermediates that certain elements, such as group 13 elements (B, Al, and Ga), enhance the selectivity of C 2+ species relative to pure Cu by facilitating CO dimerization. The theoretical work is corroborated by preliminary testing of eCO 2 R activity and selectivity of candidate dilute Cu-based alloy catalyst films prepared by electron beam evaporation in a zero-gap gas diffusion electrode-based reactor. Of all studied alloys, dilute CuAl was found to be the most active and selective toward C 2+ products like ethylene, consistent with the theoretical predictions. Finally, we attribute the improved performance of dilute CuAl alloys to more favorable dimerization reaction energetics of bound CO species relative to that on pure Cu. In a broader context, the results presented here demonstrate the power of our simulation framework in terms of rational catalyst design.

30 DIRECT ENERGY CONVERSION↗

DOE OSTI2022

Understanding the Impacts of Support–Polymer Interactions on the Dynamics of Poly(ethyleneimine) Confined in Mesoporous SBA-15

Supported amines are a promising class of CO 2 sorbents offering large uptake capacities and fast uptake rates. Additionally, among supported amines, poly(ethyleneimine) (PEI) physically impregnated in the mesopores of SBA-15 silica is widely used. Within these composite materials, the chain dynamics and morphologies of PEI strongly influence the CO 2 capture performance, yet little is known about chain and macromolecule mobility in confined pores. Here, we probe the impact of the support–PEI interactions on the dynamics and structures of PEI at the support interface and the corresponding impact on CO 2 uptake performance, which yields critical structure–property relationships. The pore walls of the support are grafted with organosilanes with different chemical end groups to differentiate interaction modes (spanning from strong attraction to repulsion) between the pore surface and PEI. Combinations of techniques, such as quasi-elastic neutron scattering (QENS), 1 H T 1 –T 2 relaxation correlation solid-state NMR, and molecular dynamics (MD) simulations, are used to comprehensively assess the physical properties of confined PEI. We hypothesized that PEI would have faster dynamics when subjected to less attractive or repulsive interactions. However, we discover that complex interfacial interactions resulted in complex structure–property relationships. Indeed, both the chain conformation of the surface-grafted chains and of the PEI around the surface influenced the chain mobility and CO 2 uptake performance. By coupling knowledge of the dynamics and distributions of PEI with CO 2 sorption performance and other characteristics, we determine that the macroscopic structures of the hybrid materials dictate the first rapid CO 2 uptake, and the rate of CO 2 sorption during the subsequent gradual uptake stage is determined by PEI chain motions that promote diffusive jumps of CO 2 through PEI-packed domains.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Leveraging Artificial Intelligence to Expedite Antibody Design and Enhance Antibody–Antigen Interactions

This perspective sheds light on the transformative impact of recent computational advancements in the field of protein therapeutics, with a particular focus on the design and development of antibodies. Cutting-edge computational methods have revolutionized our understanding of protein–protein interactions (PPIs), enhancing the efficacy of protein therapeutics in preclinical and clinical settings. Central to these advancements is the application of machine learning and deep learning, which offers unprecedented insights into the intricate mechanisms of PPIs and facilitates precise control over protein functions. Despite these advancements, the complex structural nuances of antibodies pose ongoing challenges in their design and optimization. Our review provides a comprehensive exploration of the latest deep learning approaches, including language models and diffusion techniques, and their role in surmounting these challenges. We also present a critical analysis of these methods, offering insights to drive further progress in this rapidly evolving field. The paper includes practical recommendations for the application of these computational techniques, supplemented with independent benchmark studies. These studies focus on key performance metrics such as accuracy and the ease of program execution, providing a valuable resource for researchers engaged in antibody design and development. Through this detailed perspective, we aim to contribute to the advancement of antibody design, equipping researchers with the tools and knowledge to navigate the complexities of this field.

59 BASIC BIOLOGICAL SCIENCES↗

NASA NTRS2016

Evaluation of Methods for Molecular Dynamics Simulation of Ionic Liquid Electric Double Layers

We investigate how systematically increasing the accuracy of various molecular dynamics modeling techniques influences the structure and capacitance of ionic liquid electric double layers (EDLs). The techniques probed concern long-range electrostatic interactions, electrode charging (constant charge versus constant potential conditions), and electrolyte polarizability. Our simulations are performed on a quasi-two-dimensional, or slab-like, model capacitor, which is composed of a polarizable ionic liquid electrolyte, [EMIM][BF4], interfaced between two graphite electrodes. To ensure an accurate representation of EDL differential capacitance, we derive new fluctuation formulas that resolve the differential capacitance as a function of electrode charge or electrode potential. The magnitude of differential capacitance shows sensitivity to different long-range electrostatic summation techniques, while the shape of differential capacitance is affected by charging technique and the polarizability of the electrolyte. For long-range summation techniques, errors in magnitude can be mitigated by employing two-dimensional or corrected three dimensional electrostatic summations, which lead to electric fields that conform to those of a classical electrostatic parallel plate capacitor. With respect to charging, the changes in shape are a result of ions in the Stern layer (i.e. ions at the electrode surface) having a higher electrostatic affinity to constant potential electrodes than to constant charge electrodes. For electrolyte polarizability, shape changes originate from induced dipoles that soften the interaction of Stern layer ions with the electrode. The softening is traced to ion correlations vertical to the electrode surface that induce dipoles that oppose double layer formation. In general, our analysis indicates an accuracy dependent differential capacitance profile that transitions from the characteristic camel shape with coarser representations to a more diffuse profile with finer representations.

Haskins, Justin B.↗

DOE OSTI2025

Resolution of Selectivity Steps of CO Reduction Reaction on Copper by Quantum Monte Carlo

Electrochemical reduction of carbon monoxide to valuable fuels and chemicals on copper surfaces remains a challenging area in catalysis due to a limited understanding of adsorption mechanisms and reaction pathways. Although density functional theory (DFT)-based studies have investigated these processes, their accuracy varies across different functionals. Here, in this study, we present the application of fixed-node diffusion Monte Carlo (FNDMC) to benchmark the adsorption energies of CO*, H*, and key CO reduction reaction (CORR) intermediates, COH* and CHO* on the Cu(111) surface. Our results for CO* and H* adsorption energies closely align with experimentally measured chemisorption reactions, highlighting the limitations of DFT and providing site-specific energy comparisons that are often not available experimentally. Additionally, we explore the effect of explicit solvation, demonstrating how water stabilizes the COH* over CHO*, thus suggesting a critical role of COH* in CORR. Finally, we release our high-accuracy FNDMC benchmarks for testing and developing new DFT functionals for electrocatalysis. Overall, this study underscores the potential of FNDMC for detailed surface chemistry studies and offers new insights into catalytic processes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗