I am currently reading this paper, about an effective field theory approach to Bose-Einstein condensation. I quote
We start by considering a simple non-relativistic many-body problem of spinless bosonic particles at finite temperature. In units of $\hbar=1$ and $k_B = 1$, its Euclidean action functional is given by $$ I[\psi, \psi^\dagger] = \int_0^\beta\text d \tau\int_{\mathbb{R}}\text d^3x \mathcal L $$ with the Lagrangian density $$ \mathcal L = \frac 1 2 [\psi^\dagger(x) \partial_\tau\psi(x)-\psi(x)\partial_\tau\psi^\dagger(x)] - \frac{1}{2m}\nabla\psi^\dagger(x)\cdot\nabla\psi(x) - \lambda[\psi(x)^\dagger\psi(x)]^2 $$ where $\beta = 1/T$ denotes the inverse of temperature. Here, we write $x = (\text x,\tau)$ and $\tau$ is the imaginary "time". All fields fulfill the boundary conditions as of $\psi(\text x, \tau + \beta) = \psi(\text x, \tau)$.
Can somebody explain the thoughts/intuition which lead to this Lagrangian? I do not understand what this "imaginary time" is supposed to be since it is here related to the temperature of the system.
I recognize the similarities to the Schrödinger action with dirac-delta interaction but there is something odd about the time derivatives in this theory which I do not understand.
I am also confused about the boundary condition $\psi(\text x, \tau + \beta) = \psi(\text x, \tau)$ which is probably related to my confusion about $\tau$. Any insight is highly appreciated.
My Knowledge: I am familiar with the mean-field description of a Bose-Einstein condensate using the Hartree method and the resulting Gross-Pitaevskii equations. And I am familiar with the Thomas-Fermi approximation to the Gross-Pitaevskii equations. I have also some knowledge of scalar quantum field theory and the Higgs-Mechanism. But I am really struggling to make sense of the above effective field Lagrangian.