Search NASA⌕ Search

SEARCH · Search NASA

Results for “Molecular mobility”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Molecular Dynamics Simulations of Supercritical Carbon Dioxide and Water using TraPPE and SWM4-NDP Force Fields

The increased levels of carbon dioxide (CO 2 ) emissions due to the combustion of fossil fuels and the consequential impact on global climate change have made CO 2 capture, storage, and utilization a significant area of focus for current research. In most electrochemical CO 2 applications, water is used as a proton donor due to its high availability and mobility and use as a polar solvent. Additionally, supercritical CO 2 is a promising avenue for electrochemical applications due to its unique chemical and physical properties. Consequently, understanding the interactions between water and supercritical CO 2 is of great importance for future electrochemical applications. Molecular dynamics (MD) simulation is a powerful tool that enables atomistic-resolution dynamics of molecular systems, which can complement and guide future experimental investigations. This study employed atomistic MD to study the cosolubilities, codiffusivities, and structure of supercritical CO 2 and water systems, with a polarizable water model (SWM4-NDP) and a nonpolarizable CO 2 model (TraPPE). Additionally, ab initio MD simulations were used to better understand how atomistic polarizable/nonpolarizable models compare to explicit modeling of electron densities. The polarizable water model exhibited substantial improvement in water-associated properties. In conclusion, we anticipate the development of a compatible polarizable CO 2 model to yield similar improvement, providing a pathway for realizing novel high-pressure electrochemical systems.

25 ENERGY STORAGE↗
DOE OSTI2023

Determining the influence of temperature and pressure on the structural stability in a polyurea elastomer

The structural response of amorphous elastomers to applied pressure and temperature is of significant interest given their potential for blast and impact protection applications. While understanding disordered structures is often challenging, this study focuses on the ordered correlations that are present at the atomic level even in amorphous materials. This local order was characterized in polyurea elastomers with two molecular weights up to 6 GPa and ~280 °C using in situ energy dispersive X-ray diffraction at the Advanced Photon Source via total scattering techniques. The resulting structure factors and pair distribution functions indicate that the applied pressure reduced the free volume but had little effect on covalent bonds for both molecular weights. Temperature-induced order-disorder transitions that occur readily at ambient pressures within the temperature range studied here were not observed, their onset likely delayed by the limited chain mobility associated with the free volume reduction. In conclusion, this study demonstrates that total scattering offers insights into the high-pressure structure of disordered polymers that are otherwise inaccessible using more conventional X-ray scattering techniques.

42 ENGINEERING↗
DOE OSTI2022

Hydrogen-Bond-Promoted Planar Conformation, Crystallinity, and Charge Transport in Semiconducting Diazaisoindigo Derivatives

Conformational control of π-conjugated molecules using intramolecular noncovalent bonds represents a promising strategy to tailor the solid-state molecular packing and electronic properties of these materials. Here, we report the design and synthesis of two model compounds featuring intramolecular hydrogen bonds formed between a center diazaisoindigo unit (the acceptor) and flanking indole units (the donor). Computational and experimental investigations show that these hydrogen bonds enthalpically stabilize the coplanar molecular conformation by over 10 kcal/mol. The formation of these hydrogen bonds is also slightly favorable in terms of entropy, ensuring high temperature stability of the planar conformation. In this work, thermal annealing of thin films of these compounds imparts high crystallinity and orientation in the solid state, while the non-hydrogen bond control only gave an amorphous solid. Field-effect transistor devices fabricated from these thin films exhibit hole mobilities up to 0.270 cm 2 V -1 s -1 , in contrast to the lack of measurable charge carrier mobility for the non-hydrogen bond control. This work demonstrates an efficient synthetic strategy to incorporate robust intramolecular hydrogen bonds into conjugated π-systems and elucidates the mechanism on how such hydrogen bonds promote the desired molecular conformation, solid-state packing, and electronic performances of conjugated organic materials.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Structure, Dynamics, and Hydrogen Transport in Amorphous Polymers: An Analysis of the Interplay between Free Volume Element Distribution and Local Segmental Dynamics from Molecular Dynamics Simulations

Polymers are attractive membrane materials owing to their mechanical robustness and relatively inexpensive fabrication. An important indicator of membrane performance are free volume elements (FVE): microporous void spaces created by the inefficient packing of bulky groups along the polymer chain. FVEs tend to degrade over time, as polymer chains reorganize irreversibly. While it is widely accepted that polymer flexibility has an impact on membrane transport properties, the molecular nature of this impact is still not well understood. By the establishment of a correlation between local chain dynamics and the distribution of free volume elements (FVEs), penetrant transport can be regulated more efficiently in amorphous polymer membranes. In this work, we implement all-atom molecular dynamics (MD) simulations to explore the relationship between chain dynamics and free volume in three polymers with different levels of backbone flexibility: polymethylpentene (PMP), polystyrene (PS), and HAB-6FDA thermally rearranged polymer (TRP). We construct these polymers at different temperatures and examine how temperature impacts the FVE distribution and segmental mobility. Our analysis shows that chain segments near FVEs have higher mobility compared with the atoms in the bulk; the extent of this difference increases with chain flexibility. Increasing the chain flexibility by increasing the temperature results in a broader FVE distribution. Rigid polymers such as TRP show the most robust FVE distribution and are not significantly affected by the temperature change. To capture penetrant diffusion through the polymer matrix, hydrogen is inserted and the diffusion is measured at different temperatures; hydrogen mobility is influenced by the FVE structure and overall mobility of polymer chains. Lastly, at low temperatures, hydrogen mobility is influenced by void distribution, while at high temperatures, polymer dynamics dictate hydrogen transport.

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

Effect of Molecular Conformation on Intermolecular Interactions and Photovoltaic Performances of Giant Molecule Acceptors

The molecular conformation of giant molecule acceptors (GMAs) plays a significant role in regulating the intermolecular interactions and their photovoltaic performances in organic solar cells (OSCs). For the linear GMA GT‐l, the stronger homo‐molecular interaction causes its aggregation being weakly affected by the donor, thus forming an ordered molecular stacking and proper phase separation in its blend film. The star‐shaped GMA GT‐s‐based blend film shows a dominant hetero‐molecular interaction that suppresses the aggregation of the donor and acceptor, resulting in smaller phase separation and more uniform vertical phase distribution. While for another star‐shaped GMA GTs, the weakest hetero‐molecular interaction causes its blend film to form larger phase separation. Therefore, the GT‐l based OSC with PM6 as donor shows the highest charge mobilities, the fastest charge transfer (CT) process, reduced energy loss and less charge recombination, contributing to a higher power conversion efficiency (PCE) of 19.03%. Comparatively, the PCEs of the OSCs based on GTs and GT‐s are 18.05% and 17.58% respectively. Notably, all the three GMAs based OSCs show excellent thermal stability and long‐term storage stability. This study provides a facile strategy by tuning the linking unit and its connecting mode for designing highly efficient and stable organic photovoltaic materials.

Chemistry↗
DOE OSTI2021

Multiscale Concurrent Atomistic-Continuum (CAC) modeling of multicomponent alloys

We report strengthening in complex multicomponent systems such as solid solution alloys is controlled primarily by the dynamic interactions between dislocation lines and heterogeneously distributed solute species. Modeling of extended defect length scales in such multicomponent systems becomes prohibitively expensive, motivating the development of reduced order approaches. This work explores the application of the Concurrent Atomistic-Continuum (CAC) method to model dislocation mobility in random alloys at extended length scales. By employing recently developed average-atom interatomic potentials, the average “bulk” material response in coarse-grained regions interacts with true random solute species in the atomistic-scale domain. We demonstrate that spurious stresses in domain resolution transition regions are eliminated entirely due to the CAC formulation. Simultaneously, the key details of local stress fluctuation due to randomness in the dislocation core region are captured, and fluctuating stress smoothly decays to the long-range dislocation stress field response. Dislocation mobility calculations, for line lengths over 400 nm, are computed as a function of alloy composition in the model FeNiCr system and compared to full molecular dynamics (MD). The results capture the composition-dependent trends, while reducing degrees of freedom by nearly 40%. This approach can be readily extended to any system described by an EAM potential and facilitates the study of large-scale defect dynamics in complex solute environments to support computational alloy design.

36 MATERIALS SCIENCE↗
DOE OSTI2022

High-density polarization-induced 2D electron gases in N-polar pseudomorphic undoped GaN/Al 0.85 Ga 0.15 N heterostructures on single-crystal AlN substrates

The polarization difference and band offset between Al(Ga)N and GaN induce two-dimensional (2D) free carriers in Al(Ga)N/GaN heterojunctions without any chemical doping. A high-density 2D electron gas (2DEG), analogous to the recently discovered 2D hole gas in a metal-polar structure, is predicted in a N-polar pseudomorphic GaN/Al(Ga)N heterostructure on unstrained AlN. We report the observation of such 2DEGs in N-polar undoped pseudomorphic GaN/AlGaN heterostructures on single-crystal AlN substrates by molecular beam epitaxy. With a high electron density of ~4.3 ×10 13 /cm 2 that maintains down to cryogenic temperatures and a room temperature electron mobility of ~450 cm 2 /V s, a sheet resistance as low as ~320 Ω/$\Box$ is achieved in a structure with an 8 nm GaN layer. Finally, these results indicate significant potential of AlN platform for future high-power RF electronics based on N-polar III-nitride high electron mobility transistors.

42 ENGINEERING↗
DOE OSTI2021

Heterovalent semiconductor structures and devices grown by molecular beam epitaxy

Heterovalent structures consisting of group II-VI/group III-V compound semiconductors offer attractive properties, such as a very broad range of bandgaps, large conduction band offsets, high electron and hole mobilities, and quantum-material properties such as electric-field-induced topological insulator states. These properties and characteristics are highly desirable for many electronic and optoelectronic devices as well as potential condensed-matter quantum-physics applications. Here, we provide an overview of our recent studies of the MBE growth and characterization of zincblende II-VI/III-V heterostructures as well as several novel device applications based on different sets of these materials. By combining materials with small lattice mismatch, such as ZnTe/GaSb (Δa/a~0.13%), CdTe/InSb (Δa/a~0.05%), and ZnSe/GaAs (Δa/a~0.26%), epitaxial films of excellent crystallinity were grown once the growth conditions had been optimized. Cross-sectional observations using conventional and atomic-resolution electron microscopy revealed coherent interfaces and close to defect-free heterostructures. Measurements across CdTe/InSb interfaces indicated a limited amount (~1.5 nm) of chemical intermixing. Results for ZnTe/GaSb distributed Bragg reflectors, CdTe/Mg x Cd 1-x Te double heterostructures, and CdTe/InSb two-color photodetectors are briefly presented, and the growth of a rock salt/zincblende PbTe/CdTe/InSb heterostructure is also described.

Materials Science↗
DOE OSTI2024

Enhancement of Conjugated Polymer Microstructure and Mixed-Conducting Properties via Chalcogenophene Heteroatom Substitution

Heteroatom substitution is a powerful tool to tune the intra- and intermolecular structure of conjugated polymers as well as their resulting optoelectronic and electrochemical properties. Here, a series of oligoethylene glycol bithiophene chalcogenophene polymers (p(g3T2-X)) with systematically varied furan, thiophene, selenophene, and tellurophene comonomers have been synthesized for mixed ionic-electronic conducting applications. Their microstructures have been thoroughly characterized ex situ and in situ with X-ray scattering, and their mixed conducting properties have been probed in electrochemical transistor testbeds. Chalcogenophene heteroatom choice was found to clearly dictate the polymer microstructure (crystallite dimensionality and orientation) and tune mixed conducting properties. Proceeding down Group 16, from O to Se systematically directed the molecular ordering of 2D polymer crystallites from face-on (O) to mixed (S) to edge-on (Se) orientations, with Te driving the polymer to form well-oriented edge-on 3D crystallites. Heteroatom dictated crystallite quality, and orientation tuned relative ionic transport by 2 orders of magnitude. Hole mobility (μ hole ) and mixed conducting figure of merit (μC*) were each tuned over an order of magnitude depending on heteroatom choice, with the Te-containing polymer reaching μ hole = 3.60 cm 2 V –1 s –1 and μC* = 483 F cm –1 V –1 s –1 , due to improved molecular ordering. Insights from this polymer series highlight target microstructures for enhanced mixed conduction in future conjugated polymers.

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

Increasing thermostability of the key photorespiratory enzyme glycerate 3‐kinase by structure‐based recombination

As global temperatures rise, improving crop yields will require enhancing the thermotolerance of crops. One approach for improving thermotolerance is using bioengineering to increase the thermostability of enzymes catalysing essential biological processes. Photorespiration is an essential recycling process in plants that is integral to photosynthesis and crop growth. The enzymes of photorespiration are targets for enhancing plant thermotolerance as this pathway limits carbon fixation at elevated temperatures. We explored the effects of temperature on the activity of the photorespiratory enzyme glycerate kinase (GLYK) from various organisms and the homologue from the thermophilic alga Cyanidioschyzon merolae was more thermotolerant than those from mesophilic plants, including Arabidopsis thaliana. To understand enzyme features underlying the thermotolerance of C. merolae GLYK (CmGLYK), we performed molecular dynamics simulations using AlphaFold-predicted structures, which revealed greater movement of loop regions of mesophilic plant GLYKs at higher temperatures compared to CmGLYK. Based on these simulations, hybrid proteins were produced and analysed. These hybrid enzymes contained loop regions from CmGLYK replacing the most mobile corresponding loops of AtGLYK. Two of these hybrid enzymes had enhanced thermostability, with melting temperatures increased by 6 °C. One hybrid with three grafted loops maintained higher activity at elevated temperatures. Whilst this hybrid enzyme exhibited enhanced thermostability and a similar Km for ATP compared to AtGLYK, its Km for glycerate increased threefold. This study demonstrates that molecular dynamics simulation-guided structure-based recombination offers a promising strategy for enhancing the thermostability of other plant enzymes with possible application to increasing the thermotolerance of plants under warming climates.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Engineering donor–acceptor conjugated polymers for high-performance and fast-response organic electrochemical transistors

To date, high-performance organic electrochemical transistors (OECTs) have mostly been based on polythiophene systems. Donor–acceptor (D–A) conjugated polymers are expected to be promising materials for OECTs owing to their high mobility and comparatively low crystallinity (good for ion diffusion). However, the OECT performance of D–A polymers lags far behind that of the polythiophenes. Here we synergistically engineered the backbone and side chain of a series of diketopyrrolopyrrole (DPP)-based D–A polymers and found that redox potential, molecular weight, solution processability, and film microstructures all have a severe impact on their performance. After systematic engineering, P(bgDPP-MeOT2) exhibited the best figure-of-merit (μC*) of 225 F cm -1 V -1 s -1 , amongst the highest performance of the reported D–A polymers. Besides, the DPP polymers exhibited high hole mobility of over 1.6 cm 2 V -1 s -1 , leading to fast response OECTs with a record low turn-off response time of 30 μs. The polymer also exhibited good operation stability with a current retention of 98.8% over 700 electrochemical switching cycles. This work reveals the complexity and systematicness in the development of D–A polymer based high-performance OECTs.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Toward AlGaN channel HEMTs on AlN: Polarization-induced 2DEGs in AlN/AlGaN/AlN heterostructures

Due to its high breakdown electric field, the ultra-wide bandgap semiconductor AlGaN has garnered much attention recently as a promising channel material for next-generation high electron mobility transistors (HEMTs). A comprehensive experimental study of the effects of Al composition x on the transport and structural properties is lacking. We report the charge control and transport properties of polarization-induced 2D electron gases (2DEGs) in strained AlGaN quantum well channels in molecular-beam-epitaxy-grown AlN/Al x Ga 1-x N/AlN double heterostructures by systematically varying the Al content from x = 0 (GaN) to x = 0.74, spanning energy bandgaps of the conducting HEMT channels from 3.49 to 4.9 eV measured by photoluminescence. This results in a tunable 2DEG density from 0 to 3.7 × 10 13 cm 2 . The room temperature mobilities of x ≥ 0.25 AlGaN channel HEMTs were limited by alloy disorder scattering to below 50 cm 2 /(V.s) for these 2DEG densities, leaving ample room for further heterostructure design improvements to boost mobilities. A characteristic alloy fluctuation energy of ≥1.8 eV for electron scattering in AlGaN alloy is estimated based on the temperature dependent electron transport experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Molecular shifts in dissolved organic matter along a burn severity continuum for common land cover types in the Pacific Northwest, USA

Increasing wildfire severity is of growing concern in the western United States, with consequences for the production, composition, and mobilization of dissolved organic matter (DOM) from terrestrial to aquatic systems. Our current understanding of wildfire impacted DOM (often termed pyrogenic DOM) composition is largely built from temperature-based studies that can be difficult to extrapolate to field conditions, which are often defined by ‘burn severity’, or the post-wildfire impact observed at a site. Thus, burn severity can encapsulate a broader range of fire and environmental conditions not exclusive to temperature. Biogeochemical studies that describe DOM along burn severity continuums remain limited but are needed to better link DOM composition with field conditions post-fire. Here, in this study, we addressed this need with an experimental open air burn simulation that generated chars from vegetation representative of major land cover types in the western United States. The chars were leached to simulate DOM mobilization potential. The DOM composition was characterized by ultra-high resolution mass spectrometry (HR-MS) and UV/VIS absorbance and fluorescence. Our results indicated that the shifts of DOM production and composition along a burn-severity gradient depends on the land cover type that was burned, with the degree of change dependent on the composition of the starting parent vegetation material. Fluorescence signatures indicated a strong convergence across land cover types to more aromatic DOM with increasing severity, while HR-MS indicated an increase in the production of aromatic nitrogen containing DOM with increasing severity. Results from this study enhance our ability to describe DOM composition in a framework that can be more directly related with field and remote-sensing based metrics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Interactions between molecular-scale biogeochemical processes and hyporheic exchange for understanding coupled Fe-S-C cycling in iron-rich riparian wetlands (Final Technical Report)

Riparian wetlands are dynamic interfaces that exert strong control over water quality, contaminant mobility, and greenhouse gas emissions. These systems are characterized by hyporheic exchange between oxic surface water and anoxic groundwater, which generates steep redox gradients and promotes spatially and temporally variable microbial activity. Despite growing recognition of tightly coupled iron (Fe), sulfur (S), and carbon (C) cycling in these environments, the mechanisms governing these interactions—particularly under low-sulfate freshwater conditions—remain poorly constrained. This project developed a mechanistic understanding of how hydrologic variability and microbial processes interact to control Fe–S–C cycling in iron-rich riparian wetlands. Using a multi-scale and multi-method approach integrating field observations, geochemical and spectroscopic analyses, metagenomics, and reactive transport modeling, we demonstrate that “cryptic” sulfur cycling—rapid sulfur transformations involving intermediate-valence species—plays a dominant and previously underrecognized role in freshwater wetlands. These processes persist despite low sulfate concentrations and significantly influence iron reduction, carbon mineralization, and methane dynamics. The results show that cryptic sulfur cycling enhances dissolved Fe 2+ production, regulates methane concentrations, and is strongly controlled by climate-driven hydrologic fluxes. By linking hydroclimate, subsurface flow, and biogeochemical reactions, this work provides a predictive framework for understanding how wetland systems respond to environmental change, with direct implications for water quality and carbon cycling.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Effect of the chlorine substitution position of the end-group on intermolecular interactions and photovoltaic performance of small molecule acceptors

The structure–property relationships of small molecular acceptors (SMAs) are a key issue in the molecular design of new-generation acceptor materials for further improving the device efficiencies of polymer solar cells (PSCs). Herein, three couples of SMA isomers were synthesized, based on three central fused ring units and two 1,1-dicyanomethylene-3-indanone (IC) isomer electron-withdrawing terminal units with chlorine substitutions in different positions of its benzene ring: Cl-1 with chlorine on the same side as a C=O group of IC and Cl-2 with chlorine on the same side as the CN groups of IC. Through systematical investigation, we found that the chlorine substitution position of the terminal groups has a regular and significant influence on the molecular packing and photovoltaic performance of the SMAs. The molecular packing behavior of the SMAs is closely related to and determined by the configuration of their terminal groups, no matter which central fused ring of the SMAs is used. In particular, the Cl-1-based SMAs possess a stronger crystallinity with long range ordering packing in their molecular plane direction, while the more abundant and stereoscopic π–π intermolecular interaction in the Cl-2-based SMAs promotes the molecules to form three-dimensional charge transporting channels and leads to their red-shifted absorption and higher electron mobilities. Therefore, the Cl-2-based PSCs exhibit a higher power conversion efficiency (PCE) compared to that of the Cl-1-based devices, and the best PCE of a Cl-2 SMA-based PSC reached 16.42%. These results highlight the importance of the investigation of intermolecular interactions, packing and the arrangement of the SMAs in the solid-state, which may provide direct insights for exploring the relationship between the molecular structure and property of the photovoltaic materials. Moreover, we envision that if fragments such as end groups or side chains with more diverse molecular interactions are added into the design and the subsequent synthesis of the SMAs, this may be beneficial to promoting molecular π–π accumulation and further improving the molecular order, forming suitable molecular packing and morphology in the resulting blend films, and finally affecting the efficiency of the PSCs.

14 SOLAR ENERGY↗
DOE OSTI2021

Stability and mobility of tungsten clusters on tungsten (110) surface: Ab initio and atomistic simulations

Quantifying the surface transport properties of tungsten (W) is of prime importance to understand the formation of nano-fuzz in fusion plasma-facing conditions. The stability and mobility of W adatom clusters (W n , n = 2-9) on the W(110) surface has been investigated by computer simulations, including ab initio calculations using density functional theory (DFT) and molecular statics (MS) simulations with multiple W interatomic potentials. The DFT results demonstrate that the sequential binding energy generally increases with number of W adatoms, except for the 5 th and 7 th W adatoms. The most common elemental migration steps of W n (n>2) clusters are observed to consist of monomer and dimer hops, while larger W n clusters can also diffuse by dissociation and recombination of smaller clusters. The threshold migration energy of W 9 is the highest, then followed by W 8 , W 4 , and W 6 , while W 3 , W 5 , and W 7 have similar migration energies. Compared to DFT, each interatomic potential evaluated overestimates the binding energies of W n clusters. Additionally, the embedded-atom potential developed by Juslin and Wirth adequately predicts the threshold migration energy of W n (n>2) clusters, although it predicts different underlying migration mechanisms. The results show that interaction mechanism between W adatoms controls the stability and mobility of W n clusters on the W(110) surface.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Effect of Helium Flux on Near-surface Helium Accumulation in Plasma-Exposed Tungsten

We report results of object kinetic Monte Carlo (OKMC) simulations aimed at understanding the effect of helium flux on the near-surface helium accumulation in plasma-facing tungsten, which is initially defect-free and has a W(100) surface orientation. These OKMC simulations are performed at 933 K for fluxes ranging from 10 22 to 4 × 10 25 He/m 2 s, with 100 eV helium atoms impinging on a W(100) surface up to a maximum fluence of 4×10 19 He/m 2 . In the near-surface region, helium clusters interact elastically with the free surface. The interaction is attractive and results in the drift of mobile helium clusters towards the surface as well as increased trap mutation rates. The associated kinetics and energetics of the above-mentioned processes obtained from molecular dynamics simulations are also considered. The OKMC simulations indicate that as the flux decreases, the retention of implanted helium decreases, and its depth distribution shifts to deeper below the surface in initially pristine tungsten. Furthermore, the fraction of retained helium diffusing into the bulk increases as well, so much so that for 10 22 He/m 2 s, almost all of the retained helium diffused into the bulk with minimal/negligible near-surface helium accumulation. At a given flux, with increasing fluence, the fraction of retained helium initially decreases and then starts to increase after reaching a minimum. The occurrence of the retention minimum shifts to higher fluences as the flux decreases. Although the near-surface helium accumulation spreads deeper into the material with decreasing flux and increasing fluence, the spread appears to saturate at depths between 80 and 100 nm. Finally, we present a detailed analysis of the influence of helium flux on the size and depth distribution of total helium and helium bubbles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

The mobility of small, over-pressurized helium bubbles in tungsten at 2000 K

Fusion reactor environments inevitably lead to the formation of high-pressure helium bubbles whose nucleation, growth, and diffusion strongly impact the performance of plasma-facing components. Here this research describes a diffusion mechanism of over-pressurized bubbles via a sequence of Frenkel pair nucleation, self-interstitial migration, and Frenkel pair annihilation. Molecular dynamics was used to simulate the diffusion of small bubbles in tungsten at 2000 K with helium-per-vacancy ratios in the range of 4.5 to 7. The diffusion coefficients are calculated and their dependence on helium content, number of vacancies, and number of attached self-interstitials is characterized. It is found that bubbles are most mobile when the nucleation/annihilation rates of Frenkel pairs are nearly equal and when the bubbles nucleate and annihilate a single self-interstitial. All bubbles experience a peak diffusivity, which can be as high as 10 -11 m 2 /s decreasing with bubble size. The calculated diffusion coefficients provide valuable insight into the mobility of small, high-pressure bubbles, and can be used as input parameters in mesoscale models to improve predictions of plasma-surface interactions. (LA-UR-21-21881)

36 MATERIALS SCIENCE↗