We construct the fundamental solution of time dependent linear ordinary Dirac system in terms of unknown phase functions. This construction gives approximate representation of solutions which is useful for the study of asymptotic behavior. Introducing analog of Rayleigh quotient for differential equations we generalize Hartman–Wintner asymptotic integration theorems with the error estimates for applications to the Dirac system. We also introduce the adiabatic invariants for the Dirac system, which are similar to the adiabatic invariant of Lorentz’s pendulum. Using a small parameter method it is shown that the change in the adiabatic invariants approaches zero with the power speed as a small parameter approaches zero. As another application we calculate the transition probabilities for the Dirac system. We show that for the special choice of electromagnetic field, the only transition of an electron to the positron with the opposite spin orientation is possible.