Фазовый сдвиг одномерного рассеяния (конечная яма) — нефизический?

Я рассчитываю фазовый сдвиг из одномерной потенциальной ямы. Это кажется очень простым, но я просто запутался в этом.

Пусть имеется потенциальный колодец глубины$V_0$и пространственная ширина$2a$. Рассмотрим решение рассеяния с массовым параметром$m$и энергия$E$, инцидент слева. Волновые функции слева и справа от потенциальной ямы имеют вид

$$\psi_{-\infty}(x)=Ae^{ikx}+Be^{-ikx}$$

$$\psi_{\infty}(x)=E^{ikx}$$

Где$k=\sqrt{2mE}/\hbar$,$A$находится в произвольной комплексной амплитуде приходящей волны,$B$– комплексная амплитуда отраженной волны,$E$– комплексная амплитуда прошедшей волны (конечно, заскважинного падения). Что меня беспокоит, так это фазовый сдвиг переданной волны относительно входящей волны, поэтому я хочу$E$с точки зрения$A$. Это можно рассчитать как

$$\frac{E}{A}=\frac{e^{-2ika}}{\cos(2k'a)-i\frac{\epsilon}{2}\sin(2k'a)}$$

Где (новые) параметры$k'$и$\epsilon$даны

$$k'=\frac{\sqrt{2m(E+V_0)}}{\hbar},\,\,\,\,\epsilon=\frac{k'}{k}+\frac{k}{k'}$$

Взяв квадратный модуль, можно получить правильное выражение для коэффициента передачи (я не буду писать его здесь, но это проверка безопасности).

Вот моя проблема. Фазовый сдвиг будет просто определяться фазой$E/A$. Когда я вычисляю фазу (или просто смотрю на выражение и записываю его), я получаю

$$\delta=\arctan \left(\frac{\epsilon}{2}\tan (2k'a)\right) - 2ka$$

Вот примерный график этого (удобные единицы, случайные значения параметров):

Это говорит о том, что фазовый сдвиг неограниченно увеличивается (или плавно перемещается через$90^o$и$-90^o$. Я ожидал бы, что фазовый сдвиг будет уменьшаться с увеличением энергии, в конечном итоге устанавливаясь на нуле. Я почти уверен, что я что-то напутал. Вы не могли бы мне помочь?