Quantum harmonic oscillator with harmonic perturbation

Evgn

1. The problem statement, all variables and given/known data

Let's consider a harmonic oscillator with a harmonic perturbation:
[tex]H = \frac{p^2}{2} + \frac{x^2}{2} + a \frac{x^2}{2}.[/tex]
Exact solution is known, but we want to derive it using perturbation theory. More specifically, suppose we want to obtain a series for the ground state energy [itex]E(a)[/itex]. This series is known to be a convegent one (when [itex]|a| < 1[/itex]) with sum equal to the exact answer ([itex]1/2 \sqrt{1+a}[/itex]). We want to comfirm this fact by a direct calculation.

2. Relevant equations and 3. The attempt at a solution

Well-knows recurrence relations of Rayleigh-Schrodinger theory are applicable, and I've manually checked first five orders. The calculations quickly become very cumbersome, and there seems to be no easy way to obtain a general expression for n-th order (known in advance from the exact answer).

What is the most straightforward way to obtain this series? Will the Kato's theory be of any help in this case?

Thank you in advance.

dextercioby

Are you using a ladder operator approach ? The calculations should be the simplest.

Evgn

Sure, I use them to calculate matrix elements of perturbation. This part of calculation is indeed simple. What seems to be not simple is the calculation of the n-th order of perturbation series and not just 1st and 2nd like most textbooks do.

1st order correction is [itex]A |2\rangle[/itex], 2nd order - [itex]B |2\rangle + C |4\rangle[/itex], 3rd order - [itex]D |2\rangle + F |4\rangle + G |6\rangle[/itex] and so on. The number of terms grows with n.

There should be an elegant way to do the calculation even in the framework of an "old" perturbation theory. I deliberately avoid using diagramms here.