In this paper, we study a nematic liquid crystals system in three-dimensional whole space ℝ3 and obtain the time decay rates of the higher-order spatial derivatives of the solution by the method of spectral analysis and energy estimates if the initial data belongs to L1ℝ3 additionally.