numerical integration – Instability in NIntegrate

enter image description here
I am solving the following equation with NIntegrate. However, the solution branches off at a certain time instance. When I increase the time step size, the instability occurs at larger time instant. Any solution for getting rid of the instability? NB: I have already changed AccuracyGoal and PrecisionGoal but does not work.