(741) |

where the are complex numbers. Thus, the initial state is some linear superposition of the unperturbed energy eigenstates. In the absence of the time-dependent perturbation, the time evolution of the system is given by

(742) |

Now, the probability of finding the system in state at time is

(743) |

Clearly, with , the probability of finding the system in state at time is exactly the same as the probability of finding the system in this state at the initial time . However, with , we expect to vary with time. Thus, we can write

where . Here, we have carefully separated the fast phase oscillation of the eigenkets, which depends on the unperturbed Hamiltonian, from the slow variation of the amplitudes , which depends entirely on the perturbation (i.e., is constant if ). Note that the eigenkets , appearing in Equation (744), are

Schrödinger's time evolution equation yields

It follows from Equation (744) that

(746) |

We also have

(747) |

where use has been made of the time-independence of the kets . According to Equation (745), we can equate the right-hand sides of the previous two equations to obtain

(748) |

Left-multiplication by yields

where

(750) |

and

(751) |

Here, we have made use of the standard orthonormality result, . Suppose that there are linearly independent eigenkets of the unperturbed Hamiltonian. According to Equation (749), the time variation of the coefficients , which specify the probability of finding the system in state at time , is determined by coupled first-order differential equations. Note that Equation (749) is exact--we have made no approximations at this stage. Unfortunately, we cannot generally find exact solutions to this equation, so we have to obtain approximate solutions via suitable expansions in small quantities. However, for the particularly simple case of a two-state system (i.e., ), it is actually possible to solve Equation (749) without approximation. This solution is of great practical importance.