Shuxiang Zhou

A recent experimental study on a spent uranium dioxide (UO2) fuel sample from Belgium Reactor 3 identified a unique pair structure formed by the noble metal phase (NMP) and fission gas [xenon (Xe)] precipitate. However, the fundamental mechanism behind this structure remains unclear. The present study aims to provide an understanding of the interaction between five different metal precipitates [molybdenum (Mo), ruthenium (Ru), palladium (Pd), technetium (Tc), and rhodium (Rh)] and the Xe fission gas atoms in UO2, by using density functional theory (DFT) in combination with the Hubbard U correction to compute the formation energies involved. All DFT + U calculations were performed with occupation matrix control to ensure antiferromagnetic ordering of UO2. The calculated formation and binding energies of the Xe and solid fission products in the NMP reveal that these metal precipitates form stable pair structures with Xe. Notably, the formation energy of Xe–metal pairs is lower than that of the isolated single defects in all instances, with Pd and Mo showing the most favorable binding energy, likely accounting for the observed pair structure formation.

Machine learning interatomic potentials (MLIPs) have greatly enhanced molecular dynamics (MD) simulations, achieving near-first-principles accuracy in thermal conductivity studies. In this work, we reveal that this accuracy, observed in BAs and diamond at sub-Debye temperatures, stems from an accidental error cancelation: classical statistics overestimates specific heat while underestimating phonon lifetimes, balancing out in thermal conductivity predictions. However, this balance is disrupted when isotopes are introduced, leading MLIP-based MD to significantly underpredict thermal conductivity compared to experiments and quantum statistics-based Boltzmann transport equation. This discrepancy arises not from classical statistics affecting phonon–isotope scattering rates but from its impact on the interplay between phonon–isotope and phonon–phonon scattering in the normal scattering-dominated BAs and diamond. This work underscores the limitations of MLIP-based MD for thermal conductivity studies at sub-Debye temperatures.

The ultrafast dynamics of the antiferrodistortive phase transition in perovskite SrTiO3 is monitored via time-domain Brillouin scattering. Using femtosecond optical pulses, we initiate a thermally driven tetragonal-to-cubic structural transformation and detect the crystal phase through changes in the frequency of Brillouin oscillations (BO) induced by propagating acoustic phonons. Coupling the measured BO frequency with a spatiotemporal heat diffusion model, we demonstrate that, for a sample kept in the tetragonal phase, deposition of sufficient thermal energy induces a rapid transformation of the heat-affected region to the cubic phase. The initial phase change is followed by a slower reverse cubic-to-tetragonal phase transformation occurring on a timescale of hundreds of picoseconds. We attribute this ultrafast phase transformation in the perovskite to a structural resemblance between atomic displacements of the R-point soft optic mode of the cubic phase and the tetragonal phase, both characterized by anti-phase rotation of oxygen octahedra. The structural relaxation time exhibits a strong temperature dependence consistent with the prediction of the equation of motion describing collective oxygen octahedra rotation based on the energy landscape of the phenomenological Landau theory of phase transitions. Evidence of such a fast structural transition in perovskites can open up new avenues in information processing and energy storage sectors.

In this Letter, we report that the fourth-order interatomic force constants (4th-IFCs) are significantly sensitive to the energy surface roughness of exchange-correlation (XC) functionals in density functional theory calculations. This sensitivity, which is insignificant for the second- (2nd-) and third-order (3rd-) IFCs, varies for different functionals in different materials and can cause misprediction of thermal conductivity by several times of magnitude. As a result, when calculating the 4th-IFCs using the finite difference method, the atomic displacement needs to be taken large enough to overcome the energy surface roughness, in order to accurately predict phonon lifetime and thermal conductivity. We demonstrate this phenomenon on a benchmark material (Si), a high-thermal conductivity material (BAs), and a low thermal conductivity material (NaCl). For Si, we find that the LDA, PBE, and PBEsol XC functionals are all smooth to the 2nd- and 3rd-IFCs but all rough to the 4th-IFCs. This roughness can lead to a prediction of nearly one order of magnitude lower thermal conductivity. For BAs, all three functionals are smooth to the 2nd- and 3rd-IFCs, and only the PBEsol XC functional is rough for the 4th-IFCs, which leads to a 40% underestimation of thermal conductivity. For NaCl, all functionals are smooth to the 2nd- and 3rd-IFCs but rough to the 4th-IFCs, leading to a 70% underprediction of thermal conductivity at room temperature. With these observations, we provide general guidance on the calculation of 4th-IFCs for an accurate thermal conductivity prediction.

Computing vibrational properties of crystals in the presence of complex defects often necessitates the use of (semi-)empirical potentials, which are typically not well characterized for perfect crystals. Here we explore the efficacy of a commonly used embedded-atomempirical interatomic potential for the U _x Th_1−x O₂ system, to compute phonon dispersion, lifetime, and branch specific thermal conductivity. Our approach for ThO₂ involves using lattice dynamics and the linearized Boltzmann transport equation to calculate phonon transport properties based on second and third order force constants derived from the empirical potential and from first-principles calculations. For UO₂, to circumvent the accuracy issues associated with first-principles treatments of strong electronic correlations, we compare results derived from the empirical interatomic potential to previous experimental results. It is found that the empirical potential can reasonably capture the dispersion of acoustic branches, but exhibits significant discrepancies for the optical branches, leading to overestimation of phonon lifetime and thermal conductivity. The branch specific conductivity also differs significantly with either first-principles based results (ThO₂) or experimental measurements (UO₂). These findings suggest that the empirical potential needs to be further optimized for robust prediction of thermal conductivity both in perfect crystals and in the presence of complex defects.

Broadband spin‐controlled surface plasmon polariton (SPP) launching and radiation via L‐shaped optical slot nanoantennas are proposed and demonstrated experimentally. The phase retardation and spectra overlapping between two resonant plasmon modes in the L‐shaped optical slot nanoantenna lie at the origin of this effect. SPP launching in two perpendicular directions are controlled by the spin of the incident light. Broadband directional launching of SPPs is achieved and the extinction ratio keeps larger than 7 dB within a wavelength range of 150 nm. Furthermore, the photon spin of SPP radiation via the L‐shaped optical slot nanoantennas is controlled by SPP propagation directions. These investigations provide a route for spin‐controlled nanophotonic applications.