Abstract A three species food chain model with Beddington–DeAngelis functional response is investigated. The local stability analysis is carried out and global behavior is simulated numerically for a biologically feasible choice of parameters. The persistence conditions of a food chain model are established. The bifurcation diagrams are obtained for different parameters of the model after intensive numerical simulations. The results of simulations show that the model could exhibit chaotic dynamics for realistic and biologically feasible parametric values. Finally, the effect of immigration within prey species is investigated. It is observed that adding small amount of constant immigration to prey species stabilize the system.