Time dependent perturbation theory

Time dependent perturbation theory addresses the transition between one eigenstate to another eigenstate when there is a small perturbation with respect to time such that it doesn't change the eigenfunctions.

The Hamiltonian for these problems of perturbation is, $$ H = H_0 + \lambda {H'(t)}$$ where now the perturbing Hamiltonian is time dependent. Solution for $H_0$ is, $$\psi_n(x,t) = \phi_n(x) e^{i\frac{E_n}{\hbar}t} = \phi_n e^{-i\omega_nt}$$ where $$H_0\phi_n = E^0_n\phi_n$$ and $ \omega_n = E_n/\hbar$ and the general state vector is given by, $$ \Psi = \sum_nc_n\psi_n(x,t)$$
In our problem we start with a particular eigenstate and assume that, after switching on the perturbation the state evolves to a general state always expressed as the linear combination of initial eigenstates where the time dependence is accompanied with the coefficient as $$\Psi (x,t) = \sum_n c_n(t) \phi_n(x) e^{i\omega{t}}$$  Because of this assumption, the particle now can be found in any other eigenstate and needn't to be in the same eigenstate forever. As usual, the modulus square of the coefficient gives the probability of finding the particle in any one of eigenstates, except now it changes with perturbation time. Substituting our wave function in general Schrodinger's equation (to find the time dependence of c(t)), $$i\hbar\frac{\partial\Psi}{\partial{t}} = H\Psi$$ gives (taking the inner product with $\psi_k$)$$ i\hbar\dot{c_k} = \lambda\sum_n \langle\phi_k\vert{H'}\vert\phi_n\rangle e^{i(\omega_k-\omega_n)t} f(t) c_n$$ where f(t) is the time-dependent part of the perturbing Hamiltonian H' (in variable separable form).  In the limit $\lambda\rightarrow{0}$ the coefficients are all constant in time and so we seek a solution of the form, $$ c_k(t) = c_k^0 + \lambda{c_k^1(t)} + \lambda^2c_k^2(t) + \cdots$$ with corresponding order of correction in the upper index.  With this expansion we get, $$ i\hbar \dot{c_k}^0 = 0 $$ which means $c_k^0 = constant$ whatever the value of $c_n^0$ is at the beginning of the perturbation remains a constant throughout the perturbation. In specific, if we denote the starting time $t_0$ then $c_n^0(t_0) = \delta_{nl}$ where we assumed the state is in a specific eigenstate at the start.

Similarly, the first order equation gives, $$ i\hbar\dot{c_k}^1 = \sum_n H'_{kn} f(t) c_n^0 e^{i(\omega_k-\omega_n)t}$$ substituting for $c_n^0$ we get, $$ i\hbar\dot{c_k}^1 = H'_{kl}f(t)e^{i(\omega_k-\omega_l)t}$$ finally gives the expression for first order correction as,$$c_k^1(t) = \frac{H'_{kl}}{i\hbar} \int_{t_0}^t f(t) e^{i\omega_{kl}t'} \,dt'$$ where $\omega_{kl} = \omega_k-\omega_l$ From this we can transition probability for each eigenstate by taking modulus square of the above coefficient.

Time Independent Perturbation Theory - Degenerate

Similar to the non-degenerate case, we take our Hamiltonian $$H = H_0 + H'$$ except now, the initial unperturbed Hamiltonian has degenerate states. Assuming it is q- fold degenerate, if the symmetry giving rise to the degeneracy is disturbed, it will produce distinct 'q' non degenerate states. Our motive is to calculate the energy levels of these newly produced distinct energy levels. Let us consider two different initial eigenstates $E_1^0$ and $E_{q+1}^0$ [it is equal to the zeroth order term in the expansion] where $E_1^0$ is q-fold degenerate with states $$E_1^0 = E_2^0 = .... = E_q^0$$. As we have seen in non-degenerate perturbation theory, the coefficients $$c_{jn} = \frac{H'_{jn}}{E_n^0 - E_j^0}\,\,\,j\neq{n}\,\,and\,\, j,n\leq{q}$$ will become infinite for two different energy levels which has the same energy value. This situation can be avoided only if the numerator term is zero for all $j,n \leq {q}$ (Not the diagonal terms because $j\neq{n}$). It is equivalent to diagonalizing the sub-matrix $H'_{jn}\,\,\,\,for j,n\leq{q}$. This can be achieved by transforming to the new set of eigenfunctions from the initial set of eigenfunctions.
Thus, our problem for solving the energy values of the degenerate case reduced to diagonalizing the sub-matrix of H'. The diagonal elements are the change in energy levels of the initial degenerate states due to the perturbation. 
Representing the 'q' new eigenfunctions that diagonalize the $H'_{jn}$ as $\tilde{\psi_n}$ we write, $$ \tilde{\psi_n} = \sum_{m=1}^q a_{nm} \psi_m^0$$ Since, it will diagonalize, $$ \langle\tilde{\psi_n}\vert\tilde{\psi_j}\rangle = H'_{nj}\delta_{nj} $$ So, our new set of basis functions are, $$ Basis \,= \left\{\tilde{\psi_1},\tilde{\psi_2},\cdots\tilde{\psi_q},\psi_{q+1}^0,\cdots\right\}$$The diagonal element of the sub-matrix H' is the first order correction to the energy i.e. $$E'_n = \langle\tilde{\psi_n}\vert\tilde{\psi_n}\rangle = H'_{nn} \,\,n\leq{q}$$
Corresponding eigenfunctions are determined from $$H'\tilde{\psi_n} = E'_n\tilde{\psi_n}$$ substituting for new eigenfunctions in terms of initial wavefunctions, $$H' \sum_{m=1}^q a_{nm}\psi_m^0 = E'_n \sum_{m=1}^q a_{nm} \psi_m^0$$ Taking inner product with $\psi_j^0$ we have, $$\sum_{m=1}^q \left(H'_{jm} - E'_n\delta_{jm}\right)a_nm = 0\,\,\,n,j\leq{q} $$ It could be written simply in matrix form. For non-trivial solution for {a_nm} the determinant of coefficient matrix should vanish gives the secular equation in this case as, $$determinant\,\,\vert{H'_{jm}}-E'_n\mathbb{I}\vert = 0 $$ Using these the new eigenfunctions and eigenvalues are calculated.