Molecular dynamics simulation of shock induced ejection on fused silica surface

Shock response and surface ejection behaviors of fused silica are studied by using non-equilibrium molecular dynamics combining with the Tersoff potential. First, bulk structure and Hugoniot curves of fused silica are calculated and compared with experimental results. Then, the dynamical process of surface ejection behavior is simulated under different loading velocities ranging from 3.5 to 5.0 km∕s, corresponding to shock wave velocities from 7.1 to 8.8 km∕s. The local atomistic shear strain parameter is used to describe the local plastic deformation under conditions of shock compression or releasing. Our result shows that the shear strain is localized in the bottom area of groove under the shock compression. Surface ejection is observed when the loading velocity exceeds 4.0 km∕s. Meanwhile, the temperature of the micro-jet is ∼5574.7 K, which is close to experiment measurement. Several kinds of structural defects including non-bridging oxygen are found in the bulk area of the sample after ejection.