Unitary transformation (quantum mechanics)

Recall that by the definition of a unitary matrix, ${\displaystyle U^{\dagger }U=1}$. Beginning with the Schrödinger equation,

${\displaystyle {\dot {\psi }}=-{\frac {i}{\hbar }}H\psi }$,

we can therefore insert ${\displaystyle U^{\dagger }U}$ at will. In particular, inserting it after ${\displaystyle H/\hbar }$ and also premultiplying both sides by ${\displaystyle U}$, we get

${\displaystyle U{\dot {\psi }}=-{\frac {i}{\hbar }}\left(UHU^{\dagger }\right)U\psi \quad \quad (1)}$.

Next, note that by the product rule,

${\displaystyle {\frac {\mathrm {d} }{\mathrm {d} t}}\left(U\psi \right)={\dot {U}}\psi +U{\dot {\psi }}}$.

Inserting another ${\displaystyle U^{\dagger }U}$ and rearranging, we get

${\displaystyle U{\dot {\psi }}={\frac {\mathrm {d} }{\mathrm {d} t}}{\Big (}U\psi {\Big )}-{\dot {U}}U^{\dagger }U\psi \quad \quad (2)}$.

Finally, combining (1) and (2) above results in the desired transformation:

${\displaystyle {\frac {\mathrm {d} }{\mathrm {d} t}}{\Big (}U\psi {\Big )}=-{\frac {i}{\hbar }}{\Big (}UH{U^{\dagger }}+i\hbar \,{\dot {U}}{U^{\dagger }}{\Big )}{\Big (}U\psi {\Big )}\quad \quad \left(3\right)}$.

If we adopt the notation ${\displaystyle {\breve {\psi }}:=U\psi }$ to describe the transformed wave function, the equations can be written in a clearer form. For instance, ${\displaystyle (3)}$ can be rewritten as

${\displaystyle {\frac {\mathrm {d} }{\mathrm {d} t}}{\breve {\psi }}=-{\frac {i}{\hbar }}{\breve {H}}{\breve {\psi }}\quad \quad \left(4\right)}$,

which can be rewritten in the form of the original Schrödinger equation,

${\displaystyle {\breve {H}}{\breve {\psi }}=i\hbar {\operatorname {d} \!{\breve {\psi }} \over \operatorname {d} \!t}.}$

The original wave function can be recovered as ${\displaystyle \psi =U^{\dagger }{\breve {\psi }}}$.

Unitary transformations can be seen as a generalization of the interaction (Dirac) picture. In the latter approach, a Hamiltonian is broken into a time-independent part and a time-dependent part,

${\displaystyle H(t)=H_{0}+V(t)\quad \quad (a)}$.

In this case, the Schrödinger equation becomes

${\displaystyle {\dot {\psi _{I}}}=-{\frac {i}{\hbar }}\left(e^{iH_{0}t/\hbar }Ve^{-iH_{0}t/\hbar }\right)\psi _{I}}$, with ${\displaystyle \psi _{I}=e^{iH_{0}t/\hbar }\psi }$.[4]

The correspondence to a unitary transformation can be shown by choosing ${\textstyle U(t)=\exp \left[{+iH_{0}t/\hbar }\right]}$. As a result, ${\displaystyle {U^{\dagger }}(t)=\exp \left[{-iH_{0}t}/\hbar \right].}$

Using the notation from ${\displaystyle (0)}$ above, our transformed Hamiltonian becomes

${\displaystyle {\breve {H}}=U\left[H_{0}+V(t)\right]U^{\dagger }+i\hbar {\dot {U}}U^{\dagger }\quad \quad (b)}$

First note that since ${\displaystyle U}$ is a function of ${\displaystyle H_{0}}$, the two must commute. Then

${\displaystyle UH_{0}U^{\dagger }=H_{0}}$,

which takes care of the first term in the transformation in ${\displaystyle (b)}$, i.e. ${\displaystyle {\breve {H}}=H_{0}+UV(t)U^{\dagger }+i\hbar {\dot {U}}U^{\dagger }}$. Next use the chain rule to calculate

${\displaystyle {\begin{aligned}i\hbar {\dot {U}}U^{\dagger }&=i\hbar \left({\operatorname {d} \!U \over \operatorname {d} \!t}\right)e^{-iH_{0}t/\hbar }\\&=i\hbar {\Big (}iH_{0}/\hbar {\Big )}e^{+iH_{0}t/\hbar }e^{-iH_{0}t/\hbar }\\&=i\hbar \left({iH_{0}}/\hbar \right)\\&=-H_{0},\\\end{aligned}}}$

which cancels with the other ${\displaystyle H_{0}}$. Evidently we are left with ${\displaystyle {\breve {H}}=UVU^{\dagger }}$, yielding ${\displaystyle {\dot {\psi _{I}}}=-{\frac {i}{\hbar }}UVU^{\dagger }\psi _{I}}$ as shown above.

When applying a general unitary transformation, however, it is not necessary that ${\displaystyle H(t)}$ be broken into parts, or even that ${\displaystyle U(t)}$ be a function of any part of the Hamiltonian.

Consider an atom with two states, ground ${\displaystyle |g\rangle }$ and excited ${\displaystyle |e\rangle }$. The atom has a Hamiltonian ${\displaystyle H=\hbar \omega {|{e}\rangle \langle {e}|}}$, where ${\displaystyle \omega }$ is the frequency of light associated with the g-e transition. Now suppose we illuminate the atom with a drive at frequency ${\displaystyle \omega _{d}}$ which couples the two states, and that the time-dependent driven Hamiltonian is

${\displaystyle H/\hbar =\omega |e\rangle \langle e|+\Omega \ e^{i\omega _{d}t}|g\rangle \langle e|+\Omega ^{*}\ e^{-i\omega _{d}t}|e\rangle \langle g|}$

for some complex drive strength ${\displaystyle \Omega }$. Because of the competing frequency scales (${\displaystyle \omega }$, ${\displaystyle \omega _{d}}$, and ${\displaystyle \Omega }$), it is difficult to anticipate the effect of the drive (see driven harmonic motion).

Without a drive, the phase of ${\displaystyle |e\rangle }$ would oscillate relative to ${\displaystyle |g\rangle }$. In the Bloch sphere representation of a two-state system, this corresponds to rotation around the z-axis. Conceptually, we can remove this component of the dynamics by entering a rotating frame of reference defined by the unitary transformation ${\displaystyle U=e^{i\omega t|e\rangle \langle e|}}$. Under this transformation, the Hamiltonian becomes

${\displaystyle H/\hbar \to \Omega \,e^{i(\omega _{d}-\omega )t}|g\rangle \langle e|+\Omega ^{*}\,e^{i(\omega -\omega _{d})t}|e\rangle \langle g|}$.

If the driving frequency is equal to the g-e transition's frequency, ${\displaystyle \omega _{d}=\omega }$, resonance will occur and then the equation above reduces to

${\displaystyle {\breve {H}}/\hbar =\Omega \ |g\rangle \langle e|+\Omega ^{*}\ |e\rangle \langle g|}$.

Without getting into details, we can already predict that the dynamics will involve an oscillation between the ground and excited states at frequency ${\displaystyle \Omega }$.[4]

As another limiting case, suppose the drive is far off-resonant, ${\displaystyle |\omega _{d}-\omega |\gg 0}$. We can figure out the dynamics in that case without solving the Schrödinger equation directly. Suppose the system starts in the ground state ${\displaystyle |g\rangle }$. Initially, the Hamiltonian will populate some component of ${\displaystyle |e\rangle }$. A small time later, however, it will populate roughly the same amount of ${\displaystyle |e\rangle }$ but with completely different phase. Thus the effect of an off-resonant drive will tend to cancel itself out. This can also be expressed by saying that an off-resonant drive is rapidly rotating in the frame of the atom.

These concepts are illustrated in the table below, where the sphere represents the Bloch sphere, the arrow represents the state of the atom, and the hand represents the drive.

	Lab frame	Rotating frame
Resonant drive	Resonant drive in the lab frame	Resonant drive in a frame rotating with the atom
Off-resonant drive	Off-resonant drive in the lab frame	Off-resonant drive in a frame rotating with the atom

The example above could also have been analyzed in the interaction picture. The following example, however, is more difficult to analyze without the general formulation of unitary transformations. Consider two harmonic oscillators, between which we would like to engineer a beam splitter interaction,

${\displaystyle g\,ab^{\dagger }+g^{*}\,a^{\dagger }b}$.

This was achieved experimentally with two microwave cavity resonators serving as ${\displaystyle a}$ and ${\displaystyle b}$.[5] Below, we sketch the analysis of a simplified version of this experiment.

In addition to the microwave cavities, the experiment also involved a transmon qubit, ${\displaystyle c}$, coupled to both modes. The qubit is driven simultaneously at two frequencies, ${\displaystyle \omega _{1}}$ and ${\displaystyle \omega _{2}}$, for which ${\displaystyle \omega _{1}-\omega _{2}=\omega _{a}-\omega _{b}}$.

${\displaystyle H_{\mathrm {drive} }/\hbar =\Re \left[\epsilon _{1}e^{i\omega _{1}t}+\epsilon _{2}e^{i\omega _{2}t}\right](c+c^{\dagger }).}$

In addition, there are many fourth-order terms coupling the modes, but most of them can be neglected. In this experiment, two such terms which will become important are

${\displaystyle H_{4}/\hbar =g_{4}{\Big (}e^{i(\omega _{b}-\omega _{a})t}ab^{\dagger }+{\text{h.c.}}{\Big )}c^{\dagger }c}$.

(H.c. is shorthand for the Hermitian conjugate.) We can apply a displacement transformation, ${\displaystyle U=D(-\xi _{1}e^{-i\omega _{1}t}-\xi _{2}e^{-i\omega _{2}t})}$, to mode ${\displaystyle c}$. For {{carefully chosen amplitudes, this transformation will cancel ${\displaystyle H_{\textrm {drive}}}$ while also displacing the ladder operator, ${\displaystyle c\to c+\xi _{1}e^{-i\omega _{1}t}+\xi _{2}e^{-i\omega _{2}t}}$. This leaves us with

${\displaystyle H/\hbar =g_{4}{\Big (}e^{i(\omega _{b}-\omega _{a})t}ab^{\dagger }+e^{i(\omega _{a}-\omega _{b})t}a^{\dagger }b{\big )}(c^{\dagger }+\xi _{1}^{*}e^{i\omega _{1}t}+\xi _{2}^{*}e^{i\omega _{2}t})(c+\xi _{1}e^{-i\omega _{1}t}+\xi _{2}e^{-i\omega _{2}t})}$.

Expanding this expression and dropping the rapidly rotating terms, we are left with the desired Hamiltonian,

${\displaystyle H/\hbar =g_{4}\xi _{1}^{*}\xi _{2}e^{i(\omega _{b}-\omega _{a}+\omega _{1}-\omega _{2})t}\ ab^{\dagger }+{\text{h.c.}}=g\,ab^{\dagger }+g^{*}\,a^{\dagger }b}$.