Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “effective diffusivity coefficient”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Insights into Rotational and Translational Dynamics in Mixtures of Ethylene Glycol and Choline Chloride Using Nuclear Magnetic Resonance Techniques

This work examines molecular dynamics and interactions in ethylene glycol–choline chloride (EG–ChCl) mixtures across 0–33 mol % ChCl, spanning the true eutectic region near 17–20 mol % and the commonly used 1:2 formulation. We combine pulsed-field-gradient (PFG) diffusion, fast-field-cycling (FFC) relaxometry, temperature-dependent 13 C T 1 , and nuclear Overhauser effect spectroscopy (NOESY) to disentangle local from macroscopic dynamics. PFG and FFC show that both translational and average rotational motions largely track the strong increase in viscosity with ChCl content, with ethylene glycol consistently diffusing faster than the choline cation and no global dynamical anomaly at the eutectic composition. More subtle, site-specific composition effects nevertheless emerge. The ratio of the diffusion coefficient of the hydroxyl group of choline to the diffusion coefficient of the methyl group of choline displays a shallow minimum in the 17–25 mol % region, indicating a modest change in how the hydroxyl-bearing end of choline samples the underlying translational motion relative to the methyl groups. 13 C T 1 analysis shows that rotational correlation times at 25 °C generally increase with ChCl, reflecting viscosity-coupled slowing, while the CH 2 –N α site exhibits a small but reproducible deviation from this monotonic trend near the eutectic. NOESY spectra at similar compositions reveal enhanced cross-relaxation between EG and choline protons, consistent with increased headgroup–solvent contact density rather than a wholesale structural rearrangement. Overall, our multitechnique study demonstrates that EG–ChCl dynamics are predominantly viscosity-dominated, with the eutectic region acting as a subtle dynamical crossover where specific choline segments become maximally coupled to the hydrogen-bond network. These insights refine the structure–dynamics picture of choline-chloride DESs and provide practical guidance for tuning composition in electrochemical, separation, and catalytic applications.

diffusion↗
DOE OSTI2021

Mass accommodation and gas–particle partitioning in secondary organic aerosols: dependence on diffusivity, volatility, particle-phase reactions, and penetration depth

Abstract. Mass accommodation is an essential process for gas–particle partitioning of organic compounds in secondary organic aerosols (SOA). The mass accommodation coefficient is commonly described as the probability of a gas molecule colliding with the surface to enter the particle phase. It is often applied, however, without specifying if and how deep a molecule has to penetrate beneath the surface to be regarded as being incorporated into the condensed phase (adsorption vs. absorption). While this aspect is usually not critical for liquid particles with rapid surface–bulk exchange, it can be important for viscous semi-solid or glassy solid particles to distinguish and resolve the kinetics of accommodation at the surface, transfer across the gas–particle interface, and further transport into the particle bulk. For this purpose, we introduce a novel parameter: an effective mass accommodation coefficient αeff that depends on penetration depth and is a function of surface accommodation coefficient, volatility, bulk diffusivity, and particle-phase reaction rate coefficient. Application of αeff in the traditional Fuchs–Sutugin approximation of mass-transport kinetics at the gas–particle interface yields SOA partitioning results that are consistent with a detailed kinetic multilayer model (kinetic multilayer model of gas–particle interactions in aerosols and clouds, KM-GAP; Shiraiwa et al., 2012) and two-film model solutions (Model for Simulating Aerosol Interactions and Chemistry, MOSAIC; Zaveri et al., 2014) but deviate substantially from earlier modeling approaches not considering the influence of penetration depth and related parameters. For highly viscous or semi-solid particles, we show that the effective mass accommodation coefficient remains similar to the surface accommodation coefficient in the case of low-volatility compounds, whereas it can decrease by several orders of magnitude in the case of semi-volatile compounds. Such effects can explain apparent inconsistencies between earlier studies deriving mass accommodation coefficients from experimental data or from molecular dynamics simulations. Our findings challenge the approach of traditional SOA models using the Fuchs–Sutugin approximation of mass transfer kinetics with a fixed mass accommodation coefficient, regardless of particle phase state and penetration depth. The effective mass accommodation coefficient introduced in this study provides an efficient new way of accounting for the influence of volatility, diffusivity, and particle-phase reactions on SOA partitioning in process models as well as in regional and global air quality models. While kinetic limitations may not be critical for partitioning into liquid SOA particles in the planetary boundary layer (PBL), the effects are likely important for amorphous semi-solid or glassy SOA in the free and upper troposphere (FT–UT) as well as in the PBL at low relative humidity and low temperature.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Ion Size Effects on the Thermodynamic, Kinetic, and Mechanical Properties during Ion Exchange in Solid-State Electrolytes

Ion exchange offers a pathway to impose residual compressive stresses to mitigate the electro-chemo-mechanical cracking of solid-state electrolytes such as lithium lanthanum zirconium oxide. This study uses a coupled multiscale framework (integrating density functional theory (DFT), molecular dynamics (MD), and continuum modeling) to examine how exchange ion size influences stress, diffusion, fracture toughness, and electronic properties. Larger isovalent ions (Na + , Ag + , K + ) were exchanged with Li + , with DFT confirming their preference for octahedral 96h sites and a linear relationship between ion size and chemical free expansion coefficient. MD simulations reveal stress and concentration effects on exchange ion diffusivity at elevated temperatures, with Na + and Ag + maintaining favorable mobility while K + showing concentration-dependent clustering. Continuum modeling predicts the range of fracture strength improvements and the required ion exchange concentration profile. It was shown that a 5% surface exchange concentration can induce ∼0.6 GPa of surface compressive stress using Na + and ∼1.0 GPa of surface compressive stress using Ag + . On the other hand, larger ion exchange species may penalize Li + diffusivity by increasing the activation volume and activation energy. Interestingly, Na + has a negligible penalty on Li-ion diffusivity. The room temperature Li + ion diffusivity is reduced by ∼40% with Ag + ion exchange. Electronic band structure analysis shows no size-dependent change in the bandgap, though Ag + introduces localized defect states near the valence band maximum. This study highlights ion size as a key factor in optimizing LLZO properties, offering a framework to improve the solid-state battery performance.

Jagad, Harsh D. [Brown Univ., Providence, RI (Unit↗
DOE OSTI2024

Role of Short-Range Order on Diffusion Coefficients in the Li–Mg Alloy

Li–Mg alloys are important because of their beneficial role in fostering uniform plating and stripping of lithium in all-solid-state batteries. The alloy Li x Mg 1–x forms a solid solution on the BCC crystal structure when the lithium content is greater than x ≈ 0.3. The activation barriers of lithium and magnesium exchanges with a vacancy, crucial for substitutional diffusion, are predicted to be exceptionally low and almost identical, with negligible dependence on the alloy composition. The equilibrium vacancy concentration at room temperature is predicted to be very low, and it also remains almost constant with no dependence on Mg content in the alloy (for x Li ≥ 0.5). Nevertheless, both experiments and kinetic Monte Carlo simulations indicate that the tracer diffusion coefficients decrease by almost an order of magnitude with the addition of Mg to the alloy. In this contribution, the crucial role that chemical short-range order plays in affecting the diffusion coefficients is studied. Chemical short-range order is found to increase the effective activation barrier for lithium and magnesium diffusion by making successive atomic hops with vacancies correlated.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Direct numerical simulation of a turbulent methane/air flame impinged by a sub-breakdown electric field

Here, the effects of an external electric field on a turbulent methane/air diffusion flame are analyzed in this work using direct numerical simulations. The analyzed configuration consists of a temporally evolving mixing layer of air and a mixture of methane and nitrogen at 1 atm that is impinged by an electric field in the direction normal to the mean mixing plane. The combustion and chemi-ionization reactions involved in the flow are modeled using finite-rate chemistry and a reduced reaction scheme consisting of 26 species and 134 reactions. The mass diffusion and ion-wind effects are modeled using a detailed description of the diffusion coefficients and electric mobilities based on kinetic theory. The presented calculations show that the turbulence generated within the mixing layer is mostly unaffected by the applied electric field for the configuration under exam. In fact, the electric body force is developed away from the mixing region where the flow is uniform. Conversely, the turbulence is able to introduce very high intermittency in the electrically charged species concentration and, consequently, in the electric body force. Such intermittency will constitute a challenge in the future formulation of combustion models that take into account ion-wind effects.

09 BIOMASS FUELS↗
DOE OSTI2025

Diffusion Behavior of Oversized Fission Products in bcc Fe Cladding: A First-Principles Study

Fuel-Cladding Chemical Interaction (FCCI) poses significant challenges in nuclear reactors, where fission products from nuclear fuel interact with Fe-based cladding materials, potentially compromising their structural integrity. This study investigates the diffusion behavior of oversized fission products, Pr, Nd, Ce, and La, within bcc Fe cladding using density functional theory (DFT), nudged elastic band (NEB) method, and self-consistent mean field (SCMF) theory. Our results reveal significant long-range vacancy binding energies, particularly up to the 6th nearest neighbor, with La exhibiting the strongest binding affinity, followed by Nd, Ce, and Pr. The NEB calculations indicate significant high barriers for the dissociation of 1nn vacancy-solute pairs for all fission products. The tracer diffusion coefficients of these fission products was derived in Arrhenius form. The significant trapping effect of vacancies by a very dilute amount of fission products reduces vacancy mobility, leading to an oversaturation of point defects, void nucleation, and swelling. These are critical issues for irradiated cladding materials. The tracer diffusion coefficients indicate that Nd diffuses the fastest, followed by La, Ce, and Pr. This study provides essential insights for developing advanced cladding materials and design strategies to mitigate FCCI, ultimately enhancing nuclear reactor safety and performance.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Generation of mean flows in rotating anisotropic turbulence: The case of solar near-surface shear layer

Results from helioseismology indicate that the radial gradient of the rotation rate in the near-surface shear layer (NSSL) of the Sun is independent of latitude and radius. Theoretical models using the mean-field approach have been successful in explaining this property of the NSSL, while global direct or large-eddy magnetoconvection models have so far been unable to reproduce this. We investigate the reason for this discrepancy by measuring the mean flows, Reynolds stress, and turbulent transport coefficients under conditions mimicking those in the solar NSSL. Simulations with as few ingredients as possible to generate mean flows were studied. These ingredients are inhomogeneity due to boundaries, anisotropic turbulence, and rotation. The parameters of the simulations were chosen such that they matched the weakly rotationally constrained NSSL. The simulations probe locally Cartesian patches of the star at a given depth and latitude. The depth of the patch was varied by changing the rotation rate such that the resulting Coriolis numbers covered the same range as in the NSSL. We measured the turbulent transport coefficient relevant for the nondiffusive (Λ-effect) and diffusive (turbulent viscosity) parts of the Reynolds stress and compared them with predictions of current mean-field theories. A negative radial gradient of the mean flow is generated only at the equator where meridional flows are absent. At other latitudes, the meridional flow is comparable to the mean flow corresponding to differential rotation. We also find that the meridional components of the Reynolds stress cannot be ignored. Additionally, we find that the turbulent viscosity is quenched by rotation by about 50% from the surface to the bottom of the NSSL. Our local simulations do not validate the explanation for the generation of the NSSL from mean-field theory where meridional flows and stresses are neglected. However, the rotational dependence of the turbulent viscosity in our simulations agrees well with theoretical predictions. Moreover, our results agree qualitatively with global convection simulations in that an NSSL can only be obtained near the equator.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Diffusion-Limited Kinetics in Reactive Systems

A proper representation of chemical kinetics is vital to understanding, modeling, and optimizing many important chemical processes. In liquid and surface phases, where diffusion is slow, the rate at which the reactants diffuse together limits the overall rate of many elementary reactions. Commonly, the textbook Smoluchowski theory is utilized to estimate effective rate coefficients in the liquid phase. On surfaces, modelers commonly resort to much more complex and expensive Kinetic Monte Carlo (KMC) simulations. Here, in this study, we extend the Smoluchowski model to allow the diffusing species to undergo chemical reactions and derive analytical formulas for the diffusion-limited rate coefficients for 3D, 2D, and 2D/3D interface cases. With these equations, we are able to demonstrate that when species react faster than they diffuse they can react orders of magnitude faster than predicted by Smoluchowski theory, through what we term “the reactive transport effect”. We validate the derived steady-state equations against particle Monte Carlo (PMC) simulations, KMC simulations, and non-steady-state solutions. Furthermore, using PMC and KMC simulations, we propose corrections that agree with all limits and the computed data for the 2D and 2D/3D interface steady-state equations, accounting for unique limitations in the associated derived equations. Additionally, we derive equations to handle couplings between diffusion-limited rate coefficients in reaction networks. We believe these equations should make it possible to run much more accurate mean-field simulations of liquids, surfaces, and liquid–surface interfaces accounting for diffusion limitations and the reactive transport effect.

Johnson, Matthew Sean↗
DOE OSTI2024

In-Situ and Ex-Situ Studies on the Morphology Changes of Polymer Pipeline Materials for Use in Hydrogen Gas Environments

The US natural gas infrastructure is a national asset that could be used to deliver hydrogen and hydrogen blends of natural gas as a pathway to reduce carbon emissions. The distribution system comprises nearly 50% plastic pipe composed of medium- and high-density polyethylene materials (MDPE and HDPE). While these materials perform adequately for natural gas, research on their hydrogen compatibility is essential to understand if any immediate and long-term risks are associated with hydrogen addition. The Blended Gas CRADA, a HyBlend project, has established a comprehensive test method for evaluating MDPE and HDPE of various plastic resin compositions of pipeline material in pure hydrogen and 20% hydrogen/80% methane blends. Both in-situ and ex-situ measurements were performed to capture hydrogen-induced changes in the polyethylene material's crystalline, amorphous and their interphase regions. We investigated MDPE and HDPE pipeline materials made from different polymer resin systems to evaluate the effects of hydrogen gas. The materials were characterized by their density, diffusion coefficient, free volume ratio, and degree of crystallinity. Various advanced characterization methods, including in situ NMR, ex situ XRD, ex situ DSC, and ex situ TDA, were used to analyze the effects of changes in crystalline, amorphous, and interphase regions due to gas exposure. Time-dependent post-decompression quasi-static tensile tests were conducted to explore the effects of gas exposure time on the mechanical behavior of the pipe materials. This work will highlight the time sensitivities during and after gas exposure. The correlation between gas-induced polyethylene morphology changes and the associated material performance will be addressed for the intended applications. These studies will show that polyethylene resin composition and material exposure are important factors when considering whether hydrogen gas affects pipeline materials positively or negatively.

Simmons, Kevin L.↗
DOE OSTI2021

Molecular dynamics simulation of the effect of cooling rate on the structure and properties of lithium disilicate glass

The effects of cooling rate on the structure and properties of lithium disilicate (LS2) glass are investigated using molecular dynamics (MD) computer simulations. The evolution of structural features such as pair distribution function, bond angle distribution, and Li coordination distribution are determined, and correlated with dynamic and static properties to elucidate the effects of cooling rate. The density, elastic moduli, and diffusion coefficient are found to be highly sensitive to cooling rate, whereas the Si pair distribution function and bond angle distribution are weakly affected by the cooling rate. Additionally, the changes of Si-O-Si bond angle and Li coordination number suggest the formation of Li cluster at lower cooling rates. Furthermore, by comparing results from other simulations and reported experiments, we confirm that the increase of cooling rate leads to an increase of conductivity and a decrease of density. Finally, at very high cooling rates, we find that all the atoms in LS2-glass do not relax simultaneously, but do so in two distinct configurations.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Mass Changes the Diffusion Coefficient of Particles with Ligand-Receptor Contacts in the Overdamped Limit

Inertia does not generally affect the long-time diffusion of passive overdamped particles in fluids. Yet a model starting from the Langevin equation predicts a surprising property of particles coated with ligands that bind reversibly to surface receptors: heavy particles diffuse more slowly than light ones of the same size. We show this by simulation and by deriving an analytic formula for the mass-dependent diffusion coefficient in the overdamped limit. We estimate the magnitude of this effect for a range of biophysical ligand-receptor systems, and find it is potentially observable for tailored micronscale DNA-coated colloids.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Temperature dependence of resistivity of carbon micro/nanostructures: Microscale spatial distribution with mixed metallic and semiconductive behaviors

The temperature coefficient of resistivity (θ T ) of carbon-based materials is a critical property that directly determines their electrical response upon thermal impulses. It could have metal- (positive) or semiconductor-like (negative) behavior, depending on the combined temperature dependence of electron density and electron scattering. Its distribution in space is very difficult to measure and is rarely studied. Here, for the first time, we report that carbon-based micro/nanoscale structures have a strong non-uniform spatial distribution of θ T . This distribution is probed by measuring the transient electro-thermal response of the material under extremely localized step laser heating and scanning, which magnifies the local θ T effect in the measured transient voltage evolution. For carbon microfibers (CMFs), after electrical current annealing, θ T varies from negative to positive from the sample end to the center with a magnitude change of >130% over <1 mm. This θ T sign change is confirmed by directly testing smaller segments from different regions of an annealed CMF. For micro-thick carbon nanotube bundles, θ T is found to have a relative change of >125% within a length of ~2 mm, uncovering strong metallic to semiconductive behavior change in space. In conclusion, our θ T scanning technique can be readily extended to nm-thick samples with μm scanning resolution to explore the distribution of θ T and provide a deep insight into the local electron conduction.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Effect of Mg and Ni impurities on tritium diffusion in lithium ceramics through cluster dynamics simulations

This study investigates the impact of Mg and Ni doping on tritium diffusion in LiAlO 2 and LiAl 5 O 8 ceramics, that are used in tritium-producing burnable absorber rods (TPBARs). Utilizing Centipede simulations across a broad temperature range (500 K to 1250 K), we explore the interplay between defect dynamics, cluster formation, and tritium mobility. In LiAlO 2 , Mg doping significantly enhances tritium diffusivity by increasing tritium interstitial concentrations and diffusion coefficients of key species, thereby doubling the overall tritium diffusivity. Ni doping, while shifting the dominant defect to Li vacancies, maintains high tritium mobility due to the low binding energy of Li vacancy-tritium complexes, which ensures effective tritium migration. In LiAl 5 O 8 , Mg and Ni doping results in a slight reduction in the diffusion coefficients of key species, yet the dramatic increase in tritium interstitial concentrations compensates, leading to a net small increase in tritium diffusivity. In conclusion, the findings highlight the critical role of defects in tritium transport and the effect of Mg and Ni defects on the performance of these ceramics in demanding nuclear environments.

Cluster dynamics↗
DOE OSTI2026

Grain Boundary Effects on Radiation Damage and Tritium Diffusion in Li–Al–O Ceramics from Molecular Dynamics and Experiments

Understanding irradiation damage and tritium transport in LiAlO2 ceramics is essential for their deployment in tritium-producing burnable absorber rods (TPBARs). Grain boundaries (GBs) play an important role in governing radiation response and tritium transport in LiAlO2 and its secondary phase LiAl5O8. Using molecular dynamics simulations, we investigated defect evolution and tritium diffusion during displacement cascades in single-crystal and bicrystal LiAlO2 and LiAl5O8 at 600 K. The results reveal that GBs suppress the damage in the LiAlO2 bicrystals as compared to single crystals, via interstitial emission. The effect is more pronounced for Li defects, than for Al and O. In LiAl5O8 the damage reduction due to introduction of GB is less pronounced overall. Tritium diffusion coefficients at GBs are higher by factors of 2–10 relative to diffusion inside grain bulk. The effect is particularly pronounced in LiAlO2, where S5, S17, and S25 GBs promote rapid tritium transport, whereas LiAl5O8 exhibits slower diffusion due to reduced free volume at its GBs. These findings evince a trade-off that while GBs mitigate radiation damage by absorbing interstitials, they simultaneously provide fast pathways for tritium migration which is not desirable for TPBAR application. The mechanistic insights gained here establish a foundation for microstructural design strategies to balance radiation tolerance and tritium retention in ceramic breeder materials.

Roy, Ankit↗
DOE OSTI2021

Temperature and salt controlled tuning of protein clusters

The formation of molecular assemblies in protein solutions is of strong interest both from a fundamental viewpoint and for biomedical applications. While ordered and desired protein assemblies are indispensable for some biological functions, undesired protein condensation can induce serious diseases. As a common cofactor, the presence of salt ions is essential for some biological processes involving proteins, and in aqueous suspensions of proteins can also give rise to complex phase diagrams including homogeneous solutions, large aggregates, and dissolution regimes. Here, we systematically study the cluster formation approaching the phase separation in aqueous solutions of the globular protein BSA as a function of temperature (T), the protein concentration (c p ) and the concentrations of the trivalent salts YCl 3 and LaCl 3 (c s ). As an important complement to structural, i.e. time-averaged, techniques we employ a dynamical technique that can detect clusters even when they are transient on the order of a few nanoseconds. By employing incoherent neutron spectroscopy, we unambiguously determine the short-time self-diffusion of the protein clusters depending on c p , c s and T. We determine the cluster size in terms of effective hydrodynamic radii as manifested by the cluster center-of-mass diffusion coefficients D. For both salts, we find a simple functional form D(c p , c s , T) in the parameter range explored. The calculated inter-particle attraction strength, determined from the microscopic and short-time diffusive properties of the samples, increases with salt concentration and temperature in the regime investigated and can be linked to the macroscopic behavior of the samples.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Diffusion behavior of lanthanide fission products in bcc Fe cladding: A first-principles study

Fuel-cladding chemical interaction poses significant challenges in nuclear reactors, where fission products generated from nuclear fuel interact with Fe-based cladding materials, potentially compromising their structural integrity. This study investigates the diffusion behavior of lanthanide fission products, Lanthanum (La), Cerium (Ce), Praseodymium (Pr), and Neodymium (Nd), within body-centered cubic (bcc) Fe cladding using the density functional theory, nudged elastic band method, and self-consistent mean field theory. Our results reveal significant vacancy binding energies, particularly with the 1st and 2nd nearest neighbors, which diminish beyond the 5th nearest neighbor, with La exhibiting the strongest binding affinity, followed by Nd, Ce, and Pr. The nudged elastic band calculations indicate significant high barriers for the dissociation of 1st nearest neighbor vacancy-solute pairs for all fission products. The tracer diffusion coefficients of these fission products were derived in an Arrhenius form, with a magnetic correction that accounts for the high-temperature paramagnetic state. The significant trapping effect of vacancies caused by a very dilute concentration of fission products reduces vacancy mobility, potentially leading to modifications in point defect supersaturation, void nucleation, and swelling under irradiation. These represent critical challenges for irradiated cladding materials. The tracer diffusion coefficients indicate that Nd diffuses the fastest, followed by La, Ce, and Pr. Furthermore, this study provides essential insights for understanding fission product transport in cladding materials and informs future design strategies to mitigate fuel-cladding chemical interaction, ultimately enhancing nuclear reactor safety and performance.

Diffusion↗
DOE OSTI2024

Using Active Learning to Rapidly Develop Machine Learned Diffusion Coefficients of CO 2 Conversion Reagents in Metal–Organic Frameworks

Here, we used a combined molecular dynamics/active learning (AL) approach to create machine learning models that can predict the diffusion coefficient of epichlorohydrin and chloropropene carbonate, the reactant and product of a common CO 2 cycloaddition reaction, in metal–organic frameworks (MOFs). Nanoporous MOFs are effective catalysts for the cycloaddition of CO 2 to epoxides. The diffusion rates within nanoporous catalysts can control the rate of reaction as the reactants and products must diffuse to the active sites within the MOF and then out of the nanoporous material for reusability. However, the diffusion process is routinely ignored when searching for new materials in catalytic applications. Here we verified improvement during the AL process by consistently tracking metrics on the same groups of MOFs to ensure consistency. Metal identity was found to have little impact on diffusion rates, while structural features like pore limiting diameter act as a threshold where a minimum value is needed for high diffusion rates. We identified the MOFs with the highest epichlorohydrin and chloropropene carbonate diffusion coefficients which can be used for further studies of reaction energetics.

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

Underlying mechanism of structural transformation between GaSb and GaAs response to intense electronic excitation

Ion irradiation of semiconductors has emerged as a promising approach for fabricating self-organized nanosystems with high atomic precision, despite often being accompanied by undesirable phenomena. Exploring the mechanisms underlying structural transformations is crucial for assessing nanostructure array types under complex irradiation environments. By quantitatively calculating the thermodynamically driven processes and analyzing the impact of intrinsic structural parameters, distinct structural transformations in response to intense electronic excitation are systematically investigated in gallium antimonide (GaSb) and gallium arsenide (GaAs) systems. In high-energy regimes, the nanofibers layer of GaSb exhibits intriguing structural discrepancy, characterized by partial nanofibers with coherent boundaries, interspersed nanopores accompanied by antisite defects and Ga precipitates, distinguishing to a series of discontinuous latent tracks that emerged within cylindrical trajectories in GaAs. Furthermore, significant diffusion behaviors of the nanohillocks are discovered in GaAs, with higher average roughness than GaSb, driven by the gradient stress distribution influenced by the free-surface effects. The deposition energy for melting phase formation, Gibbs free energy, and Ga diffusion coefficients contribute to the distinctive structural features, evidencing relatively stable morphological configurations and higher irradiation resistance in GaAs. Consequently, special optoelectronic properties associated with structural discrepancies facilitate the design and optimization of material functionalities by irradiation technologies.

36 MATERIALS SCIENCE↗