From the Gross-Pitaevskii equation \begin{equation}i\hbar\frac{\partial\psi}{\partial t}=\left(-\frac{\hbar^2}{2m}\nabla^2+V+g|\psi|^2\right)\psi\end{equation} using the variational relation \begin{equation}i\hbar\frac{\partial\psi}{\partial t}=\frac{\partial\varepsilon}{\partial \psi^*}\end{equation} we find the energy density \begin{equation}\varepsilon=\frac{\hbar^2}{2m}|\nabla\psi|^2+V|\psi|^2+\frac{g}{2}|\psi|^4\end{equation} The energy would be $E=\int d^3r \varepsilon$ and this is a prime integral of the motion, meaning it is a conserved quantity.

My questions are:

1) How do we get the variational relation?

2)How can we prove that $E$ is a conserved quantity?

  • 1
    $\begingroup$ 1) You may read the chapter 7 of Pethick and Smith's book on BEC. 2) $\epsilon$ is energy density which is not conserve, only the total energy $\endgroup$ – unsym Jul 9 '14 at 23:24
  • $\begingroup$ 1) Even if I read that book already I didn't remember that. Thank you. 2) Of course you are right, my mistake; I will amend the question $\endgroup$ – Semola Jul 9 '14 at 23:37
  • $\begingroup$ Do you know a nice way to show the conservation of $ E$ different from brute force computation? $\endgroup$ – Semola Jul 9 '14 at 23:41
  • $\begingroup$ After an integration by parts, $\epsilon$ is almost the operator that gives the time-evolution, and one could perhaps try to copy the proof of the Ehrenfest theorem. $\endgroup$ – Robin Ekman Jul 10 '14 at 0:35
  • $\begingroup$ I think it won't work because the exponential would depend on $|\psi|^2$ as well. $\endgroup$ – Semola Jul 11 '14 at 13:40

The answer to the second question is actually quite straightforward:

by computing $\partial_t E$ and using what given in the Gross-Pitaevskii for $\dot{\psi}$ and $\dot{\psi^*}$ one can check that all the terms cancel out so that $\partial_tE=0$.

To derive the expression for the energy one could also start considering a lagrangian giving the Gross-Pitaevskii when thrown inside the Euler-Lagrange machinery and derive the Hamiltonian in the usual way, so that one can go around the variational relation.

However it would be still interesting to know where that variational relation is coming from.

Does anyone have an idea? Also I am sure there must be a more elegant way to show that the energy is conserved, something less dumb than forced term by term computation.

  • $\begingroup$ Do you think Dirac-Frenkel variational principle (which can be derived as a saddle point approximation for the coherent state Feynman path integral) with the wave functional ansatz $$|\{\psi(\bar{r},t)\}>=e^{\int_{V}d^3_{}\bar{r}\psi(\bar{r},t)\Psi_{}^{\dagger}(\bar{r})}|\ \text{vac}>$$ will lead to the Gross-Pitaevskii equation? $\endgroup$ – Sunyam Oct 6 '17 at 15:16

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.