Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “multiscale structure”

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

Denoising diffusion algorithm for inverse design of microstructures with fine-tuned nonlinear material properties

Here we introduce a denoising diffusion algorithm to discover microstructures with nonlinear fine-tuned properties. Denoising diffusion probabilistic models are generative models that use diffusion-based dynamics to gradually denoise images and generate realistic synthetic samples. By learning the reverse of a Markov diffusion process, we design an artificial intelligence to efficiently manipulate the topology of microstructures to generate a massive number of prototypes that exhibit constitutive responses sufficiently close to designated nonlinear constitutive behaviors. To identify the subset of microcstructures with sufficiently precise fine-tuned properties, a convolutional neural network surrogate is trained to replace high-fidelity finite element simulations to filter out prototypes outside the admissible range. Results of this study indicate that the denoising diffusion process is capable of creating microstructures of fine-tuned nonlinear material properties within the latent space of the training data. More importantly, this denoising diffusion algorithm can be easily extended to incorporate additional topological and geometric modifications by introducing high-dimensional structures embedded in the latent space. Numerical experiments are conducted on the open-source mechanical MNIST data set (Lejeune, 2020). Consequently, this algorithm is not only capable of performing inverse design of nonlinear effective media, but also learns the nonlinear structure–property map to quantitatively understand the multiscale interplay among the geometry, topology, and their effective macroscopic properties.

42 ENGINEERING↗

Light scattering by V 4 O 7 film across the metal–insulator transition

The experimental study of the angle-resolved hemispherical light scattering by V 4 O 7 film within a broad temperature range across metal–insulator transition reveals complex structural reorganization of the film deposited on the c-cut sapphire crystal. The bidirectional scattering distribution function and the surface autocorrelation function were obtained from scattering data to visualize statistics of the spatially resolved contributions of optical inhomogeneities in normal and lateral directions to the surface. The measurements reveal an anisotropic surface roughness distribution due to the twinned domain structure, with significant anisotropy changes across the phase transition. The V 4 O 7 film deposited on sapphire leads to a polydomain structure, minimizing elastic strain energy with distinct multiscale distributions of surface domains. Near T c , the material shows the lowest roughness but the highest lateral disorder of the surface.

36 MATERIALS SCIENCE↗

Development of Composite Photocatalyst Materials that are Highly Selective for Solar Hydrogen Production and their Evaluation in Z-Scheme Reactor Designs

The key technology gap preventing a vertically stacked dual-bed particle suspension reactor from achieving the DOE MYRD&D ultimate cost target for H 2 production remains the lack of materials in particle form factor that exhibit ≥10% solar-to-H 2 energy conversion (STH) efficiency as a suspension. Therefore, our project goals centered around strategies to increase the STH efficiency by enhancing photophysical properties of perovskite oxide particles including increased visible-light absorption, increased selectivity for electrocatalysis of the H 2 evolution reaction (HER) and the O 2 evolution reaction (OER) through development of ultrathin oxide coatings, correlating composition and structure to function, and improving understanding of multiscale transport and kinetic processes.

08 HYDROGEN↗

Addressing Critical Problems in Materials Science Through Multiscale and Multimode Characterization (Project 1); Characterization and Optimization of Novel Triple-Conducting Oxide Materials for Energy Applications (Project 2) (CRADA Final Report)

PROJECT 1: Address critical problems in materials science and simultaneously advance the state-of-the-art in multiscale and multimode characterization using the combined advanced analytical capabilities and expertise of Colorado School of Mines (CSM) and the National Renewable Energy Laboratory (NREL). The primary effort of the Phase I of this CRADA is to establish the International Center for Multiscale Characterization using shared resources at both NREL and CSM. Phase II will focus on capability development and marketing, choosing candidate materials science issues in the areas of structure imaging, chemical composition mapping, and correlating properties and performance of materials for impact in energy-related, environmental and critical materials areas. The CRADA will be modified to include specific topics of concern in materials science to industry member partners. Advanced analytical capabilities and expertise at CSM and NREL will be used to advance materials understanding and performance through characterization of multiscale phenomena including structural imaging, chemical composition mapping, and other techniques correlating properties and performance of materials. PROJECT 2: As part of the International Center for Materials Characterization, work under Modification #1 will be led by Colorado School of Mines (CSM), working in collaboration with NREL staff to mentor and advise CSM postdoctoral researchers on set up of diffusion annealing experiments. The purpose of the modification is to provide for NREL staff to mentor and advise CSM postdoctoral researchers on set up of diffusion annealing experiments, including mentoring and advising the CSM-NREL team on proper Secondary Ion Mass Spectrometry (SIMS) data analysis as needed. SIMS measurements of 10-20 samples will be performed at NREL during the project duration.

08 HYDROGEN↗

Machine learning predictions for local electronic properties of disordered correlated electron systems

We present a scalable machine learning (ML) model to predict local electronic properties such as on-site electron number and double occupation for disordered correlated electron systems. Our approach is based on the locality principle, or the nearsightedness nature, of many-electron systems, which means local electronic properties depend mainly on the immediate environment. A ML model is developed to encode this complex dependence of local quantities on the neighborhood. We demonstrate our approach using the square-lattice Anderson-Hubbard model, which is a paradigmatic system for studying the interplay between Mott transition and Anderson localization. We develop a lattice descriptor based on the group-theoretical method to represent the on-site random potentials within a finite region. The resultant feature variables are used as input to a multilayer fully connected neural network, which is trained from data sets of variational Monte Carlo (VMC) simulations on small systems. We show that the ML predictions agree reasonably well with the VMC data. Our work underscores the promising potential of ML methods for multiscale modeling of correlated electron systems.

36 MATERIALS SCIENCE↗

Artificial intelligence techniques for integrative structural biology of intrinsically disordered proteins

Here, we outline recent developments in artificial intelligence (AI) and machine learning (ML) techniques for integrative structural biology of intrinsically disordered proteins (IDP) ensembles. IDPs challenge the traditional protein structure-function paradigm by adapting their conformations in response to specific binding partners leading them to mediate diverse, and often complex cellular functions such as biological signaling, self-organization and compartmentalization. Obtaining mechanistic insights into their function can therefore be challenging for traditional structural determination techniques. Often, scientists have to rely on piecemeal evidence drawn from diverse experimental techniques to characterize their functional mechanisms. Multiscale simulations can help bridge critical knowledge gaps about IDP structure-function relationships - however, these techniques also face challenges in resolving emergent phenomena within IDP conformational ensembles. We posit that scalable statistical inference techniques can effectively integrate information gleaned from multiple experimental techniques as well as from simulations, thus providing access to atomistic details of these emergent phenomena.

59 BASIC BIOLOGICAL SCIENCES↗

Laser Powder Directed Energy Deposition of Steels for Nuclear Applications

This comprehensive investigation examines the structure–property relationships in two nuclear alloy systems—Alloy 709 (A709) austenitic stainless steel and Grade 92 (G-92) ferritic/martensitic (F/M) steel—manufactured via directed energy deposition (DED) for sodium-cooled fast reactor applications. This study establishes the fundamental mechanisms for controlling microstructures for optimizing the mechanical performance of additively manufactured nuclear materials through systematic heat treatment optimization and multiscale characterization. As-deposited A709 steel develops a complex multiscale strengthening architecture consisting of a fine cellular solidification structure with diameter of 2-3 µm within10–50 µm grains, elevated dislocation densities from rapid thermal cycling, and grain boundary precipitates that activate concurrent Hall–Petch, dislocation, and precipitation hardening mechanisms to achieve exceptional properties [yield strength (YS): 603 MPa, ultimate tensile strength (UTS): 844 MPa, Vickers hardness: 220 HV] that achieve a 44% superior strength compared to that of the wrought material. Heat treatments produce different results. Solution annealing (SA) dissolves the cellular structure and reduces the hardness to 190 HV. Precipitation treatment (PT) keeps the cellular structure but adds carbides, allowing the hardness to reach 205 HV. The best approach combines both treatments (SA+PT) and creates uniform precipitate distributions with M 23 C 6 carbides at the grain boundaries and MX carbonitrides in the matrix, achieving a hardness of 195 HV. However, directional differences persist, with a 12%–15% strength variation between orientations due to the inherited layered microstructural architecture that survives aggressive heat treatment. While tensile testing at 550°C demonstrates 40%–50% thermal softening with dynamic strain aging, DED A709 steel still maintains a 71% higher YS than that of the wrought material. Ion irradiation studies (100–400 dpa) of DED A709 steel reveal progressive radiation damage with increasing void density and radiation-induced segregation causing nickel enrichment and chromium depletion, which will ultimately compromise mechanical properties. As-deposited G-92 exhibits exceptional strength (UTS: 1650–1700 MPa, 430 HV) through a complex microstructure containing both ferrite and martensite phases, a high geometrically necessary dislocation (GND) density (17.04×10 14 /m 2 ), and fine carbides. Heat treatments create distinct changes. Normalizing produces fresh martensite with the highest hardness (460 HV) and an increased GND density (20.23×10 14 /m 2 ). Tempering develops dual precipitation systems and reduces the hardness to 290 HV. The optimal approach uses sequential normalizing plus tempering, achieving balanced properties with the lowest hardness (250 HV) and a reduced GND density (11.01×10 14 /m 2 ). A processing-dependent anisotropy is observed: horizontal specimens achieve superior ductile behavior, while vertical specimens exhibit brittle failure. A tempering heat treatment successfully mitigates this anisotropic behavior by transforming the hard martensitic as-deposited structure into tempered martensite enabling both horizontal and vertical specimens to exhibit similar stress–strain characteristics with visible necking behavior. Remarkably, testing at 550°C reveals a reversal in the anisotropy, where as-deposited specimens achieve near isotropy with superior thermal stability (a 15%–20% strength reduction), while tempered specimens develop an orientation dependence with a 25%–30% strength reduction. Both alloy systems demonstrate that DED processing creates specimens with a superior strength through refined microstructural features, though with distinct strengthening mechanisms—austenitic through cellular structures and precipitates versus F/M through phase transformations and precipitates. Heat treatment optimization requires alloy-specific approaches, with A709 benefiting from controlled precipitation while G-92 requires careful phase transformation control. The results show that DED manufacturing can produce nuclear materials with exceptional performance, but directional effects and temperature-dependent behavior must be carefully considered for reactor component design and qualification.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Statistical Properties of Magnetic Structures and Energy Dissipation during Turbulent Reconnection in the Earth's Magnetotail

Here, we present the first statistical study of magnetic structures and associated energy dissipation observed during a single period of turbulent magnetic reconnection, by using the in situ measurements of the Magnetospheric Multiscale mission in the Earth's magnetotail on 26 July 2017. The structures are selected by identifying a bipolar signature in the magnetic field and categorized as plasmoids or current sheets via an automated algorithm which examines current density and plasma flow. The size of the plasmoids forms a decaying exponential distribution ranging from subelectron up to ion scales. The presence of substantial number of current sheets is consistent with a physical picture of dynamic production and merging of plasmoids during turbulent reconnection. The magnetic structures are locations of significant energy dissipation via electric field parallel to the local magnetic field, while dissipation via perpendicular electric field dominates outside of the structures. Significant energy also returns from particles to fields.

79 ASTRONOMY AND ASTROPHYSICS↗

Statistical properties of magnetic structures and energy dissipation during turbulent reconnection in the Earth's magnetotail

We present the first statistical study of magnetic structures and associated energy dissipation observed during a single period of turbulent magnetic reconnection, by using the in situ measurements of the Magnetospheric Multiscale mission in the Earth's magnetotail on 26 July 2017. The structures are selected by identifying a bipolar signature in the magnetic field and categorized as plasmoids or current sheets via an automated algorithm which examines current density and plasma flow. The size of the plasmoids forms a decaying exponential distribution ranging from subelectron up to ion scales. The presence of substantial number of current sheets is consistent with a physical picture of dynamic production and merging of plasmoids during turbulent reconnection. The magnetic structures are locations of significant energy dissipation via electric field parallel to the local magnetic field, while dissipation via perpendicular electric field dominates outside of the structures. Significant energy also returns from particles to fields.

Magnetotail↗

Multimodal Characterization of Crystal Structure and Formation in Rubrene Thin Films Reveals Erasure of Orientational Discontinuities

Abstract Multimodal multiscale characterization provide opportunities to study organic semiconducting thin films with multiple length scales, across multiple platforms, to elucidate crystallization mechanisms of the various microstructures that impact functionality. With polarized scanning transmission X‐ray and 4D‐scanning transmission electron microscopy, hybrid crystalline structures in rubrene thin films in which large crystalline domains surround a common nucleus and transition to a spherulite morphology at larger radii is observed. These high‐resolution techniques reveal how azimuthal orientational discontinuities at smaller radii are erased as spherulite morphology takes hold. In situ crystallization in the films with optical microscopy is also captured, discovering the importance of considering the initial temperature increase of a film during thermal annealing over the crystallization timescale. This kinetic information of the radial crystallization rate and of corresponding film heating kinetics is used to estimate the temperature at which the larger crystalline regions transition into a spherulite. By combining the results obtained from the different characterization modes, it is learned that thermal conditions can sensitively affect the crystallization of rubrene and other organic thin films. The observations suggest opportunities for more complex temperature‐dependent processing to maximize hybrid structures’ functionality in organic thin films and demonstrate that multimodal studies deepen the understanding of structure‐function dynamics.

Tan, Jenna A.↗

Ternary Potassium‐Bismuth‐Telluride Intermetallic Support Promotes Electrochemical Stability in Potassium Metal Anodes

We employed accumulative roll bonding to fabricate self‐standing metallurgical composite of in situ formed alkaline potassium‐bismuth‐telluride intermetallic K 2 (Bi 2/6 Te 3/6 Vac 1/6 ) embedded in potassium metal. This newly discovered thermodynamically stable potassiophilic crystal, termed “KBT”, is fcc antifluorite with K 2 Te archetype. Symmetric cells achieve 880 h of cycling at 0.5 mA cm −2 and 0.5 mAh cm −2 . Potassium metal battery (KMB) with Prussian blue (PB) cathode in carbonate electrolyte retains 80% capacity after 200 cycles at 1C. In ether‐based electrolyte with organic cathode, it achieves 80% retention after 900 cycles at 2C. Combined synchrotron X‐ray nano‐tomography, cryogenic‐focused ion beam microscopy (Cryo‐FIB) and sputter‐down X‐ray photoelectron spectroscopy (XPS) demonstrate uniform electrodeposits, versus baseline of potassium filaments intermixed with pores and coarse SEI. Binary K 3 Bi‐K and K 2 Te‐K intermetallic supports also provide improved electrochemical performance, albeit to lesser extent. Multiscale simulation provides insight into role of support structure in adatom energetics, film nucleation, early‐stage SEI morphology and interfacial stability.

36 MATERIALS SCIENCE↗

Dynamic Bubbling Balanced Proactive CO 2 Capture and Reduction on a Triple-Phase Interface Nanoporous Electrocatalyst

The formation and preservation of the active phase of the catalysts at the triple-phase interface during CO 2 capture and reduction is essential for improving the conversion efficiency of CO 2 electroreduction toward value-added chemicals and fuels under operational conditions. Designing such ideal catalysts that can mitigate parasitic hydrogen generation and prevent active phase degradation during the CO 2 reduction reaction (CO 2 RR), however, remains a significant challenge. Herein, we developed an interfacial engineering strategy to build a new SnO x catalyst by invoking multiscale approaches. This catalyst features a hierarchically nanoporous structure coated with an organic F-monolayer that modifies the triple-phase interface in aqueous electrolytes, substantially reducing competing hydrogen generation (less than 5%) and enhancing CO 2 RR selectivity (~90%). This rationally designed triple-phase interface overcomes the issue of limited CO 2 solubility in aqueous electrolytes via proactive CO 2 capture and reduction. Concurrently, we utilized pulsed square-wave potentials to dynamically recover the active phase for the CO 2 RR to regulate the production of C1 products such as formate and carbon monoxide (CO). This protocol ensures profoundly enhanced CO 2 RR selectivity (~90%) compared with constant potential (~70%) applied at -0.8 V (V vs RHE). We further achieved a mechanistic understanding of the CO 2 capture and reduction processes under pulsed square-wave potentials via in situ Raman spectroscopy, thereby observing the potential-dependent intensity of Raman vibrational modes of the active phase and CO 2 RR intermediates. Finally, this work will inspire material design strategies by leveraging triple-phase interface engineering for emerging electrochemical processes, as technology moves toward electrification and decarbonization.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Observation of a super-tetrahedral cluster of acetonitrile-solvated dodecaborate dianion via dihydrogen bonding

We launched a combined negative ion photoelectron spectroscopy and multiscale theoretical investigation on the geometric and electronic structures of a series of acetonitrile-solvated dodecaborate clusters, i.e., B 12 H 12 2- ·nCH 3 CN (n = 1–4). The electron binding energies of B 12 H 12 2- ·nCH 3 CN are observed to increase with cluster size, suggesting their enhanced electronic stability. B3LYP-D3(BJ)/ma-def2-TZVP geometry optimizations indicate each acetonitrile molecule binds to B 12 H 12 2- via a threefold dihydrogen bond (DHB) B3–H3 ⋮⋮⋮ H3C–CN unit, in which three adjacent nucleophilic H atoms in B 12 H 12 2- interact with the three methyl hydrogens of acetonitrile. The structural evolution from n = 1 to 4 can be rationalized by the surface charge redistributions through the restrained electrostatic potential analysis. Notably, a super-tetrahedral cluster of B 12 H 12 2- solvated by four acetonitrile molecules with 12 DHBs is observed. The post-Hartree–Fock domain-based local pair natural orbital- coupled cluster singles, doubles, and perturbative triples [DLPNO-CCSD(T)] calculated vertical detachment energies agree well with the experimental measurements, confirming the identified isomers as the most stable ones. Furthermore, the nature and strength of the intermolecular interactions between B 12 H 12 2- and CH 3 CN are revealed by the quantum theory of atoms-in-molecules and the energy decomposition analysis. Ab initio molecular dynamics simulations are conducted at various temperatures to reveal the great kinetic and thermodynamic stabilities of the selected B 12 H 12 2- ·CH 3 CN cluster. The binding motif in B 12 H 12 2- ·CH 3 CN is largely retained for the whole halogenated series B 12 X 12 2- ·CH 3 CN (X = F–I). This study provides a molecular-level understanding of structural evolution for acetonitrile-solvated dodecaborate clusters and a fresh view by examining acetonitrile as a real hydrogen bond (HB) donor to form strong HB interactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Insights into Supported Subnanometer Catalysts Exposed to CO via Machine-Learning-Enabled Multiscale Modeling

Subnanometer catalysts offer high noble metal utilization and superior performance for several reactions. However, understanding their structures and properties on an atomic scale under working conditions is challenging due to the large configurational space. Here, we introduce an efficient multiscale framework to predict their stability exposed to an adsorbate. The framework integrates a comprehensive toolset including density functional theory (DFT) calculations, cluster expansion, machine learning, and structure optimization. The end-to-end machine-learning workflow guides DFT data generation and enables significant computational acceleration. We demonstrate the approach for CO-adsorbed Pdn (n = 1–55) clusters on CeO 2 (111). Simulation results reveal that CO can facilitate restructuring by stabilizing smaller planar structures and bilayer structures of specific intermediate sizes, consistent with experimental reports. Metal–support interactions, preferential CO adsorption, and metal nuclearity and structure control catalyst stability. As a result, the framework allows automatic discovery of stable catalyst structures and a systematic strategy to exploit properties in the subnanometer scale.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Selenium Migration and Local Structures in Cu‐Doped CdSeTe Solar Cells after Aging

Selenium grading plays a critical role in state-of-the-art Cadmium Telluride photovoltaic cells by enhancing long-wavelength absorption and extending minority carrier lifetimes —key to enabling the current performance record of 23.08%. However, very little is understood about selenium motion. In this study, a comprehensive, multimodal, and multiscale approach is employed to investigate Se migration and local structural changes in copper (Cu)-doped CdSeTe solar cells subjected to accelerated stress. X-ray fluorescence (XRF) microscopy shows unexpected levels of Se diffusion after 500 h under heat (75°C) and light (0.8 suns, 80 mW/cm 2 ), suggesting the coexistence of fast and slow diffusion channels even at low temperatures, with unexpectedly low activation energies (<0.85 eV). X-ray Absorption Near Edge Structure (XANES) analysis indicates a preferential migration of Se atoms to anionic lattice sites and a reduction in Se-Cl co-passivation at Te-terminated dislocation cores. Furthermore, these findings point to a reconfiguration of Se local environments and highlight the potential role of extended structural defects in enabling Se transport at low temperatures. Additionally, XANES results suggest that the presence of metallic Cu across the absorber layer may contribute to back-contact degradation and reduced hole density in both fresh and aged devices.

14 SOLAR ENERGY↗

Efficient Multiscale Optoelectronic Prediction for Conjugated Polymers

Conjugated polymers represent a high-potential material class due to the tunability of their optoelectronic properties via straightforward processing protocols. However, to in silico tailor these optoelectronic properties, one must compute the conformationally dependent electronic structure over mesoscopic length and time scales. This task represents a challenging multiscale computational problem in which quantum chemistry and atomistic or coarse-grained molecular dynamics must be integrated. Recently, we introduced a methodology, called electronic coarse graining (ECG), that utilizes artificial neural networks to compute ab initio quality electronic structures using only the reduced degrees of freedom of coarse-grained models. Here, we adapt ECG to variable-molecular weight conjugated oligomers by casting the problem in the framework of sequence detection. A machine-learning (ML) architecture employing one-dimensional convolutional neural networks and long short-term memory networks is utilized to compute ground-state orbital energies, charge density distributions, and optical spectra solely from the coarse-grained model's configurational degrees of freedom. Robust molecular weight transferability for ECG is established via a A-ML approach that leverages model electronic Hamiltonians for ground and excited-state property determination. The accuracy and transferability of the ECG methodology opens the door for scalable optoelectronic property prediction in conjugated polymers directly from coarse-grained degrees of freedom.

36 MATERIALS SCIENCE↗

Subtle Molecular Changes Largely Modulate Chiral Helical Assemblies of Achiral Conjugated Polymers by Tuning Solution-State Aggregation

Understanding the solution-state aggregate structure and the consequent hierarchical assembly of conjugated polymers is crucial for controlling multiscale morphologies during solid thin-film deposition and the resultant electronic properties. However, it remains challenging to comprehend detailed solution aggregate structures of conjugated polymers, let alone their chiral assembly due to the complex aggregation behavior. Herein, we present solution-state aggregate structures and their impact on hierarchical chiral helical assembly using an achiral diketopyrrolopyrrole-quaterthiophene (DPP-T4) copolymer and its two close structural analogues wherein the bithiophene is functionalized with methyl groups (DPP-T2M2) or fluorine atoms (DPP-T2F2). Combining in-depth small-angle X-ray scattering analysis with various microscopic solution imaging techniques, we find distinct aggregate in each DPP solution: (i) semicrystalline 1D fiber aggregates of DPP-T2F2 with a strongly bound internal structure, (ii) semicrystalline 1D fiber aggregates of DPP-T2M2 with a weakly bound internal structure, and (iii) highly crystalline 2D sheet aggregates of DPP-T4. These nanoscopic aggregates develop into lyotropic chiral helical liquid crystal (LC) mesophases at high solution concentrations. Intriguingly, the dimensionality of solution aggregates largely modulates hierarchical chiral helical pitches across nanoscopic to micrometer scales, with the more rigid 2D sheet aggregate of DPP-T4 creating much larger pitch length than the more flexible 1D fiber aggregates. Combining relatively small helical pitch with long-range order, the striped twist-bent mesophase of DPP-T2F2 composed of highly ordered, more rigid 1D fiber aggregate exhibits an anisotropic dissymmetry factor (g-factor) as high as 0.09. This study can be a prominent addition to our knowledge on a solution-state hierarchical assembly of conjugated polymers and, in particular, chiral helical assembly of achiral organic semiconductors that can catalyze an emerging field of chiral (opto)electronics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Coverage-Dependent Adsorption of Hydrogen on Fe(100): Determining Catalytically Relevant Surface Structures via Lattice Gas Models

Hydrogen adatoms are a critical surface species for several reactions catalyzed by Fe surfaces such as Fischer–Tropsch synthesis, ammonia synthesis, and the hydrodeoxygenation of biomass-derived oxygenates. Parameterizing the energetics for H/Fe in terms of both coverage and configuration space can significantly aid in the development of multiscale models as well as provide atomic level insight into the dominant surface structures present under realistic reaction conditions. Here, we construct a lattice gas model for H/Fe(100), where the lateral interactions are determined from first-principles using density functional theory. Using 950 symmetrically unique H/Fe(100) configurations, we generate a cluster expansion with a predictive accuracy in terms of surface energy of 3.8 meV/site over a coverage range from 0 to 3 monolayers. Ten electronic ground state structures are identified from this thorough scan (including the structures at 0 and 3 monolayers), which were subsequently used to generate ab initio phase diagrams under a range of temperatures and pressures. Under reaction conditions typical of Fischer–Tropsch synthesis, ammonia synthesis, and biomass oxygenate hydrodeoxygenation, we find that the 1.0 monolayer structure is dominant. Furthermore, examination of the total H–H lateral interactions for the H/Fe(100) electronic ground state structures shows that H/Fe(100) can be accurately modeled via a mean-field ideal lattice gas model for coverages less than 1.0 monolayers. Altogether, this work enables the incorporation of H–H lateral interactions on Fe(100) into multiscale models, via either mean-field or site-dependent techniques, and provides atomic insight into the catalytically relevant H/Fe(100) structures for a range of heterogeneous reactions.

09 BIOMASS FUELS↗