Motivated by the recent experimental progress in the Λc decay that contains a neutron in the final state, we analyze the semileptonic decay Λc→nℓνℓ in the framework of QCD sum rules. The transition form factors are analytically computed using three-point correlation functions and the Cutkosky cutting rules, which can be extrapolated into the physical region by employing the z-series parametrization. The branching fractions of Λc→ne+νe and Λc→nμ+νμ are estimated to be (0.281±0.056)% and (0.275±0.055)%, respectively. Furthermore, we calculate as well the relevant decay asymmetry observables sensitive to new physics beyond the standard model. The numerical results of semileptonic decays Λc→Λℓνℓ are also given and confronted to the latest experimental data.