Imagine an area of space far from a black body with temperature $T$ that has emitted perfect black body radiation. So that an area of space has random photons.
Can we construct a "useful" wavefunction that could represent this? To get an exact wavefunction would presumable require us to know the exact entanglement between each pair of photons. But can we make an approximation? Let's ignore the polarisation so we approximate the electromagnetic field with a scalar field , $\phi$
My guess would be something like:
$$|\psi\rangle= \exp\left(\int\int\phi(x) e^{ix.k}\left( \frac{ k^3}{e^{k/T} - 1} \right) dx^3 dk^3\right) |0\rangle $$
using Plank's formula for black body radiation. And I'm using the exponential as a way to sum over 1, 2, 3,.. photons.
Is this anywhere close?