A question with a time dependent Hamiltonian

dEdt
Messages
286
Reaction score
2
To cut to the chase, I have to solve for the evolution of a two-state system where the system's state at time t satisfies the equation
\mathrm{i}\hbar\left( <br /> \begin{array}{cc} <br /> \dot{c}_1(t)\\ <br /> \dot{c}_2(t)<br /> \end{array} <br /> \right)=\left(<br /> \begin{array}{cc} <br /> 0 &amp; \gamma \mathrm{e}^{\mathrm{i}(\omega - \omega_{21})t}\\ <br /> \gamma \mathrm{e}^{-\mathrm{i}(\omega - \omega_{21})t} &amp; 0 <br /> \end{array} <br /> \right)\left(<br /> \begin{array}{cc} <br /> c_1(t)\\ <br /> c_2(t)<br /> \end{array} <br /> \right)
with initial conditions c_1(0)=1 and c_2(0)=0.
I decided to solve this by diagonalization. The eigenvectors of the Hamiltonian are
\frac{1}{\sqrt{2}}\left(<br /> \begin{array}{cc} <br /> \mathrm{e}^{\mathrm{i}(\omega - \omega_{21})t}\\ <br /> 1<br /> \end{array} <br /> \right) and \frac{1}{\sqrt{2}}\left(<br /> \begin{array}{cc} <br /> \mathrm{e}^{\mathrm{i}(\omega - \omega_{21})t}\\ <br /> -1<br /> \end{array} <br /> \right), with eigenvalues \gamma and -\gamma, respectively.

The change-of-basis matrix then is just
U^\dagger = \frac{1}{\sqrt{2}}\left( <br /> \begin{array}{cc} <br /> \mathrm{e}^{-\mathrm{i}(\omega - \omega_{21})t} &amp; 1\\ <br /> \mathrm{e}^{-\mathrm{i}(\omega - \omega_{21})t} &amp; -1<br /> \end{array} <br /> \right). So in the new basis, the state vector becomes
\frac{1}{\sqrt{2}}\left( <br /> \begin{array}{cc} <br /> \mathrm{e}^{-\mathrm{i}(\omega - \omega_{21})t} &amp; 1\\ <br /> \mathrm{e}^{-\mathrm{i}(\omega - \omega_{21})t} &amp; -1<br /> \end{array} <br /> \right)\left(<br /> \begin{array}{cc} <br /> c_1(t)\\ <br /> c_2(t)<br /> \end{array} <br /> \right):=\left(<br /> \begin{array}{cc} <br /> d_1(t)\\ <br /> d_2(t)<br /> \end{array} <br /> \right).

Therefore, we have
\mathrm{i}\hbar\left(<br /> \begin{array}{cc} <br /> \dot{d}_1(t)\\ <br /> \dot{d}_2(t)<br /> \end{array} <br /> \right)=\left(<br /> \begin{array}{cc} <br /> \gamma &amp; 0\\ <br /> 0 &amp; -\gamma<br /> \end{array} <br /> \right)\left(<br /> \begin{array}{cc} <br /> d_1(t)\\ <br /> d_2(t)<br /> \end{array} <br /> \right)
which can be solved easily. Once the d's are found, the c's are easily computed.

Problem is, this doesn't give the right answer. I suspect it has something to do with using a time-dependent change of basis matrix, but I'm not sure.
 
Physics news on Phys.org
I think your suspicions are right.

From \mathrm{i}\hbar\left( <br /> \begin{array}{cc} <br /> \dot{c}_1(t)\\ <br /> \dot{c}_2(t)<br /> \end{array} <br /> \right)=\left(<br /> \begin{array}{cc} <br /> 0 &amp; \gamma \mathrm{e}^{\mathrm{i}(\omega - \omega_{21})t}\\ <br /> \gamma \mathrm{e}^{-\mathrm{i}(\omega - \omega_{21})t} &amp; 0 <br /> \end{array} <br /> \right)\left(<br /> \begin{array}{cc} <br /> c_1(t)\\ <br /> c_2(t)<br /> \end{array} <br /> \right)

You can write
\mathrm{i}\hbar U^\dagger \left( <br /> \begin{array}{cc} <br /> \dot{c}_1(t)\\ <br /> \dot{c}_2(t)<br /> \end{array} <br /> \right)= U^\dagger \left(<br /> \begin{array}{cc} <br /> 0 &amp; \gamma \mathrm{e}^{\mathrm{i}(\omega - \omega_{21})t}\\ <br /> \gamma \mathrm{e}^{-\mathrm{i}(\omega - \omega_{21})t} &amp; 0 <br /> \end{array} <br /> \right) U U^\dagger \left(<br /> \begin{array}{cc} <br /> c_1(t)\\ <br /> c_2(t)<br /> \end{array} <br /> \right) = \left(<br /> \begin{array}{cc} <br /> \gamma &amp; 0\\ <br /> 0 &amp; -\gamma<br /> \end{array} <br /> \right)\left(<br /> \begin{array}{cc} <br /> d_1(t)\\ <br /> d_2(t)<br /> \end{array} <br /> \right)

But the left hand side cannot be equated to ##\mathrm{i}\hbar\left(
\begin{array}{cc}
\dot{d}_1(t)\\
\dot{d}_2(t)
\end{array}
\right)## due to the time dependence of ##U^\dagger##:

##\left(
\begin{array}{cc}
\dot{d}_1(t)\\
\dot{d}_2(t)
\end{array}
\right) = \frac{d}{dt} \left(U^\dagger \left(
\begin{array}{cc}
c_1(t)\\
c_2(t)
\end{array}
\right)\right) \neq U^\dagger \left(
\begin{array}{cc}
\dot{c}_1(t)\\
\dot{c}_2(t)
\end{array}
\right)##
 
you just need to write the two eqns represented by them,differentiate the second one and put the values from the first one.That's all.
 
TSny said:
But the left hand side cannot be equated to ##\mathrm{i}\hbar\left(
\begin{array}{cc}
\dot{d}_1(t)\\
\dot{d}_2(t)
\end{array}
\right)##

I understand your reasoning, and I can't find a flaw with it, but the conclusion is really bothersome. It's as if the state doesn't obey the Schrodinger equation if its matrix representation uses energy eigenkets as a basis. But how can that be right?

andrien said:
you just need to write the two eqns represented by them,differentiate the second one and put the values from the first one.That's all.
I'll have to do this I guess, but I'd like to see why my approach doesn't work.
 
Anyone?
 
dEdt said:
I understand your reasoning, and I can't find a flaw with it, but the conclusion is really bothersome. It's as if the state doesn't obey the Schrodinger equation if its matrix representation uses energy eigenkets as a basis. But how can that be right?I'll have to do this I guess, but I'd like to see why my approach doesn't work.

TSny showed you why it doesn't work. I'm not sure exactly what you mean when you say "It's as if the state doesn't obey the schrodinger equation if its matrix representation uses energy eigenkets as a basis". In what way does it appear to you to not obey the Schrodinger equation?

Perhaps the issue is that your "change of basis" - your unitary transformation - is really a change of reference frame. You are moving into a rotating frame. There are no true energy eigenkets for this system because the energy is not conserved in time (though it is conserved on average). Perhaps one can say this is similar to the case in classical mechanics: when you move into a rotating frame in classical mechanics you generate "ficticious forces" into Newton's equations of motion. Similarly, your change of reference frame here is generating what look to you like new terms into Schrodinger's equation.

I don't have a quantum book handy right now, but this is (or is pretty similar to) a fairly standard example of a spin in a rotating magnetic field. A book on NMR that treats the problem quantum mechanically should have a discussion, though I don't know if it would exactly address your philosophical problem.
 
Mute said:
Perhaps the issue is that your "change of basis" - your unitary transformation - is really a change of reference frame. You are moving into a rotating frame. There are no true energy eigenkets for this system because the energy is not conserved in time (though it is conserved on average).

Thanks, now it makes sense.
 
Back
Top