Many PDEs involving fractional Laplacian are naturally set in unbounded domains with underlying solutions decay very slowly, subject to certain power laws. Their numerical solutions are under-explored. In this talk we will develop spectral Galerkin/collocation methods using rational basis (or modified mapped Gegenbauer functions) for such models in unbounded domains.