SEARCH · Engineering Papers
Results for “energy density functional”
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.
Decarbonizing nitrogen fertilizer production via the electrochemical nitrogen oxidation reaction
Nitric acid is an important commodity chemical with extensive applications in both agricultural and industrial sectors. However, current production methods involve a combination of the Haber–Bosch and Ostwald processes, which are both energy and carbon emission intensive. The electrochemical nitrogen oxidation reaction (NOR) to produce nitric acid or nitrates shows great potential as an environmentally friendly method for producing fertilizers under mild conditions. The key to progress in this field lies in understanding the fundamental mechanistic insights and establishing robust experimental methods, which is essential for the efficient design and synthesis of electrocatalysts for NOR. Additionally, poor gas mass transport in conventional electrochemical reactors at present lead to lower NOR activity, thereby limiting the progress in this field. In this work, we present a synergistic computational and experimental approach to map out the potential chemical and electrochemical steps and determine the energetics on PtO 2 catalyst to gain mechanistic insights into NOR. Here, this study marks the first attempt to perform NOR in a vapor-fed reactor designed using advanced (additive) manufacturing. The vapor-fed reactor significantly improved the N 2 mass transport to the catalyst, allowing us to report the highest rate for nitrate production to date at 3.3 μmol cm -2 h -1 at 2.01 V vs RHE.
Effect of Sn microalloying on the nucleation of L1 2 Al 3 Zr precipitates in a dilute aluminum-zirconium alloy
While L1 2 -Al 3 Zr nanoprecipitates provide a balance between strengthening and good electrical conductivity, the precipitation of L1 2 -Al 3 Zr in aluminum requires aggressive heat treatments. An improved age-hardening response was observed during isochronal aging of an Al-0.24Zr (wt%) alloy when microalloyed with Sn. A new mechanism termed Low melting point Element-Assisted Nucleation (LEAN) is proposed to explain the lower temperature nucleation of L1 2 -Al 3 Zr precipitates observed in this alloy based on the addition of a low melting point element, such as Sn. Characterization verified the first-principles density functional theory prediction that Zr and Sn atoms cluster during homogenization owing to the favorable binding energy of Zr-Sn-vacancy triplets. Direct microstructural observations revealed these clusters form Sn nanoprecipitates that assist the nucleation of L1 2 -Al 3 Zr at 200°C, where L1 2 -Al 3 Zr precipitation is not expected due to the low diffusivity of Zr atoms in Al. At higher temperatures (≳350°C), the acceleration of L1 2 -Al 3 Zr precipitation is driven by faster Zr diffusion in Al with Sn microalloying and the nuclei formed via the LEAN mechanism. In conclusion, this combination of mechanisms explains the improvement in age hardening through L1 2 -Al 3 Zr precipitation with Sn microalloying.
Activation of dislocations in Mg with solute Y
Mg-Y cast alloy shows excellent ductility (elongation to failure > 15%) compared with pure Mg and commercial Mg cast alloys. By monitoring the microstructure evolution during an in situ tensile test of a Mg-2.5 wt%Y alloy, we identify the activation of prismatic slip, which is rare in Mg. Synchrotron X-ray micro-beam Laue diffraction (μ-Laue) and transmission electron microscopy revealed the morphology of prismatic slip bands and individual dislocations. Density functional theory and molecular dynamics calculations indicate that solute Y can significantly reduce the stacking fault energy (SFE) along direction on prismatic plane in Mg lattice and thus facilitate the nucleation of dislocations during deformation. The presence of free dislocations in the Mg lattice can also lead to nucleation of {10–12} twins even under unfavorable geometric conditions.
Influence of External Conditions on the Black-to-Yellow Phase Transition of CsPbI 3 Based on First-Principles Calculations: Pressure and Moisture
All-inorganic CsPbI 3 perovskite has emerged as a promising candidate for next-generation solar cells. However, owing to the poor structural stability of black phase CsPbI 3 perovskite at room temperature, it spontaneously transforms to the yellow, photoinactive, nonperovskite phase at room temperature, limiting further development of CsPbI 3 perovskite solar cells. Here, to better understand the mechanism driving such undesirable phase transformation, we examine the thermodynamics and kinetics associated with γ-to-δ phase transformation using density functional theory calculation. A solid-state nudged elastic band method was employed to find minimum energy pathways assuming a concerted solid-state phase transformation between these two phases. The gas molecules representing atmospheric (H 2 O, O 2 , and N 2 ) and inert (Ar) ambient environments, as well as applied external pressure, were considered to examine the influence of external conditions on the black-to-yellow phase transition. Our calculation reveals that, with increasing pressure, the thermodynamic driving force for converting the γ-phase to the δ-phase increases, whereas a larger activation energy must be overcome for the phase transition to occur. We also investigate the moisture-induced phase transformation with an H 2 O molecule and its dissociated species (H + /OH – ) and demonstrate that the reaction energy barriers can be significantly lowered in the presence of H 2 O or OH – . On the other hand, other nonpolar species (O 2 , N 2 , and Ar) have a negligible effect on the phase transformation kinetics, while they might reduce the thermodynamic driving force for the phase transformation and suppress the undesirable phase transformation. Our theoretical prediction could support recent observations that the γ-to-δ phase transformation occurs rapidly and catalytically in the presence of moisture, whereas in dry argon and dry oxygen atmospheres, the γ-CsPbI 3 remains stable.
Impact of SO 2 on NiFe Nanoparticle Exsolution and Dissolution from LaFe 0.9 Ni 0.1 O 3 Perovskite Oxides
Ni-doped LaFeO 3 perovskite oxide is a promising cathode material for solid oxide electrolysis cells (SOECs) designed for CO 2 /H 2 O coelectrolysis. Here, the performance of LaFe 0.9 Ni 0.1 O 3 is being investigated under real-world conditions that include exposure to acid gases, such as SO 2 , relevant to SOEC operation. Experiments show that LaFe 0.9 Ni 0.1 O 3 exsolves NiFe nanoparticles, along with the formation of surface SO 4 2– and SO 3 2– after being exposed to 200 ppm of SO 2 . This suggests that the ionic diffusion of Ni 3+ and Fe 3+ between the bulk and the surface remains unaffected throughout the exsolution–dissolution–exsolution cycle. Thermochemical water splitting has been employed as a probe reaction to evaluate the catalytic properties of the exsolved NiFe nanoparticles. These nanoparticles demonstrated improved hydrogen production compared to bare perovskite oxide substrates. However, after exposure to SO 2 , the formation of Fe-rich NiFe nanoparticles led to poor thermocatalytic performance and rapid deactivation of the perovskite at elevated temperatures. Density functional theory (DFT) analysis was utilized to validate the experimental findings, indicating a significantly negative reaction energy for water splitting over exsolved Fe, as well as stronger binding of SO 2 to Fe than to Ni. Computational analysis further suggests that the presence of surface sulfate promotes the formation of Fe-rich NiFe nanoparticles, aligning with the experimental results. Overall, this study clarifies how SO 2 affects the structure of SOEC perovskite oxide candidate materials. Future engineering efforts should focus on enhancing nanoparticle exsolution and sulfur resistance, which is crucial for improving the hydrogen production capacity of La-based perovskite oxides for electro- and thermocatalytic water splitting in real environments containing acid gases.
Nuclear–Electronic Orbital QM/MM Approach: Geometry Optimizations and Molecular Dynamics
Hybrid quantum mechanical/molecular mechanical (QM/MM) methods allow simulations of chemical reactions in atomistic solvent and heterogeneous environments such as proteins. Herein, the nuclear–electronic orbital (NEO) QM/MM approach is introduced to enable the quantization of specified nuclei, typically protons, in the QM region using a method such as NEO-density functional theory (NEO-DFT). This approach includes proton delocalization, polarization, anharmonicity, and zero-point energy in geometry optimizations and dynamics. Expressions for the energies and analytical gradients associated with the NEO-QM/MM method, as well as the previously developed polarizable continuum model (NEO-PCM), are provided. Geometry optimizations of small organic molecules hydrogen bonded to water in either dielectric continuum solvent or explicit atomistic solvent illustrate that aqueous solvation can strengthen hydrogen-bonding interactions for the systems studied, as indicated by shorter intermolecular distances at the hydrogen-bond interface. We then performed a real-time direct dynamics simulation of a phenol molecule in explicit water using the NEO-QM/MM method. Furthermore, these developments and initial examples provide the foundation for future studies of nuclear–electronic quantum dynamics in complex chemical and biological environments.
From Molecules to Solids: A vdW-DF-C09 Case Study of the Mercury Dihalides
The mercury dihalides show a remarkable diversity in the structural preferences in their minimum energy structure types, spanning molecular to strongly bound ionic solids. A challenge in the development of density functional methods for extended systems is to arrive at strategies that serve equally well such a broad range of bonding modes or structural preferences. The chemical bonding and the stabilities of mercury dihalides and the general utility and reliability of the van der Waals density functional with C09 exchange (vdW-DF-C09) in predicting or describing the energetics and structural preferences in these metal dihalides is examined. Here, we show that, in contrast with the uncorrected generalized gradient approximation of the Perdew-Burke-Erzenhoff (PBE) exchange-correlation functional, qualitative and quantitative patterns in the bonding of the mercury dihalide solids are well reproduced with vdW-DF-C09 for the full series of HgX 2 systems for X = F, Cl, Br, and I. The possible existence of a low-temperature cotunnite polymorph for HgF 2 and PbF 2 is posited.
On Proton Defects and the Phase Transformation of NiO 2
One approach to high-capacity Li-ion battery cathodes is to increase the Ni content. Unfortunately, Ni-rich materials undergo a phase transition and volume collapse at a high state of charge (SOC), which degrades battery performance. At a high SOC, NiO 2 undergoes a phase transition from the O 3 to the O 1 phase. The O 1 phase appears to be more thermodynamically stable than the O 3 phase. However, density functional theory (DFT) often predicts higher stability for O 3 , depending on the DFT approximation. The energy difference between the O 1 and O 3 phases is on the order of 10 meV/atom. Furthermore, we find that proton defects are more stable in the O 1 phase than in the O 3 phase by about 90 meV per proton. Defects in NiO 2 energetically favor the O 3 -O 1 phase transition observed in a high SOC in Li $x$ NiO 2 .
Predicting Device Parameters for Dye-Sensitized Solar Cells from Electronic Structure Calculations to Reproduce Experiment.
Given that improvements to the power-conversion efficiency (PCE) of dye-sensitized solar cells (DSSCs) have slowed in recent years, a means to accurately predict the device parameters yielded by trial dyes in silico, without having to synthesize them, would be extremely valuable to speed up the design process. Currently, the best-performing methods of calculating device parameters rely on a set of experimentally determined kinetic coefficients. In practice, it is very difficult to measure these kinetic parameters accurately, limiting the overall accuracy of such predictive methods. This work proposes a model to obtain key parameters such as J(SC), V-OC, and PCE using only the results from density functional theory (DFT) and time-dependent DFT calculations, noting that rates of electron-transfer steps are ultimately linked to the electronic structure of the dye center dot center dot center dot TiO2 working electrode. Six organic DSSC dyes from dissimilar chemical classes (L0, L1, L2, WS-2, WS-92, and C281) were chosen to demonstrate the power of this approach. Their a priori known experimentally determined device performance metrics served to validate our predictions. The greatest absolute error in our predicted PCE values was 0.36% relative to the experiment, while the greatest fractional error was 0.042. This indicates that the proposed model offers a dramatic improvement on previous predictive methods for DSSC device parameters, both in accuracy and in consistency. Moreover, such a predictive model has great potential to be applied to other photovoltaic applications, further enabling the design of novel, highly efficient photoactive materials.
Heteroatom-Doped Graphenes as Actively Interacting 2D Encapsulation Media for Mg-Based Hydrogen Storage
Nanoencapsulation using graphene derivatives enables the facile fabrication of two-dimensional (2D) nanocomposites with unique microstructures and has been generally applied to many fields of energy materials. Particularly, metal hydrides such as MgH 2 encapsulated by graphene derivatives have emerged as a promising hybrid material for overcoming the disadvantageous properties of Mg-based hydrogen storage. Although the behavior of the graphene–Mg nanoencapsulation interface has been studied for many composite materials, the direct modification of graphene with nonmetal foreign elements for changing the interfacial behavior has been limitedly reported. In this regard, using B-doped graphene and N-doped graphene as nanoencapsulation media for tuning the interfacial behavior of graphene derivative–Mg nanoparticles, we present altered hydrogen storage kinetics of heteroatom-doped (B and N) graphene–Mg composites. The effect of heteroatom doping is studied in terms of bonding configurations and heteroatom doping concentrations. The enhancement in hydrogen uptake was observed for all of the heteroatom-doped graphene–Mg nanocomposites. On the other hand, a few samples exhibit significantly low activation energy at the early stage of desorption, which can be related to the facilitated nucleus formation. Density functional theory calculation indicates that B-doping and N-doping accelerate hydrogen absorption kinetics in different ways, aiding charge transfer and inducing surface deformation of Mg nanoparticles, respectively. Their effects can be augmented in the presence of structural defects on graphene, such as vacancies, pores, or graphene edges. These results demonstrate that hydrogen storage kinetics of Mg-based systems can be altered by utilizing heteroatom-doped graphene oxide derivatives as 2D nanoencapsulation media, suggesting that the addition of a nonmetal doping element can also be applied to Mg-based hydrogen storage by modifying the nanoencapsulation interface without forming Mg alloy phases.
Mechanistic Insights into the Electrochemical Oxidation of 5-Hydroxymethylfurfural on a Thin-Film Ni Anode
The electrochemical oxidation of alcohols is being explored as a favorable substitute for the oxygen evolution reaction owing to its capability to generate high-value products and lower overpotentials. Herein, we present a systematic investigation into the electrochemical oxidation of 5-hydroxymethylfurfural (HMF), a model biomass platform chemical, on a thin-film nickel catalyst, aiming to investigate the underlying reaction mechanism and shed light on the role of the catalyst’s microenvironment and phase on activity and product selectivity. Utilizing a combined experimental and computational approach, we demonstrate that NiOOH is the active phase for HMF oxidation. Additionally, we find a substantial impact of the electrochemical environment, particularly the electrolyte pH, on the reaction. Under highly alkaline conditions (pH = 13), higher activity for HMF oxidation is observed, accompanied by an increased selectivity toward 2,5-furandicarboxylic acid (FDCA) production. Conversely, a less alkaline environment (pH = 11) results in diminished HMF oxidation activity and a higher preference for the partial oxidation product 2,5-diformylfuran (DFF). Mechanistic insights from DFT studies reveal that geminal diols that are present under highly alkaline conditions undergo hydride transfer via HMFCA, while a shift to an alkoxide route occurs at a lower pH, favoring the DFF pathway. Hydride transfer energetics are also strongly affected by the surface Ni oxidation state. Finally, this integrated approach, bridging experimental and computational insights, provides a general framework for investigating the electrochemical oxidation of aldehydes and alcohols, thereby advancing rational design strategies in electrocatalysts for alcohol electro-oxidation reactions.
Why Are 5-Thioglycopyranosyl Donors More Axially Selective than their Glycopyranosyl Counterparts? A Low and Variable Temperature NMR Spectroscopy and Computational Study
5-Thioglycopyranosyl donors differ in reactivity and selectivity from simple glycopyranosyl donors. An extensive study has been conducted on the nature and stability of the reactive intermediates generated on the activation of per-O-acetyl and per-O-methyl 5-thioglucopyranosyl donors and the corresponding glucopyranosyl donors. Variable temperature nuclear magnetic resonance (NMR) studies with per-O-methylated or per-O-acetyl glycosyl sulfoxides and trichloroacetimidates on activation with trifluoromethanesulfonic anhydride or trimethylsilyl triflate are reported. These show that following initial adduct formation with the promoter conversion of the 5-thioglucopyranosyl donors to the 5-thioglucopyranosyl triflates requires higher temperatures than that of the glucopyranosyl donors to the glucopyranosyl triflates. It is demonstrated that neighboring group participation is a less important phenomenon for the peracetylated thioglucosyl donors than for the peracetylated glucosyl donors. A simple thiocarbenium ion was generated by protonation of 2,3-dihydro-4H-thiopyran at low temperature and characterized by NMR spectroscopy. However, the corresponding 5-thioglucopyranosyl thenium ions were not observed in any of the NMR studies of the 5-thiopyranosyl donors: the electron-withdrawing C–O bonds around the thiopyranoside core discourage thiocarbenium ion formation, just as they discourage oxocarbenium ion formation. Density functional theory (DFT) calculations reveal the tetrahydrothiopyranyl thiocarbenium ion to be approximately 2.5 kcal/mol lower in energy than the corresponding tetrahydropyranyl oxocarbenium ion relative to the corresponding covalent triflates. However, the computations reveal a 5.8 kcal/mol activation barrier for conversion of the tetrahydrothiopyranyl triflate to the thiocarbenium ion, while formation of the oxocarbenium ion–triflate ion pair from tetrahydropyranyl triflate requires only 2.6 kcal·mol –1 . Overall, the greater axial selectivity of 5-thioglycopyranosyl donors compared to analogous glycopyranosyl donors derives from (i) the lower kinetic reactivity necessitating higher reaction temperatures, (ii) the greater stability of the thiocarbenium ion over the oxocarbenium ion facilitating equilibration under thermodynamic conditions, (iii) the greater magnitude of the anomeric effect in the 5-thiosugars, and (iv) decreased neighboring group participation in the per-esterified 5-thiosugars.
Thermodynamics and electronic structure of adsorbed and intercalated plumbene in graphene/hexagonal SiC heterostructures
Abstract Graphene-covered hexagonal SiC substrates have been frequently discussed to be appropriate starting points for epitaxial overlayers of Xenes, such as plumbene, or even their deposition as intercalates between graphene and SiC. Here, we investigate, within density functional theory, the plumbene deposition for various layer orderings and substrate terminations. By means of total energy studies we demonstrate the favorization of the intercalation versus the epitaxy for both C-terminated and Si-terminated 4H-SiC substrates. These results are explained in terms of chemical bonding and by means of layer-resolved projected band structures. Our results are compared with available experimental findings.
Data-driven design of novel halide perovskite alloys
The great tunability of the properties of halide perovskites presents new opportunities for optoelectronic applications as well as significant challenges associated with exploring combinatorial chemical spaces. Here, in this work, we develop a framework powered by high-throughput computations and machine learning for the design and prediction of mixed cation halide perovskite alloys. In a chemical space of ABX 3 perovskites with a selected set of options for A, B, and X species, pseudo-cubic structures with B-site mixing are simulated using density functional theory (DFT) and several properties are computed, including stability, lattice constant, band gap, vacancy formation energy, refractive index, and optical absorption spectrum, using both semi-local and hybrid functionals. Neural networks (NN) are used to train predictive models for every property using tabulated elemental properties of A, B, and X site species as descriptors. Starting from a DFT dataset of 229 points, we use the trained NN models to predict the structural, energetic, electronic, and optical properties of an extensive dataset of 17955 compounds, and perform high-throughput screening in terms of stability, band gap, and defect tolerance, to obtain 392 promising compounds that are ranked as potential absorbers according to their photovoltaic figure of merit. Compositional trends in the screened set of attractive mixed cation halide perovskites are revealed and additional computations are performed on selected compounds. The data-driven design framework developed here is promising for designing novel mixed compositions and can be extended to a wider perovskite chemical space in terms of A, B, and X species, different kinds of mixing at the A, B, or X sites, non-cubic phases, and other properties of interest.
Transition-state correlations for predicting thermochemistry of adsorbates and surface reactions
Estimating thermochemical properties from linear correlations may provide a pathway to circumvent expensive density functional theory (DFT) calculations for quantities such as pre-exponentials and temperature corrections to DFT energies. Here, in this work, we construct thermochemical scaling relations between C 1 –C 6 n-alkanes in the gas phase and adsorbed alkyl chains extending from several transition metal surfaces, and examine changes in the slope and fit between metals and adsorption sites. We subsequently add –OH, –NH 2 , C$=$O, and C$=$C functional groups to the C 1 –C 6 molecules and demonstrate strong linear correlations for thermochemistry across all species. We broaden the correlations to incorporate transition states of C 1 –C 6 n-alkane dehydrogenation reactions, where thermochemistry for computationally prohibitive transition-state calculations can be quickly assessed. Additionally, we rationalize the linearity of thermochemical correlations based on the composition of the homologous series and theoretical assessments. As an application of the correlations, we estimate pre-exponentials for elementary surface reactions of ethane and propane hydrogenolysis on Ru(0001), which is of relevance to plastic hydrogenolysis. Depending on kinetically important steps, entropic contributions may be necessary to include in certain reaction mechanisms; in contrasting examples, entropies are found to be relatively insignificant for ethane hydrogenolysis but pertinent for propane hydrogenolysis.
Wettability of ultra-small pores of carbon electrodes by size-asymmetric ionic fluids
Recently, we studied the phase behavior of ionic fluids under confinement using the classical density functional theory within the framework of the restricted primitive model. Additionally, the theoretical results indicate that narrowing the pore size may lead to a drastic reduction in the electric double layer capacitance, while increasing the surface electrical potential would improve the ionic accessibility of micropores. In this work, we extend the theoretical investigation to systems containing size-asymmetric electrolytes that may exhibit a vapor-liquid like phase transition in the bulk phase. The effects of pore size and surface electric potential on the phase diagram and microscopic structures of the confined electrolytes were studied over a broad range of parameters. We found that decreasing the pore size or increasing the surface potential could destabilize the liquid phase in micropores, and capillary evaporation could occur regardless of the size asymmetry between cations and anions. Compared to that in a symmetric ionic system, the vapor-liquid phase separation is more likely to take place as the size asymmetry becomes more pronounced. The phase transition would alter the “accessibility” of ions to micropores and lead to coexisting micropores with different surface charge densities as identified by Monte Carlo simulation.