Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density functional theory development”

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 1,135 records · Page 63
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Bending as a control knob for the electronic and optical properties of phosphorene nanoribbons

We have assessed mechanical bending as a powerful controlling tool for the electronic structure and optical properties of phosphorene nanoribbons. We use state-of-the-art density functional approximations in our work. The overall performance of the recently developed meta-generalized gradient approximation (meta-GGA) mTASK [Phys. Rev. Materials 5, 063803 (2021)] functional establishes the method as a useful alternative to the screened hybrid HSE06 for better band gaps of phosphorene nanoribbons. We present a detailed analysis to interpret the optical absorption of bent phosphorene nanoribbons using the GW-Bethe-Salpeter approximation. Finally, we demonstrate the important role of the unoccupied in-gap state and conclude that this in-gap state in armchair nanoribbons introduced by bending can significantly affect the properties of low-energy excitons and add some useful opportunities for applications in optoelectronic devices.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Nuclear charge densities in spherical and deformed nuclei: Toward precise calculations of charge radii

Precise measurements of atomic transitions affected by electron-nucleus hyperfine interactions offer sensitivity to explore basic properties of the atomic nucleus and study fundamental symmetries, including the search for new physics beyond the standard model of particle physics. Specifically, such measurements, augmented by atomic and nuclear calculations, will enable extraction of the higher-order radial moments of the charge-density distribution in spherical and deformed nuclei. The new data impose higher precision requirements on a theoretical description. The nuclear charge density is composed of the proton point distribution folded with the nucleonic charge distributions. The latter induce subtle relativistic corrections due to the coupling of nucleon magnetic moments with the nuclear spin-orbit density. Additional corrections come from the effect of center-of-mass projection. We assess the precision of nuclear charge density calculations by studying the behavior of relativistic and center-of-mass motion corrections to the second and fourth charge radial moments. Special attention has been paid to the magnetic spin-orbit density associated with the local variations of the spin-orbit current. The calculations for semimagic and open-shell nuclei are performed in the framework of self-consistent mean-field theory using quantified energy density functionals and density-dependent pairing forces. We used the general expression for the spin-orbit form factor that is valid for spherical and deformed nuclei. We studied the impact of various correction terms on the charge radii, fourth radial moments, diffraction radii, and surface thickness of spherical and deformed nuclei. The spin-orbit corrections to charge radial moments and surface thickness show strong shell fluctuations which can make an appreciable effect when aiming at high-precision predictions of isotopic shifts. The inclusion of relativistic and center-of-mass corrections impacts the quality of energy density functionals optimized to charge radii data. To establish reliable constraints on the existence of new forces from isotope shift measurements, precise calculations of nuclear charge densities of deformed nuclei are needed. The proper inclusion of the spin-orbit charge density and other correction terms is essential when aiming at extraction of subtle effects which become particularly visible in isotopic trends. It is also important when developing high-quality nuclear energy density functionals optimized using heterogeneous datasets involving absolute charge radii, differential charge radii, and charge form factor properties deduced from electron-scattering data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Gas-phase ion-molecule interactions in a collision reaction cell with triple quadrupole-inductively coupled plasma mass spectrometry: Investigations with N 2 O as the reaction gas

Nitrous oxide (N 2 O) was used as a reaction gas to investigate the gas phase ion-molecule interactions using the Agilent 8900 QQQ-ICP-MS. A multi-element standard containing 45 elements with masses ranging from 9 to 208 u was measured in the presence and absence of N 2 O. The main product ion species observed were oxides and nitrides. Comparison of the N 2 O reaction results with similar measurements conducted with O 2 revealed that N 2 O was more effective at forming oxides in general: the elements Cd and Pb were shown to produce oxides with N 2 O where the reaction did not occur with O 2 . Nitrous oxide was also shown to produce a significant amount of nitride species in a few cases. The general reactivity was shown to be consistent with density functional theory (DFT)-predicted reaction enthalpies, such that all predicted exothermic reactions produced product ions at levels at least 1% of the unreacted ion. Our results show that reaction enthalpy is a reasonable predictor of reactivity with N2O on the timescales of the interactions in non-thermal ICP-MS/MS systems. Our work demonstrates the utility of two relatively new platforms (commercial elemental ICP-MS/MS and EMSL Arrows interface to the NWChem program suite), which allows for the study of a large number of elements within a short period. While DFT with the basis sets utilized here is not the most accurate computational method, it is also not computationally expensive and is shown to be suitable for predicting gas phase reactivity in the QQQ-ICP-MS for the majority of ions studied. Here, the ease and rapidity of data collection and DFT calculations has the potential to be very impactful for the identification of targeted reaction chemistries to be leveraged for analytical method development, such as for the inline separation of isobaric interferences from analytes of interest.

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

The Intermetallic Reactivity Database: Compiling Chemical Pressure and Electronic Metrics toward Materials Design and Discovery

Here, the advent of high-throughput Density Functional Theory (DFT) calculations has supported the creation of large databases containing the quantitative output necessary for constructing theoretical phase diagrams and predicting physical properties. In this Article, we present a complementary resource, the Intermetallic Reactivity Database (IRD), focused on the chemical bonding features of solid-state structures and indicators of potential structural transformations. Each IRD entry augments common features, such as band structures and density of states (DOS) distributions, with chemically motivated information including DFT-Chemical Pressure (CP) schemes and visualizable representations of the atomic charges. Together, these data types enable the rationalization and prediction of potential structural phenomena encountered in intermetallic chemistry, as we illustrate with four examples: the origins of the Y 2 Ni 2 Mg structure in terms of CP features of its parent structures, the anticipation of intergrowth phases from the net atomic CPs collected in Al-containing binary phases, the correlation between trends in the CP schemes of CaCu 5 -type phases and experimentally observed structural variations, and finally, the development of theoretical methodology with the testing of a streamlined method generating DFT-CP schemes. Altogether, these examples highlight how the IRD supports the creation of models of structural chemistry that extend beyond the bounds of its entries.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Light-Induced Transformation of Virus-Like Particles on TiO 2

Titanium dioxide (TiO 2 ) shows significant potential as a self-cleaning material to inactivate severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and prevent virus transmission. This study provides insights into the impact of UV-A light on the photocatalytic inactivation of adsorbed SARS-CoV-2 virus-like particles (VLPs) on a TiO 2 surface at the molecular and atomic levels. X-ray photoelectron spectroscopy, combined with density functional theory calculations, reveals that spike proteins can adsorb on TiO 2 predominantly via their amine and amide functional groups in their amino acids blocks. We employ atomic force microscopy and grazing-incidence small-angle X-ray scattering (GISAXS) to investigate the molecular-scale morphological changes during the inactivation of VLPs on TiO 2 under light irradiation. Notably, in situ measurements reveal photoinduced morphological changes of VLPs, resulting in increased particle diameters. These results suggest that the denaturation of structural proteins induced by UV irradiation and oxidation of the virus structure through photocatalytic reactions can take place on the TiO 2 surface. The in situ GISAXS measurements under an N 2 atmosphere reveal that the virus morphology remains intact under UV light. This provides evidence that the presence of both oxygen and UV light is necessary to initiate photocatalytic reactions on the surface and subsequently inactivate the adsorbed viruses. The chemical insights into the virus inactivation process obtained in this study contribute significantly to the development of solid materials for the inactivation of enveloped viruses.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Chemistry and Interfacial Structure Promoting Quasi-van der Waals Epitaxial Growth of WS 2 Nanosheets on Sapphire for Prospective Application in Field-Effect Transistors

How do chemical and structural modifications to the supporting crystal surface affect the subsequent van der Waals (vdW) or quasi(Q)-vdW epitaxial growth of 2D nanocrystals? Developing an atomic-scale picture of such an interfacial system is crucial for understanding its impact on the physical and chemical properties of the supported 2D materials. The elucidation of the interfacial structure and chemistry needed to promote the Q-vdW epitaxial growth of 2D tungsten disulfide (WS 2 ) nanocrystals contributes to the growth mechanism understanding, thus pushing forward the integration of such atomically thin semiconductors toward real field-effect transistor applications. In addition to an atomic-force microscopy top view, we showcase a combination of X-ray techniques for a top-to-bottom investigation of the complexities of the buried interface structures. Furthermore, this approach uses X-ray photoelectron spectroscopy, X-ray standing wave excited X-ray fluorescence, and crystal truncation rod scattering to produce a highly resolved chemical-state-specific 3D atomic map for the extended interface structure of WS 2 /α-Al 2 O 3 (001). Employing these detailed analysis methods, along with density functional theory to further refine the picoscale structure, we demonstrate how two different types of interface engineering during the pregrowth stage lead to significant differences in the chemical and structural modifications to the terminal surface of c-face sapphire, which in turn leads to substantial differences in the submonolayer growth of supported WS 2 2D nanocrystals in terms of lateral domain sizes, epitaxial registry, vdW gaps, and stability.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Machine learning enabled discovery of application dependent design principles for two-dimensional materials

Abstract The unique electronic and mechanical properties of two-dimensional (2D) materials make them promising next-generation candidates for a variety of applications. Large-scale searches for high-performing 2D materials are limited to calculating descriptors with computationally demanding first-principles density functional theory. In this work, we alleviate this issue by extending and generalizing crystal graph convolutional neural networks to systems with planar periodicity and train an ensemble of models to predict thermodynamic, mechanical and electronic properties. We carry out a screening of nearly 45,000 structures for two separate applications: mechanical strength and photovoltaics. By collecting statistics of the screened candidates, we investigate structural and compositional design principles that impact the properties of the structures surveyed. Our approach recovers some well-accepted design rules: hybrid organic-inorganic perovskites with lead and tin tend to be good candidates for solar cell applications and titanium based MXenes usually have high stiffness coefficients. Interestingly, other members of the group 4 elements also contribute to increasing the mechanical strength of MXenes. For all-inorganic perovskites, we discover some compositions that have not been deeply studied in the field of photovoltaics and thus open up paths for further investigation. We open-source the code-base to spur further development in this space.

Venturi, Victor (ORCID:0000000156727549)↗
DOE OSTI2024

Multi-catalytic active site biochar-based catalysts for glucose isomerized to fructose: Experiments and density functional theory study

In this study, this work provides an innovative method for preparing different isomerization catalysts by impregnating different proportions of MgCl 2 and AlCl 3 and combining different K compounds on cellulose-derived biochar, followed by pyrolysis. Results show MgO and Al(OH) 3 existing in 4 Mg- 1 Al-C catalyst can obtain better catalytic effect on glucose isomerization than the singe of Al presenting in 0 Mg- 1 Al-C catalyst. Moreover, the synergism effects of the multi-catalytic active sites such as β-, γ -Al(OH) 3 , KCl, MgO, and K 4 H 2 (CO 3 ) 3 in Mg-Al-KHCO 3 -C catalyst can further lead to an increase in glucose isomerization, compared to the 4 Mg- 1 Al-C catalyst. The X-ray diffraction results present that the value of O/Al in Mg-Al-KHCO 3 -C catalyst is as high as 13.38, which provides many unsaturated acidic catalysis sites and benefits the glucose isomerization. Simultaneously, the TPD results reveal that the main active sites (MgO, Al(OH) 3 , and K 4 H 2 (CO 3 ) 3 ) in Mg-Al-KHCO 3 -C catalyst can provide weakly acidic and basic sites and avoid strongly acidic and basic sites to excessively attack the glucose. Based on the DFT analysis, the results indicate that the MgO has a great effect on the ring-opening reaction to form acyclic glucose, while Al(OH) 3+ has a great effect on promoting acyclic glucose hydrogen transfer isomerized to form fructose. Compared to other carbon-based metal catalysts, the prepared Mg-Al-KHCO 3 -C has excellent catalytic performance, which gives a higher fructose yield (38.7%) and selectivity (87.72%) and glucose conversion (44.12%) at 100 °C in 30 min. In this study, we develop a highly efficient Mg-Al-K-biochar catalyst for glucose isomerization and provide an efficient method for cellulose valorization.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Thermodynamics of monoclinic and tetragonal hafnium dioxide (HfO 2 ) at ambient pressure

Accurate and precise thermodynamic models of oxide compounds and their phases are important for calculating the phase stability of oxide materials. We develop and use a coupled quantum mechanical and molecular dynamics approach to create thermodynamic models of hafnia (HfO2) polymorphs from 0 K to 3000 K at ambient pressure. The approach is based on the quasi-harmonic approximation below the Debye temperature and on ab-initio molecular dynamics calculations above the Debye temperature to predict constant pressure heat capacities (C-p). A Bayesian model provides interpolated values between these regimes. As a case study, we develop thermodynamic models of monoclinic and tetragonal HfO2 polymorphs. The predicted heat capacities are in excellent agreement with experiment, and the predicted temperature of the monoclinic to tetragonal phase transition (2173 K) is in good agreement with the experimental value (2078K). These results provide a comprehensive and accurate thermodynamic model of the monoclinic and tetragonal phases of hafnia on a broad range of temperatures and can serve as input to CALPHAD assessment of multi-component hafnia-based phase diagrams.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Ab Initio modeling of Near-Edge EELS spectra for chemisorbed molecules

Abstract Electron energy loss spectroscopy (EELS) has recently been applied to probe chemisorbed molecules on metal nanostructures, but a fundamental understanding of the correlation between these spectra and the electronic structures of the adsorbates has been limited. We report here on the insights afforded by time-dependent density functional theory to decipher the energy loss near edge structure (ELNES) of EELS spectra associated with chemisorption. These first-principles calculations simulate the ELNES-EELS spectra for chemisorbed CO on various facets of Au and Pt. Computational predictions of key signatures such as the ‘red shift’ and reductions in the peak intensity of the 2 π * and 6 σ * peaks, as compared to free CO in the gas phase, are validated in comparison to experimentally collected EELS spectra. These signatures are revealed to arise from changes in the electronic structure in terms of unoccupied density of states associated with the chemisorption process. They are consistent with a Blyholder model that incorporates donation and back-donation of electrons. They are also characteristic of the chemisorption process, such as the choice of metal, site of adsorption and the coverage and distribution of adsorbates. Our simulations thus provide guidelines for the use of ELNES-EELS to characterize the atomic structure and adsorption property of nanostructured surfaces and facilitate the development of advanced nanomaterials for catalytic applications.

Materials Science↗
DOE OSTI2025

Bayesian prior construction for uncertainty quantification in first-principles statistical mechanics

First-principles statistical mechanics enables the prediction of thermodynamic and kinetic properties of materials, but is computationally expensive. Many approaches require surrogate models to calculate energies within Monte Carlo or molecular dynamics simulations. Inexpensive surrogates such as cluster expansions enable otherwise intractable calculations by interpolating data from higher accuracy methods, such as Density Functional Theory (DFT). Surrogate models introduce uncertainty into downstream calculations, in addition to any uncertainty inherent to DFT calculations. Bayesian frameworks address this by quantifying uncertainty and incorporating expert knowledge through priors. However, constructing effective priors remains challenging. This work introduces and describes practical strategies for building Bayesian cluster expansions, focusing on basis truncation, hyperparameter selection, and ground state replication. We analyze multiple basis truncation schemes, compare cross-validation to the evidence-approximation for hyperparameter optimization, and provide methods to find and enforce ground-state-preserving models through priors. Additionally, we compare the uncertainties between different approximations to DFT (LDA, PBE, SCAN) against the uncertainty introduced with the use of cluster expansion surrogate models. These approaches are demonstrated on the BCC Li x Mg 1-x and Li x Al 1-x alloys, which are both of interest for solid-state Li batteries. Our results provide guidelines for constructing and utilizing Bayesian cluster expansions, thereby improving the transparency of materials modeling. Furthermore, the approaches and insights developed in this work can be transferred to a wide range of cluster expansion surrogate models, including the atomic cluster expansion and related machine-learned interatomic potential architectures.

Alloy theory↗
DOE OSTI2020

First-principles study of magnetism and electric field effects in 2D systems

This review article provides a bird's-eye view of what first-principles based methods can contribute to next-generation device design and simulation. After a brief overview of methods and capabilities in the area, the authors focus on published work by their group since 2015 and current work on CrI3. The authors introduce both single- and dual-gate models in the framework of density functional theory and the constrained random phase approximation in estimating the Hubbard U for 2D systems vs their 3D counterparts. A wide range of systems, including graphene-based heterogeneous systems, transition metal dichalcogenides, and topological insulators, and a rich array of physical phenomena, including the macroscopic origin of polarization, field effects on magnetic order, interface state resonance induced peak in transmission coefficients, spin filtration, etc., are covered. For CrI3, the authors present their new results on bilayer systems such as the interplay between stacking and magnetic order, pressure dependence, and electric field induced magnetic phase transitions. The authors find that a bare bilayer CrI3, graphene|bilayer CrI3|graphene, hexagonal boron nitride (h-BN)|bilayer CrI3|h-BN, and h-BN|bilayer CrI3|graphene all have a different response at high field, while at small field, the difference is small except for graphene|bilayer CrI3|graphene. The authors conclude with discussion of some ongoing work and work planned in the near future, with the inclusion of further method development and applications.

Cheng, Hai-Ping (ORCID:0000000159901725)↗
DOE OSTI2021

Machine learning elastic constants of multi-component alloys

The present manuscript explores application of machine learning methods for determining elastic constants and other derived mechanical properties of multi-component alloys. Here, a number of machine learning models, including linear regression, neural network and random forest based models, are trained and tested on a dataset of binary alloys generated using density functional theory (DFT) calculations and spanning over a large number of elemental species in the periodic table. Starting with a wide range of simple and easily accessible compositionally-averaged elemental features, a correlation-based feature selection strategy was used to systematically down-select a set of most relevant features towards the prediction of the elasticity tensor components. The true predictive performance and the associated uncertainties of the models were established by testing on unseen data and bootstrapping, respectively. A single and pair-wise feature partial dependence analysis was performed to visualize the average property trends in the multi-dimensional feature space in order to further understand the achieved predictive performance. The utility of the trained model is further demonstrated by obtaining sufficiently accurate yet highly efficient approximations for bulk modulus, Young’s modulus, shear modulus and Poisson’s ratio for alloys beyond the binary space (i.e., two-component alloys) on which the model was originally trained. More importantly, we test and validate the predictive performance of the developed model directly against the experimentally measured elastic constants of technologically relevant multi-component alloys (such as, Ni- and Ti-based alloys). Finally, utility of such a data-enabled route is demonstrated by predicting the possible range of various elastic properties for vast composition space available within the five component Ni-Cr-Fe-Mo-W alloy system in a high-throughput manner.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Understanding Inlet Concentration Effects on the Electrocatalytic Conversion of CO 2 to Formic Acid in Gas-Fed Electrolyzers

The electrochemical CO 2 reduction reaction (CO2RR) to produce value-added products remains a developing technology for utilizing waste CO 2 streams. Most device-level CO2RR studies use pure CO 2 gas feeds; however, the effect of dilute CO 2 on the electrolyzer performance is an important consideration for large-scale electrolyzer operation, single-pass conversion, and real-world CO 2 source utilization. This work investigates the effect that the CO 2 concentration has on the performance of formic acid (HCOOH) producing tin oxide (SnO 2 ) and bismuth oxide (Bi 2 O 3 ) catalysts in an electrolyzer device setting. Surprisingly, SnO2 demonstrated an approximately 20% increase in HCOOH selectivity (Faradaic efficiency) when the CO 2 concentration decreased from 100 to 20%. In contrast, Bi 2 O 3 consistently demonstrated high selectivity toward HCOOH across the same CO 2 concentration range. The effects of the CO 2 concentration on selectivity were further investigated with half-cell experiments and in situ Raman spectroscopy, which revealed dynamic changes in the cathodic overpotential and chemical state of the catalyst that depended on the CO 2 concentration. Density functional theory calculations showed how changes in the surface oxidation state of Sn, varying from fully oxidized SnO 2 to metallic Sn(0), affect the thermodynamic barriers of the three main observed products: HCOOH, CO, and H 2 . Our results indicate that dilute CO 2 concentrations required larger cathodic overpotentials to sustain a fixed current density, which, in turn, pushed the Sn-based catalyst toward a more reduced surface that was favorable to HCOOH formation. On the other hand, the Bi-based catalyst remained in a metallic state at CO2RR-relevant potentials and demonstrated a consistent product selectivity regardless of CO 2 concentration. These findings highlight how varying the CO 2 inlet gas concentrations affects the chemical state of catalysts and the resulting performance metrics.

42 ENGINEERING↗
DOE OSTI2023

Rational Design and Reticulation of Infinite qbe Rod Secondary Building Units into Metal–Organic Frameworks through a Global Desymmetrization Approach for Inverse C 3 H 8 /C 3 H 6 Separation

Abstract The development of reticular chemistry has enabled the construction of a large array of metal–organic frameworks (MOFs) with diverse net topologies and functions. However, dominating this class of materials are those built from discrete/finite secondary building units (SBUs), yet the designed synthesis of frameworks involving infinite rod‐shaped SBUs remain underdeveloped. Here, by virtue of a global linker desymmetrization approach, we successfully targeted a novel Cu‐MOF (Cu‐ASY) incorporating infinite Cu‐carboxylate rod SBUs with its structure determined by micro electron diffraction (MicroED) crystallography. Interestingly, the rod SBU can be simplified as a unique cylindric sphere packingqbetubule made of [4 3 .6 2 ] tiles, which further connect the tritopic linkers to give a newly discovered 3,5‐connectedgfcnet. Cu‐ASY is a permanent ultramicroporous material featuring 1D channels with highly inert surfaces and shows a preferential adsorption of propane (C 3 H 8 ) over propene (C 3 H 6 ). The efficiency of C 3 H 8 selective Cu‐ASY is validated by multicycle breakthrough experiments, giving C 3 H 6 productivity of 2.2 L/kg. Density functional theory (DFT) calculations reveal that C 3 H 8 molecules form multiple C−H⋅⋅⋅π and atypical C−H⋅⋅⋅ H−C van der Waals interactions with the inner nonpolar surfaces. This work therefore highlights the linker desymmetrization as an encouraging and intriguing strategy for achieving unique MOF structures and properties.

Chemistry↗
DOE OSTI2023

Assessment of exchange-correlation functionals on oxygen vacancy formation energies of metal oxides

Oxygen vacancies are a common defect in oxides and play vital roles in many technological applications such as oxygen separation, catalytic reactors, solid oxide fuel cells, and solar thermochemical water splitting. The oxygen vacancy formation energy is directly related to the oxide reduction enthalpy and the vacancy concentration and, therefore, is a key quantity underlying these applications. The complexity of measuring oxygen vacancy formation energies experimentally suggests the utility of calculations based on density functional theory (DFT). However, the calculated results are strongly dependent on the exchange-correlation functionals and other parameters such as the Hubbard U . In this work, we compare the performance of the meta-GGA strongly constrained and appropriately normed (SCAN) functional and the commonly used semilocal generalized gradient approximation (GGA) functionals with experimental values of structural parameters, band gaps, magnetic structures, and oxygen vacancy formation energies for six representative oxides, i.e., CaMnO 3 , SrMnO 3 , LaMnO 3 , YMnO 3 , LaFeO 3 , and CeO 2 . Further, our results show that SCAN usually has better agreement with the experimental lattice constants and band gaps and larger magnetic moments than the commonly used GGA functionals. Although SCAN overestimates the oxygen vacancy formation energies of transition metal oxides and therefore requires unusually large Hubbard U values to reproduce the experimental reduction enthalpies, it does predict the correct oxygen vacancy formation energy of CeO 2 , which challenges the commonly used GGA and hybrid functionals. Our results underscore the challenges that exist in describing these complex oxides by DFT and may shed light on developing more accurate exchange-correlation functionals for oxygen vacancy formation energy calculations.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Modulating Operational Conditions to Mitigate Deactivation in Formate Dehydrogenation on Pd Phases

In heterogeneous catalysis, poisoning by surface-bound intermediates poses a major barrier to sustained catalyst performance in (de)hydrogenation reactions. Formate/bicarbonate systems, as liquid organic hydrogen carriers (LOHCs), offer a CO 2 -integrated, low-temperature pathway for hydrogen storage and release, making them attractive for circular energy applications. However, their lower hydrogen density and susceptibility to catalyst deactivation limit their competitiveness compared to conventional LOHCs like methylcyclohexane. Here, this study investigates the mechanistic origins of formate (HCOO – ) dehydrogenation and associated deactivation on Pd interfaces. Using density functional theory (DFT) simulations, we show that under thermocatalytic conditions, strongly bound formate accumulates on the catalyst surface (Pd(111)), blocking active sites, raising activation barriers, and leading to progressive performance loss. Because formate adsorption involves charge transfer, we exploit its sensitivity to electronic structure by modulating the electrochemical potential of the catalyst. Our results reveal that hydrogen transfer from water and formate exhibits opposing potential dependencies, providing insights into the opposing driving forces behind both catalytic activity and poisoning. To further probe this phenomenon, we examine electrochemically induced phase transitions in Pd, focusing on PdO(100) and PdH(110), which are stable under oxidizing and reducing potentials, respectively, and demonstrate enhanced dehydrogenation activity between −0.4 and 0.2 V vs standard hydrogen electrode (SHE). Complementary thermal treatments help decouple kinetic and thermodynamic contributions to intermediate binding. These findings underscore the critical role of the catalyst phase and external stimuli in dictating poison-active site interactions and highlight phase engineering as a promising strategy to mitigate deactivation. This work offers mechanistic insights and design principles for developing more resilient and efficient catalysts for LOHC applications under realistic operating conditions.

Pd phase↗
DOE OSTI2021

Improving and Understanding the Hydrogen Evolving Activity of a Cobalt Dithiolene Metal–Organic Framework

Despite the promising previous reports on the development of electrocatalytic dithiolene-based metal–organic frameworks (MOFs) for the hydrogen evolution reaction (HER), these materials often display poor reproducibility of the HER performance because of their poor bulk properties upon integration with electrode materials. In this work, we demonstrate an in-depth investigation of the electrocatalytic HER activity of a cobalt 2,3,6,7,10,11-triphenylenehexathiolate (CoTHT) MOF. To enhance the durability and charge transport properties of the constructed CoTHT/electrode architecture, CoTHT is deposited as an ink composite (1) composed of Nafion and carbon black. We leverage here the well-established use of catalyst inks in the literature to increase adhesion of the catalyst to the electrode surface and to improve the overall electrical conductivity of the integrated catalyst/electrode. The utilization of the composite 1 leads to a significant improvement in the overpotential (η) to reach a current density of 10 mA/cm 2 (η = 143 mV) compared to prior reports, resulting in the most active MOF-based electrocatalyst for the HER that contains only earth-abundant elements. Extensive density functional theory (DFT) calculations were applied to understand the structure of CoTHT and the mechanistic pathways of the HER. The computational results suggest that an AB stacking geometry is energetically favorable, where one layer is slipped by 1.6 Å relative to the neighboring one along the a and b vectors. Additionally, the DFT calculations indicate that the catalytic cycle likely involves a Volmer discharge step to generate a cobalt hydride, followed by a Heyrovsky step to form a cobalt–H 2 intermediate, and finally the dissociation of H 2 .

DFT calculations↗