Correlating Optical and Structural Properties of CO on Transition Metal Surfaces
Not Available
SEARCH · Engineering Papers
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.
Not Available
Aminopolymers are attractive sorbents for CO 2 direct air capture applications due to their high density of amine groups, which can readily react with atmospheric levels of CO 2 to form chemisorbed species. The identity of these chemisorbed species and the functional groups that form upon oxidative degradation depends on both material properties and processing conditions, forming a variety of carbonyl-type sites such as ammonium carbamates, bicarbonates, carbonates, carbamic acids, ureas, and amides. 13 C solid-state nuclear magnetic resonance (NMR) is often used to help elucidate the identity of these reacted species, but it is challenging due to the narrow chemical shift range of carbonyl sites. Herein, we demonstrate the application of a two-dimensional (2D) chemical shift anisotropy (CSA) recoupling pulse sequence (ROCSA) to obtain CSA tensor values at each isotropic chemical shift, overcoming limitations of isotropic peak resolution. CSA tensor values describe the local chemical environment and can readily differentiate between the chemisorbed and degradation products. To aid identification, we also developed a k-nearest neighbor (kNN) classification model to distinguish the functional groups via their CSA tensor parameters. This methodology was demonstrated on poly(ethylenimine) in γ-Al 2 O 3 exposed to CO 2 and showed that the chemisorbed products are ammonium carbamate and a mixed carbamate–carbamic acid species. The sample was analyzed again after desorption at 100 °C inducing mild degradation, and the remaining products were strongly bound carbamate and urea species. In conclusion, the combination of 2D CSA measurements coupled with a kNN classification model enhances the ability to accurately identify chemisorbed or degradation products in complex carbon capture materials.
Ammonia (NH 3 ) is one of the most important industrial chemicals. The conventional NH 3 synthesis method-the Haber–Bosch process-converts atmospheric nitrogen (N 2 ) into NH 3 using H 2 with an iron (Fe) catalyst. However, this process requires high pressures (100–200 atm) and temperatures (700–800 K) near thermal equilibrium. Recently, Fe-based nanocatalysts have been reported to produce promising NH 3 yields under atmospheric pressures and temperature-modulated nonequilibrium conditions. Understanding the mechanism of nonequilibrium catalysis with programmed temperature variation could help to optimize this fully electrified and less energy-intensive process. Although reactive molecular dynamics (RMD) simulations can be a useful tool to model nonequilibrium catalytic processes, they require the development of accurate force fields (i.e., interatomic potentials). Here, we present a machine-learned (ML) force field within the Deep Potential MD (DPMD) framework, trained using periodic density functional theory (DFT) calculations, to model NH 3 synthesis on Fe catalysts with various surface adsorbates such as *N, *H, *N 2 , *H 2 , *NH, *NH 2 , and *NH 3 . Here, we generated the DFT data from static models of elementary reactions on the most stable (110) surface of body-centered cubic Fe, which then were augmented by data from constant number of particles–volume–temperature (NVT) DFT-MD trajectories at various temperatures. Finally, we utilized the fully optimized ML force field to investigate reaction dynamics at an Fe(110) surface at linearly increasing temperatures using NVT-DPMD simulations. Our simulations indicate that pulsed temperature ramping could prove favorable for NH3 synthesis. For example, we conducted ramping under multiple sets of conditions: (i) from 900 to 1200 K over periods of 0.1–0.3 ns for Fe surfaces precovered with N or NH along with H; and (ii) from 300 to 600 K over 0.1–0.3 ns for Fe surfaces precovered with NH 3 . While our simulations so far are limited to short time scales (very rapid heating), these observations shed light on the mechanism of the high NH 3 synthesis rate achieved in a novel temperature-modulated nonequilibrium catalytic reactor using pulsed heating and cooling.
In van der Waals materials, coupling between adjacent layers is weak, and consequently interlayer interactions are weakly screened. This opens the possibility to profoundly modify the electronic structure, e.g., by applying electric fields or with adsorbates. Here, in this study, we show for the case of the topologically trivial semimetal 1T′-MoTe 2 that potassium dosing at room temperature significantly transforms its band structure. With a combination of angle-resolved photoemission spectroscopy, scanning tunneling microscopy, X-ray photoemission spectroscopy, and density functional theory we show that (i) for small concentrations of K, 1T′-MoTe 2 undergoes a Lifshitz transition with the electronic structure shifting rigidly, and (ii) for larger K concentrations 1T′-MoTe 2 undergoes significant band structure transformation. Our results demonstrate that the origin of this electronic structure change stems from alkali metal intercalation.
Here, in this study, we investigated the solid residual poly(vinyl alcohol) (PVA) and TiO 2 after nonthermal air, CO 2 , and N 2 plasma depolymerization in the absence and presence of TiO 2 . Scanning electron microscopy studies showed the absence of highly viscous tar and carbonaceous residues on the surfaces of PVA and TiO 2 . In the absence of TiO 2 , PVA particles exhibited micron-sized holes on their surfaces, whereas in the presence of TiO 2 , the surface roughness of PVA particles was observed at the submicron scale. These observations suggest that TiO 2 facilitates the even distribution of nonthermal plasma at a submicron scale, leading to a more uniform depolymerization of PVA surfaces. Raman, Fourier transform infrared spectroscopy, and X-ray absorption spectroscopy showed that (i) the surface of residual PVA contains mainly ketone functional groups and less C-H bonds than the pristine PVA and (ii) further confirmed the absence of highly viscous tar and carbonaceous residues on both used TiO 2 and residual PVA. The nuclear magnetic resonance and mass spectroscopy suggested that the PVA is growing back to poly polyvinyl acetate by the esterification reaction, and the ethers are produced by the acetal reaction between PVA and aldehyde. The transmission electron microscopy and X-ray diffraction analysis indicated no major crystal structural change of the TiO 2 catalyst after the plasma reactions. This study demonstrates that nonthermal plasma-assisted depolymerization is a viable alternative to thermal depolymerization, offering the unique advantage of converting polymer wastes into gaseous small organic molecules without generating recalcitrant viscous tar and carbonaceous residues on the surfaces of the polymer and TiO 2 catalysts.
Carbon capture and storage (CCS) technologies, along with CO 2 capture and conversion methods, have emerged as crucial research areas to address rising CO 2 emissions. In this study, we seek to understand the mechanistic role of amines in enabling lower-energy pathways for CO 2 conversion. Our research focuses on the development and analysis of dual-functional materials (DFMs) engineered for the reactive capture and conversion (RCC) of CO 2 into methane, utilizing Ru catalysts grafted with amine groups. We employ Density Functional Theory (DFT) calculations using methylamine as a model amine to investigate the impact of amine groups on CO 2 methanation on a Ru(0001) surface, both in the presence and absence of amine groups. The amine ligand alters the carbon coordination environment, promoting direct C–O dissociation and potentially destabilizing the CO* adsorbate, thereby reducing the risk of CO poisoning. Additionally, we observe a preference for hydrogenation, although it becomes more energetically uphill in the amine-bound scenario. Our experiments, however, report similar CO 2 conversion and CH 4 production rates over the synthesized catalysts “Ru/TiO 2 ” and the amine (N-(2-aminoethyl)-3-aminoproplytrimethoxysilane (“diaminosilane”)) deposited catalyst “Diamine−Ru/TiO 2 ”. By constructing comparative reaction-free energy diagrams and performing microkinetic modeling (MKM) simulations, we link our theoretical findings with experimentally observed CO 2 uptake, conversion, and methane production rates. A microkinetic model was employed to investigate the anomaly, showing reduced amine–carbon complex coverage and increased CO 2 coverage at all temperatures. The MKM simulations consistently confirmed these trends. In conclusion, this comprehensive approach offers key insights into the role of the amine-CO 2 bond in methanation, highlighting a pathway toward lower-energy, more efficient CO 2 capture and conversion processes.
Not Available
Computation of optical properties using conventional time-dependent density functional theory (TD-DFT) is time-consuming and memory-intensive. In this study, we investigate the accuracy and efficiency of the density functional tight binding (DFTB) framework with newly optimized Slater–Koster (SK) parameters for modeling the structural, electronic properties, and absorption spectra of silver and gold nanoclusters and their alloys. Our investigation of the ground state (GS) properties demonstrates that the newly developed GS-SK parameters enable DFTB to closely approximate DFT-calculated bond lengths for octahedron, tetrahedron, icosahedra, and truncated octahedron with sizes Ag n /Au n (n = 19, 20, 38, 55), nanoclusters and Ag 20 /Au 20 nanoalloys, with a maximum deviation of approximately 0.15 Å. Formation energy results indicate that the GS-SK parameters can closely estimate changes in formation energies with alloy composition, and the comparison of electronic structures for Ag 20 , Au 20 , and AgAu alloy nanoclusters using the DFTB approximation reveals good agreement in the projected density of states (DOS) profiles and energy levels. A second set of SK parameters, ES-SK, has been developed to describe excited state (ES) properties, including the absorption spectra of silver octahedron Ag 19 , tetrahedral Ag n (n = 20, 56, 84), truncated octahedron Ag 38 , and icosahedra Ag 55 closed-shell clusters and their gold and alloy counterparts over a broad range of alloy compositions. This parametrization uses TD-DFTB calculations and fine-tunes the d and p eigenvalues by comparing them to reference absorption spectra from first-principles TD-DFT. This enables the generation of absorption spectra that closely match the reference spectra when plasmon excitation is dominant, as demonstrated by studying the plasmonic properties of icosahedral Ag n and Au n (n = 309 and 561) nanoparticles. This includes the rapid loss in plasmon quality when Au partially replaces Ag in alloy clusters. Furthermore, these results provide a foundation for addressing computational bottlenecks in plasmonics and with new prospects for applications in the quantum plasmonics for bimetallic alloys.
The coarsening of the Ni particles in the hydrogen electrode of solid oxide cells (SOCs) is an important degradation mechanism. Here, in this paper, density-functional theory and kinetic Monte Carlo methods are used to explore our recent hypothesis that the surface diffusion of NiH may cause faster Ni coarsening in electrolysis cell mode under an overpotential. Using both methods, the diffusion constant or diffusivity of NiH on Ni (111) is determined as the product of the surface coverage and single-molecule diffusivity for the first time considering all possible diffusion paths. It is then determined versus overpotential at the triple-phase boundary of the hydrogen electrode assuming a typical operating temperature of the SOC. Under a significant overpotential, the diffusivity of NiH is found to be sufficiently large to support the above hypothesis that NiH may promote Ni coarsening. However, based on the adsorption configurations identified, the dissociation and reformation of NiH on Ni (111) could occur. Thus, more work is needed to develop a model of Ni coarsening considering both molecular and dissociated forms of NiH.
Non-orthogonal configuration interaction with fragment calculations are presented for a number of compounds that show singlet fission properties: (i) four perylene-diimide derivatives, (ii) crystalline pentacene and its (B,N)-substituted variant, and (iii) a regular and a distorted stack of three indolonaphthyridine molecules. The electronic couplings between the singlet excitonic states (S 1 ) and the singlet-coupled double triplet (T 1 T 1 ), the so-called singlet fission coupling, were computed from ensembles with two and three molecules, and except for some small deviations when charge transfer states were included, results are virtually the same. Ensembles of three molecules were used to study the mechanisms of triplet separation, double triplet diffusion, and singlet and triplet exciton diffusion. The calculations show that apart from the standard mechanism for the generation of two uncoupled triplet states (S 1 → T 1 T 1 → T 1 ...T 1 ), there are two other possible pathways: the direct generation from the singlet excitonic state (S 1 → T 1 ...T 1 ) and the process in which the excitonic state evolves in a superposition of T 1 T 1 and T 1 ...T 1 states. Furthermore, the electronic coupling for triplet diffusion is in general much smaller than for singlet diffusion.
Here, we demonstrate the ability to grow epitaxial thin films of manganese–chromium nitride (Mn x Cr 1–x N) across the entire composition range on magnesium oxide (MgO) (1 0 0) surfaces by molecular beam epitaxy. By independently controlling the ratio of the atomic fluxes of manganese and chromium, the compositional dependence of structural and electrical transport properties is determined. Soft X-ray absorption shows that the valence states of Cr( 3+ ) and Mn( 3+ ) remain identical among these MnCrN epitaxial films. The metallic transport behavior of binary manganese nitride (MnN) and chromium nitride (CrN) is suppressed in the ternary alloys. Changes in room temperature electrical resistivity of up to 100 times and 5 orders of magnitude at cryogenic temperatures are observed. This work expands the variety of nitride compounds available for epitaxial integration and could provide additional flexibility for the development of metal–semiconductor–insulator junctions for spintronic devices.
OsO 4 and RuO 4 are molecular oxides with unique tetrameric structures and rare +8 oxidation states. Accurately modeling their properties remains challenging for density functional theory (DFT) due to weak intertetramer interactions, which standard functionals fail to capture. Here, in this work, we show that the van der Waals (vdW)-corrected density functional (vdW-DF-optB86b) provides structural parameters that are much closer to experimental values than the standard generalized gradient approximation, with volume predictions that fall within the experimentally observed range. Phonon band structure analysis shows that the inclusion of vdW interactions stabilizes soft phonon modes, highlighting the importance of dispersion corrections for accurate predictions of lattice dynamics. Experimental measurements of the phonon density of states for OsO 4 , obtained via inelastic neutron scattering, demonstrate good agreement with our vdW-DF-optB86b calculations. These results validate OsO 4 and RuO 4 as valuable benchmarks for structural and vibrational calculations via vdW-corrected DFT methods and offer insights for studying the broader class of sparse molecular materials.
The CO 2 methanation reaction (aka the Sabatier reaction) is one of the important carbon utilization technologies to reduce atmospheric greenhouse gases. In this study, kinetic modeling of the Sabatier reaction proceeding via the CO route has been presented by considering weak and strong chemisorption concepts rooted in electronic theory. Based on the single-route mechanism of the Sabatier reaction, the dependence of the reaction rate on the position of the catalyst Fermi level with respect to the antibonding states of the reagent (CO 2 ) and the product (CH 4 ), is demonstrated with the classic Sabatier volcano-shape. Furthermore, the model was applied to explain experimental results of reaction rates for Ni- and Ru-based catalytic systems and to understand origin of enhancement of turnover frequencies established for dynamic catalysis.
Surface-enhanced Raman scattering (SERS) intensity of two-dimensional (2D) materials critically depends on the resonant conditions and factors such as the substrate interferences and molecule adsorption fluctuations, making comprehensive investigation, understanding, and optimization of 2D materials-assisted SERS challenging. Here, in this work, the wavelength-dependent SERS of van der Waals structures of 2D materials is systematically investigated, focusing on the intrinsic frequency-dependent Raman tensors by first-principles method while ruling out other extrinsic factors in experiments. Distinct enhancement profiles are found for different 2D materials, among which MoS2 and graphene exhibit remarkably strong and broadband enhancement effects. For stacked multilayers and heterostructures of 2D materials, the calculated SERS addresses the significance of the first contact monolayer effect. Based on the above theory, the van der Waals sandwich structures are proposed and investigated as the SERS substrates, verifying a further significantly enhanced SERS performance. This resonant first-principles study demonstrates a comprehensive and analytical way to explore and promote the SERS of van der Waals structures.
Porphyrin-based metal–organic frameworks (MOFs) offer a unique platform for building porous donor–acceptor networks that exhibit long-lived charge separation and transport upon incorporation of electron acceptor guest species. Here, in this work, porphyrin-based MOFs, PCN-222(H 2 ) and PCN-222(Zn), synthesized as nanoparticle suspensions, are successfully infiltrated with fullerene acceptor molecules, C 60 or PC 61 BM, in both polar and nonpolar solvent environments. The location and relative binding strength of these guest species are evaluated through a combination of N 2 physisorption measurements, photoluminescence quenching, and UV–vis absorption titration experiments. Semiempirical tight binding calculations are used to screen potential locations of the fullerene guest within the MOF pores, and hybrid density functional theory (DFT)-computed interaction energies confirm the energetically favorable positions. The fundamental photophysics of these donor–acceptor host–guest combinations are probed using ultrafast transient absorption spectroscopy. Sub-picosecond electron transfer involving initial exciplex population is observed, with slow charge recombination lifetimes on the order of τ ~1 ns for all systems in both dimethylformamide and 1,4-dioxane. Charge recombination occurs through population of fullerene and/or framework porphyrin triplet states depending on the porphyrin metalation status. The photophysics of the fullerene-loaded MOFs are discussed in the context of relevant porphyrin–fullerene donor–acceptor molecules to highlight the unique role of the framework environment in dictating photoinduced electron transfer and decay pathways.
Not Available
Shock-induced detonation of insensitive high explosives based on 1,3,5-triamino-2,4,6-trinitrobenzene starts with formation of hot spots at microstructural defects but has eluded atomistic modeling treatment at micron length scales. To this end, we performed multimicron scale all-atom molecular dynamics (MD) simulations of hot spots that form during the collapse of cylindrical pores with diameters between 10 and 300 nm. Our MD simulations show that hot spots formed at pores larger than 20 nm exhibit temperature fields with scale-invariant features for sizes up to at least 300 nm. Through a continuum-based grain-scale modeling framework, we span and extend beyond the size scales currently accessible to MD and find that hot spot scale invariance is a general feature that arises when the mechanical strength is insensitive to strain rate. Finally, our results demonstrate the applicability of all-atom MD to simulate the complicated dynamical evolution of micron-sized systems and bolster confidence in insights from MD simulations of materials that exhibit strength with negligible rate dependence over the relevant intervals.
Catalytic hydrodeoxygenation (HDO) following catalytic fast pyrolysis (CFP) offers an approach to convert the vapor-phase product of biomass pyrolysis to a stable bio-oil product by reducing the oxygen content. Fundamental insights into the HDO of carboxylic acids, which are a corrosive and acidic CFP product, on promising catalyst materials, such as Pt/TiO 2 , are needed to inform the design of multifunctional HDO catalysts with improved carbon efficiency. In this contribution, density functional theory (DFT) calculations were used to assess the role of Pt-metal and Pt-TiO 2 -interface sites on acetic acid HDO (AA-HDO), and to determine the effect of interfacial oxygen vacancies at the Pt-TiO 2 interface, by calculating the reaction energetics for key AA-HDO surface intermediates and elementary steps on each site type. Pt-metal sites, modeled via Pt(111), preferred to form undesired decarboxylation products (CH 4 and CO 2 ), whereas Pt-TiO 2 -interface sites, modeled via an anatase-supported Pt nanowire, favored the formation of desired deoxygenation products (acetaldehyde and ethane). Interfacial-vacancy sites lowered the activation energy barrier for the first C-O bond-scission step in AA-HDO, predicted to be the rate-limiting step for AA-HDO at the Pt-TiO 2 interface in the absence of a vacancy. These atomistic insights reveal the importance of metal-metal oxide interface sites in AA-HDO selectivity and can be used to inform the rational design of improved HDO catalysts.