- #1
Opi_Phys
- 4
- 0
Hello,
I tried to consider one-dimensional anharmonic (cubic) oscillator in second quantization formalism. And I found that energy difference between the first excited state and the ground state (as well as energy difference between the second excited stated and the the first one) as function of number of used basis functions show no convergence on this number. Should it be? The more basis functions I use the less is energy difference between the first excited and ground states. Should it be?
In more details, I start form Hamiltonian in x-p representation:
H = p^2/2 + 1/2(k + a*x)*x^2, where "a" is anharmonicity.
Then I use transition from x and p operator to creation and annihilation operators and obtain the following Hamiltonian:
H = epsilon1 * ap * am + V + Delta * (am + ap) * (am + ap) * (am + ap),
where ap - creation operator, and am - annihilation operator.
epsilon, V, Delta - some constants.
Then I rewrite Hamiltonian in the eigen states of Harmonic oscillators. It means that I multiply previous Hamiltonian (from left and right side) on unitary operator:
H_new = I * H_old * I,
where I = cket(0)*bra(0) + cket(1)*bra(1) + ... + cket(n)*bra(n).
Then new Hamiltonian I write in matrix form:
sum_{i,j=1}^{n} cket(i) * H(i,j) * bra(j)
And finally I diagonalize matrix H(i,j) to find energy levels of the system.
And I see no convergence of result as a function of "n"...
I tried to consider one-dimensional anharmonic (cubic) oscillator in second quantization formalism. And I found that energy difference between the first excited state and the ground state (as well as energy difference between the second excited stated and the the first one) as function of number of used basis functions show no convergence on this number. Should it be? The more basis functions I use the less is energy difference between the first excited and ground states. Should it be?
In more details, I start form Hamiltonian in x-p representation:
H = p^2/2 + 1/2(k + a*x)*x^2, where "a" is anharmonicity.
Then I use transition from x and p operator to creation and annihilation operators and obtain the following Hamiltonian:
H = epsilon1 * ap * am + V + Delta * (am + ap) * (am + ap) * (am + ap),
where ap - creation operator, and am - annihilation operator.
epsilon, V, Delta - some constants.
Then I rewrite Hamiltonian in the eigen states of Harmonic oscillators. It means that I multiply previous Hamiltonian (from left and right side) on unitary operator:
H_new = I * H_old * I,
where I = cket(0)*bra(0) + cket(1)*bra(1) + ... + cket(n)*bra(n).
Then new Hamiltonian I write in matrix form:
sum_{i,j=1}^{n} cket(i) * H(i,j) * bra(j)
And finally I diagonalize matrix H(i,j) to find energy levels of the system.
And I see no convergence of result as a function of "n"...