The quantum version of a non-linear oscillator, previouly analyzed at the classical level, is studied. This is a problem of quantization of a system with position-dependent mass of the form m = (1 + λx 2 ) −1 and with a λ-dependent nonpolynomial rational potential. This λ-dependent system can be considered as a deformation of the harmonic oscillator in the sense that for λ → 0 all the characteristics of the linear oscillator are recovered. Firstly, the λ-dependent Schrödinger equation is exactly solved as a Sturm-Liouville problem and the λ-dependent eigenenergies and eigenfunctions are obtained for both λ > 0 and λ < 0. The λ-dependent wave functions appear as related with a family of orthogonal polynomials that can be considered as λ-deformations of the standard Hermite polynomials. In the second part, the λ-dependent Schrödinger equation is solved by using the Schrödinger factorization method, the theory of intertwined Hamiltonians and the property of shape invariance as an approach. Finally, the new family of orthogonal polynomials is studied. We prove the existence of a λ-dependent Rodrigues formula, a generating function and λ-dependent recursion relations between polynomials of different orders.