<p style='text-indent:20px;'>The aim of this work is to continue the analysis, started in [<xref ref-type="bibr" rid="b10">10</xref>], of the dynamics of a point-mass particle <inline-formula><tex-math id="M1">\begin{document}$ P $\end{document}</tex-math></inline-formula> moving in a galaxy with an harmonic biaxial core, in whose center sits a Keplerian attractive center (e.g. a Black Hole). Accordingly, the plane <inline-formula><tex-math id="M2">\begin{document}$ \mathbb{R}^2 $\end{document}</tex-math></inline-formula> is divided into two complementary domains, depending on whether the gravitational effects of the galaxy's mass distribution or of the Black Hole prevail. Thus, solutions alternate arcs of Keplerian hyperbolæ with harmonic ellipses; at the interface, the trajectory is refracted according to Snell's law. The model was introduced in [<xref ref-type="bibr" rid="b11">11</xref>], in view of applications to astrodynamics.</p><p style='text-indent:20px;'>In this paper we address the general issue of periodic and quasi-periodic orbits and associated caustics when the domain is a perturbation of the circle, taking advantage of KAM and Aubry-Mather theories.</p>