Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “surface energetics”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Insights on keto-hydroperoxide formation from O 2 addition to the beta-tetrahydrofuran radical

Cyclic ethers provide an interesting case study of low-temperature oxidation chemistry, especially in relevance to biofuels. A recent experimental study revealed new questions regarding the low-temperature oxidation mechanism of tetrahydrofuran (THF) concerning the formation of keto-hydroperoxides. In particular, keto-hydroperoxides originating from the THF-β radical were not captured accurately by current literature models, motivating this work to calculate the energetics of the first and second O 2 -addition pathways for THF radicals. Electronic structure calculations at the CCSD(T)/cc-pV∞Z//M06-2X/cc-pVTZ level of theory were used to generate potential energy surfaces for the α-C and β-C THF radicals and subsequent pathways to the formation of the keto-hydroperoxide isomers. These are the first theoretical calculations of the second O 2 -addition radical pathways for the THF-β radical. Results from the theorical work provided further insight into the low-temperature oxidation of THF. This included identifying the pathways most likely to form the keto-hydroperoxide isomers observed in prior experimental work; and detecting that the shortcomings in prior models are likely due to uncertainties in R + THF abstraction reaction rates. Here, these conclusions will motivate future work for accurate THF kinetic model development.

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

Structural chemistry of penta- and hexanitrato thorium($\tiny{IV}$) complexes isolated using N–H donors

In this study, a series of fifteen tetravalent thorium phases were prepared. The compounds were isolated from acidic aqueous nitrate solutions using protonated nitrogen heterocycles of varying hydrogen-bond donation strength. Structural analysis via single crystal X-ray diffraction showed that the structures are built from pentanitrato, [Th(NO 3 ) 5 (H 2 O) 2 ] 1- , and hexanitrato, [Th(NO 3 ) 6 ] 2- , molecular units, with the latter being far more prevalent in the solid state. The vibrational properties of the compounds were examined using Raman and IR spectroscopy; the spectra are dominated by stretches characteristic of nitrate and the organic ions. The relative energetics of nitrate complexation was examined using electronic structure theory. These results confirmed that there are clear thermodynamic sinks for the penta- and hexanitrato structural units that were observed experimentally. Additionally, electrostatic surface potentials (ESPs) were calculated in an effort to better understand the counterion stabilization of the complexes. The ESP surfaces showed that the position of the water and nitrate molecules and the coordination geometry of the metal complex had a clear effect on the polarizability of the two structural motifs. Despite limited speciation of the Th–nitrate structural units, the compounds exhibit rich supramolecular chemistry resulting from hydrogen bonding of the Th complexes with the organic N–H donors and π–π stacking interactions from the protonated N-heterocycles.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Stabilized hole-selective layer for high-performance inverted p-i-n perovskite solar cells

P-i-n geometry perovskite solar cells (PSCs) offer simplified fabrication, greater amenability to charge extraction layers, and low-temperature processing over n-i-p counterparts. Self-assembled monolayers (SAMs) can enhance the performance of p-i-n PSCs but ultrathin SAMs can be thermally unstable. Here we report a thermally robust hole-selective layer comprised of nickel oxide (NiO x ) nanoparticle film with a surface-anchored (4-(3,11-dimethoxy-7H-dibenzo[c,g]carbazol-7-yl)butyl)phosphonic acid (MeO-4PADBC) SAM that can improve and stabilize the NiO x /perovskite interface. The energetic alignment and favorable contact and binding between NiO x /MeO-4PADBC and perovskite reduced the voltage deficit of PSCs with various perovskite compositions and led to strong interface toughening effects under thermal stress. The resulting 1.53–electron-volt devices achieved 25.6% certified power conversion efficiency and maintained >90% of their initial efficiency after continuously operating at 65 degrees Celsius for 1200 hours under 1-sun illumination.

14 SOLAR ENERGY↗
DOE OSTI2022

Sulfur-enhanced dynamics of coinage metal(111) surfaces: Step edges versus terraces as locations for metal-sulfur complex formation

The propensity of trace amounts of sulfur adsorbed on coinage metal(111) surfaces to dramatically enhance surface dynamics has been demonstrated by STM observations of accelerated 2D island decay for Cu and Ag. It is generally accepted that this enhancement is due to the formation of adsorbed metal-sulfur complexes, which facilitate surface mass transport of the metal. These complexes were originally proposed to form on terraces following the extraction of metal atoms from step edges and subsequent combination with sulfur on the terraces. However, even when thermodynamically feasible, this mechanism may not be kinetically viable for some complexes due to limited coupling of the complex concentration to the surface diffusion flux of metal atoms. Focusing on the case of Cu, we assess various scenarios where complexes are formed either on terraces or instead directly at step edges, the latter being a new paradigm. In this work, a new pathway is proposed for the formation on terraces. A rich variety of structures incorporating S at step edges exist, which could provide a viable source for complexes, at least from a thermodynamic perspective. However, it is necessary to also assess the activation barrier for complex formation and detachment from step edges. This is facilitated by the nudged-elastic-band analysis of the minimum energy path for this process utilizing machine-learning derived potentials based on density functional theory energetics for the metal-sulfur system.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Surface enrichment dictates block copolymer orientation

Orientation of block copolymer (BCP) morphology in thin films is critical to applications as nanostructured coatings. Despite being well-studied, the ability to control BCP orientation across all possible block constituents remains challenging. Here, in this study, we deploy coarse-grained molecular dynamics simulations to study diblock copolymer ordering in thin films, focusing on chain makeup, substrate surface energy, and surface tension disparity between the two constituent blocks. We explore the multi-dimensional parameter space of ordering using a machine-learning approach, where an autonomous loop using a Gaussian process (GP) control algorithm iteratively selects high-value simulations to compute. The GP kernel was engineered to capture known symmetries. The trained GP model serves as both a complete map of system response, and a robust means of extracting material knowledge. We demonstrate that the vertical orientation of BCP phases depends on several counter-balancing energetic contributions, including entropic and enthalpic material enrichment at interfaces, distortion of morphological objects through the film depth, and of course interfacial energies. BCP lamellae are found more resistant to these effects, and thus more robustly form vertical orientations across a broad range of conditions; while BCP cylinders are found to be highly sensitive to surface tension disparity.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2021

On the structure sensitivity of and CO coverage effects on formic acid decomposition on Pd surfaces

Using density functional theory calculations, the Pd-catalyzed vapor-phase formic acid decomposition was studied, with a focus on the structure sensitivity and CO coverage effects. A comprehensive reaction network was developed on both the (111) and (100) facets of Pd, at CO coverages of 0 and 5/9 monolayer (ML). Pd(100) was determined to be more reactive than Pd(111) at both CO coverages. The introduction of 5/9 ML CO decreased the activity of both facets significantly, due to destabilization of the surface intermediates and transition states on the CO-decorated surfaces. Three reaction pathways were explored on the clean surfaces: the formate (HCOO) pathway, the carboxyl (COOH) pathway leading to the formation of CO 2 , and the COOH pathway leading to the formation of CO (COOH→CO). Based on the DFT-derived energetics alone, it appears that all three pathways contribute to the reaction on clean Pd, whereas the presence of 5/9 ML of CO inhibits the HCOO pathway on both facets and favors the COOH→CO pathway on the (111) facet, but the COOH→CO 2 one on the (100) facet. Moreover, at high CO coverages, alternative spectator CO-assisted adsorbate decomposition pathways were discovered, which could potentially play a role in formic acid decomposition on Pd catalysts under realistic reaction conditions.

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

Simulations and Experiments Reveal Effect of Nanopores on Helium Diffusion in Quartz

The diffusion properties of noble gases in minerals are widely used to reconstruct the thermal histories of rocks. In this work, we combine density functional theory (DFT) calculations with laboratory experiments to investigate controls on helium diffusion in quartz. DFT calculations for perfect α-quartz predict substantially lower activation energies and frequency factors for helium diffusion than observed in laboratory experiments, especially in the [001] direction. These results imply that no helium could be retained in quartz at Earth surface temperatures, which conflicts with observations of partial cosmogenic 3 He retention over geologic time scales. Here, we implement a model of helium diffusion in α-quartz modulated by nanopore defects that disrupt energetically favorable diffusion pathways. In this model, we find that laboratory-determined diffusivities can be most closely reproduced when a helium atom encounters ~70 nanopore sites per million interstitial sites. The results of our model indicate that diffusion of helium in natural quartz, like other noble gases in other minerals, can be significantly modulated by extended defects.

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

Convergent ab initio analysis of the multi-channel HOBr + H reaction

High-level potential energy surfaces for three reactions of hypobromous acid with atomic hydrogen were computed at the CCSDTQ/CBS//CCSDT(Q)/complete basis set level of theory. Focal point analysis was utilized to extrapolate energies and gradients for energetics and optimizations, respectively. The H attack at Br and subsequent Br–O cleavage were found to proceed barrierlessly. The slightly submerged transition state lies −0.2 kcal mol−1 lower in energy than the reactants and produces OH and HBr. The two other studied reaction paths are the radical substitution to produce H2O and Br with a 4.0 kcal mol−1 barrier and the abstraction at hydrogen to produce BrO and H2 with an 11.2 kcal mol−1 barrier. The final product energies lie −37.2, −67.9, and −7.3 kcal mol−1 lower in energy than reactants, HOBr + H, for the sets of products OH + HBr, H2O + Br, and H2 + BrO, respectively. Additive corrections computed for the final energetics, particularly the zero-point vibrational energies and spin–orbit corrections, significantly impacted the final stationary point energies, with corrections up to 6.2 kcal mol−1.

Chemistry↗
DOE OSTI2026

Cluster Dynamics Simulations of Intra-Granular Fission Gas Bubble Size and Pressure Evolution in UO 2

Fission gases such as xenon (Xe) play a critical role in determining the behavior and response of nuclear fuel. Given that Xe has little solubility in UO 2 , it accumulates and forms bubbles, which significantly impact fuel performance. Intra- and inter-granular bubble nucleation and growth can lead to fuel swelling, and once bubbles interconnect at grain boundaries, fission gas can be released into the plenum. At low temperatures, limited uranium vacancy mobility can restrict swelling, therefore causing the bubbles to become highly pressurized. Consequently, this can induce micro-cracking, promote fission gas release (increasing the likelihood of cladding failure), and even lead to fuel pulverization under accident conditions such as a loss of coolant accident. As bubble evolution is strongly influenced by local temperature and fission rate, markedly different behavior occurs across the radial profile of the fuel pellet. Capturing the mechanisms that underpin bubble evolution is therefore important to predict these behaviors in the fuel. Previous models describing important mechanisms informed by lower length scale simulations have been developed under the NEAMS program. These can describe the evolution of a single bubble type (i.e., single value for radius and pressure) at each position in the pellet, for instance using the Centipede cluster dynamic code. However, in reality, a full distribution in bubble sizes and pressures exists within the microstructure at a given position in the pellet. To address this the cluster dynamics code Xolotl, which can predict Xe and vacancy phase space (i.e., bubble distributions) for intra-granular bubbles, has been used before. Prior work benchmarked the Xolotl code against the Centipede cluster dynamics code to ensure compatibility and to verify that mobile defect properties are adequately transferred between the two codes, along with some physics improvements. In this work, we go further by introducing a physics-based set of improvements that will allow us to accurately predict bubble size distributions and internal bubble pressures under representative UO 2 irradiation conditions. The improvements include (i) coupling bubble-defect reaction energies to a virial equation of state (EOS), (ii) including a bubble surface tension contribution, (iii) incorporating radiation-induced re-solution of Xe and vacancies, (iv) enabling pressure-driven dislocation loop punching through an effective emission of interstitial clusters informed by interstitial loop energetics, (v) accounting for radiation induced athermal diffusion of Xe, and (vi) implementing a Booth-type grain boundary sink representation for all mobile defects and defect clusters. After these modifications, we observe good agreement of Xolotl fission gas bubble size and concentration predictions with legacy experimental measurements. Additionally, it allows the distribution of Xe bubble pressures and radius to also be predicted and compared to data produced through the Advanced Fuels Campaign (AFC) program. Here, we have done this by running simulations under conditions similar to the AFC post-irradiation examination (PIE) samples irradiated at North Anna 2 light water reactor (LWR). Our results shows excellent agreement with these experimental measurements.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Acetylene Semi-Hydrogenation on a Perovskite Oxyhydride Surface: Insights from First Principles and Microkinetic Modeling

Although perovskite oxyhydrides (POHs) have shown promise as a support for selective hydrogenation, it is unknown if they themselves can catalyze alkyne semi-hydrogenation. Here, we use first-principles density functional theory, coupled with microkinetic modeling, to investigate acetylene semi-hydrogenation on a prototypical POH, cubic BaTiO 2.5 H 0.5 (BTOH). Two different mechanisms are examined on a representative surface of BTOH under the reaction conditions: both are based on the Horiuti–Polanyi mechanism, but the way of H 2 dissociation is different. In mechanism 1, a lattice hydride H atom and then a surface adsorbed H atom sequentially hydrogenate the adsorbed acetylene. In mechanism 2, two lattice hydride H atoms from the BTOH sequentially hydrogenate the adsorbed acetylene. In both mechanisms, the H atoms are replenished from gas phase H 2 dissociation. Using density functional theory (DFT), we have calculated pathways for these two mechanisms in hydrogenation of acetylene to ethylene, as well as further to ethane. Microkinetic modeling based on the DFT energetics indicates that at 523 K, hydrogenation to ethylene and ethane via mechanism 1 occurs at a rate up to two orders of magnitude faster than via mechanism 2. A selectivity analysis for the temperature range of 373–673 K shows that the product observed is essentially only ethylene within the more active mechanism 1. This is because at the steady state, there is a significant amount of surface anion vacancies that facilitate heterolytic H 2 dissociation and help stabilize a vinyl intermediate. In conclusion, the present findings suggest that POHs are capable of selective hydrogenation, thanks to their rich hydrogen surface-redox chemistry.

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

Enabling Particulate Materials Processing Science for High-Consequence, Small-Lot Precision Manufacturing

This Laboratory Directed Research and Development project developed and applied closely coupled experimental and computational tools to investigate powder compaction across multiple length scales. The primary motivation for this work is to provide connections between powder feedstock characteristics, processing conditions, and powder pellet properties in the context of powder-based energetic components manufacturing. We have focused our efforts on multicrystalline cellulose, a molecular crystalline surrogate material that is mechanically similar to several energetic materials of interest, but provides several advantages for fundamental investigations. We report extensive experimental characterization ranging in length scale from nanometers to macroscopic, bulk behavior. Experiments included nanoindentation of well-controlled, micron-scale pillar geometries milled into the surface of individual particles, single-particle crushing experiments, in-situ optical and computed tomography imaging of the compaction of multiple particles in different geometries, and bulk powder compaction. In order to capture the large plastic deformation and fracture of particles in computational models, we have advanced two distinct meshfree Lagrangian simulation techniques: 1.) bonded particle methods, which extend existing discrete element method capabilities in the Sandia-developed , open-source LAMMPS code to capture particle deformation and fracture and 2.) extensions of peridynamics for application to mesoscale powder compaction, including a novel material model that includes plasticity and creep. We have demonstrated both methods for simulations of single-particle crushing as well as mesoscale multi-particle compaction, with favorable comparisons to experimental data. We have used small-scale, mechanical characterization data to inform material models, and in-situ imaging of mesoscale particle structures to provide initial conditions for simulations. Both mesostructure porosity characteristics and overall stress-strain behavior were found to be in good agreement between simulations and experiments. We have thus demonstrated a novel multi-scale, closely coupled experimental and computational approach to the study of powder compaction. This enables a wide range of possible investigations into feedstock-process-structure relationships in powder-based materials, with immediate applications in energetic component manufacturing, as well as other particle-based components and processes.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Effects of 2D filler on rheology of additive manufacturing polymers: Simulation and experiment on polyetherketoneketone-mica composites

Polymeric composites have attracted increasing interest in fundamental studies and industrial applications due to their often-superior characteristics over conventional pure polymers. In this work, we studied the relationship between the rheological properties, critical in processing methods such as additive manufacturing, and the structure of a polymer-2D filler system by experiment and simulation methods. Considering the high performance polymer polyetherketoneketone (PEKK) with 2D filler (mica) composite as an example, the effects of filler size, filler content, surface chemistry and shear rate were systematically investigated. Here, it is demonstrated that the interaction between the polymer and 2D filler brings two different consequences on tuning the viscosity of polymer-filler composite. The stronger interaction increases the shear stress of the polymer-2D filler composite system, leading to higher viscosity. On the other hand, the stronger interaction can promote more absorption of polymer chains on the surface, which makes the local polymer network disentangled and decreases the continuous viscosity from the polymer network. This work provides a fundamental understanding on the relationship between macroscopic rheological properties and the structural properties and energetics at molecular level of polymer-2D filler composite.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Investigating Combinations of Alkali Metal Oxides and Hydrogenation Catalysts for Reactive Capture of CO2 to Useful C1 Products

As the world endures environmental crises associated with climate change and the rise of atmospheric CO2 concentrations from anthropogenic CO2 emission, carbon capture and utilization (CCU) technologies are increasingly necessary. Reactive carbon capture (RCC) technologies, in which capture and conversion of CO2 occur in a single reactor, are more energetically and economically attractive by avoiding the need to purify, compress, and transport the captured CO2. To this end, dual function materials (DFMs) - composed of sorbents and catalysts co-dispersed on the same high surface area carrier - have been developed. The sorbent component allows for selective capture of CO2 from a gas stream and the catalyst component subsequently performs the in-situ conversion of the adsorbed CO2 upon introduction of a reactive gas (typically H2). The end product of the most established DFM, comprised of Ru and/or Ni with an alkaline sorbent, is methane via the CO2 methanation reaction. While renewable methane would be an excellent transition fuel, fossil methane is inexpensive (averaging $2.57/MMBTU in pre-pandemic 2019) and the economics of renewable methane utilization are noncompetitive. This requires the design and investigation of DFMs that enable CO2 capture and conversion to more valuable and more useful C1 products like CO or methanol (average price of methanol was $20.61/MMBTU in pre-pandemic 2019). These products can then be further upgraded to high energy density synthetic fuels, and related carbonaceous products, for more sustainable alternatives in industries that are difficult to decarbonize, specifically heavy duty vehicles and aviation. Herein, we report various sorbent + catalyst combinations to achieve the production of useful C1 products through reactive capture of CO2.

carbon capture↗
DOE OSTI2023

Reduced Energetic Disorders in Dion–Jacobson Perovskites for Efficient and Spectral Stable Blue LEDs

Metal halide perovskites have witnessed great success in green, red, and near-infrared light-emitting diodes (LEDs), yet blue LEDs still lag behind. Reducing undesired energetic disorders – broad n -phases and halide segregation – is considered as the most critical strategy to further improve the performances. Here, the study reports a newly designed and synthesized di-ammonium ligand with rigid π -conjugated rings and additional methyl groups to construct Dion–Jacobson (DJ) structure. Augmented coordination from the extra ammonium site and increased effective bulkiness from methyl groups lead to better distribution control over conventional mono-ammonium ligands. This enhances the radiative recombination of blue emissions in the film with homogeneous energy landscape and improved surface morphology, as evidenced by a series of imaging and mapping techniques. As a result, it demonstrates DJ perovskite LEDs (PeLEDs) with peak external quantum efficiencies of ≈4% at 484 nm and ≈11% at 494 nm, which are among the top reported for pure DJ phase-based PeLEDs in the corresponding wavelength regions. The results deepen the understanding of regulating energetic disorders in perovskite materials via molecular engineering.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Correlation of Band Bending and Ionic Losses in 1.68 eV Wide Band Gap Perovskite Solar Cells

Abstract Perovskite solar cells (PSCs) are promising for high‐efficiency tandem applications, but their long‐term stability, particularly due to ion migration, remains a challenge. Despite progress in stabilizing PSCs, they still fall short compared to mature technologies like silicon. This study explores how different piperazinium salt treatments using iodide, chloride, tosylate, and bistriflimide anions affect the energetics, carrier dynamics, and stability of 1.68 eV bandgap PSCs. Chloride‐based treatments achieved the highest power conversion efficiency (21.5%) and open‐circuit voltage (1.28 V), correlating with stronger band bending and n‐type character at the surface. At the same time, they showed reduced long‐term stability due to increased ionic losses. Tosylate‐treated devices offered the best balance, retaining 96.4% efficiency after 1000 h (ISOS‐LC‐1I). These findings suggest that targeted surface treatments can enhance both efficiency and stability in PSCs.

14 SOLAR ENERGY↗
DOE OSTI2020

From electricity to fuels: Descriptors for C 1 selectivity in electrochemical CO 2 reduction

Electrochemical reduction of carbon dioxide (CO 2 RR) over transition metals follows a complex reaction network. Here, we combine observations from experimental literature with a theoretical analysis of energetics to rationalize that not all intermediates in the reduction of CO 2 are formed through direct protonation steps. We derive a selectivity map for two-electron products (carbon monoxide (CO) and formate) on pure metal surfaces using only the CO and OH binding energies as descriptors. For the pure metals that are selective towards CO formation, the variation of the CO binding energy is sufficient to further subdivide the map into domains that predominantly form H 2 , CO, and more reduced products. Our analysis rationalizes experimentally observed product distributions in CO 2 RR across pure metal systems. Overall, we highlight the need for additional material screening descriptors for CO 2 R and the importance of considering competition from the elementary steps of the hydrogen evolution reaction.

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

The predicted shapes of voids and Xe bubbles in UO 2

Morphology is a fundamental attribute when investigating voids and bubbles in UO 2 . Here we use molecular dynamics and Monte Carlo simulations to predict the lowest energy shapes for voids and bubbles in UO 2 . The energies of the {100}, {110} and {111} surfaces have been calculated and used to predict the equilibrium void shape from Wulff construction. This equilibrium shape is compared to low energy faceted voids exhibiting different relative proportions of each family of terminating surfaces. It is found that the equilibrium Wulff shape does not represent the lowest energy morphology for nm void sizes at temperatures between 300 K and 1200 K. Furthermore, the lowest energy faceted voids are slightly more energetically favourable than spherical voids, and as Xe is added, and bubble pressure increases, the faceted morphology becomes even more favourable than the spherical shape.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Ideal internal kink stability in presence of plasma flow and neoclassical toroidal viscosity due to energetic particles

Abstract The influence of energetic particles (EPs) on the ideal internal kink mode, in rotating tokamak plasmas, is numerically investigated by simultaneously solving MHD-kinetic hybrid equations together with a toroidal momentum balance equation utilizing the MARS-Q code (Liu et al 2013 Phys. Plasmas 20 042503). The neoclassical toroidal viscous (NTV) torque, induced by precessional drift resonances of trapped energetic particles, acts as the momentum sink term to damp the plasma flow. Quasi-linear initial value simulations show local reduction of the flow amplitude and enhancement of the flow shear near the q = 1 rational surface ( q is the safety factor) due to EP induced NTV. Both effects in turn destabilize the internal kink mode. These numerical findings are robust against the initial linear stability of internal kink, the initial plasma flow profile, as well as the equilibrium distribution model for EPs.

Physics↗