The following treatment is fairly common in literature, and often referred as time-dependent perturbation theory in a more advanced form.
Model
We assume a one-dimensional quantum harmonic oscillator of massm and chargee. The expression for the potential energy of this system is this of the harmonic oscillator. The total wavefunction is denoted by Ψ, and the spatial part of the wave function is ψ. As we deal with stationary states, the total wave function is a solution of the Schrödinger equation and reads with eigenvalue. The probability of transition from the fundamental level labelled 0 to a level labelled 1 under an electromagnetic stimulation is analysed below.
A two level model
For this situation, we write the total wave function as a linear combination for a two-levels system: The coefficients c0,1 are time-dependent. They represent the proportion of the state in the total wave function with time, thus they represent the probability of the wave-function to fall in one of the two state when an observer will collapse the wave function. As we deal with a two-level system, we have the normalisation relation :
Perturbation
The electromagnetic stimulation will be a uniform electric field, oscillating with a frequency ω. This is very similar to the semi-classical analysis of the behaviour of an atom or a molecule under a polarized electromagnetic plane wave. Thus, potential energy will be the sum of the unperturbed potential and of the perturbation and reads:
From the Schrödinger equation to ''c''1 time-dependence
On the right part, the total hamiltonian is the sum of the unperturbed hamiltonian and the external perturbation. This allows to substitute the eigenvalues of the stationary states in the total hamiltonian. Thus we write: Using the Schrödinger equation above, we end up with
Extract the ''c''1(''t'') time dependence
We use now the bra–ket notation to avoid cumbersome integrals. This reads : Then we multiply by and end up with the following The two different levels are orthogonal, so. Also we are working with normalized wave functions, so. Finally, This latter equation expresses the time variation of c1 with time. This is the crux of our calculation, since by then, we can deduce exactly its expression from the differential equation we obtained.
Solving the time-dependent differential equation
There is no proper way in general to evaluate, unless we have a precise knowledge of the two unperturbed wave function, that is to say unless we can solve the non-perturbed Schrödinger equation. In the case of the harmonic potential, the wave functions solutions of the one-dimensional quantum harmonic oscillator are known as Hermite polynomials.
We made several assumptions to get to the final result. First we suppose that c1 = 0, because at time t = 0, the interaction of the field with the matter did not start. That impose for the total wave function to be normalized that c0 = 1. We use these conditions, and we can write, at t = 0: Again, in this non-relativistic picture, we remove the time dependence outside. The quantity is called the transition moment integral. Its dimensions are · and SI units A·s·m. It can be measured experimentally, or calculated analytically if one know the expression of the spatial wave function for both the energy levels. It can be the case if we deal with an harmonic oscillator like it's the case here. We will not it : as the transition moment from the level 0 to the level 1. Finally, we end with
Solving the first order differential equation
The remaining task is to integrate this expression to obtain c1. However, we must recall from the previous approximations we made, we are here at time t = 0. So the solution we obtain from integration will be only valid as long as |c0|2 is still very close to 1, that is to say for very short time after the perturbation began to act. We suppose that the time dependent perturbation has the following form, to make the computation easier. This is a scalar quantity, as we assumed from the beginning a scalar charged particle and a one dimensional electric field. So we have to integrate the following expression : We can write and doing the variable change we obtain the correct form of the Fourier transform :
Using the Fourier transform
where is the rectangular function. We notice from the previous equation that c1 is the Fourier transform of the product of a cosine with a square of width t'. From then, the formalism of Fourier transforms will make the work easier. We have Where sinc is the cardinal sinus function in its normalized form. The convolution with the Dirac distribution will translate the term on the left of the sign. We obtain finally
Interpretation
The probability of a transition is given in general for a multi-level system by the following expression:
Final result
The probability to fall in the 1 state corresponds to. This is really easy to compute from all the tedious calculation we made previously. We observe in the equation that has a very simple expression. Indeed, the phase factor, varying with t, disappears naturally. So we obtain the expression
Conclusion
We made the hypothesis that the stimulation was a complex exponential. However a true electric field is real valued. A further analysis should take it in account. Also, we always assume that t is very small. We should keep it in mind before to conclude.
