“…A different choice of the parameters results in a different correction to the original potential profile V (q), and the above Hamiltonian always has the following form: Hence, the typical investigations into (quasi-)exact solvability of PDM quantum systems consist in finding simultaneously a pair of an effective potential U(q) and a mass function m(q) for which the PDM Hamiltonian (1.2) admits (a number of) exact eigenfunctions in closed form. Up to now, two different methods have been frequently employed, namely, coordinate transformations including point canonical transformations [1,2,5,6,7,8,9,12,16,18,20,25], and supersymmetric methods [3,4,7,8,11,13,14,19,20,21,22,24,25]. The latter approaches were also applied to many-body PDM quantum systems [26].…”