Yongfeng Zhang

Molten salts are promising for various energy applications including fuel and solar cells and nuclear energy. These applications face a common challenge: corrosion of structural materials by impurities such as H₂O. This work employs ab-initio molecular dynamics simulations to study H₂O induced corrosion of FeCr alloys in molten NaF and NaCl salts. H₂O is found highly stable in both salts, with infrequent, reversible dissociation into OH⁻ and H⁺ along with HF or HCl formation. The dissociation tendency correlates positively with the electronegativity and negatively with the size of halogen atoms. Accordingly, H₂O reaches the salt/metal interface as a molecule before reacting with metal. Reduction of H⁺ is found to occur without simultaneous oxidation of specific metal atoms such as Cr, suggesting sequential instead of the commonly proposed concurrent reduction and oxidation. The reduced H atoms prefer to stay at the interface and may re-enter NaF but not NaCl, highlighting the influence of salt chemistry.

Uranium dioxide (UO₂) is the primary fuel material that is used in current nuclear reactors. As one of the most fundamental material parameters, grain boundary (GB) energy strongly influences many fuel properties, and the influences depend on the characters and properties of individual GBs. Using molecular dynamics simulations, a high throughput survey of GB energy in UO₂ was carried out for the purpose of elucidating the roles of GB geometry such as misorientation and inclination, as well as the bonding nature of UO₂, in affecting GB energy. GB energies in CeO₂ were calculated as well for comparison with UO₂ to investigate the generality of GB energy anisotropy in fluorite phase oxides. The results show significant GB energy anisotropy in both UO₂ and CeO₂ that is associated with the cubic symmetry of the fluorite structure. More interestingly, the GB anisotropy is found to be dependent not only on the crystal structure but also the ionic bonding. As such, the GB energy anisotropy in fluorite oxides has significant differences compared with that in fcc metals. The data obtained and the increased knowledge on GB anisotropy will facilitate GB engineering for nuclear fuels with improved properties.

Despite their exceptional mechanical and corrosion properties, duplex stainless steels (DSS) have not found widespread use in high-temperature applications due to concerns over thermal aging and embrittlement at elevated operational temperatures (> 300 °C). The present study investigated the effect of thermal aging time on the electrochemical properties of lean and standard grade DSS that are exposed to a range of pressurized water reactors containing LiOH and H₃BO₃. The results indicated that the electrolyte chemistry plays a significant role in the corrosion behavior of the DSS alloys. Corrosion resistance decreased with thermal aging time for all DSS alloys; however, standard grade DSS (2205 and 2209-w) alloys showed better corrosion resistance than lean grades (2003, 2101, 2101-w). The presence of dissolved oxygen in the electrolytes resulted in a significant increase in corrosion rate for the DSS alloys, but it did not affect the general trend of corrosion rates with aging time. All DSS alloys became vulnerable to pitting corrosion due to chloride addition, but the pitting resistance decreased with increasing thermal aging time. Increased boron B content resulted in degradation of corrosion resistance of the DSS alloys, while minor changes in pH did not show a significant change in corrosion resistance. Mechanical and metallurgical characterization coupled with electrochemical characterization of the DSS alloys gave a comprehensive insight into the effects of thermal aging on the electrochemical response of the DSS.

Graphical abstract

Materials performance can be significantly degraded due to bubble generation. In this work, the bubble growth process is elaborated in Cu by atomistic modeling to bridge the gap of experimental observations. Upon continuous He implantation, bubble growth is accommodated first by nucleation of dislocation network from bubble surface, then formation of dissociated prismatic dislocation loop (DPDL), and final DPDL emission in $$\langle 110\rangle$$ ⟨ 110 ⟩ directions. As the DPDL is found capable of collecting He atoms, this process is likely to assist the formation of self-organized bubble superlattice, which has been reported from experiments. Moreover, the pressurized bubble in solid state manifests the shape of an imperfect octahedron, built by Cu $$\{111\}$$ { 111 } surfaces, consistent with experiments. These atomistic details integrating experimental work fill the gap of mechanistic understanding of athermal bubble growth in Cu. Importantly, by associating with nanoindentation testings, DPDL punching by bubble growth arguably applies to various FCC (face-centered cubic) metals such as Au, Ag, Ni, and Al.

The ion exchange and point defect models are two prominent models describing the role of anions, such as chlorides, in the degradation of passive oxide films. Here the thermodynamic feasibility of critical steps of Cl-induced degradation of a hydroxylated α-Cr₂O₃ (0001) surface, as proposed by these two models, are studied. Both models begin with Cl substitution of surface OH and H₂O, which becomes less favorable with increasing Cl coverage. The initial stages of Cl-induced breakdown of the α-Cr₂O₃ depend on Cl coverage and the presence of O vacancy near the surface as follows: (1) neither Cl insertion (supporting the ion exchange model) nor Cr vacancy formation (supporting the point defect model) is feasible at low Cl coverages except in the presence of O vacancies near the surface, where Cl insertion is thermodynamically feasible even at low coverages, (2) in the absence of O vacancies, Cr vacancy formation becomes feasible from 10/12 ML onwards whereas Cl insertion by exchange with subsurface OH only becomes feasible at full coverage. This implies that at higher coverages Cl-induced degradation first initiatesthrough a vacancy formation mechanism, but both insertion and vacancy formation would be feasible at full coverage.

The ion exchange and point defect models are two prominent models describing the role of anions, such as chlorides, in the degradation of passive oxide films. Here the thermodynamic feasibility of critical steps of Cl-induced degradation of a hydroxylated α-Cr₂O₃ (0001) surface, as proposed by these two models, are studied. Both models begin with Cl substitution of surface OH and H₂O, which becomes less favorable with increasing Cl coverage. The initial stages of Cl-induced breakdown of the α-Cr₂O₃ depend on Cl coverage and the presence of O vacancy near the surface as follows: (1) neither Cl insertion (supporting the ion exchange model) nor Cr vacancy formation (supporting the point defect model) is feasible at low Cl coverages except in the presence of O vacancies near the surface, where Cl insertion is thermodynamically feasible even at low coverages, (2) in the absence of O vacancies, Cr vacancy formation becomes feasible from 10/12 ML onwards whereas Cl insertion by exchange with subsurface OH only becomes feasible at full coverage. This implies that at higher coverages Cl-induced degradation first initiatesthrough a vacancy formation mechanism, but both insertion and vacancy formation would be feasible at full coverage.

Ordering and self-organization are critical in determining the dynamics of reaction-diffusion systems. Here we show a unique pattern formation mechanism, dictated by the coupling of thermodynamic instability and kinetic anisotropy. Intrinsically different from the physical origin of Turing instability and patterning, the ordered patterns we obtained are caused by the interplay of the instability from uphill diffusion, the symmetry breaking from anisotropic diffusion, and the reactions. To understand the formation of the void/gas bubble superlattices in crystals under irradiation, we establish a general theoretical framework to predict the symmetry selection of superlattice structures associated with anisotropic diffusion. Through analytical study and phase field simulations, we found that the symmetry of a superlattice is determined by the coupling of diffusion anisotropy and the reaction rate, which indicates a new type of bifurcation phenomenon. Our discovery suggests a means for designing target experiments to tailor different microstructural patterns.

The generation and motion of crystalline defects during plastic deformation are critical processes that determine the mechanical properties of a crystal. The types of defect generated are not only related to the symmetry of a crystal but also associated with the symmetry-breaking process during deformation. Proposed here is a new mathematical framework to capture the intrinsic coupling between crystal symmetry and deformation-induced symmetry breaking. Using a combination of group theory and graph theory, a general approach is demonstrated for the systematic determination of the types of crystalline defect induced by plastic deformation, through the construction of a crystal deformation group and a deformation pathway graph. The types of defect generated in the deformation of a face-centered cubic crystal are analyzed through the deformation pathway graph and compared with experimental observations.

Systematic ab initio calculations based on density functional theory have been performed to gain fundamental understanding of the interactions between noble gas atoms (He, Ne, Ar and Kr) and bcc transition metals in groups 5B (V, Nb and Ta), 6B (Cr, Mo and W) and 8B (Fe).

Nano-structured superlattices may have novel physical properties and irradiation is a powerful mean to drive their self-organization. However, the formation mechanism of superlattice under irradiation is still open for debate. Here we use atomic kinetic Monte Carlo simulations in conjunction with a theoretical analysis to understand and predict the self-organization of nano-void superlattices under irradiation, which have been observed in various types of materials for more than 40 years but yet to be well understood. The superlattice is found to be a result of spontaneous precipitation of voids from the matrix, a process similar to phase separation in regular solid solution, with the symmetry dictated by anisotropic materials properties such as one-dimensional interstitial atom diffusion. This discovery challenges the widely accepted empirical rule of the coherency between the superlattice and host matrix crystal lattice. The atomic scale perspective has enabled a new theoretical analysis to successfully predict the superlattice parameters, which are in good agreement with independent experiments. The theory developed in this work can provide guidelines for designing target experiments to tailor desired microstructure under irradiation. It may also be generalized for situations beyond irradiation, such as spontaneous phase separation with reaction.

This report presents an accelerated kinetic Monte Carlo (KMC) method to compute the diffusivity of hydrogen in hcp metals and alloys, considering both thermally activated hopping and quantum tunneling. The acceleration is achieved by replacing regular KMC jumps in trapping energy basins formed by neighboring tetrahedral interstitial sites, with analytical solutions for basin exiting time and probability. Parameterized by density functional theory (DFT) calculations, the accelerated KMC method is shown to be capable of efficiently calculating hydrogen diffusivity in α-Zr and Zircaloy, without altering the kinetics of long-range diffusion. Above room temperature, hydrogen diffusion in α-Zr and Zircaloy is dominated by thermal hopping, with negligible contribution from quantum tunneling. The diffusivity predicted by this DFT + KMC approach agrees well with that from previous independent experiments and theories, without using any data fitting. The diffusivity along <c> is found to be slightly higher than that along <a>, with the anisotropy saturated at about 1.20 at high temperatures, resolving contradictory results in previous experiments. Demonstrated using hydrogen diffusion in α-Zr, the same method can be extended for on-lattice diffusion in hcp metals, or systems with similar trapping basins.