We have two spectra: R050209 and R060124 . Samples are: Molybdenite from Boss Mountain mine, Hendrix Lake, British Columbia, Canada (Source: Lloyd Twaites, Owner: RRUFF), Measured Chemistry: Mo1.00S2.00, and Molybdenite from Crown Point mine, near Chelan, Washington, USA, (Source: University of Arizona Mineral Museum 11966, Owner: RRUFF), Measured Chemistry: Mo1.00S2.00. Samples were unoriented and the instrument was Thermo Almega XR 532nm @ 100% of 150mW. In the lower part of the Figure , we can see the two main peaks of molybdenite. In the Application Note by Angela Flack, 2022, the band are given as E2g (actually, E 1 2g) band at 385 cm -1 and A1g band at 405 cm -1 . "The E2g band is due to in-plane vibrations, and the A1g band attributed to out-of-plane vibrations" . Note the asymmetry of the peaks. We could use several fitted Lorentzians and Gaussians or q-Gaussians to decompose the peaks, but here we want to test the use of asymmetric q-Gaussians. Let us start from two q-BWF functions. To show the decomposition, let us use Fityk software . Fityk is a curve fitting and data analysis application, free and open source. It can perform a Levenberg-Marquardt algorithm, among other best fit algorithms. The function of merit is the weighted sum of squared residuals (WSSR), also called chi square. Fityk has several built-in functions but allows to execute a script with user-defined functions (see in the Appendix how we can define in Fityk q-Gaussian, split-q-Gaussian and q-BWF functions).