D.34 The adi­a­batic the­o­rem

Con­sider the Schrö­din­ger equa­tion

\begin{displaymath}
{\rm i}\hbar \frac{\partial\Psi}{\partial t} = H \Psi
\end{displaymath}

If the Hamil­ton­ian is in­de­pen­dent of time, the so­lu­tion can be writ­ten in terms of the Hamil­ton­ian en­ergy eigen­val­ues $E_{\vec n}$ and eigen­func­tions $\psi_{\vec n}$ as

\begin{displaymath}
\Psi = \sum_{\vec n}c_{\vec n}(0) e^{{\rm i}\theta_{\vec n}...
...vec n}
\qquad
\theta_{\vec n}= - \frac{1}{\hbar} E_{\vec n}t
\end{displaymath}

Here ${\vec n}$ stands for the quan­tum num­bers of the eigen­func­tions and the $c_{\vec n}(0)$ are ar­bi­trary con­stants.

How­ever, the Hamil­ton­ian varies with time for the sys­tems of in­ter­est here. Still, at any given time its eigen­func­tions form a com­plete set. So it is still pos­si­ble to write the wave func­tion as a sum of them, say like

\begin{displaymath}
\Psi = \sum_{\vec n}\bar c_{\vec n}e^{{\rm i}\theta_{\vec n...
... \theta_{\vec n}= -\frac{1}{\hbar} \int E_{\vec n}{ \rm d}t %
\end{displaymath} (D.18)

But the co­ef­fi­cients $\bar{c}_{\vec n}$ can no longer be as­sumed to be con­stant like the $c_{\vec n}(0)$. They may be dif­fer­ent at dif­fer­ent times.

To get an equa­tion for their vari­a­tion, plug the ex­pres­sion for $\Psi$ in the Schrö­din­ger equa­tion. That gives:

\begin{displaymath}
{\rm i}\hbar\sum_{\vec n}\bar c_{\vec n}^{ \prime}e^{{\rm ...
...\vec n}\bar c_{\vec n}e^{{\rm i}\theta_{\vec n}} \psi_{\vec n}
\end{displaymath}

where the primes in­di­cate time de­riv­a­tives. The mid­dle sum in the left hand side and the right hand side can­cel against each other since by de­f­i­n­i­tion $\psi_{\vec n}$ is an eigen­func­tion of the Hamil­ton­ian with eigen­value $E_{\vec n}$. For the re­main­ing two sums, take an in­ner prod­uct with an ar­bi­trary eigen­func­tion $\langle\psi_{\underline{\vec n}}\vert$:

\begin{displaymath}
{\rm i}\hbar \bar c_{\underline{\vec n}}^{ \prime}e^{{\rm ...
...angle\psi_{\underline{\vec n}}\vert \psi_{\vec n}'\rangle
= 0
\end{displaymath}

In the first sum only the term ${\vec n}$ $\vphantom0\raisebox{1.5pt}{$=$}$ $\underline{\vec n}$ sur­vived be­cause of the or­tho­nor­mal­ity of the eigen­func­tions. Di­vide by ${\rm i}{\hbar}e^{{\rm i}\theta_{\underline{\vec n}}}$ and re­arrange to get
\begin{displaymath}
\bar c_{\underline{\vec n}}^{ \prime}
= - \sum_{\vec n}e^...
...underline{\vec n}}\vert\psi_{\vec n}'\rangle \bar c_{\vec n} %
\end{displaymath} (D.19)

This is still ex­act.

How­ever, the pur­pose of the cur­rent de­riva­tion is to ad­dress the adi­a­batic ap­prox­i­ma­tion. The adi­a­batic ap­prox­i­ma­tion as­sumes that the en­tire evo­lu­tion takes place very slowly over a large time in­ter­val $T$. For such an evo­lu­tion, it helps to con­sider all quan­ti­ties to be func­tions of the scaled time vari­able $t$$\raisebox{.5pt}{$/$}$$T$. Vari­ables change by a fi­nite amount when $t$ changes by a fi­nite frac­tion of $T$, so when $t$$\raisebox{.5pt}{$/$}$$T$ changes by a fi­nite amount. This im­plies that the time de­riv­a­tives of the slowly vary­ing quan­ti­ties are nor­mally small, of or­der 1/$T$.

Con­sider now first the case that there is no de­gen­er­acy, in other words, that there is only one eigen­func­tion $\psi_{\vec n}$ for each en­ergy $E_{\vec n}$. If the Hamil­ton­ian changes slowly and reg­u­larly in time, then so do the en­ergy eigen­val­ues and eigen­func­tions. In par­tic­u­lar, the time de­riv­a­tives of the eigen­func­tions in (D.19) are small of or­der 1/$T$. It then fol­lows from the en­tire equa­tion that the time de­riv­a­tives of the co­ef­fi­cients are small of or­der 1$\raisebox{.5pt}{$/$}$$T$ too.

(Re­call that the square mag­ni­tudes of the co­ef­fi­cients give the prob­a­bil­ity for the cor­re­spond­ing en­ergy. So the mag­ni­tude of the co­ef­fi­cients is bounded by 1. Also, for sim­plic­ity it will be as­sumed that the num­ber of eigen­func­tions in the sys­tem is fi­nite. Oth­er­wise the sums over ${\vec n}$ might ex­plode. This book rou­tinely as­sumes that it is good enough to ap­prox­i­mate an in­fi­nite sys­tem by a large-enough fi­nite one. That makes life a lot eas­ier, not just here but also in other de­riva­tions like {D.18}.)

It is con­ve­nient to split up the sum in (D.19):

\begin{displaymath}
\bar c_{\underline{\vec n}}^{ \prime}
= - \langle\psi_{\u...
...underline{\vec n}}\vert\psi_{\vec n}'\rangle \bar c_{\vec n} %
\end{displaymath} (D.20)

Un­der the stated con­di­tions, the fi­nal sum can be ig­nored.

How­ever, that is not be­cause it is small due to the time de­riv­a­tive in it, as one ref­er­ence claims. While the time de­riv­a­tive of $\psi_{\vec n}$ is in­deed small of or­der 1/$T$, it acts over a time that is large of or­der $T$. The sum can be ig­nored be­cause of the ex­po­nen­tial in it. As the de­f­i­n­i­tion of $\theta_{\vec n}$ shows, it varies on the nor­mal time scale, rather than on the long time scale $T$. There­fore it os­cil­lates many times on the long time scale; that causes op­po­site val­ues of the ex­po­nen­tial to largely can­cel each other.

To show that more pre­cisely, note that the for­mal so­lu­tion of the full equa­tion (D.20) is, [41, 19.2]:

\begin{displaymath}
\bar c_{\underline{\vec n}}(t) = e^{{\rm i}\gamma_{\underli...
...i_{\underline{\vec n}}\vert\psi_{\underline{\vec n}}'\rangle %
\end{displaymath} (D.21)

To check this so­lu­tion, you can just plug it in. Note in do­ing so that the in­te­grands are taken to be func­tions of $\bar{t}$, not $t$.

All the in­te­grals are neg­li­gi­bly small be­cause of the rapid vari­a­tion of the first ex­po­nen­tial in them. To ver­ify that, rewrite them a bit and then per­form an in­te­gra­tion by parts:

\begin{eqnarray*}
\lefteqn{
\int_{\bar t=0}^t - \frac{{\rm i}}{\hbar}(E_{\vec ...
...}(E_{\vec n}-E_{\underline{\vec n}})}
\right)'
{ \rm d}\bar t
\end{eqnarray*}

The first term in the right hand side is small of or­der 1$\raisebox{.5pt}{$/$}$$T$ be­cause the time de­riv­a­tive of the wave func­tion is. The in­te­grand in the sec­ond term is small of or­der 1$\raisebox{.5pt}{$/$}$$T^2$ be­cause of the two time de­riv­a­tives. So in­te­grated over an or­der $T$ time range, it is small of or­der 1$\raisebox{.5pt}{$/$}$$T$ like the first term. It fol­lows that the in­te­grals in (D.21) be­come zero in the limit $T\to\infty$.

And that means that in the adi­a­batic ap­prox­i­ma­tion

\begin{displaymath}
\bar c_{{\vec n}} = c_{{\vec n}}(0) e^{{\rm i}\gamma_{{\vec...
...t \langle\psi_{{\vec n}}\vert\psi_{{\vec n}}'\rangle{ \rm d}t
\end{displaymath}

The un­der­bar used to keep $\underline{\vec n}$ and ${\vec n}$ apart is no longer needed here since only one set of quan­tum num­bers ap­pears. This ex­pres­sion for the co­ef­fi­cients can be plugged in (D.18) to find the wave func­tion $\Psi$. The con­stants $c_{{\vec n}}(0)$ de­pend on the ini­tial con­di­tion for $\Psi$. (They also de­pend on the choice of in­te­gra­tion con­stants for $\theta_{{\vec n}}$ and $\gamma_{{\vec n}}$, but nor­mally you take the phases zero at the ini­tial time).

Note that $\gamma_{{\vec n}}$ is real. To ver­ify that, dif­fer­en­ti­ate the nor­mal­iza­tion re­quire­ment to get

\begin{displaymath}
\langle\psi_{\vec n}\vert\psi_{\vec n}\rangle = 1
\quad\Lo...
...}\rangle +
\langle\psi_{\vec n}\vert\psi_{\vec n}'\rangle = 0
\end{displaymath}

So the sum of the in­ner prod­uct plus its com­plex con­ju­gate are zero. That makes it purely imag­i­nary, so $\gamma_{\vec n}$ is real.

Since both $\gamma_{{\vec n}}$ and $\theta_{{\vec n}}$ are real, it fol­lows that the mag­ni­tudes of the co­ef­fi­cients of the eigen­func­tions do not change in time. In par­tic­u­lar, if the sys­tem starts out in a sin­gle eigen­func­tion, then it stays in that eigen­func­tion.

So far it has been as­sumed that there is no de­gen­er­acy, at least not for the con­sid­ered state. How­ever it is no prob­lem if at a fi­nite num­ber of times, the en­ergy of the con­sid­ered state crosses some other en­ergy. For ex­am­ple, con­sider a three-di­men­sion­al har­monic os­cil­la­tor with three time vary­ing spring stiff­nesses. When­ever any two stiff­nesses be­come equal, there is sig­nif­i­cant de­gen­er­acy. De­spite that, the given adi­a­batic so­lu­tion still ap­plies. (This does as­sume that you have cho­sen the eigen­func­tions to change smoothly through de­gen­er­acy, as per­tur­ba­tion the­ory says you can, {D.79}.)

To ver­ify that the so­lu­tion is in­deed still valid, cut out a time in­ter­val of size $\delta{T}$ around each cross­ing time. Here $\delta$ is some num­ber still to be cho­sen. The parts of the in­te­grals in (D.21) out­side of these in­ter­vals have mag­ni­tudes $\varepsilon(T,\delta)$ that be­come zero when $T\to\infty$ for the same rea­sons as be­fore. The parts of the in­te­grals cor­re­spond­ing to the in­ter­vals can be es­ti­mated as no more than some fi­nite mul­ti­ple of $\delta$. The rea­son is that the in­te­grands are of or­der 1$\raisebox{.5pt}{$/$}$$T$ and they are in­te­grated over ranges of size $\delta{T}$. All to­gether, that is enough to show that the com­plete in­te­grals are less than say 1%; just take $\delta$ small enough that the in­ter­vals con­tribute no more than 0.5% and then take $T$ large enough that the re­main­ing in­te­gra­tion range con­tributes no more than 0.5% too. Since you can play the same game for 0.1%, 0.01% or any ar­bi­trar­ily small amount, the con­clu­sion is that for in­fi­nite $T$, the con­tri­bu­tion of the in­te­grals be­comes zero. So in the limit $T\to\infty$, the adi­a­batic so­lu­tion ap­plies.

Things change if some en­ergy lev­els are per­ma­nently de­gen­er­ate. Con­sider an har­monic os­cil­la­tor for which at least two spring stiff­nesses are per­ma­nently equal. In that case, you need to solve for all co­ef­fi­cients at a given en­ergy level $E_{\underline{\vec n}}$ to­gether. To fig­ure out how to do that, you will need to con­sult a book on math­e­mat­ics that cov­ers sys­tems of or­di­nary dif­fer­en­tial equa­tions. In par­tic­u­lar, the co­ef­fi­cient $\bar{c}_{\underline{\vec n}}$ in (D.21) gets re­placed by a vec­tor of co­ef­fi­cients with the same en­ergy. The scalar $\gamma_{\underline{\vec n}}$ be­comes a ma­trix with in­dices rang­ing over the set of co­ef­fi­cients in the vec­tor. Also, $e^{{\rm i}\gamma_{\underline{\vec n}}}$ gets re­placed by a fun­da­men­tal so­lu­tion ma­trix, a ma­trix con­sist­ing of in­de­pen­dent so­lu­tion vec­tors. And $e^{-{\rm i}\gamma_{\underline{\vec n}}}$ is the in­verse ma­trix. The sum no longer in­cludes any of the co­ef­fi­cients of the con­sid­ered en­ergy.

More re­cent de­riva­tions al­low the spec­trum to be con­tin­u­ous, in which case the nonzero en­ergy gaps $E_{\underline{\vec n}}-E_{\vec n}$ can no longer be as­sumed to be larger than some nonzero amount. And un­for­tu­nately, as­sum­ing the sys­tem to be ap­prox­i­mated by a fi­nite one helps only par­tially here; an ac­cu­rate ap­prox­i­ma­tion will pro­duce very closely spaced en­er­gies. Such prob­lems are well out­side the scope of this book.