As we can see from our Fig. that the largest intensity is at 668.5 cm -1 . Let us start fitting within its range as in the Figure . We will find the q-Gaussian shoulder at 630.8 cm -1 . ). For the deconvolution, three q-Gaussians have been used (values of the q-parameters are given in the figure). The misfit is proposed in the lower part of the plot. On the right, the same fit is shown with log scale for y-axis (semi log scale). Data and q-Gaussians are given as functions of integers n (equally spaced points used in fitting), for the x-axis which is representing the Raman shift. A convenient scale is used for the y-axis (intensity axis). The fitting calculation is obtained by minimizing the sum of the squares of the deviations (sum from n=1 to n=656). The frequencies of the q-Gaussian centres are also given in the image on the left.