Endohedral fullerenes, or endofullerenes, are chemical systems of fullerene cages encapsulating single atoms or small molecules. These species provide an interesting challenge of Potential Energy Surface determination as examples of non-covalently bonded, bound systems. While the majority of studies focus on C60 as the encapsulating cage, introducing some anisotropy by using a different fullerene, e.g., C70 can unveil a double well potential along the unique axis. By approximating the potential as a pairwise Lennard-Jones (LJ) summation over the fixed C cage atoms, the parameter space of the Hamiltonian includes three tunable variables: (M, ɛ, σ) representing the mass of the trapped species, the LJ energy, and length scales respectively. Fixing the mass and allowing the others to vary can imitate the potentials of endohedral species trapped in more elongated fullerenes. We choose to explore the LJ parameter space of an endohedral atom in C70 with ɛ ∈ [20, 150 cm−1], and σ ∈ [2.85, 3.05 Å]. As the barrier height and positions of these wells vary between [1, 264 cm−1] and [0.35, 0.85 Å] respectively, using a 3D direct product basis of 1D harmonic oscillator (HO) wavefunctions centred at the origin where there is a local maximum is unphysical. Instead we propose the use of a non-orthogonal basis set, using 1D HO wavefunctions centred in each minimum and compare this to other choices. The ground state energy of the X@C70 is tracked across the LJ parameter space, along with its corresponding nuclear translational wavefunctions. A classification of the wavefunction characteristics, namely the prolateness and “peanut-likeness” based on its statistical moments is also proposed. Excited states of longer fullerenes are assigned quantum numbers, and the fundamental transitions of Ne@C70 are tracked across the parameter space.