X.D. Yuan

Thermal conductivity of vitreous silica from molecular dynamics simulations: The effects of force field, heat flux and system size

The thermal conductivity of vitreous silica is computed using the direct method in molecular dynamics simulations with three sets of empirical force fields, including the BKS, Teter, and ReaxFF, to investigate their performance in thermal characterization. Various heat flux and system sizes are used in the simulations to evaluate the statistical uncertainty and the finite-size effect. While all these potentials can reproduce realistic silica structures, the ReaxFF provides better agreement with experiments at 300 K than the BKS and Teter, which is due to its improved description of low-frequency vibrations. Increasing the heat flux and cross-sectional area tends to reduce the calculated standard deviation induced by thermal fluctuations, thus contributing to more accurate thermal conductivity predictions.

UV-induced modification of fused silica: Insights from ReaxFF-based molecular dynamics simulations

Atomic structural modification and defect processes of fused silica resulting from UV-laser irradiation are studied by a combination of molecular dynamics (MD) simulations and the Reactive Force Field (ReaxFF). Bond state transitions by laser excitation are modeled as the result of localized recoils during energy deposition. Computations of pair distribution functions and bond angle distributions of the irradiated structure reveal that fused silica undergoes significant changes in terms of Si-O, Si-Si pair distances and Si-O-Si bond angles, which are attributed to the formation of silicon and oxygen coordination defects. It is found that nonbridging oxygen is responsible for the decreased Si-O bond length, while laser-induced five-coordinated silicon leads to small Si-O-Si bond angles in 2-membered rings.