Introduction

Exercise 12.5 from Modern Condensed Matter Physics (Girvin and Yang, 2019) asks to construct a Gaussian wave packet in the lowest Landau level in the Landau gauge, such that it is localized as closely as possible around some point R(Rx,Ry)\mbf R\ceq\p{R_x,R_y}.

Actually, we can prove that the smallest wave packet is a Gaussian wave packet. Here is the derivation.

The problem

First, for readers who are not familiar with the Landau levels, here is a brief introduction. For an electron confined in the xyxy plane under a magnetic field B=Bz^\mbf B=B\bhat z, its Hamiltonian is H=12me(px2+(pyeBcx)2)H=\fr1{2m_e}\p{p_x^2+\p{p_y-\fr{eB}cx}^2} under the Landau gauge A=Bxy^\mbf A=Bx\bhat y. Its eigenstates in the position representation are ψnk ⁣(x,y)=eikyHn ⁣(xlkl)e(xkl2)2/2l2\fc{\psi_{nk}}{x,y}=\e^{\i ky}\fc{H_n}{\fr xl-kl} \e^{-\p{x-kl^2}^2/2l^2} labeled by nNn\in\bN and kRk\in\bR, where HnH_n is the Hermite polynomial of degree nn and lc/eBl\ceq\sqrt{\hbar c/eB}. States with the same nn are degenerate in energy (En=(n+1/2)eB/mecE_n=\p{n+1/2}\hbar eB/m_ec) and make up the nnth Landau level. The Landau level with n=0n=0 is called the lowest Landau level.

The problem, now, is this optimization problem: minakΨ|x2+y2|ΨstΨ|Ψ=1,Ψ|x|Ψ=Rx,Ψ|y|Ψ=Ry\begin{align*} \min_{a_k}\quad&\mel{\Psi}{x^2+y^2}{\Psi}\\ \st\quad&\braket{\Psi}{\Psi}=1,\\ &\mel{\Psi}{x}{\Psi}=R_x,\\ &\mel{\Psi}{y}{\Psi}=R_y \end{align*} (optimizing x2+y2\a{x^2+y^2} is equivalent to optimizing σx2+σy2\sgm_x^2+\sgm_y^2 because x\a x and y\a y are both fixed), where Ψ\ket\Psi is defined as the state whose position representation is Ψ ⁣(x,y)=dkakeikye(xkl2)2/2l2.\fc\Psi{x,y}=\int\d k\,a_k\e^{\i ky}\e^{-\p{x-kl^2}^2/2l^2}.

The solution

Consider the moment-generating function M ⁣(u,v)Ψ|eux+vy|Ψ=dxdyeux+vydkakeikye12l2(xkl2)2dkakeikye12l2(xkl2)2=dkdkakakdxeux12l2(xkl2)212l2(xkl2)2dyexp ⁣(vy+i(kk)y)2πδ ⁣(kkiv)=2πdkakak+ivdxexp ⁣(ux12l2(xkl2)212l2(x(k+iv)l2)2)lπexp ⁣(14l2(4ku+u2+2iuv+v2))=2π3/2lexp ⁣(14l2(u2+2iuv+v2))dkakak+ivekl2u=2π3/2ldkak(ak+kl2aku+iakv+14l2(1+2k2l2)aku2+14(l2ak2ak)v2+i2l2(ak+2kak)uv+),\begin{align*} \fc M{u,v}&\ceq\mel{\Psi}{\e^{ux+vy}}{\Psi}\\ &=\iint\d x\d y\,\e^{ux+vy} \int\d k\,a_k^*\e^{-\i ky}\e^{-\fr1{2l^2}\p{x-kl^2}^2} \int\d k'\,a_{k'}\e^{\i k'y}\e^{-\fr1{2l^2}\p{x-k'l^2}^2}\\ &=\iint\d k\d k'\,a_k^*a_{k'}\int\d x\,\e^{ ux-\fr1{2l^2}\p{x-kl^2}^2-\fr1{2l^2}\p{x-k'l^2}^2 }\underbrace{\int\d y\,\fc\exp{vy+\i\p{k'-k}y}}_{2\pi\fc\dlt{k'-k-\i v}}\\ &=2\pi\int\d k\,a_k^*a_{k+\i v}\underbrace{\int\d x\,\fc\exp{ ux-\fr1{2l^2}\p{x-kl^2}^2-\fr1{2l^2}\p{x-\p{k+\i v}l^2}^2 }}_{l\sqrt\pi\fc\exp{\fr14l^2\p{4ku+u^2+2\i uv+v^2}}}\\ &=2\pi^{3/2}l\fc\exp{\fr14l^2\p{u^2+2\i uv+v^2}} \int\d k\,a_k^*a_{k+\i v}\e^{kl^2u}\\ &=2\pi^{3/2}l\int\d k\,a_k^*\left( a_k+kl^2a_ku+\i a_k'v+\fr14l^2\p{1+2k^2l^2}a_ku^2 \right.\\&\qquad\qquad\qquad\qquad\left. {}+\fr14\p{l^2a_k-2a_k''}v^2 +\fr\i2l^2\p{a_k+2ka_k'}uv+\cdots \right), \end{align*} where akdak/dka_k'\ceq\d a_k/\d k and akd2ak/dk2a_k''\ceq\d^2a_k/\d k^2. On the other hand, we have M ⁣(u,v)=Ψ|1+ux+uy+12u2x2+12v2y2+uvxy+|Ψ.\fc M{u,v}=\mel{\Psi}{1+ux+uy+\fr12u^2x^2+\fr12v^2y^2+uvxy+\cdots}{\Psi}. Compare the expansion coefficients, and we have Ψ|Ψ=2π3/2ldkakak,Ψ|x|Ψ=2π3/2l3dkakkak,Ψ|y|Ψ=2iπ3/2ldkakak,Ψ|x2|Ψ=12π3/2l3dkak(1+2k2l2)ak,Ψ|y2|Ψ=12π3/2ldkak(l22ak)ak.\begin{align*} \braket{\Psi}{\Psi}&=2\pi^{3/2}l\int\d k\,a_k^*a_k,\\ \mel{\Psi}{x}{\Psi}&=2\pi^{3/2}l^3\int\d k\,a_k^*ka_k,\\ \mel{\Psi}{y}{\Psi}&=2\i\pi^{3/2}l\int\d k\,a_k^*a_k',\\ \mel{\Psi}{x^2}{\Psi}&=\fr12\pi^{3/2}l^3\int\d k\,a_k^*\p{1+2k^2l^2}a_k,\\ \mel{\Psi}{y^2}{\Psi}&=\fr12\pi^{3/2}l\int\d k\,a_k^*\p{l^2-2a_k''}a_k. \end{align*}

Define φ ⁣(k)ak2π3/2l\fc\vphi k\ceq a_k\sqrt{2\pi^{3/2}l}. Define fictitious position and momentum operators acting on φ\vphi as Ξφ:kkφ ⁣(k),Πφ:kiφ ⁣(k).\Xi\vphi:k\mapsto k\fc\vphi k,\quad \Pi\vphi:k\mapsto-\i\fc{\vphi'}k. Using the constraints of the original optimization problem and abusing the bra–ket notation on φ\vphi, we have φ|φ=1,φ|Ξ|φ=Rxl2,φ|Π|φ=Ry.\braket{\vphi}{\vphi}=1,\quad\mel\vphi\Xi\vphi=\fr{R_x}{l^2},\quad \mel\vphi\Pi\vphi=-R_y. The objective function then becomes Ψ|x2+y2|Ψ=12l2+φ|H|φ,\mel{\Psi}{x^2+y^2}{\Psi}=\fr12l^2+\mel{\vphi}{\mcal H}{\vphi}, where HΠ2/2+l4Ξ2/2\mcal H\ceq \Pi^2/2+l^4\Xi^2/2 is a fictitious Hamiltonian, which is the Hamiltonian of a harmonic oscillator with mass 11 and angular frequency ωl2\omg\ceq l^2.

The optimization problem can now be re-stated in terms of φ\ket\vphi as minφφ|H|φstφ|φ=1,φ|Ξ|φ=Rx/ω,φ|Π|φ=Ry.\begin{align*} \min_{\ket\vphi}\quad&\mel\vphi{\mcal H}{\vphi}\\ \st\quad&\braket\vphi\vphi=1,\quad\mel\vphi\Xi\vphi=R_x/\omg,\quad\mel\vphi\Pi\vphi=-R_y. \end{align*} Physically, this means that we want to find the state of a harmonic oscillator with the given expectation values of position and momentum and the lowest energy. To find it, we can use Hisenberg’s uncertainty principle: H=12Π2+12ω2Ξ2=12(Π2+σΠ2)+12ω2(Ξ2+σΞ2)=12σΠ2+12ω2σΞ2+12Ry2+12Rx2ωσΠσΞ+12R212ω+12R2.\begin{align*} \a{\mcal H}&=\fr12\a{\Pi^2}+\fr12\omg^2\a{\Xi^2}\\ &=\fr12\p{\a\Pi^2+\sgm_\Pi^2}+\fr12\omg^2\p{\a{\Xi^2}+\sgm_\Xi^2}\\ &=\fr12\sgm_\Pi^2+\fr12\omg^2\sgm_\Xi^2+\fr12R_y^2+\fr12 R_x^2\\ &\ge\omg\sgm_\Pi\sgm_\Xi+\fr12R^2 \ge\fr12\omg+\fr12R^2. \end{align*} The equality in the first “\ge” is achieved when σΠ=ωσΞ\sgm_\Pi=\omg\sgm_\Xi, and that in the second “\ge” is achieved when the uncertainty principle is saturated. As we know from quantum mechanics, the coherent state of a harmonic oscillator satisfies both conditions. The wavefunction of this state is φ ⁣(k)=(ωπ)1/4exp ⁣(12ω(kRxω)2iRyk).\fc\vphi k=\p{\fr\omg\pi}^{1/4} \fc\exp{-\fr12\omg\p{k-\fr{R_x}{\omg}}^2-\i R_yk}. Express the final result in terms of aka_k: ak=12πeikRye12l2(Rxkl2)2.a_k=\fr1{\sqrt2\pi}\e^{-\i kR_y}\e^{-\fr1{2l^2}\p{R_x-kl^2}^2}. We may work out the integral to get the wave function of the wave packet: Ψ ⁣(x,y)=12πlexp ⁣(14l2((xRx)2+(yRy)22i(x+Rx)(yRy))).\fc{\Psi}{x,y}=\fr1{\sqrt{2\pi}l}\fc\exp{-\fr1{4l^2}\p{ \p{x-R_x}^2+\p{y-R_y}^2-2\i\p{x+R_x}\p{y-R_y} }}. This is a Gaussian wave packet centered at R\mbf R with covariance matrix Diag ⁣(l2,l2)\opc{Diag}{l^2,l^2}.

Further problems

The optimal wave packet is indeed Gaussian. This makes me curious about whether this is a coincidence or not.

Another thing worth noting is that this result is actually the Dirac delta wave function peaking at R\mbf R projected into the lowest Landau level. This was actually my first idea to solve the problem. I was like: well, isn’t the Dirac delta the smallest possible wave packet by all means? If the basis is complete, I can surely combine them into a Dirac delta, and it would be very easy to work out aka_k in this case. Then, I was like: nah, merely a single Landau level is not complete, so I cannot do that anyway. I then did not even bother to proceed with this approach and went on to trying other methods. It turns out that this approach is actually correct—at least it gives the same result as the correct approach.