13.07.2015 Views

PhD thesis in English

PhD thesis in English

PhD thesis in English

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

|2. Diagonalization of Transition Amplitudes10 5 0.001 0.01 0.1| E 0(p) (∆, L, t) - E 0exact10 -5 110 -1010 -1510 -20∆ = 0.05∆ = 0.10∆ = 0.2010 -25Figure 2.12: Deviations |E (p)k(∆, L, t) − Eexact k | as a function of t for k = 0, for themodified Pöschl-Teller potential. Energy eigenvalues are obta<strong>in</strong>ed us<strong>in</strong>g effectiveaction levels p = 1, 3, 5, 7, 9, 11, 13 and L = 5, with the parameters χ = 0.5, λ = 15.5,∆ = 0.02. Dashed l<strong>in</strong>es <strong>in</strong> correspond to the discretization error (2.19).taction (which does not depend on δ) near x = 0. As can be seen, this does not affectthe obta<strong>in</strong>ed numerical results.Table 2.3(top) gives the obta<strong>in</strong>ed energy spectra for the modified Pöschl-Tellerpotential with the parameters χ = 0.5, λ = 15.5. If necessary, the precision ofobta<strong>in</strong>ed energy levels can be further <strong>in</strong>creased by appropriately chang<strong>in</strong>g the discretizationparameters. Contrary to the situation for anharmonic oscillator, whererelative error of numerically calculated low-ly<strong>in</strong>g energy levels did not change significantly,here we see that the <strong>in</strong>crease <strong>in</strong> the error is substantial. This is causedby the fact that this potential has only a small f<strong>in</strong>ite set of discrete bound states, soenergy levels k ∼ 10 correspond to the very top of the discrete spectrum. In practicalapplications such pathological situations are not encountered, but as we can see,even this can be dealt with by the proper choice of discretization parameters. Thequality of numerically calculated eigenfunctions is assessed <strong>in</strong> Table 2.3(bottom),where we give a symmetric matrix of scalar products 〈ψ k |ψlexact 〉 of numerically calculatedand analytic eigenfunctions. As we can see, the overlap between analyticand numeric eigenfunctions is excellent, and they are orthogonal with high precision,which is preserved even for higher energy levels. We have also verified that forparameters given <strong>in</strong> the caption of Table 2.3 and with the discretization step of theorder ∆ = 10 −3 eigenfunctions of all bound states can be accurately reproduced.46

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!