I am considering a system of two interacting scalar fields: ψ, and ϕ. The Lagrangian is given by:
L[ψ]=12∂μψ∂μψ+12∂μχ∂μχ−mψ22ψ2−mχ22χ2−σ24χ2ψ2
ψ is initially in a near vacuum state, while χ is initially supposed to be in a thermal state. For simplifying the calculations, the thermal state is taken as a state of definite occupation numbers, following the Bose-Einstein distribution. The interaction between χ and ψ can be assumed to be small, and we can treat χ as an unchanging background for ψ. The exact definition of "small" is not given here, as it will be clear later that it is also part of my question.
I want to calculate, perturbatively, the excitation of ψ as function of time due to its interaction with the thermal bath χ. Since I need to know the dynamics of ψ at an arbitrary time, the final state cannot be assumed to be an equilibrium or free state, so the Gell-Mann-Low theorem does not apply. The problem is I know neither a general method to calculate the effect of χ on ψ, nor an approximating method that still captures that effect. I am satisfied to know the least general method that gives a non-zero excitation, preferably using canonical approach. (I have a limited knowledge of QFT of non-equilibrium systems and finite temperature field theory, so I want to avoid them if possible.)
It would be most ideal for me to know what excitation effects we can capture considering different simplifying assumptions (for example, based on different levels of smallness of the coupling between the fields), and (for each simplifying case) the specific methods we can use to account for these effects. I know a complete answer is very unlikely, so a partial one would be very much appreciated still.
One canonical approach I tried is given below. Is there a different approach, more sophisticated perhaps, but still written in canonical/elementary language, that at least gives a non-zero ψ particle production effect?
I understand that canonical approach must eventually fail. But for my purposes, accuracy is not very important. I am ready to decrease the complexity/generality of the problem, by adding more and more assumptions, perhaps until the canonical approach is adequate. Of course the more generality is preserved, the happier I am. So another way to put it is: What is the most complicated case that we can treat using canonical language?
(One approach that may produce result and can be treated using canonical language is of quasiparticle, but I don't know how to implement it. At some point the concept of particle is no longer valid, but if the coupling is weak enough, then χ's effect on ψ might be sufficiently described by thinking of ψ as a collection of dressed particles - quasiparticles.)
----------------------------------------------------------------------------------------------------------------------------------------
[The failed try] We derive the Euler-Lagrange equation out of the lagrangian above,
(◻+λ2|ϕ0|2)ψ=−σ22χ2ψ
Then, take the Fourier transform
d2ψk(t)dt2+ω2(t)ψk(t)=Jk(t)
Where
ψk(t) and
Jk(t) are the Fourier transforms of
ψ(x) and
−σ22χ2ψ respectively.
and then I define the following "Hamiltonian"
"Hk"=12p2k+12ω2kψ2k−Jk(t)qk
with
qk=ψk, and
pk=˙ψk. We assume that we can still define the annihilation operator as
ak(t)=√ωk2[qk(t)+iωkpk(t)]
The evolution of this operator is given by the following Heisenberg equation
dakdt+iωkak=i√2ωkJk(t)
So
ak(t)=a0ke−iωkt+i√2ωk∫t0dt′eiωk(t′−t)Jk(t′)=a0ke−iωkt+Ak(t)
where
a0k is the evolving unperturbed annihilation operator, and we denoted the second term as
Ak(t).
The plan is take the expectation value of the time varying occupation number operator Nk=a†kak. Since I do not know how to construct a state in the interacting theory, the best I can try is to sandwich the occupation number operator by a state |s⟩ obtained by taking the direct product of a free state of ψ with definite occupation numbers and a free thermal state of χ:
⟨s|a†k(t)ak(t)|s⟩≈⟨s|(a0k)†a0k|s⟩+⟨s|(a0k)†Akeiωkt+A†ka0ke−iωkt|s⟩=αk+2Re[⟨s|(a0k)†Akeiωkt|s⟩]=αk+2Re[i√2ωk∫t0dt′eiωkt′.−σ22∫d3p∫d3q⟨s|(a0k)†χqχp−qψk−p|s⟩]
Using ψk(t)=1√2[akv∗k(t)+a†−kvk(t)] and χk(t)=1√2[bkV∗k(t)+b†−kVk(t)], we have
⟨s|(a0k)†χqχp−qψk−p|s⟩eiωkt=1√8⟨s|(a0k)†[a0k−pv∗k−p+(a0p−k)†vk−p][bqV∗q+b†−qVq][bp−qV∗p−q+b†q−pVp−q]|s⟩
=1√8⟨s|[(a0k)†a0k−pv∗k−p][bqb†q−pVp−qV∗q+b†−qbp−qV∗p−qVq]|s⟩
=1√8[αkv∗kδ(3)(p)][(βq+1)V∗qVqδ(3)(p)+βqVqV∗qδ(3)(p)]
where
αk and
βk are, respectively, the
ψ particle number and
χ particle number number of
k mode. Therefore,
⟨s|(a0k)†Ak|s⟩eiωkt=−iσ28√ωk∫t0dt′eiωkt′[αk√ωke−iωkt′]∫d3q[(2βq+1)1ωχq]=imaginary
where we used
vk=eiωkt′√ωk. Therefore
⟨Nk(t)⟩≈αk
We see that this approach produces a trivial result: no χ particle is produced. This is to be expected, since the ψk operators only act on the ψ part of the Hilbert space, and the χk operators only act on the χ parts. So, there is no mixing among the modes of ψ and χ.