Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Energy landscapes”

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 OSTI2024

Vacancy diffusion barrier spectrum and diffusion correlation in multicomponent alloys

Vacancy diffusion serves a crucial role in many important kinetic behaviors and properties of multicomponent alloys. Essential questions, however, persist regarding how chemical complexity affects diffusion and what unique characteristics, if any, set these alloys apart from traditional metals. Using neural network kinetics model, we study vacancy diffusion in NbMoTa alloy across a broad temperature range (2600 to 800 K). Unlike pure metals, the two key diffusion parameters—diffusion correlation factor f and activation energy ΔG m —are not constant in alloys, but instead substantially decrease with decreasing temperature. This temperature dependence arises from a reduced number of active vacancy jump pathways at lower temperatures, leading to more correlated diffusion. Upon examining vacancy diffusion throughout the entire compositional space of the Nb-Mo-Ta system, we discover that the slowest vacancy diffusion surprisingly occurs in the non-equimolar region, rather than the equimolar concentration where the configurational entropy is highest. The diffusion barrier spectrum, characterizing the diffusion energy landscape, is an intrinsic material characteristic, which controls both f and ΔG m and, thereby, the diffusivity. Lastly, we find that the vacancy diffusion rate drops noticeably in the presence of local chemical order in the NbMoTa system, particularly for MoTa alloys with long-range B2 order.

Defects↗
DOE OSTI2025

On the distribution of stacking faults at dissociated medium-angle grain boundaries: Crystallographic geometry and metastability

Grain boundaries in FCC metals with low stacking-fault energy can form in dissociated configurations of stacking faults. Perhaps the most studied example is 9R stacking at boundaries near Σ3{112}, where the distribution of stacking faults is related to the emission of Shockley partial dislocations. Here, we combine atomic-resolution electron microscopy, atomistic simulations, and dislocation theory to demonstrate that boundaries vicinal to Σ33a support the stabilization of dissociated 9R stacking within a narrow range of inclinations. This boundary is interesting since its misorientation (20.05°) lies in the medium-angle regime, just past the upper misorientation limit for low-angle boundaries, motivating questions for how best to describe it in terms of dislocations. Our HAADF-STEM observations of thin film bicrystals, supported by atomistic modeling, reveal that this inclination dependence arises from specific geometric constraints on the arrangement of Shockley partial dislocations at the interface. Quantification of stacking-fault distributions across multiple boundaries indicates that the density and spacing of faults closely follow the ideal 9R motif, with subtle variations reflecting the complex energy landscape of these boundaries. Through energetic analysis, we establish the presence of competing metastable states enabled by variations in stacking sequences, emphasizing the significant role of crystallographic geometry. We generalize our analysis as a function of misorientation, showing how 9R at the Σ33a boundary is related to previous observations and calculations of HCP at a 29.7° boundary. This study provides a crystallographically grounded framework connecting dislocation structures, stacking-fault distributions, and metastability at grain boundaries in FCC metals.

Atomistic modeling↗
DOE OSTI2024

Colloidal phenomena reflect the interplay between interfacial solution structure, interparticle forces, and dynamical response

The physical and chemical properties of colloidal systems are central to both natural and synthetic processes that enable a wide range of energy and environmental applications. These properties are an emergent outcome of the structure that colloidal systems adopt in response to the interplay of electrical, van der Waals, hydration, and steric forces coupled to the inherent fluctuations of the system. Here, in this study, we review the nature of those forces and show how variations in their relative importance lead to distinct free energy landscapes and consequent dynamic behaviors, illustrated for colloidal dispersions, nanocrystal superlattices, and mesocrystlline solids formed through oriented attachment. We end by discussing the challenges to future progress arising from the need to seamless couple from molecular to particle to ensemble scales.

42 ENGINEERING↗
DOE OSTI2022

CASM — A software package for first-principles based study of multicomponent crystalline solids

CASM is a software package that enables first-principles based studies of crystalline materials. It has been designed to treat coupled chemical, mechanical, vibrational, and magnetic degrees of freedom to determine ground state and finite temperature properties of crystals. The symmetry of the underlying parent crystal structure is used to enumerate perturbations of the parent crystal structure and generate derivative structures which can be input to first-principles calculations in order to explore the ground state energy landscape. CASM algorithmically constructs cluster expansions that fully couple discrete and continuous degrees of freedom, and generates highly efficient code to evaluate the cluster expansion basis functions. Widely used machine learning methods are integrated for fitting expansion coefficients to first-principle calculations. The fully parameterized cluster expansions can be combined with (kinetic) Monte Carlo methods to calculate finite temperature thermodynamic and kinetic properties. CASM Alloy Manager identifies distinct parent crystal structures in an alloy system, creating individual projects for each, and integrating the results. The integrated infrastructure facilitates the linkage between first-principles statistical mechanics predictions with higher length scale computational methods, such as phase field simulations. In conclusion, CASM projects are designed to be easy to share in repositories, to re-use, and to extend to include additional chemical species or types of degrees of freedom.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Phase-field modeling of solid-state metathesis reactions with the charge neutrality constraint

In this work, we present a phase-field model that captures the evolution of ionic concentrations and phase fractions during solid-state metathesis (SSM) reactions where diffusion limits the rate of transformation. The evolution of the mole fraction of each ion is obtained via governing equations that describe the reduction of a free energy, which includes an energy landscape with local minima located at compositions corresponding to stable products. We utilized two Lagrange multipliers to impose constraints of electroneutrality as well as on the sum of mole fractions, which were then eliminated to derive set of two partial differential equations that describe the dynamics of the mole fraction evolution. From these governing equations, the expressions for effective mobilities for the cations and the anions were obtained. We first study the effect of mobilities of ions on the reaction kinetics, using a simple model considering the ions with an identical absolute value of charge numbers. The simulation results show that the overall characteristic mobility, defined as the sum of the two effective ionic mobilities, provides an excellent measure of the rate at which reaction progresses and that the ratio of the effective mobilities of the anions and the cations signifies the manner by which the reaction progresses. We then generalize the model to consider ions with different charge numbers and tuned the mobility of ions based on their diffusion coefficients reported in the literature and experimental data from a thin-film experiment for the synthesis of FeS 2 to demonstrate the capability of the model to predict the phase evolution during SSM reactions. In particular, the simulation predicts nonplanar phase evolution, which is recently observed in thin-film reactions for the synthesis of FeS 2 via transmission electron microscopy. The approach can serve as a basis for models for phase transformations in other multiphase ionic mixtures, such as in all-solid-state batteries and in ionic liquids.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Efficient Screening of Bi–Metallic Electrocatalysts for Glycerol Valorization

Glycerol is a byproduct of biodiesel production and, as such, it is of limited economic value. By means of electrooxidation, glycerol can be used as a feedstock for scalable hydrogen production, in addition to conversion to value-added products. The development of novel and efficient catalytic electrode materials for the anodic side of the reaction is a key towards a hydrogen-based energy economy. In the present study, a computational screening protocol combining DFT, scaling relations, and microkinetic modeling allows for a rational selection of novel catalysts that can deliver efficient glycerol electrooxidation, low cost of production, and environmental sustainability. Activity and chemical selectivity towards hydrogen production on pure metal catalysts is discussed in terms of volcano-shaped plots. We find that the selectivity in the glycerol oxidation reaction is influenced by a different energy landscape when in the presence of water and best classified by a comparison of O—H and C—H bond-breaking barriers. In addition, we screened 3570 bi-metallic catalysts in the AB (L1 0 ) and A 3 B (L1 2 ) ordered structures for activity, stability, price, and toxicity. By filtering based on the criteria for toxicity, resistance to oxidation, miscibility, and price, we have identified 5 L1 0 structured catalysts (AgPd, AuPd, PtSb, CuPt, and AgPt) and 20 L1 2 catalysts (Ga 3 Ta, In 3 Ta, Ir 3 W, Ir 3 Mo, Cu 3 Pt, Ir 3 Ta, Ir 3 Re, Pd 3 Bi, Pd 3 Cu, Pd 3 W, Pd 3 Co, Pd 3 Sn, Pd 3 Mo, Pd 3 Ag, Pd 3 Ga, Pd 3 Ta, Au 3 Ru, Pd 3 In, Au 3 Ir, and Pd 3 Au) that are all predicted to show high activity. We also identify an additional 37 L1 0 and 92 L1 2 structured electrocatalysts with an anticipated medium-high activity.

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

Metal dissolution from first principles: Potential-dependent kinetics and charge transfer

Metal dissolution and its inverse process are integral to both corrosion and electrodeposition; however, many mechanistic details regarding the dissolution process are challenging to decipher. These include how ion dissolution kinetics and charge transfer are influenced by the competition between metal and solvent interactions under an electrode potential. Here, in this work, we introduce a computational framework based on density functional theory with grand-canonical treatment of electrons to directly predict the potential energy landscape for metal dissolution at a constant potential. Using aluminum as an example, we demonstrate that dissolution kinetics is governed by competing kinetics between two physical processes associated with metal–metal bond breaking and ion-migration within the electrical double layer, respectively. We identify a kinetic transition between regimes dominated by each of these processes and show that this transition depends on the operating electrode potential, among other key factors. It is further found that kinetics and thermodynamics of these processes can be described with a simple, one-parameter Marcus-theory-type model. Beyond offering new understanding of charge transfer during dissolution, our simulation protocol provides a recipe for directly predicting other important quantities in electrochemical reactions from first principles that are difficult to measure, such as the symmetry factor.

36 MATERIALS SCIENCE↗
DOE OSTI2022

High temperature stability and transport characteristics of hydrogen in alumina via multiscale computation

Here, the impact of hydrogen charge states on the stability and transport characteristics of hydrogen interstitials in alumina polymorphs is evaluated by multiscale computational methods including density functional theory (DFT), ab initio molecular dynamics (AIMD) and machine learned force fields. Thermodynamic calculations show that the protonic H i +1 interstitial is the most stable defect species for most values of the electronic bandgap in both and amorphous alumina (Al 2 O 3 ). Further, active learned Gaussian approximation potentials (GAP) were developed using AIMD data to study temperature dependent long time proton diffusion in alumina. Diffusivity calculations from GAP-MD simulations are found to be comparable with of the AIMD data, while being ~340 times faster and scalable to larger systems. Comparisons with diffusivity values for other interstitial charge states (H i 0 and H i -1 ) and published experimental literature indicate that H i +1 diffusion is the likely mechanism of hydrogen transport. A good agreement is obtained between H i +1 diffusivity calculated in α-Al 2 O 3 from DFT: 5.05 10 -3 exp(-0.81 eV/k B /T) cm 2 /s and reported experiment: 9.7X10 -4 exp(-0.83 eV/k B /T) cm 2 /s. Computationally and experimentally calculated energy barriers (0.81 and 0.83 eV respectively) only differ by 2.5%. Similarly, the pre-exponential diffusion coefficients only differ by 0.5 orders of magnitude. Moreover, the diffusivity of H i +1 in amorphous Al 2 O 3 in the 1000–2000 K range is calculated to be 2.53X10 -2 exp(-0.89 eV/k B /T), just one order of magnitude higher than the corresponding value in α-Al 2 O 3 . This suggests that local structural disorder does not significantly affect the energy landscape and diffusion behavior of H i +1 in Al 2 O 3 . Overall, these results show promise for the application of alumina polymorphs as hydrogen permeation barriers.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Economic assessment of seismic monitoring for underground hydrogen storage

Underground hydrogen storage (UHS) plays a key role in the energy landscape. However, like other subsurface engineering technologies, UHS may cause leakage into the groundwater or atmosphere and possibly induce local seismicity. To reduce these risks, seismic monitoring could be a viable technique to track the UHS plume, detect leakages, and locate induced seismicity events. Seismic monitoring has been proposed to safely monitor UHS, but research in this area is still new and requires field studies. Lab and theoretical studies have demonstrated the validity of seismic monitoring for UHS. Therefore, it is imperative to analyze the economic feasibility of seismic monitoring for UHS. Hence, we develop a cost model and open-source Python code for seismic monitoring that considers types of seismometers, comprehensive operational scenarios, detection thresholds, and long-term leakage monitoring. A case study is further provided to validate the cost model on reservoir simulations of UHS. We find that the levelized cost for a 10-year operating UHS site will range on the order of ∼0.003 $\$$/kg. The methods developed in this study could also be applied to the monitoring of groundwater, gas, and/or wastewater injection.

08 HYDROGEN↗
DOE OSTI2025

Machine learning guided search for energetically favorable metal borocarbide ternary compounds

In this work, we employ machine-learning (ML) combined with first principles calculations to efficiently search for the energetically favorable metal borocarbide (M-B-C) ternary compounds with M being the group 1–3 metal elements. Using a crystal graph convolutional neural network (CGCNN) ML approach followed by first-principles calculations, we predicted 47 energetically favorable stable and metastable ternary Na-B-C, Ca-B-C, and La-B-C ternary compounds with their decomposition energy (E d ) below or within 100 meV/atom from the currently known convex hulls. Phonon spectra and electronic structures of the 14 energetically favorable stable structures are also investigated by first-principles calculations. By substituting the metal atoms in the 29 energetically favorable non-equivalent template structures of Na (Ca, La)-B-C with other group 1–3 elements in the periodic table, we further obtain 22 stable structures and 52 metastable structures (E d ≤100 meV/atom with respect to the known convex hulls) for Li-B-C, K-B-C, Rb-B-C, Mg-B-C, Sr-B-C, Ba-B-C, Sc-B-C and Y-B-C ternary compounds. New convex hulls including our newly predicted stable ternary structures and the known stable structures are constructed for the M-B-C systems. The results obtained by our ML guided first-principles calculations enrich our knowledge in the structure and energy landscape of metal borocarbide ternary compounds and provide useful guidance for further experimental synthesis and discovery.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Energetics of water expulsion from intervening space between two particles during aggregation

Solvent expulsion away from an intervening region between two approaching particles plays important roles in particle aggregation yet remains poorly understood. Here, in this work, we use metadynamics molecular simulations to study the free energy landscape of removing water molecules from gibbsite and pyrophyllite slit pores representing the confined spaces between two approaching particles. For gibbsite, removing water from the intervening region is both entropically and enthalpically unfavorable. The closer the particles approach each other, the harder it is to expel water molecules. For pyrophyllite, water expulsion is spontaneous, which is different from the gibbsite system. A smaller pore makes the water removal more favorable. When water is being drained from the intervening region, single chains of water molecules are observed in gibbsite pore, while in pyrophyllite pore water cluster is usually observed. Water-gibbsite hydrogen bonds help stabilize water chains, while water forms clusters in pyrophyllite pore to maximize the number of hydrogen bonds among themselves. This work provides the first assessment into the energetics and structure of water being drained from the intervening region between two approaching particles during oriented attachment and aggregation.

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

Twin nucleation and growth in hexagonal close-packed metals: The role of slip-mediated plasticity on twin embryo formation and evolution

Twinning is a key deformation mechanism in hexagonal close-packed (hcp) metals, which are typified by a lack of easily activated slip systems that can accommodate a general state of loading. Pragmatically, the nucleation and evolution of twin domains occurs concomitantly with slip, such that the eventual twin network is conditioned by both external and internal stresses resulting, among others, from the evolution of dislocations. However our understanding of the interplay between dislocations, and twin nucleation and stability remains limited. This work focuses on elucidating the influence of dislocation-mediated plasticity on the formation, growth, and stability of twin embryos. First, a new mesoscale spectral crystal plasticity-twinning framework is developed and used to quantify the change in the free energy landscape following the nucleation {$10\bar{1}2$} of twins in Mg. Representative twin morphologies are modeled under conditions of limited and profuse slip activity, which are emulative of small- and bulk-scale samples, respectively. Then, the driving traction profiles around twin embryos are investigated via a sharp interface approach to obtain insights into how concurrent slip-mediated plasticity can influence the growth/stabilization of nanometric twins. The driving traction profiles are further utilized to determine the stability of twin embryos post loading. The initial dislocation density in the samples, and within the different domains (i.e., twin vs. parent), is seen to have a significant effect on the twin nucleation stress. Namely, the activation of high levels of concomitant plasticity in the parent grain is seen to significantly drive the formation of nanometric twin nuclei at stresses as low as 250 MPa. Further, the sharp interface analysis reveals that profuse plasticity in the parent grain simultaneously alters the forward and back stresses, such that the magnitude/polarity of the driving tractions become increasingly favorable for nanometric twin growth when slip is active. Finally, prior plasticity in the parent grain is seen to result in favorable driving tractions for nanometric twin growth, even at applied stresses as low as ~ 100 MPa. In conclusion, these results are in stark contrast to a case without any dislocations, wherein applied stresses as high as ~ 500 MPa are necessary to grow the twin domains.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Small uranium and oxygen interstitial clusters in UO 2 : An empirical potential study

The aggregation of defect clusters in UO 2 under irradiations creates dislocation loops and cavities, which in turn contribute to swelling of the fuels, leading to detrimental pellet-clad mechanical interactions. The kinetics of the defect clusters, however, have received less attention and are not well understood. As such, using two interatomic potentials, the Morelon potential and a many-body model by Cooper, Rushton, and Grimes (CRG), the energetics and kinetics of small uranium interstitial clusters, with and without bound oxygen interstitials, are investigated with molecular dynamics and accelerated molecular dynamics simulations. For the single bound anti-Schottky interstitial cluster (U i 2O i ), both potentials find the interstitialcy migration mechanism to be a low barrier migration mechanism. However, for the CRG potential, an oxygen Frenkel pair assisted interstitial cluster migration mechanism with a lower migration barrier is identified during accelerated molecular dynamics simulations. For the Morelon potential, the interstitial clusters have geometrically compact structures, and small interstitial clusters (two and three anti-Schottky clusters) are more mobile than the single bound anti-Schottky interstitial cluster; for the CRG potential, the interstitial clusters have geometrically non-compact structures, yielding a more complex energy landscape. The rapid migration of interstitials and interstitial clusters are important for the irradiation response of UO 2 and correlates well with the very low barriers reported in a few historic experimental studies.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2023

Up up down down left right left right B A Start for the catalytic hackers of programmable materials

Catalysts have advanced over the last century to accelerate and control reactions based on static active sites. More effective catalysis can be achieved using catalysts that change with time over the course of a reaction, providing a dynamic free energy landscape that is tailored to each step in the reaction sequence. Here, the catalyst is modulated via an input program that directs the surface to change physically or electronically with time, providing information regarding the extent and duration of change optimized for specific combinations of chemistry and programmable catalyst surfaces. While in its infancy, programmable catalysis is advancing with parallel efforts to establish fundamental principles of dynamic catalysts, design of programmable materials, and strategies to design input programs that will control catalysis for faster and more selective reactions.

catalyst↗
DOE OSTI2020

Formation and energetics of amorphous rare earth (RE) carbonates in the RE 2 O 3 –CO 2 –H 2 O system

Amorphous materials are of interest in the genesis of crystalline solids. They are often the first to precipitate from supersaturated solutions and are in the nanometer size range, transforming into a nanocrystalline metastable phase, and then to thermodynamically stable phases. One must seek knowledge of the composition, structure, thermodynamics, kinetics, and reaction pathways that relate the amorphous and crystalline phases. The subject of this study is the amorphous phases in the system RE 2 O 3 -CO 2 –H 2 O (RE- La, Nd, Dy,Yb). They were synthesized via direct precipitation and urea hydrolysis, and characterized by powder XRD, DTA with mass spectrometric analysis of the evolved gases. Phases have the non-stoichiometric composition RE 2 O 3 ·xCO 2 ·yH 2 O (1 < x < 3, 3 < y <7), which is different from that of crystalline simple carbonates and hydroxycarbonates. Thus we consider them to be amorphous precursors, rather than amorphous carbonates. High temperature oxide melt solution calorimetry in molten sodium molybdate solvent was used to derive their enthalpies of formation from oxides and elements. Increase of energetic stability per mole of RE 2 O 3 compounds occurs in the order: RE 2 O 3 → 2RE(OH)3→ 2REO(OH) → RE 2 O 2 CO 3 → Amorphous precursor → 2(REOHCO 3 ) → RE 2 (CO 3 ) 3 ·yH 2 O. We calculated the enthalpies of possible transformations of amorphous precursors to crystalline phases and conclude that thermodynamics determines the dependence of the crystallizing products on the temperature and partial pressure of CO 2 and H 2 O. Amorphous precursors are clearly intermediate in the synthesis, and ternary phases containing both H 2 O and CO 2 compete with each other in terms of thermodynamic stability. The energy landscape obtained here will allow one to directly synthesize specific products and control their functionality.

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

Suppression of Disorder in Benzamide and Thiobenzamide Crystals by Fluorine Substitution

Disorder is a common feature of molecular crystals that complicates determination of structures and can potentially affect electric and mechanical properties. Suppression of disorder is observed in otherwise severely disordered benzamide and thiobenzamide crystals by substituting hydrogen with fluorine in the ortho-position of the phenyl ring. Fluorine occupancies of 20–30% are sufficient to suppress disorder without changing the packing motif. Crystal structure prediction calculations reveal a much denser lattice energy landscape for benzamide compared to 2-fluorobenzamide, suggesting that fluorine substitution makes disorder less likely.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Lattice Disorder and Oxygen Migration Pathways in Pyrochlore and Defect-Fluorite Oxides

Atomic-scale disorder plays an important role in the chemical and physical properties of oxide materials. The structural flexibility of pyrochlore-type oxides allows for crystal-chemical engineering of these properties. Compositional modification can push pyrochlore oxides toward a disordered defect-fluorite structure with anion Frenkel pair defects that facilitate oxygen migration. The local structure of the long-range average cubic defect-fluorite was recently claimed to consist of randomly arranged orthorhombic weberite-type domains. Here, we show, using low-temperature neutron total-scattering experiments, that this is not the case for Zr-rich defect-fluorites. By analyzing data from the pyrochlore/defect-fluorite Y 2 Sn 2–x Zr x O 7 series using a combination of neutron pair distribution function and big-box modelling, we have differentiated and quantified the relationship between anion sub-lattice disorder and Frenkel defects. These details directly influence the energy landscape for oxygen migration and are crucial for simulations and design of new materials with improved properties.

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

Predicting and Synthesizing Interface Stabilized 2D Layers

The compound (Pb 2 MnSe 3 ) 0.6 VSe 2 was predicted to be kinetically stable based on density functional theory (DFT) calculations on an island of Pb 2 MnSe 3 between layers of VSe 2 . This approach provides a high degree of freedom by not forcing interlayer lattice match, making it ideal to investigate the likelihood of formation of new incommensurate layer misfit structures. The free space around the island is critical, as it allows atoms to diffuse and hence exploring the local energy landscape around the initial configuration. (Pb 2 MnSe 3 ) 0.6 VSe 2 was synthesized via a near diffusionless reaction from precursors where a repeating sequence of elemental layers matches the local composition and layer sequence of the predicted compound. The VSe 2 layer consists of a Se–V–Se trilayer with octahedral coordination of the V atoms. The Pb 2 MnSe 3 layer consists of three rock-salt-like planes, with a MnSe layer between the planes of PbSe. The center MnSe plane stabilizes the puckering of the outer PbSe layers. Electrical properties indicate that (Pb 2 Mn 1 Se 3 ) 0.6 VSe 2 undergoes a charge density wave transition at ~100 K and orders ferromagnetically at 35 K. Overall, the combination of theory and experiment enables a faster convergence to new heterostructures than either approach in isolation.

36 MATERIALS SCIENCE↗