Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “microstructure simulation”

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 1,261 records · Page 70

State of the art, gaps, and prospects in fusion materials theory and modelling

Advancing the theory and simulation of materials for fusion applications remains a key component of global roadmaps aimed at delivering much-needed fusion power. Especially as the drive for commercial application increases, prototypes must be designed against radiation damage before the relevant experimental data can be collected and cost reductions that are possible by testing materials in silico become even more important. Here, we summarise the state of the art as it emerged during the 7 th Fusion Materials Theory & Modelling Workshop that took place in 2024, with the aim to highlight present gaps and future directions for the fusion materials modelling community. Of particular interest were the effects of transmutations, chemical complexity with the development of novel alloys and interatomic potentials, advancements in modelling high-dose microstructures, comparison with experimental data and multiscale models for structural assessment relying on high-performance computing and virtual reality.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Sensitivity of Thermal Predictions to Uncertain Surface Tension Data in Laser Additive Manufacturing

To understand the process-microstructure relationships in additive manufacturing (AM), it is necessary to predict the solidification characteristics in the melt pool. This study investigates the influence of Marangoni driven fluid flow on the predicted melt pool geometry and solidification conditions using a continuum finite volume model. A calibrated laser absorptivity was determined by comparing the model predictions (neglecting fluid flow) against melt pool dimensions obtained from single laser melt experiments on a nickel super alloy 625 (IN625) plate. Using this calibrated efficiency, predicted melt pool geometries agree well with experiments across a range of process conditions. When fluid mechanics is considered, a surface tension gradient recommended for IN625 tends to overpredict the influence of convective heat transfer, but the use of an intermediate value reported from experimental measurements of a similar nickel super alloy produces excellent experimental agreement. Despite its significant effect on the melt pool geometry predictions, fluid flow was found to have a small effect on the predicted solidification conditions compared to processing conditions. We find that this result suggests that under certain circumstances, a model only considering conductive heat transfer is sufficient for approximating process-microstructure relationships in laser AM. Extending the model to multiple laser passes further showed that fluid flow also has a small effect on the solidification conditions compared to the transient variations in the process. Furthermore, limitations of the current model and areas of improvement, including uncertainties associated with the phenomenological model inputs are discussed.

36 MATERIALS SCIENCE↗

Structural uniformity and compositional homogeneity of solid-phase alloyed rod

Solid-phase processes have emerged as an alternative to fusion-based alloying to avoid coarse microstructures, undesirable phase formation, and high energy consumption. However, achieving uniform distribution of alloying elements during friction-based processing remains challenging due to highly heterogeneous thermomechanical conditions. This work evaluates the structural uniformity and compositional homogeneity of Al–Cu–Zn alloyed rods produced by friction extrusion (FE) and establishes the role of the rotational speed to feed rate ratio (N/V) on alloying effectiveness. A systematic matrix of FE experiments was conducted at constant extrusion ratio with N/V values ranging from 3.7 to 300. Compositional uniformity was assessed along the rod length (ICP-OES), in three dimensions (X-ray computed tomography), and at the microscale (SEM–EDS), supported by a gray-level co-occurrence matrix (GLCM)–based homogeneity metric. Smoothed particle hydrodynamics (SPH) simulations were used to reveal material flow and thermomechanical fields. Results show that N/V = 100 produces a high-shear mixing zone that eliminates the unmixed core and enables near-full dissolution and dispersion of Cu and Zn. At lower N/V, a laminar flow region persists at the rod center, causing segregation and large composition gradients. The combined experimental–computational analysis provides mechanistic insight into the transition from fragmented particle dispersion to thermomechanically assisted metallurgical mixing. This study establishes processing–structure relationships for solid-phase alloying and provides guidance for achieving homogenized compositions comparable to wrought alloys via rapid, scalable FE processing.

Aluminum↗

Stochastic 3D reconstruction of cracked polycrystalline NMC particles using 2D SEM data

Li-ion battery performance is strongly influenced by the 3D microstructure of its cathode particles. Cracks within these particles develop during calendaring and cycling, reducing connectivity but increasing reactive surface, making their impact on battery performance complex. Understanding these contradictory effects requires a quantitative link between particle morphology and battery performance. However, informative 3D imaging techniques are time-consuming, costly and rarely available, such that analyses often have to rely on 2D image data. This paper presents a novel stereological approach for generating virtual 3D cathode particles exhibiting crack networks that are statistically equivalent to those observed in 2D sections of experimentally measured particles. Consequently, 2D image data suffices for deriving a full 3D characterization of cracked cathodes particles. Such virtually generated 3D particles could serve as geometry input for spatially resolved electro-chemo-mechanical simulations to enhance our understanding of structure-property relationships of cathodes in Li-ion batteries.

36 MATERIALS SCIENCE↗

A Combined Computational, Experimental, and Technology Development Approach to In-Space Laser Manufacturing Maturation at NASA Marshall Space Flight Center

In-space manufacturing (ISM) is emerging as a field vital to continued access and capabilities in the space environment. NASA Marshall Space Flight Center (MSFC) is advancing the frontier of in-space laser manufacturing (ISLM) techniques through work initially focused on maturing laser beam welding (LBW) and laser forming (LF) for use in space. Such techniques proffer the ability to assemble and join structures in space from sheet metal or other stock – extant satellites, in situ resource utilization of Lunar regolith, etc. – by forming to desired shapes and then joining via in-space welding (ISW). ISLM processes are useful for assembly, joining, modification, and repair of structures in free space and on the Lunar surface such as large observatories, antennas, trusses, blast/thermal/radiation shields, pressure vessels, and more. However, these techniques are not yet qualified & certified (Q&C) for regular application in space. It would be prohibitively expensive, laborious, and time-consuming to perform Q&C via traditional experimental approaches as data collection & experimentation in space is resource-intensive. As such, benchmark experiments and focused, properly instrumented technology demonstration efforts in space can collect sufficient data that – when combined with verified computational models in an integrated computational materials engineering (ICME) approach – can validate ICME tools capable of translating more readily obtained ground data to in-space, in situ, computationally informed Q&C of ISLM techniques. Several ISLM projects at MSFC are obtaining the data required to validate ICME tools through both ground and flight experiments. A parabolic flight experiment of LBW under vacuum is manifested for August 2024, including both microgravity and Lunar gravity profiles. This collaboration with the Ohio State University is investigating common aerospace alloys such as 316L stainless steel, 2219 aluminum alloy, and Ti64 titanium alloy. In situ data collection includes videography, thermography, and reference thermocouples to build a thermal model of the welds. This will elucidate the relevant physics when combined with post-flight microstructural examination and mechanical testing. MSFC is also progressing towards a suborbital flight experiment of LBW under vacuum, which could provide reams of data on ISW during sustained, high-quality reduced gravity. The effect of combined thermal (cryogenic and high-temperature) and vacuum exposure on both LBW (NASA-funded) and LF (DARPA-funded) is being investigated through ground experiments. In addition to the copious data collected during these ground experiments, ruggedization of LBW hardware will also be pursued. The datasets from these experiments will be used to validate computational models which will inform future ISLM efforts in an ICME framework. A variety of techniques across lengths scales, from CALPHAD-driven thermodynamics & kinetics to phase field modeling of solidification to kinetic Monte Carlo simulations of grain evolution at the mesoscale, will be employed to accelerate the infusion and eventual Q&C of LBW and LF for use in space. The development of data-driven surrogate models to bridge ground to flight experiments and thereby reduce the need for resource-intensive experiments in space will also be investigated. These ICME techniques, surrogate models, and datasets from ground testing can also be employed to advance manufacturing in terrestrial environments.

in-space welding↗

Intermediate-phase engineering via dimethylammonium cation additive for stable perovskite solar cells

Achieving the long-term stability of perovskite solar cells is arguably the most important challenge required to enable widespread commercialization. Understanding the perovskite crystallization process and its direct impact on device stability is critical to achieving this goal. The commonly employed dimethyl-formamide/dimethyl-sulfoxide solvent preparation method results in a poor crystal quality and microstructure of the polycrystalline perovskite films. In this work, we introduce a high-temperature dimethyl-sulfoxide-free processing method that utilizes dimethylammonium chloride as an additive to control the perovskite intermediate precursor phases. By controlling the crystallization sequence, we tune the grain size, texturing, orientation (corner-up versus face-up) and crystallinity of the formamidinium (FA)/caesium (FA) y Cs 1-y Pb(I x Br 1-x )3 perovskite system. A population of encapsulated devices showed improved operational stability, with a median T80 lifetime (the time over which the device power conversion efficiency decreases to 80% of its initial value) for the steady-state power conversion efficiency of 1,190 hours, and a champion device showed a T80 of 1,410 hours, under simulated sunlight at 65 degrees C in air, under open-circuit conditions. In conclusion, this work highlights the importance of material quality in achieving the long-term operational stability of perovskite optoelectronic devices.

14 SOLAR ENERGY↗

Electrical and Thermal Performance of Ga₂O₃–Al₂O₃–Diamond Super-Junction Schottky Barrier Diodes

The design space of Ga 2 O 3 -based devices is severely constrained due to its low thermal conductivity and absence of viable p-type dopants. In this work, we discuss the limits of operation of a novel Ga 2 O 3 –Al 2 O 3 –diamond-based super-junction device concept, which can alleviate the constraints associated with Ga 2 O 3 -based devices. The improvements achieved using the proposed device concept are demonstrated through electrical and thermal simulations of Ga 2 O 3 –Al 2 O 3 –diamond-based super-junction Schottky barrier diodes (SJ-SBDs) and non-punch-through or conventional Schottky barrier diodes (NP-SBDs). The SJ-SBD enables operation below the R ON -breakdown voltage limit of Ga 2 O 3 NP-SBD, enabling >4 kV blocking voltage at R ON of 1–3 mΩ cm 2 . Here, the maximum switching frequency of SJ-SBD may be only a few kHz, as it is limited by the activation energy of acceptors (0.39 eV) in the diamond. Crucially, compared with NP-SBD, the use of diamond also results in ~60% reduction in temperature rise during static power dissipation. Polycrystalline diamond (PCD) properties depend on detailed microstructure and benefits compared to ideal Ga 2 O 3</sub NP-SBD arise for diamond critical fields ≥6 MV/cm and thermal conductivities as low as 50–150 W/(m • K).

36 MATERIALS SCIENCE↗

Microstructural Changes and Pb Mobility During the Zircon to Reidite Transformation: Implications for Planetary Impact Chronology

Impact events modify and leave behind a complex history of rock metamorphism on terrestrial planets. Evidence for an impact event may be recorded in physical changes to minerals, such as mineral deformation and formation of high P-T polymorphs, but also in the form of chemical fingerprints, such as enhanced elemental diffusion and isotopic mixing. Here we explore laboratory shock-induced physical and chemical changes to zircon and feldspar, the former of which is of interest because its trace elements abundances and isotope ratios are used extensively in geochemistry and geochronology. To this end, a granular mixture of Bishop Tuff sanidine and Kuehl Lake zircon, both with well characterized Pb isotope compositions, was prepared and then shocked via a flat plate accelerator. The peak pressure of the experiment, as calculated by the impedance matching method, was ~24 GPa although a broader range of P-T conditions is anticipated due to starting sample porosity. Unshocked and shocked materials were characterized via scanning electron microscopy (SEM), electron backscatter diffraction (EBSD), and Raman spectroscopy. These methods show that the starting zircon material had abundant metamict regions, and the conversion of the feldspar to glass in the post-shock material. Analyses of the shocked product also yielded multiple occurrences of the high-pressure ZrSiO 4 polymorph reidite, with some domains up to 300 μm across. The possibility of U-Pb system disturbance was evaluated via laser ablation-inductively coupled plasma-mass spectrometry (LA-ICP-MS) and secondary ion mass spectrometry (SIMS). The isotopic data reveal that disturbance of the U-Pb geochronometer in the reidite was minimal (<2% for the main U-Pb geochronometers). To better constrain the P-T conditions during the shock experiment, we complement impedance matching pressure calculations with iSALE2D impact simulations. The simulated results yield a range of P-T conditions experienced during the experiment and show that much of the sample may have reached >30 GPa, which is consistent with formation of reidite. In the recovered shocked material, we identified lamellae of reidite, some of which interlock with zircon lamellae. Reidite {112} twins were identified, which we interpret to have formed to reduce stress between the crystal structure of the host zircon and reidite. These two findings support the interpretation that shear transformation enabled the transition of zircon to reidite. The size and presence of reidite found here indicate that this phase is probably common in impact-shocked crustal rocks that experienced ~25 to ~35 GPa, especially when the target material has porosity. Additionally, shock loading of the zircon and transformation to reidite at these pressures in porous materials is unlikely to significantly disturb the U-Pb system in zircon and that the reidite inherits the primary U and Pb elemental and isotopic ratios from the zircon.

Bishop Tuff↗

Implementation and experimental validation of nonlocal damage in a large-strain elasto-viscoplastic FFT-based framework for predicting ductile fracture in 3D polycrystalline materials

Ductile materials, such as metal alloys, can undergo substantial deformation before failure. Additionally, these materials are usually of polycrystalline composition and exhibit strongly anisotropic behavior at small length scales. Previously developed fast Fourier transform (FFT)-based models can model ductile fracture of isotropic materials or the elastic–plastic behavior of anisotropic polycrystalline materials; however, there remains a need to couple both capabilities. This work extends a large-strain FFT-based crystal plasticity model to simulate ductile fracture of polycrystalline materials. In this work, a triaxiality-based continuum damage mechanics (CDM) formulation is incorporated into a large-strain elasto-viscoplastic FFT (LS-EVPFFT) framework. The CDM formulation is augmented with an integral-based nonlocal regularization approach that correctly handles gas-phase material necessary to model unconstrained surfaces. To validate the damage-enabled LS-EVPFFT framework, mesoscale copper tensile coupons were machined using microwire electrical discharge machining and experimentally characterized using electron backscatter diffraction. In-situ optical digital image correlation was performed during uniaxial testing to provide a side-by-side comparison of the experimental and computational strain fields and stress–strain responses. The damage-enabled LS-EVPFFT framework can simulate the complete macroscopic stress–strain response of ductile polycrystals to failure. The model reproduces necking behavior that qualitatively agrees with experimental observations. By leveraging the relatively low computational cost of the damage-enabled LS-EVPFFT framework, the framework presented here allows the ductile fracture response of 3D polycrystalline materials to be tractably predicted.

36 MATERIALS SCIENCE↗

Nucleation of Dislocations in 3.9 nm Nanocrystals at High Pressure

As circuitry approaches single nanometer length scales, it has become important to predict the stability of single nanometer-sized metals. The behavior of metals at larger scales can be predicted based on the behavior of dislocations, but it is unclear if dislocations can form and be sustained at single nanometer dimensions. Here, we report the formation of dislocations within individual 3.9 nm Au nanocrystals under nonhydrostatic pressure in a diamond anvil cell. We used a combination of x-ray diffraction, optical absorbance spectroscopy, and molecular dynamics simulation to characterize the defects that are formed, which were found to be surface-nucleated partial dislocations. These results indicate that dislocations are still active at single nanometer length scales and can lead to permanent plasticity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Probing Surface Plasmon Dynamics in Periodic Nanostructures through Ultrafast Electron Microscopy

Surface plasmon polaritons (SPPs) can be manipulated to localize and guide light in subwavelength distances, enabling them to find applications in a wide range of areas, from sensing to quantum computing. Among several methods of SPP excitation, periodic arrays of nano- and microstructures are of particular interest, as they enable engineering SPP properties through structural parameters. Here, in this study, using the photon-induced near-field electron microscopy (PINEM) technique, we investigated the mode formation, coupling, interference, and decay of SPPs in square and hexagonal arrays of circular nanoholes under both visible and near-infrared excitation. Polarization-resolved analysis revealed the key factors governing SPP localization and interference patterns, showing that the periodicity and symmetry of the array primarily determine the SPP interference patterns and their orientation, while pump polarization mainly modulates their intensity. Time-resolved PINEM measurements demonstrated the spatial dependence of the SPP temporal characteristics. In addition, cathodoluminescence (CL) spectroscopy was employed to examine the intrinsic plasmonic characteristics of the structure. Finite difference time domain (FDTD) simulations showed strong agreement with both PINEM and CL measurements on the spatial and spectral behavior of SPPs. Understanding the spatiotemporal dynamics of SPPs on nanostructures beyond the diffraction limit is crucial for optimizing plasmonic structures for advanced photonic and quantum technologies.

Plasmonics↗

N-Doping improves charge transport and morphology in the organic non-fullerene acceptor O-IDTBR

Molecular doping has been shown to improve the performance of various organic (opto)electronic devices. When compared to p-doped systems, research into n-doped organic small-molecules is relatively limited, primarily due to the lack of suitable dopants and the often encountered unfavourable microstructural effects. These factors have prevented the use of n-doping in a wider range of existing materials, such as non-fullerene acceptors (NFAs), that have already shown great promise for a range of (opto)electronic applications. Here, we show that several different molecular n-dopants, namely [1,2-b:2',1'-d]benzo[i][2.5]benzodiazocine potassium triflate adduct (DMBI-BDZC), tetra-n-butylammonium fluoride (TBAF) and 4-(2,3-dihydro-1,3-dimethyl-1H-benzimidazol-2-yl)-N,N-dimethylbenzenamine (N-DMBI), can be used to n-dope the molecular semiconductor O-IDTBR, a promising NFA, and increase the electron field-effect mobility to >1 cm 2 V -1 s -1 . By combining complementary experimental techniques with computer simulations of doping and charge carrier dynamics, we show that improved charge transport arises from synergistic effects of n-type doping and morphological changes. Specifically, a new, previously unreported dopant-induced packing orientation results in one of the highest electron mobility values reported to-date for an NFA molecule. Overall, this work highlights the importance of dopant–semiconductor interactions and their impact on morphology, showing that dopant-induced molecular packing motifs may be generic and a key element of the charge transport enhancement observed in doped organics.

36 MATERIALS SCIENCE↗

Effects of low earth orbit on the optical performance of multi-layer enhanced high reflectance mirrors

Two mirror designs developed for space applications were flown along with a standard mid-infrared design on the leading and trailing edges of the Long Duration Exposure Facility (LDEF). Preliminary observations of induced changes in optical performance of ZnS-coated mirrors and impact-related microstructural and microchemical effects are described in the proceedings of the First LDEF Post-Retrieval Symposium. In this paper, effects of the induced environment and meteoroid/debris impacts on mirror performance are described in more detail. Also, an analysis of reflectance spectra using the results of Auger and secondary ion mass spectroscopy (SIMS) profiling measurements are used to identify an optical-degradation mechanism for the ZnS-coated mirrors. Structural damage associated with a high-velocity impact on a (Si/Al2O3)-coated mirror was imaged optically and with scanning electron and atomic force microscopy (SEM and AFM). Scanning Auger and SIMS analysis provided chemical mapping of selected impact sites. The impact data suggest design and fabrication modifications for obtaining improved mechanical performance using a design variation identified in preflight laboratory simulations. Auger surface profile and SIMS imaging data verified the conclusion that secondary impacts are the source of contamination associated with the dendrites grown on the leading-edge ZnS-coated test samples. It was also found that dendrites can be grown in the laboratory by irradiating contaminated sites on a trailing-edge ZnS-coated sample with a rastered electron beam. These results suggest a mechanism for dendrite growth.

Donovan, Terence↗

Sensitivity of UO 2 fuel performance to microstructural evolutions driven by dilute additives

Use of dilute additives to nuclear fuel is being considered to increase the security of commercial fuel management through traceability of fabricated fuel elements. Taggants, as additives are denoted when included for traceability purposes, may also improve fuel performance, as demonstrated in Cr-containing uranium dioxide as described in the literature, and they may also improve fuel safety. In fact, studies have shown that some additives affect fuel material properties such as grain size and density after sintering. Given the possible range of elements that could be used as additives, the impact of such fuel property variations on the fuel’s thermomechanical behavior becomes relevant. These effects can be evaluated through a sensitivity study of standard fuel models to analyze changes in these properties using a fuel performance code. In this work, the BISON code is being used to investigate these effects through a 2D axisymmetric model of smeared UO 2 fuel pellets and ZIRLO® cladding under realistic pressurized water reactor core irradiation conditions. Here, randomly sampled densities and grain sizes within specified ranges are used as input parameters in the simulations, and several fuel model-related outputs are evaluated. The thermomechanical response of the cladding is also addressed in this study. The simultaneous variation of both input parameters offers a more comprehensive path to identify key sensitivities. Outputs explored include temperature, fission gas release, creep, and radial stress. Results show that although most of these outputs are sensitive to grain size to a certain extent, density mainly affects fuel temperature and elastic strain. Furthermore, sensitivities can vary depending on the radial position within the fuel pellet.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Enhanced Radiation Damage Tolerance of Amorphous Interphase and Grain Boundary Complexions in Cu-Ta

Amorphous interfacial complexions are particularly resistant to radiation damage and have been primarily studied in alloys with good glass-forming ability, yet recent reports suggest that these features can form even in immiscible alloys such as Cu-Ta under irradiation. In this work, the mechanisms of damage production and annihilation due to primary knock-on atom collisions are investigated for amorphous interphase and grain boundaries in a Cu-Ta alloy using atomistic simulations. Amorphous complexions, in particular amorphous interphase complexions that separate Cu and Ta grains, result in less residual defect damage than their ordered counterparts. Stemming from the nanophase chemical separation in this alloy, the amorphous complexions exhibit a highly heterogeneous distribution of atomic excess volume, as compared to a good glass former like Cu-Zr. Complexion thickness, a tunable structural descriptor, plays a vital role in damage resistance. Thicker interfacial films are more damage-tolerant because they alter the defect production rate due to differences in intrinsic displacement threshold energies during the collision cascade. Overall, the findings of this work highlight the importance of interfacial engineering in enhancing the properties of materials operating in radiation-prone environments and the promise of amorphous complexions as particularly radiation damage-tolerant microstructural features.

36 MATERIALS SCIENCE↗

Hard superellipse phases: particle shape anisotropy & curvature

In this work, we report computer simulations of two-dimensional convex hard superellipse particle phases vs. particle shape parameters including aspect ratio, corner curvature, and sidewall curvature. Shapes investigated include disks, ellipses, squares, rectangles, and rhombuses, as well as, shapes with non-uniform curvature including rounded squares, rounded rectangles, and rounded rhombuses. Using measures of orientational order, order parameters, and a novel stretched bond orientational order parameter, we systematically identify particle shape properties that determine liquid crystal and crystalline phases including their coarse boundaries and symmetry. We observe phases including isotropic, nematic, tetratic, plastic crystals, square crystals, and hexagonal crystals (including stretched variants). Our results catalog known benchmark shapes, but include new shapes that also interpolate between known shapes. Our results indicate design rules for particle shapes that determine two-dimensional liquid, liquid crystalline, and crystalline microstructures that can be realized via particle assembly.

42 ENGINEERING↗

Myco-Architecture off Planet: Growing Surface Structures at Destination

Our work focused on filling major select key technical knowledge gaps, and technical aspects to be addressed in a Phase II proposal were identified such as the secretion of glues or plastics by the fungi or bacteria to form biocomposites. Specifically, we raised the technical TRL by assessing the growth of in-house mycelial-producing strains on potential food substrates and analyzed the advantages and disadvantages to their use off planet. Growth of the mycelia on sawdust and powdered nutrients including dried, powdered algae, was used as a baseline. We measured growth as a function of temperature, tested the relevant material properties of the mycelia products produced, and measured the mass of input and output volume for production per volume of material. Mycotecture was produced in a bag simulating the proposed mission implementation. Post-production such as heat treatment (as is done with terrestrial mycotecture)was assessed.The material properties of the dried, frozen and baked mycelial outputs included density, strength,thermal insulation and flame retardation. Tensile/compression testing on biocomposite material specimens was performed to obtain three-dimensional elastic constants, as well as to examine deformation and fracture behavior. These mechanical tests were performed on both dormant and activated samples to understand the evolution of the structural material. In addition, experimental data obtained from the mechanical testing was used to build a failure prediction model that accounts for material anisotropy. This aided in future structural design through a quantitative understanding of the mechanical limits of the material. To examine the texture, defects, fracture surfaces of the material, high-resolution microstructural imaging was utilized before and after the mechanical tests. This microstructural analysis informed us about the macro-structural behavior and influenced structural design. Thermal decomposition analysis was also performed tounderstand the thermal limits.

Biocomposites↗

Electronic effects in radiation-induced collision cascades in nickel

The accurate treatment of electronic effects in multi-million-atom simulations of radiation-induced collision cascades is crucial for reliable predictions of primary radiation damage. In this work, we explore the fidelity of a recently developed two-temperature molecular dynamics model implementing an electron density-dependent coupling of electronic and atomic subsystems for cascade simulations in nickel. We show that the parameter-free model realistically captures the instantaneous energy losses during all stages of the highly nonequilibrium cascade process. Our simulations predict two distinct coupling regimes, corresponding to the rapid energy losses through electronic stopping in the early stages of the cascade and to the slow equilibration through the electron-phonon coupling mechanism in the later stages, without the use of separate models or coupling terms. The intermediate stage of the cascade dynamics displays a complex energy transfer between the subsystems, which cannot be interpreted by comparison to either electronic stopping or electron-phonon coupling theories. We therefore compare the predicted atomic mixing, which is sensitive to the energy losses during the intermediate cascade stage, with experimental ion beam mixing measurements. We find good agreement with the experiments, validating the coupling model for the intermediate stage of the cascade. Predictions of final defect numbers and cluster sizes are found in line with the results from conventional electronic stopping-based methods, while significantly reducing the theoretical uncertainty in the outcomes of conventional models stemming from arbitrary choices of thresholds for different coupling terms. Our results represent a notable improvement in cascade damage predictions in nickel, providing validation of the electron density-dependent coupling model for radiation damage simulations in general. The results lead us to propose an interpretation of the electronic energy losses in the intermediate regime of velocities, where we find an effectively nonlinear dissipation.

Crystal defects↗