Search NASA⌕ Search

SEARCH · Search NASA

Results for “Mixing thermodynamics”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Thermodynamic modeling of CsF with LiF-NaF-KF for molten fluoride-fueled reactors

Gibbs energy models were developed to describe the thermochemical behavior of CsF in molten FLiNaK (46.5LiF-11.5NaF-42KF mol%), a proposed molten salt reactor (MSR) fuel solvent and coolant, as cesium is of concern due to its high radiotoxicity and volatility. Initially, it was necessary to obtain a more accurate Gibbs energy function for CsF which required fitting parameters to reported vapor pressures over condensed phase CsF. The pseudo-binary systems CsF-LiF, CsF-NaF and CsF-KF were then evaluated utilizing phase equilibria and enthalpy of mixing (Δ mix H) values, together with original differential scanning calorimetry (DSC) measurements performed for the CsF-LiF and CsF-KF systems. The CsF-LiF-NaF, CsF-LiF-KF and CsF-NaF-KF pseudo-ternary system representations were obtained by interpolation of the constituent pseudo-binary systems, with DSC measurements performed for the CsF-LiF-NaF system to corroborate the calculated liquidus temperature. Ultimately, the pseudo-ternary systems were interpolated to obtain Gibbs energy models for the pseudo-quaternary CsF-LiF-NaF-KF system, supported by DSC measurements at low CsF compositions (1–10 mol%), yielding computed equilibria and cesium-containing vapor pressures that compare favorably with reported values. In conclusion, the Molten Salt Thermal Properties Database – Thermochemical (MSTDB-TC) was subsequently expanded to include these Gibbs energy models allowing description of the thermochemical behavior of the CsF-LiF-NaF-KF system.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Stabilization of Miscible Aqueous Phases via Diffusion‐Controlled Multifunctional Nanoparticle‐Ligand Complexation

Liquid‐in‐liquid structuring by harnessing miscible aqueous domains can be achieved by inducing thermodynamically defined phase separation using aqueous solutions of incompatible solutes, where immiscibility is dictated by the concentrations of the solutes. For instance, in aqueous two‐phase systems (ATPSs), the interfacial tension between the two different aqueous solutions is so small that it stabilizes the phase‐separated domains, but not their shape, thus falling short of designing tunable, robust structures. Here, we introduce a diffusion‐controlled strategy that enables liquid‐in‐liquid compartmentalization absent a thermodynamically defined interface or bulk liquid‐liquid phase separation using a barrier formed in situ at the initial contact boundary of an aqueous dispersion of multifunctional nanoparticles and a solution of multifunctional ligands in pure water or a water/alcohol mixture. The nanoparticles and ligands are entropically driven to disperse uniformly but interact at the initial contact boundary more rapidly than bulk mixing, forming a layer that establishes a barrier akin to a semi‐permeable membrane and allows the formation of tubular compartments. A wealth of opportunities is enabled by the proposed technique, where examples in separations, soft electronic devices, and ultralight porous electromagnetic shielding materials are shown.

cryogels↗
DOE OSTI2022

Comparison of Na 2 SO 4 , K 2 SO 4 and Na 2 SO 4 -K 2 SO 4 deposit induced hot corrosion of a β-NiAl coating

Here, hot corrosion behaviour of a β-NiAl coating induced by Na 2 SO 4 , K 2 SO 4 , and Na 2 SO 4 -K 2 SO 4 deposits was investigated at 700 °C in air + 300 ppm SO 2 . To provide an insight into the corrosion mechanisms, phase equilibria of the Na 2 SO 4 -K 2 SO 4 -NiO-SO 3 system were calculated using an in-house developed thermodynamic database. The attack rates and morphologies could be correlated with the deposit compositions and their reactions with oxides. K 2 SO 4 and Na 2 SO 4 -K 2 SO 4 deposits are more corrosive than Na 2 SO 4 . This is attributed to a reduction of the melting temperature of the deposit and the formation of Ni-containing mixed sulphate phases in the presence of K 2 SO 4 .

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

Methods and apparatus to facilitate alkali metal transport during battery cycling, and batteries incorporating same

An anode includes a mixed ionic-electronic conductor (MIEC) with an open pore structure. The open pore structure includes open pores to facilitate motion of an alkali metal into and/or out of the MIEC. The open pore structure thus provides open space to relieve the stresses generated by the alkali metal when charging/discharging a battery. The MIEC is formed from a material that is thermodynamically and electrochemically stable against the alkali metal to prevent the formation of solid-electrolyte interphase (SEI) debris and the formation of dead alkali metal. The MIEC may also be passive (the MIEC does not store or release alkali metal). In one example, the open pore structure may be an array of substantially aligned tubules with a width less than about 300 nm, a wall thickness between about 1 nm to about 30 nm, and a height of at least 10 um arranged as a honeycomb.

Li, Ju↗
DOE OSTI2021

Adsorbate Free Energies from DFT-Derived Translational Energy Landscapes

Adsorption free energies are fundamental to surface chemistry and catalysis. Standard models combine some assumed analytical form of the translational potential energy surface, often parameterized against density functional theory (DFT) calculations, with an analytical expression for the resultant translational densities of states (DOS), free energy, and entropy. Here we compare the performance of such models against numerical evaluations of the DOS and thermodynamic functions derived from solutions to the translational Schrodinger equation. We compare results for a translational potential energy surface (PES) derived from nudged eleastic band calculations with those obtained from adsorbate rastering across a series of monatomic (O, S, C, N, H) and polyatomic (NH x ) adsorbates on (100) Pt and Au facets. We find that analytical models as commonly parameterized have mixed performance for describing the translational PES and that the consequences for computed free energies are modest but potentially significant in microkinetic models. Here, numerical solutions are possible for modest to no additional computational cost over analytical models and thus should be considered when reliable free energy estimates are needed or translational potential energy surfaces are available.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Thermodynamic controls on the synthesis, structure and reactivity of materials for energy (Final Report)

The major objectives have been: a) to advance and use unique calorimetric capabilities to determine the energetics of oxide and other materials having possible energy applications, with an emphasis on new materials for batteries, fuel cells and solar energy conversion, b) to understand, in terms of both macroscopic energetics and microscopic structure and bonding, the interplay of defect chemistry, oxidation-reduction, and size effects at the bulk and nanoscale in determining materials properties. This project emphasizes mixed ionic-electronic conduction materials containing transition metals of variable valence, especially perovskites and lithium containing phases, with M = Mn, Fe, Co, Ni. The thermochemistry of new families of complex halide, chalcogenide, nitride and hybrid (organic-inorganic) perovskite materials were explored. The major and unique technique employed is high temperature oxide melt solution calorimetry, with careful attention given to materials synthesis and characterization. The thermochemical studies will be coupled with structural studies, especially to characterize order-disorder phenomena by in situ high temperature diffraction techniques.

14 SOLAR ENERGY↗
DOE OSTI2025

Insights into the Oxidation Mechanism of Vivianite to Metavivianite from First-Principles Calculations

Vivianite, a hydrous ferrous iron-phosphate mineral (Fe 3 (PO 4 ) 2 ·8H 2 O), readily oxidizes in contact with air yielding the less hydrous mixed-valent iron-phosphate mineral metavivianite. This topotactic transformation nominally occurs by oxidative dehydrogenation in which outgoing electrons from the iron sublattice are charge compensated by hydrolysis of structural water. However, the details of this internal charge balancing mechanism that allows the structure to remain electrostatically stable remain unknown. Here, in this study, we use density functional theory (DFT) calculations and ab initio thermodynamics (AIT) to evaluate the energetics of this process in terms of hydrogen release as a function of environmental variables such as partial pressure and temperature. The results show a thermodynamic driving force for vivianite phase transformation as oxidation progresses that is triggered by its rigid structure that has a limited accommodation for hydrogen vacancies. In contrast, metavivianite has a more flexible lattice and hydrogen bond network that stabilizes these hydrogen defects. Metavivianite is shown to be a stable intermediate for 66% residual Fe 2+ down to 33% Fe 2+ , below which other phases/structures such as santabarbarite should be more thermodynamically favorable. Our study provides a basis for experimental tests of our mechanistic findings and helps fill a basic knowledge gap about the solid-state process that defines how vivianite interacts with its surrounding environment.

Sassi, Michel [Pacific Northwest National Laborato↗
DOE OSTI2019

Coexistence of static and dynamic magnetism in the Kitaev spin liquid material Cu 2 IrO 3

Searching for a Kitaev spin liquid phase motivated intense research on the honeycomb iridate materials. However, access to a spin liquid ground state has been hindered by magnetic ordering. Cu 2 IrO 3 is a new honeycomb iridate without thermodynamic signatures of a long-range order. Here, we use muon spin relaxation to uncover the magnetic ground state of Cu 2 IrO 3 . We find a two-component depolarization with slow- and fast-relaxation rates corresponding to distinct regions with dynamic and static magnetism coexisting in the material. X-ray absorption spectroscopy and first-principles calculations identify a mixed copper valence as the origin of this behavior. Furthermore, our results suggest that a minority of Cu 2+ ions nucleate regions of static magnetism, whereas the majority of Cu + /Ir 4+ on the honeycomb lattice give rise to a Kitaev spin liquid.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Local cation order and ferrimagnetism in compositionally complex spinel ferrites

We present an exploration of a family of compositionally complex cubic spinel ferrites featuring combinations of Mg, Fe, Co, Ni, Cu, Mn, and Zn cations, systematically investigating the average and local atomic structures, chemical short-range order, magnetic spin configurations, and magnetic properties. All compositions result in ferrimagnetic average structures with extremely similar local bonding environments; however, the samples display varying degrees of cation inversion and, therefore, differing apparent bulk magnetization. Additionally, first-order reversal curve analysis of the magnetic reversal behavior indicates varying degrees of magnetic ordering and interactions, including potentially local frustration. Finally, reverse Monte Carlo modeling of the spin orientation demonstrates a relationship between the degree of cation inversion and the spin collinearity. Collectively, these observations correlate with differences in synthesis procedures. This work provides a framework for understanding magnetic behavior reported for “high-entropy spinels,” revealing many are likely compositionally complex oxides with differing degrees of chemical short-range order—not meeting the community established criteria for high or medium entropy compounds. Moreover, this work highlights the importance of reporting complete sample processing histories and investigating local to long-range atomic arrangements when evaluating potential entropic mixing effects and assumed property correlations in high entropy materials.

36 MATERIALS SCIENCE↗
DOE OSTI2021

High-pressure phase diagram of beryllium from ab initio free-energy calculations

In this report we use first-principles molecular dynamics simulations coupled with the thermodynamic integration method to study the hexagonal close-packed (hcp) to body-centered cubic (bcc) transition and melting of beryllium up to a pressure of 1600 GPa. We derive the melting line by equating solid and liquid Gibbs free energies and represent it by a Simon-Glatzel fit T m = 1564 K [1 + P/(15.6032 GPa)] 0.383 , which is in good agreement with previous two-phase simulations <6000 K. We also derive the hcp-bcc solid-solid phase boundary and show that the quasiharmonic approximation underestimates the stability of the hcp structure, predicting lower transition pressures between hcp and bcc phases. Our results are consistent with the stability regime predicted by the phonon quasiparticle method. We also predict that the hcp-bcc-liquid triple point is located at 164.7 GPa and 4314 K. In addition, we compute the shock Hugoniot curve and show that it is in good agreement with experiments, intersecting our derived melting curve at ~235 GPa and 4900 K. Finally, we make predictions for future ramp compression experiments. Starting with an isentropic compression of the liquid, we predict the path to intersect the melting line at low pressure and temperature, then to continue along the melting line over a large temperature interval of 7000 K as the sample remains in the mixed solid-liquid state before it enters the solid phase.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Vortex Structure of Head Bubble in Convective Cloud Starting Plume

Abstract A developing cumulus cloud (Cu) was modeled, and dynamic, thermodynamic, and microphysical properties of an ascending head bubble reproducing the upper part of a developing Cu were investigated. The data for analysis are taken from 10-m-resolution LES of trade wind Cu under BOMEX conditions. The detection of a rising bubble is carried out using wavelet filtering of the velocity fields and microphysical fields, while a low-frequency signal of the filtering is associated with the convective-scale structure of cloud. We substantiate and discuss the representation of the bubble as a vortex ring, and estimate the parameters of this vortex ring. The simplest Hill’s vortex was chosen as a model of a vortex ring inside cloud. Analytical approximations of the radial profiles of the vertical velocity and of conservative quantities (such as total water mixing ratio and liquid water potential temperature inside and outside the bubble) are obtained. The spatial structure of these quantities is investigated using analytical expressions. Analytical models for spatial distributions of liquid water content (LWC) and adiabatic fraction (AF) are also designed and analyzed. The results demonstrate the existence of a cloud core with high values of LWC and AF up to the height of 1800 m. The horizontally averaged value of the adiabatic fraction, calculated analytically using the Hill’s vortex concept, is evaluated as 0.39, which is the typical AF value in the upper parts of such Cu. The vertical profiles of different important quantities characterizing cloud structure are presented. The analysis performed in this study allows us to conclude that a rising vortex ring plays the dominating role in formation of the thermodynamic and microphysical structure of developing Cu. Significance Statement 1) Dynamic and thermodynamic fields of a developing cumulus cloud simulated by high-resolution LES with spectral bin microphysics are separated into convective and turbulent components by means of the wavelet technique. 2) The analysis of convective component of the cloud revealed the existence of a vortex ring at the developing stage of the cloud and evaluate its parameters. 3) The analysis performed in this study allows us to conclude that a rising vortex ring plays the important role in formation of the thermodynamic and microphysical structure of developing Cu. 4) The study provides a novel insight into the cloud–environment interaction. 5) The approximating equations describing the vortex ring can be usefully applied for developing new schemes of convective parameterization.

Meteorology & Atmospheric Sciences↗
DOE OSTI2025

Atomic-scale understanding of oxide growth and dissolution kinetics of Ni-Cr alloys

Aqueous corrosion of metals is governed by formation and dissolution of a passivating, multi-component surface oxide. Unfortunately, a detailed atomistic description is challenging due to the compositional complexity and the need to consider multiple kinetic factors simultaneously. To this end, we combine experiments with a first-principles-derived, multiscale computational framework that transcends thermodynamic descriptions to explicitly simulate the kinetic evolution of surface oxides of Ni-Cr alloys as a function of composition, temperature, pH, and applied voltage. In the absence of pitting, we identify three distinct voltage regimes, which are kinetically dominated by oxide growth, dissolution, and competitive dissolution and reprecipitation. Evolving compositional gradients and oxide thickness are revealed, including a transition between a metastable Ni-Cr mixed oxide and a thick, porous Ni-dominated oxide. Beyond elucidating the underlying physics, we highlight the need for competing kinetics in models to properly predict the transition from passivation to corrosion. Our results provide a key step towards co-design of alloy composition alongside environmental conditions for sustainable use across a variety of critical energy and infrastructure applications.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Solventless, Ambient-Pressure Production of Bio-Based Lactones Over Earth-Abundant, Mixed Metal Oxide Catalysts for Circular Polyesters

Transitioning to a circular plastics economy will require use of renewable feedstocks, energy-efficient processes, and closed-loop recyclable polymers, such as polyesters. A key challenge lies in sustainably sourcing monomers used to make recyclable polyesters. This work presents a catalytic platform utilizing earth-abundant Cu(x)Ca(1-x)O mixed metal oxides for the oxidative dehydrocyclization of bio-based diols to lactones, which are advantaged for energy-efficient ring-opening polymerization. Operating below 200 degrees C, at ambient pressure, and without solvent, the process uses air as the sole oxidant, achieving high yields of lactones across a broad substrate scope of C4-8 diols in the liquid phase. The oxidative dehydrocyclization reaction is thermodynamically downhill due to water formation and energy-efficient compared to incumbent, non-redox pathways utilized in fossil carbon-based industrial processes for lactone production. Mechanistic studies reveal facile redox cycling of Cu2+-O(Ca2+)-Cu2+ interfacial sites unique to the developed catalyst. Techno-economic analysis and life cycle assessment estimate 40% lower energy demand and 15% lower GHG intensity per mass of butyrolactone produced compared to the fossil carbon-based route. Liquid-phase oxidative dehydrocyclization offers a promising approach for scalable lactone production from renewable, bio-based diols to enable circular polyesters.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A generalized and adaptable tensor-contraction-based cluster expansion formalism for multicomponent solids

Density functional theory (DFT)-based simulations of materials have first-principles accuracy, but are very computationally expensive. For simulating various properties of multi-component alloys, the cluster expansion (CE) technique has served as the standard workaround to improve computational efficiency. However, the standard CE technique is difficult to extend to exotic and/or low-symmetry lattices, often implemented via iteration over particular cluster types, which must be enumerated per lattice structure. In this work, we introduce the tensor cluster expansion (TCE), implemented in the open-source code tce-lib, which maps correlation functions to mixed tensor contractions, eliminating the need to iterate over cluster types and additionally making the calculation of correlation functions well-suited for massively parallel architectures like GPUs. We show that local interaction energies are an immediate consequence of the TCE formalism, yielding nearly $\mathcal{O}$(1) energy difference calculations. We then use this formalism to fit CE models for the TaW and CoNiCrFeMn systems, and use these models to respectively compute the enthalpy of mixing curve and Cowley short-range order parameters, showing excellent agreement with ground truth data.

Cluster expansion↗
DOE OSTI2024

Exploiting Saturation Regimes and Surface Effects to Tune Composite Design: Single Platelet Nanocomposites of Peptoid Nanosheets and CaCO 3

Mineral-polymer composites found in nature exhibit exceptional structural properties essential to their function, and transferring these attributes to the synthetic design of functional materials holds promise across various sectors. Biomimetic fabrication of nanocomposites introduces new pathways for advanced material design and explores biomineralization strategies. This study presents a novel approach for producing single platelet nanocomposites composed of CaCO 3 and biomimetic peptoid (N-substituted glycines) polymers, akin to the bricks found in the brick-and-mortar structure of nacre, the inner layer of certain mollusc shells. The significant aspect of the proposed strategy is the use of organic peptoid nanosheets as the scaffolds for brick formation, along with their controlled mineralization in solution. Here, we employ the B28 peptoid nanosheet as a scaffold, which readily forms free-floating zwitterionic bilayers in aqueous solution. The peptoid nanosheets were mineralized under consistent initial conditions (σ calcite = 1.2, pH 9.00), with variations in mixing conditions and supersaturation profiles over time aimed at controlling the final product. Nanosheets were mineralized in both feedback control experiments, where supersaturation was continuously replenished by titrant addition and in batch experiments without a feedback loop. Complete coverage of the nanosheet surface by amorphous calcium carbonate was achieved under specific conditions with feedback control mineralization, whereas vaterite was the primary CaCO 3 phase observed after batch experiments. Thermodynamic calculations suggest that time-dependent supersaturation profiles as well as the spatial distribution of supersaturation are effective controls for tuning the mineralization extent and product. We anticipate that the control strategies outlined in this work can serve as a foundation for the advanced and scalable fabrication of nanocomposites as building blocks for nacre-mimetic and functional materials.

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

Disorder-induced magnetoelastic behaviors of MnTexSbyBi1-x-y alloys

This dataset contains input and output files from density functional theory (DFT) simulations used to study the disorder-induced magnetoelastic behaviors of MnTexSbyBi1-x-y (0 ≤ x + y ≤ 1) alloys and their binary end members MnTe, MnSb, and MnBi. The alloys adopt the hexagonal NiAs-type (nickeline) structure and span ternary (MnTexSb1-x, MnTexBi1-x, MnBixSb1-x), and quaternary compositions across the full MnTe–MnSb–MnBi composition triangle. For each alloy composition, the dataset provides DFT calculations in three magnetic configurations: A-type antiferromagnetic (AFM), C-type AFM, and ferromagnetic (FM). Every magnetic configuration folder contains the fully relaxed crystal structure (CONTCAR), VASP input parameters (INCAR), and the main VASP output file (OUTCAR), from which total electronic energies, Mn magnetic moments, lattice parameters, and percent volume changes between magnetic states are extracted. These data are used to construct compositional phase diagrams, evaluate thermodynamic stability (formability), and map magnetoelastic responses across the alloy space. For A-type AFM and FM configurations, additional data are provided as follows: (i) FORCE_CONSTANTS and thermal_properties.yaml files at the top level of A-type_AFM/ and FM/ folders — present only for compositions marked with an asterisk (*) in Table I of the main text. These are derived from Phonopy finite-displacement calculations on full disordered 128-atom supercells and provide vibrational free energy, entropy (Svib)contribution from explicit disorder calculations. (Table I of the associated main manuscript) (ii) A VCA/ subfolder within A-type_AFM/ and FM/, containing FORCE_CONSTANTS and thermal_properties.yaml from Virtual Crystal Approximation phonon calculations (without spin-orbit coupling). VCA data are available for all compositions and are used to estimate vibrational contributions to the Gibbs free energy across the full composition space. (iii) A SOC/ subfolder containing CONTCAR, INCAR, and OUTCAR from spin-orbit coupling calculations, providing relativistic corrections to electronic energies and lattice parameters (Tables S2–S3 of the SM, and Table I of the main manuscript). (iv) A SOC/VCA/ subfolder containing FORCE_CONSTANTS and thermal_properties.yaml from VCA phonon calculations performed within the SOC framework, combining relativistic and vibrational thermodynamic corrections. The computed properties are used to map the AFM–FM magnetic crossover near MnTe0.75Sb0.25, demonstrate disorder- and spin-induced phonon broadening, identify a semiconductor-to-metal crossover, and quantify the pronounced magnetoelastic volume response near the magnetic phase boundary.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Hole Trapping in Halide Perovskites Induces Phase Segregation

Metal halide perovskites have garnered a great deal of attention for their applications in photovoltaics, LEDs, and radiation detection. The ease of solution processing high-quality perovskite semiconductors with large absorption coefficients and tolerance to native defects is decidedly attractive. Additionally, the ability to precisely tune the band gap of halide perovskites through compositional alloying of the halide ion is of particular interest for a range of applications, especially for tandem solar cells. However, under steady state light irradiation, an initially homogeneous mixed halide perovskite (MHP) will form local domains that are rich in one halide ion (e.g., Br or I). This light-induced phase segregation in MHPs forms iodide-rich domains that act as charge carrier traps and lowers the efficiency of perovskite-based devices. Thus, phase segregation poses a serious challenge to the implementation of MHPs in real-world device settings. Interestingly, when a phase segregated MHP film is placed in the dark, entropic driving forces become dominant and the segregated perovskite remixes and returns to its initially homogeneous state. Several key mechanistic details of phase segregation have been elucidated over the years. However, there are still aspects of halide segregation that are not clear, and there is ongoing debate in the literature as to what are the key factors that contribute to the mechanism. This Account discusses recent results that point to the specific role of hole trapping in phase segregation. Interestingly, generation of holes through above-band-gap excitation or through electrochemical injection increases ion migration and leads to phase segregation. The thermodynamic and redox properties of halide perovskites provide a strong driving force for hole trapping and oxidation of iodide species in MHPs. However, mobile halide species within the perovskite lattice take time to migrate and generate halide-rich domains. When in contact with a nonpolar solvent, the migration of iodine species is further extended to expulsion of iodine from the perovskite film. Thus, the mobility of halides and their susceptibility to hole-induced oxidation play a crucial role in determining the long-term stability of metal halide perovskites. Strategies to gain kinetic control over ion migration to slow phase segregation are needed to overcome these hurdles and achieve stable mixed halide perovskites. Modification of the perovskite composition through introduction of different cations or halide ions, or introduction of low-dimensional perovskite phases may suppress phase segregation. Furthermore, in achieving stability and improving the efficiency of perovskite solar cells and light emitting devices with minimal impacts, suppression of segregation remains the key factor.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Changes in dislocation punching behavior due to hydrogen-seeded helium bubble growth in tungsten

The accumulation of gas atoms in tungsten is a topic of long-standing interest to the plasma-facing materials community due the metal's use as a divertor material in some tokamak fusion reactors. The nucleation and growth of He/H gas bubbles (along with their isotopes) can result from impinging fluxes of these gases which give rise to damage at the W divertor surface. The inclusion of He or H in W has been studied extensively by the community, finding that He bubbles modify the surface through periodic dislocation punching and bursting mechanisms while H bubbles impact the metal through plastic-strain induced material failure. However, the mechanisms which are present during the combined flux of both He and H is not well-studied atomistically. Motivated by this, an atomistic modeling study is conducted using molecular dynamics to assess the behavior of mixed concentration He:H bubbles in W. Here we find that the introduction of H into a growing He bubble results in a dramatic change in the nature and presence of dislocation loops which are typically generated via dislocation punching in over-pressurized He bubbles. Most notably, at high H concentrations, there is a switchover in energetic favorability from glissile 1/2<111> dislocations to sessile <100> dislocations. This thermodynamic crossover could imply a significant reduction in W surface morphology changes than with pure He bubbles and, additionally, we show this to have implications on the trapping of H in the bubbles and their associated dislocations.

36 MATERIALS SCIENCE↗