Abstract
In this paper we consider a tumor-immune system interaction model with immune response delay, in which a nonmonotonic function is used to describe immune response to the tumor burden and a time delay is used to represent the time for the immune system to respond and take effect. It is shown that the model may have one, two or three tumor equilibria, respectively, under different conditions. Time delay can only affect the stability of the low tumor equilibrium, and local Hopf bifurcation occurs when the time delay passes through a critical value. The direction and stability of the bifurcating periodic solutions are also determined. Moreover, the global existence of periodic solutions is established by using a global Hopf bifurcation theorem. We also observe the existence of relaxation oscillations and complex oscillating patterns driven by the time delay. Numerical simulations are presented to illustrate the theoretical results.