Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Step”

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 757 records · Page 42
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

High resolution, high conductivity components by aerosol jet printing

Highly conductive electrical traces formed over mechanical steps or on non-planar surfaces with linewidths of 10 to 100 μm and a method for forming such electrical traces are disclosed. The method employs two steps, with the first step using an aerosol jet printing (AJP) process to form thin electrical traces that serve as the seed layers for the second step. The first step preferably employs multiple passes with the AJP to create multiple seed sub-layers with improved continuity and conductivity. In the second step, the seed layers are subjected to an electrodeposition process that forms the bulk of the thickness of the electrical traces. The electrodeposition process may include one, two, or three phases at corresponding low or high biases, with low biases providing denser, more highly conductive plating sub-layers, while high biases provide plating sub-layers having better gap bridging properties.

Lavin, Judith Maria↗
DOE OSTI2020

Distribution Development for Residual Inventory at the New York West Valley Site - 20400

The New York State Energy Research and Development Authority (NYSERDA) is the owner of the Western New York Nuclear Service Center (WNYNSC), a 1,351 ha site located approximately 48 km south of Buffalo, New York. In 1962, Nuclear Fuel Services, Inc. (NFS) entered into Agreements with the Atomic Energy Commission and New York State to construct the first commercial reprocessing plant of nuclear fuel in the United States. NFS, a private company, built and operated the spent fuel reprocessing plant and waste disposal facilities, processing 640 Mg of spent nuclear fuel from 1966 to 1972 under an Atomic Energy Commission license. Nuclear fuel reprocessing operations ended in 1972 and never reopened, leaving behind radioactive and chemical wastes. Operations led to contamination in a number of facilities and locations. Some of that contamination has migrated from waste disposal zones to other layers, formations, and features on and off the WNYNSC. Phase I decommissioning activities are ongoing and involve the removal of a number of areas and structures that have been associated with contamination. The purpose of this work is to outline the approach for characterizing contamination not associated with disposed wastes, contaminated structures, or specific releases. In this work, the term, residual radiological activity, is used to describe environmental contamination that exists subsequent to the completion of Phase I decommissioning activities, that is not associated with disposed wastes, contaminated structures, or specific releases. Contamination from the Site was quantified relative to data that characterize the concentrations of radionuclides that exist in background. Background concentrations are those present in the area but having no influence from Site related activities. The existence of residual radiological activity that is elevated relative to background has the potential to contribute to future risks to human health and the environment. As a consequence, the residual inventory information is used to inform the West Valley Probabilistic Performance Assessment (PPA) model to characterize potential future risks to human health and the environment. The centralized West Valley Data Management System (DMS) was the source of information for the data assembled in this analysis. The DMS is a fairly large compilation consisting of thousands of records from investigation studies, with sample dates ranging from 1990 to present. Samples from monitoring wells, boreholes, geoprobe studies, surface water, surface soils, storm water outfalls, ventilation stack filters, plant and animal tissues, and more are included in the DMS. Results are typically reported in units of activity per unit volume. For the purpose of the analyses presented here, all results were converted into consistent units of pCi per unit volume. Since 1990, data have been collected from various locations across the WNYNSC at different times with varying frequency over the course of several decades. As a consequence, a number of potential issues can arise with respect to the assembly of a dataset that is deemed adequate for the characterization of residual radiological activity. These issues were assessed and resolved to the extent possible through careful consideration of the properties of the distributions. The intent was to use data which characterize the current state of the Site. Radionuclides can be designated to one of several groups depending on their origin. In this work the groups considered were 1) Naturally Occurring Radioactive Material (NORM), 2) fallout, and 3) Other (including power plant, medical research, etc). This grouping is a useful construct with respect to the interpretation of fixed laboratory results. For example, NORM radionuclides that exist within a decay chain should have approximately equivalent distributions of concentrations if they are representative of background conditions. Insights such as these can be used as a check to identify sample results that need to be further investigated or omitted due to issues associated with reported results from fixed laboratory analyses. This type of analysis provided a foundation for the assessment of the adequacy of sample results for use in subsequent components of an assessment. The general process for the assessment of residual radiological contamination at the Site consists of a sequence of several steps. First, for each analyte, several statistical tests were performed to assess the weight of evidence against the null hypothesis that the mean of the distribution of concentrations was equal to zero. If the mean of the distribution of concentrations for a given radionuclide was not found to be greater than zero, then it was removed from consideration as a component of the residual radiological contamination. If there was significant evidence to reject the hypothesis of the mean being equal to zero, the second step was to compare the distribution of the data from the Site to that of the corresponding background. A suite of tests was used to compare the distributions of the site and background data. The results of these tests were collectively used to determine if site data are elevated relative to background. The third step was to develop distributions using a Bayesian framework to characterize the distribution of mean of the increment present above background for each of the radionuclides. The Bayesian model implemented allowed for the comparison of site-specific records to background concentrations to better approximate contamination attributed to the Site. A final screening step was employed for radionuclides that exceed background. This screening step compared 95% upper confidence limits (UCLs) from the increment distribution developed in the previous step to the risk screening levels. This approach yields a list of analytes that were determined to be elevated relative to background.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

Cloud Process Coupling and Time Integration in the E3SM Atmosphere Model

Abstract In this study, we find significant sensitivity to the choice of time step for the Energy Exascale Earth System Model's atmospheric component, leading to large decreases in the magnitude of cloud forcing when the time step is reduced to 10 s. Reducing the time step size for the microphysics increases precipitation, leading to a drying of the atmosphere and an increase in surface evaporation. This effect is amplified when the microphysics is substepped together with other cloud physics processes. Coupling the model's dynamics and physics more frequently reduces cloud fraction at lower altitudes, while producing more cloud liquid at higher altitudes. Reducing the deep convection time step also reduces low cloud mass and cloud fraction. Together, these results suggest that cloud physics in a global circulation model can depend strongly on time step and, in particular, on the frequency with which cloud‐related processes are coupled with each other and with the model dynamics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Mineral Deposition on the Rough Walls of a Fracture

Modeling carbonate growth in fractures and pores is important for understanding carbon sequestration in the environment or when supersaturated solutions are injected into rocks. Here, we study the simple but nontrivial problem of calcite growth on fractures with rough walls of the same mineral using kinetic Monte Carlo simulations of attachment and detachment of molecules and scaling approaches. First, we consider wedge-shaped fracture walls whose upper terraces are in the same low-energy planes and show that the valleys are slowly filled by the propagation of parallel monolayer steps in the wedge sides. The growth ceases when the walls reach these low-energy configurations so that a gap between the walls may not be filled. Second, we consider fracture walls with equally separated monolayer steps (vicinal surfaces with roughness below 1 nm) and show that growth by step propagation will eventually clog the fracture gap. In both cases, scaling approaches predict the times to attain the final configurations as a function of the initial geometry and the step-propagation velocity, which is set by the saturation index. The same reasoning applied to a random wall geometry shows that step propagation leads to lateral filling of surface valleys until the wall reaches the low-energy crystalline plane that has the smallest initial density of molecules. Thus, the final configurations of the fracture walls are much more sensitive to the crystallography than to the roughness or the local curvature. The framework developed here may be used to determine those configurations, the times to reach them, and the mass of deposited mineral. Effects of transport limitations are discussed when the fracture gap is significantly narrowed.

calcite↗
DOE OSTI2022

Isolation, Characterization, and Depolymerization of L -Cysteine Substituted Eucalyptus Lignin

Lignin condensation reactions are hard to avoid or control during separation, which is a deterrent to lignin isolation and post-conversation, especially for the full utilization of lignocelluloses. Selective protection of β-aryl ether linkages in the isolation process is crucial to lignin valorization. Herein, a two-step acid/alkali separation method assisted with l-cysteine for eucalyptus lignin separation is developed, and the isolated L -cysteine lignins (LCLs) are comprehensively characterized by 2D NMR, 31 P NMR, thioacidolysis, etc. Compared to the two-step control treatment, a much higher β-O-4 content is preserved without reducing the separation efficiency assisted by L -cysteine, which is also significantly higher than alkali lignin and kraft lignin. The results of hydrogenolysis show that LCLs generate a much higher monomer yield than that of control sample. Structural analysis of LCLs suggests that lignin condensation reaction, to some extent, is suppressed by adding L -cysteine during the two-step acid/alkali separation. Further, mechanistic studies using dimeric model compound reveals that L -cysteine may be the α-carbon protective agent in the two-step separation. The role of L -cysteine in the two-step lignin isolation method provides novel insights to the selective fractionation of lignin from biomass, especially for the full valorization of lignocellulosic biomass.

09 BIOMASS FUELS↗
DOE OSTI2021

On the Structural Transformation of Ni/BaH 2 During a N 2 -H 2 Chemical Looping Process for Ammonia Synthesis: A Joint In Situ Inelastic Neutron Scattering and First-Principles Simulation Study

The demand for decarbonizing the ammonia industry by using renewable energy has invoked increasing research interests into catalyst development for effective N 2 reduction under mild conditions. Hydride-based materials are among some of the emerging catalysts for ammonia synthesis at ambient pressure and low temperatures (< 673 K). A recent chemical looping process based on Ni/BaH 2 showed the most promise as it can realize ammonia production at a temperature as low as 373 K and under ambient pressure. However, the chemical transformation of the hydride catalyst at the molecular level remains unclear in this process. Here, we report detailed in situ neutron spectroscopy and diffraction investigations along with first-principles simulations on the structural transformation of Ni/BaH 2 during the nitridation and hydrogenation steps in the chemical looping process for ammonia synthesis. It was shown that a ball-milling process of the starting Ni/BaH 2 could significantly decrease the size of BaH 2 and increase the density of defects, thus potentially enhancing the reactivity of the hydride. The evolution from BaH 2 to barium imide (BaNH) was evidenced in the inelastic neutron scattering (INS) and neutron diffraction results during the N 2 reaction step. During the hydrogenation study, in addition to the recovery of BaH 2 , a possible intermediate species, N-deficient barium imide, was also detected. In comparing the N 2 and H 2 reaction steps, the neutron results indicate that the hydrogenation step appears more difficult than the nitridation step, confirming the facile N 2 fixation property of Ni/BaH 2 catalyst in ammonia synthesis.

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

Two pathways to understanding electron transfer in reaction centers from photosynthetic bacteria: A comparison of Rhodobacter sphaeroides and Rhodobacter capsulatus mutants

The rates, yields, mechanisms and directionality of electron transfer (ET) are explored in twelve pairs of Rhodobacter (R.) sphaeroides and R. capsulatus mutant RCs designed to defeat ET from the excited primary donor (P*) to the A-side cofactors and re-direct ET to the normally inactive mirror-image B-side cofactors. In general, the R. sphaeroides variants have larger P + H B - yields (up to similar to ~90%) than their R. capsulatus analogs (up to similar to ~60%), where H B is the B-side bacteriopheophytin. Substitution of Tyr for Phe at L-polypeptide position L181 near B B primarily increases the contribution of fast P* → P + B B - → P + H B - two-step ET, where B-B is the "bridging" B-side bacteriochlorophyll. The second step (similar to ~6-8 ps) is slower than the first (similar to ~3-4 ps), unlike A-side two-step ET (P* → P + B A - → P + H A - ) where the second step (similar to ~1 ps) is faster than the first (similar to ~3-4 ps) in the native RC. Substitutions near H B , at L185 (Leu, Trp or Arg) and at M-polypeptide site M133/131 (Thr, Val or Glu), strongly affect the contribution of slower (20-50 ps) P* → P + H B - one-step superexchange ET. Here, both ET mechanisms are effective in directing electrons "the wrong way" to H B and both compete with internal conversion of P* to the ground state (similar to ~200 ps) and ET to the A-side cofactors. Collectively, the work demonstrates cooperative amino-acid control of rates, yields and mechanisms of ET in bacterial RCs and how A- vs. B-side charge separation can be tuned in both species.

Charge Separation↗
DOE OSTI2026

Advancing Ethanol-to-Jet cost Effectiveness via direct conversion to n -Butene-Rich olefins and Co-Product Valorization

Ethanol is a promising feedstock for sustainable aviation fuel production; however, conventional routes face significant energy and cost challenges, particularly due to the ethanol dehydration step to ethylene. Here, this study leverages breakthrough experimental data to perform comprehensive techno-economic and life-cycle assessments of an innovative ethanol-to-jet process. The process employs a single-step catalytic conversion, enabled by multifunctional Cu-ZrO 2 /SBA-16 catalyst, to directly upgrade ethanol into a mixed olefin stream rich in n-butene. The single-step conversion eliminates the costly ethanol dehydration step in the conventional process. High selectivity toward n-butene offers key advantages: it simplifies downstream oligomerization into jet-range hydrocarbons and enables the co-production of renewable n-butene alongside sustainable aviation fuel. The analysis estimates a minimum fuel selling price as low as $\$$2.50 per gallon, whether using corn ethanol or cellulosic ethanol from corn stover. Life cycle CO 2 equivalent emissions are projected to be as low as 10.6 g CO 2 eq/MJ sustainable aviation fuel, representing over 70% reduction compared to conventional petroleum-based jet fuel. This one-step ethanol upgrading approach not only facilitates SAF and n-butene co-production but also provides operational flexibility. The ability to tailor product outputs allows the ethanol-to-jet process to adapt to varying feedstocks, incentive programs, and market dynamics, ultimately enhancing the economic viability of sustainable aviation fuel production.

Xu, Yiling [Pacific Northwest National Laboratory ↗
DOE OSTI2022

Krylov subspace recycling for evolving structures

Krylov subspace recycling is a powerful tool when solving a long series of large, sparse linear systems that change only slowly over time. In PDE constrained shape optimization, these series appear naturally, as typically hundreds or thousands of optimization steps are needed with only small changes in the geometry. In this setting, however, applying Krylov subspace recycling can be a difficult task. As the geometry evolves, in general, so does the finite element mesh defined on or representing this geometry, including the numbers of nodes and elements and element connectivity. This is especially the case if re-meshing techniques are used. As a result, the number of algebraic degrees of freedom in the system changes, and in general the linear system matrices resulting from the finite element discretization change size from one optimization step to the next. Changes in the mesh connectivity also lead to structural changes in the matrices. In the case of re-meshing, even if the geometry changes only a little, the corresponding mesh might differ substantially from the previous one. Obviously, this prevents any straightforward mapping of the approximate invariant subspace of the linear system matrix (the focus of recycling in this work) from one optimization step to the next; similar problems arise for other selected subspaces. In this paper, we present an algorithm to map an approximate invariant subspace of the linear system matrix for the previous optimization step to an approximate invariant subspace of the linear system matrix for the current optimization step, for general meshes. This is achieved by exploiting the map from coefficient vectors to finite element functions on the mesh, combined with interpolation or approximation of functions on the finite element mesh. We demonstrate the effectiveness of our approach numerically with several proof of concept studies for a specific meshing technique.

42 ENGINEERING↗
DOE OSTI2024

Synthesis of results for Brine Availability Test in Salt (BATS) DECOVALEX-2023 Task E

In the 2023 phase of the international collaborative DECOVALEX modeling project, Task E focused on understanding thermal, hydrological, and mechanical (THM) processes related to predicting brine migration in the excavation damaged zone around a heated excavation in salt. Salt is attractive as a disposal medium for radioactive waste because it is self-healing and is essentially impermeable and non-porous in the far field. Investigation of the short-term, near-field behavior is important for radioactive waste disposal because this early period strongly controls the amount of inflowing brine. Brine leads to corrosion of waste forms and waste packages, and possible dissolution of radionuclides with brine transport being a potential transport vector to the accessible environment.The Task was divided into steps. Step 0 included matching unheated brine inflow data from boreholes at the Waste Isolation Pilot Plant (WIPP) and matching temperature observations during a Brine Availability Test in Salt (BATS) heater test. Step 1 included validation of models against a thermo-poroelastic analytical solution, and two-phase flow around an excavation. Finally, Step 2 required all the individual components covered in steps 0 and 1 to come together to match observed brine inflow behavior during the same BATS heater test.There were a range of approaches from the teams, from mechanistic to prescriptive. Given the uncertainties in the problem, some teams used one- or two-dimensional models of the processes, while other teams included more geometrical complexity in three-dimensional models. Task E was a learning experience for the teams involved, and feedback from the modeling teams has led to changes in follow-on BATS experiments at WIPP. The primary Task E lessons learned were the impact of hydrologic initialization methods (wetting up vs. drying down), the difference between confined and unconfined thermal expansion, and the large changes in permeability associated with heating and cooling.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

Asymptotic-preserving gyrokinetic implicit particle-orbit integrator for arbitrary electromagnetic fields

We extend the asymptotic preserving and energy conserving time integrator for charged-particle motion developed in Ricketson and Chacón (2020) to include finite Larmor-radius (FLR) effects in the presence of electric-field length-scales comparable to the particle gyro-radius (the gyro-kinetic limit). We introduce two modifications to the earlier scheme. The first is the explicit gyro-averaging of the electric field at the half time-step, along with an analogous modification to the current deposition, which we show preserves total energy conservation in implicit PIC schemes. The number of gyrophase samples is chosen adaptively, ensuring proper averaging for large timesteps and the recovery of full-orbit dynamics in the small time-step limit. The second modification is an alternating large and small time-step strategy that ensures the particle trajectory samples gyrophases evenly. We show that this strategy relaxes the time-step restrictions on the scheme, allowing even larger speed-ups than previously achievable. We demonstrate the new method with several single-particle motion tests in a variety of electromagnetic field configurations featuring gyro-scale variation in the electric field. Finally, the results demonstrate the advertised ability to capture FLR effects accurately even when significantly stepping over the gyration time-scale.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Origin of strain softening in a nanograined Al alloy

Here, one-step and multi-step nanoindentation experiments were conducted to investigate the mechanical response of an Al-Mg-Li alloy. Opposite to the strain hardening effect observed in the coarse-grained (CG) alloy, the nanograined (NG) alloy exhibits strain softening. Compared with multi-step nanoindentation experiments, the ultralow dislocation density and the grain boundary (GB) segregation promote a higher stress level during one-step nanoindentation experiments. Residual perfect dislocations induced by a previous loading and unloading during multi-step nanoindentation experiments lower the stress required for further plastic deformation. In addition, the storage of partial dislocations also serves to render plastic flow to commence at lower stresses. This work provides new insights for the fabrication and property optimization of NG alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Low Temperature Activation of Methane on Metal-Oxides and Complex Interfaces: Insights from Surface Science

The abundance of inexpensive, natural gas has transformed the energy landscape, whereby revealing new possibilities for sustainable chemical technologies or impacting those that have relied on traditional fossil fuels. The primary component, methane, is underutilized and wastefully exhausted, leading to anthropogenic global warming. Historically, the manipulation of methane remained “ clavis aurea ,” an insurmountable yet rewarding challenge and thus the focus of intense research. This is primarily due to an inability to dissociate C–H bonds in methane selectively, which requires a high energy penalty and is an essential prerequisite for the direct conversion of methane into a large set of value-added products. The discovery of such processes would promise an energy gainful use of natural gas benefiting several essential chemical processes associated with C1 chemistry. This first C–H bond dissociation step of the methane molecule appears in numerous catalytic mechanisms as the rate-determining step or most essential barrier sequence for all subsequent steps that follow in the production of C–C, C–O, or C x –H y –O z bonds found in value added products. A main goal is to catalytically reduce the energy barrier for the first C–H bond dissociation to be able to achieve the activation of methane at low or moderate temperatures. As such there is great value in understanding the fundamental nature of the active sites responsible for bond breaking or formation and thus be able to facilitate better control of this chemistry, leading to the development of new technologies for fuel production and chemical conversion. Surface science studies offer enhanced perspectives for a careful manipulation of bonds over the last layer atoms of catalyst surfaces, an essential factor for the design of atomically precise catalysts and unravelling of the reaction mechanism. With the advent of new surface imaging, spectroscopy, and in situ tools, it has been possible to decipher the surface chemistry of complex materials systems and further our understanding of atomic active sites on the surfaces of metals, oxides, and carbides or metal–oxide and metal–carbide interfaces. The once considered near impossible step of C–H bond activation is now observed at low temperatures with high propensity over a collection of oxide, metal–oxide, and metal–carbide systems in a conventional or inverse configuration (oxide or carbide on metal). The enabling of C–H activation at low temperature has opened interesting possibilities for the specific production of chemicals such as methanol directly from methane, a step toward facile synthesis of liquid fuels. We highlight the most recent of these results and present the key aspects of active site configurations engineered from surface science studies which enable such a simple reactive event through careful manipulation of the last surface layer of atoms found in the catalyst structure. New concepts which help in the activation and conversion of methane are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Equipartition and the Temperature of Maximum Density of TIP4P/2005 Water

Here, we simulate TIP4P/2005 water in the temperature range of 257 to 318 K with time-steps δ = 0.25, 0.50, 1.00, 2.00, and 4.00 fs. The density–temperature behavior obtained using 0.25 or 0.50 fs is in excellent agreement with each other but differs from those obtained using time steps that have been shown earlier to lead to a breakdown of equipartition. For δt = 0.25 or 0.50 fs, the temperature of maximum density (TMD) is 277.15 K and the density value is in close agreement with experiments. For δt = 1.00 fs, the TMD is 277.15 K, but the density value is shifted higher. For the other time steps considered here, the TMD is shifted to progressively lower values for longer time steps, a trend that holds for different thermostat/barostat combinations. Enhancing the water–water dispersion interaction, as has been recommended for simulating disordered proteins in TIP4P/2005, degrades the description of the liquid–vapor phase envelope. We present a simple physically transparent explanation that highlights the separation of the time scales between translational and rotational motion. We also develop a metric, χ, that we term the equipartition anomaly, to detect equipartition violations in simulations that include molecules that are treated as rigid objects. Calculating χ is shown to be straightforward and sensitive to equipartition violations. A key takeaway from this study is that using sufficiently short time steps (≤0.5 fs) to preserve equipartition is essential for obtaining meaningful liquid water properties and for producing reliable simulation data, as correct ensemble sampling is fundamental to ensure reproducibility across codes and simulation algorithms.

Asthagiri, Dilipkumar N. [Oak Ridge National Labor↗
DOE OSTI2024

Composition-Dependent Reconstructions of Titanium Oxide Clusters and Cu(111) Support via Cluster-Adatom Interactions

Mass-selected cluster deposition and scanning tunneling microscopy (STM) were used to investigate the surface structure and thermal stability of Ti 3 O y (y = 5, 6) clusters on a Cu(111) surface. STM imaging shows that the substoichiometric Ti 3 O 5 clusters aggregate into large dendritic assemblies at step edges, while the stoichiometric Ti 3 O 6 clusters are highly dispersed as single or small groupings of clusters. Cluster deposition also results in extensive reconstruction of the Cu(111) surface. The Ti 3 O 5 /Cu(111) surface exhibits one Cu-layer deep pits that are ringed by cluster assemblies that nucleate and grow at the perimeter step edges. By comparison, the Ti 3 O 6 /Cu(111) surface consists of irregular step edges and small one Cu-layer islands that are decorated with clusters. The restructuring of the Cu(111) surface is attributed to strong Cu-cluster interactions that release Cu adatoms from the step edges. This mechanism is supported by density functional theory (DFT) calculations, which predict that the Cu adatom attachment to the clusters is energetically favorable. Annealing the surfaces at elevated temperatures (400-700 K) results in higher densities of irregular steps, pits and isolated Cu islands as well as increased cluster aggregation. In conclusion, this work highlights adatom-cluster interactions as another manifestation of strong cluster-support interactions that can significantly alter the nature of the oxide-metal interface.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2020

Insights into the Interaction of Redox Active Organic Molecules and Solvents with the Pristine and Defective Graphene Surfaces from Density Functional Theory

A systematic study is reported on the interaction of two representative redox active organic molecules and two solvent molecules with pristine and defective graphene surfaces as a model of an electrode surface of a redox flow battery (RFB). The redox active molecules include a catholyte, 2,5-di-tert-butyl-1,4-dimethoxybenzene (DDB), and an anolyte molecule, benzothiadiazole (BTZ), and the solvent molecules include acetonitrile (MeCN) and ethylene carbonate (EC). The graphene defects investigated include a single vacancy, double vacancy, zigzag step edge, and armchair step edge. Computations suggest that the interactions of all molecules with a pristine graphene surface are relatively weak (0.2 to 0.8 eV) and dominated by van der Waals effects; therefore, these molecules are chemically stable upon interacting with pristine nondefective graphene. The BTZ, DDB, MeCN, and EC molecules interact strongly (1.5 to 5.5 eV) with the single vacancy and zigzag step edges of graphene that leads to the possible decomposition of the molecules with strength of interaction in the order of MeCN > BTZ > DDB > EC and MeCN > EC > BTZ > DDB, respectively. Calculations show that the BTZ, DDB, MeCN, and EC molecules interact less strongly (0.2 to 1.4 eV) to the double vacancy and armchair step edge than to the single vacancy and zigzag step edge. The binding energies of the molecules were significantly reduced when interacting with the passivated defects, suggesting that the passivation of the defects could help prevent unwanted chemical interactions between the neutral molecules and the electrode surface. In all organic RFBs, one of the crucial bottlenecks is the stability of the constituent molecules with the electrode material, and this study provides insights into the chemical interaction of selected candidate species with a model carbon electrode.

Howard, Jason D.↗
DOE OSTI2024

Unraveling the CO Oxidation Mechanism over Highly Dispersed Pt Single Atom on Anatase TiO 2 (101)

Catalysts with noble metals deposited as single atoms on metal oxide supports have recently been studied extensively due to their maximized metal utilization and potential for performing difficult chemical conversions owing to their unique electronic properties. Understanding of the reaction mechanisms on supported single-metal atoms is still limited but is highly important for designing more efficient catalysts. In this study, we report the complexity of the CO oxidation reaction mechanism on Pt single atoms supported on anatase TiO 2 (Pt SA /a-TiO 2 ) by coupling density functional theory (DFT) calculations and microkinetic analysis with kinetic measurements, in situ/operando infrared, and X-ray absorption spectroscopies. Starting from the adsorbed Pt SA occupying an O vacancy induced by reductive pretreatment, we show that CO oxidation follows a complex mechanism consisting of initiation steps to reorganize the active site and multibranch reactive cycles, with the Pt SA /a-TiO 2 catalyst not returning to its initial configuration. The initiation step consists of CO and O 2 adsorption healing the O vacancy, followed by CO oxidation using gas-phase CO to form Pt(CO). The reactive cycle alternates O 2 adsorption and dissociation to oxidize the catalyst to Pt(O)(O)(CO) and branching pathways of competing Langmuir–Hinshelwood (LH)- or Eley–Rideal (ER)-type CO oxidation steps to reduce it again to Pt(CO). In situ/operando infrared experiments, including cryogenic CO adsorption and isotopic CO exchange, confirm the combined involvement of strongly adsorbed CO and gas-phase CO in an Eley–Rideal step along the reaction cycle. Microkinetic modeling shows that Pt single atoms are present in a mixture of Pt(CO), Pt(CO)(O 2 ), Pt(O)(CO)(O 2 ), and Pt(CO)(CO 3 ) structures as the main intermediates during steady-state CO oxidation, all having the C–O vibrational stretch close to the experimentally observed value of 2115 cm –1 . Microkinetic modeling also shows that the fractional orders of CO and O 2 measured experimentally originate from multiple steps with a high degree of rate control and not from a simple competitive adsorption. The results demonstrate the complex reaction pathways that even CO oxidation on a simple single-atom system can follow, providing mechanistic insights for designing efficient Pt-based single-atom catalysts. Here, we further show that microkinetic modeling results are sensitive to changes in energies of intermediate and transition states within errors of density functional theory, which can ultimately lead to incorrect conclusions regarding the reaction pathways and most abundant reaction intermediates if not accounted for by experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Revealing a Hidden Intermediate of Rotatory Catalysis with X-ray Crystallography and Molecular Simulations

The mechanism of rotatory catalysis in ATPhydrolyzing molecular motors remains an unresolved puzzle in biological energy transfer. Notwithstanding the wealth of available biochemical and structural information inferred from years of experiments, knowledge on how the coupling between the chemical and mechanical steps within motors enforces directional rotatory movements remains fragmentary. Even more contentious is to pinpoint the rate-limiting step of a multistep rotation process. Here, using vacuolar or V 1 -type hexameric ATPase as an exemplary rotational motor, we present a model of the complete 4-step conformational cycle involved in rotatory catalysis. First, using Xray crystallography, a new intermediate or “dwell” is identified, which enables the release of an inorganic phosphate (or P i ) after ATP hydrolysis. Using molecular dynamics simulations, this new dwell is placed in a sequence with three other crystal structures to derive a putative cyclic rotation path. Free-energy simulations are employed to estimate the rate of the hexameric protein transformations and delineate allosteric effects that allow new reactant ATP entry only after hydrolysis product exit. An analysis of transfer entropy brings to light how the side-chain-level interactions transcend into larger-scale reorganizations, highlighting the role of the ubiquitous arginine-finger residues in coupling chemical and mechanical information. An inspection of all known rates encompassing the 4-step rotation mechanism implicates the overcoming of the ADP interactions with V 1 -ATPase to be the rate-limiting step of motor action.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗