Elastic properties are an important issue in explaining the behavior of seisms and to ascertain the mineralogical composition of the Earth’s shells through which seismic waves pass. Computational methods can yield an additional, detailed, free-of-heterogeneity model knowledge of the mineral series of interest. Therefore, a computational study on the influence of the interlayer cation in the muscovite-paragonite (Ms-Pg) series on the crystal cell, internal geometry, and the elastic properties was made to shed light on the mineralogical, geophysical, and geochemical properties of the series. These properties have been calculated by means of Density Functional Theory (SIESTA2.0.2 code). The crystal structure and internal geometry agreed with the range of experimental values in the literature. In general, elastic stiffness constants (EC) agreed with the known experimental values. ECs of different interlayer cation configurations for the middle concentration sample showed very similar values, except for