Amplitude effects on seismic velocities: How low can we go?

α-quartz is one of the most important SiO polymorphs because it is the basis of very common minerals, especially for seabed materials with geoscientific importance. The elastic characterization of these materials is particularly relevant when the properties governing phonon and sound propagation are involved. These studies are especially interesting for oil exploration purposes. Recently, we published a new method that constitutes to the best of our knowledge the first attempt to recreate longitudinal and transversal perturbations in a simulation box to observe their propagation through the crystal by means of a set of descriptors [D. Melgar et al., J. Phys. Chem. C 122, 3006-3013 (2018)]. The agreement with the experimental S- and P-wave velocities was rather excellent. Thus, an effort has been undertaken to deepen the particularities of this new methodology. Here, bearing in mind this encouraging initial methodology-development progress, we deepen our knowledge of the particularities of this new methodology in presenting a systematic investigation of the implementation of the perturbation source. This includes new ways of creating the perturbation, as well as analyzing the possible effects the perturbation amplitude could have on the resultant velocities. In addition, different force fields were tested to describe the interatomic interactions. The lack of dependence of the seismic velocities on the way the perturbation is created and the perturbation amplitude, and the good agreement with the experimental results are the main reasons that allow the definition of this new methodology as robust and reliable. These qualities are consolidated by the physical behavior of the calculated velocities in the presence of vacancies and under stress. The development of this method opens up a new line of research of calculating seismic velocities for geophysically relevant materials in a systematic way, with full control not only on the sample features (composition, porosity, vacancies, stress, etc.) but also on the particularities of perturbation itself, as well as determining optimal system-response metrics.