Let me ask this question again to hopefully get an answer.
Consider a free scalar field $\phi$ on a curved spacetime. The way we define the vacuum is by decomposing the field in terms of mode functions and associate to them creation and annihilation operators $a_k$, $a_k^{\dagger}$ via the inner product defined on the space of solutions of the KG equation (namely, we define $a_k:=\langle f_k, \phi \rangle$, similarly for $a^{\dagger}$) and this definition ensures that the $a_k$, $a_k^{\dagger}$ operators in the field expansion are time-independent (follows from the definition of the product). Then the vacuum is that state $|0\rangle$ s.t. $\forall k$, $a_k|0\rangle=0$.
In a dynamical background (i.e. a non-static spacetime) the Hamiltonian of the system is time dependent, thus we can let $a_k$ and $a_k^{\dagger}$ evolve (non trivially) in the Heisenberg picture, $a_k(t)$.
My question then becomes: Consider an asymptotically flat (in the past and in the future) spacetime with a field as above, see Birrell & Davis for an example, and two sets of mode functions $\{f_k\}$ and $\{g_k\}$ s.t. the first one in the past and the second one in the future, are asymptotically Minkowskian. Call $a_k$ and $b_k$ the relative creation operators.
In light of the Heisenberg evolution of these operators $a_k(t)$,$b_k(t)$, can I interpret $a_k$ as the creation operator of the $\{f_k\}$ mode functions in the infinite past and $b_k$ as the creation operator of the $\{g_k\}$ mode functions in the infinite future? A friend of mine even suggested that $a_k(+\infty)=b_k$ and $b_k(-\infty)=a_k$, is it right?
For context: I was trying to understand operationally the expectation value $${}_{in}\langle 0|N_{out}|0\rangle_{in},$$ where $N_{out}$ is the number operator for the $b_k$ operators and $|0\rangle_{in}$ is the $a_k$ vacuum.
My reasoning is, assuming that the answer to the above question is positive: the inertial observer in the past (i.e. mode functions $\{f_k\}$) prepares the state of the field in $|0\rangle_{in}$ (that is, $a_k(-\infty)|0\rangle_{in}=a_k|0\rangle_{in}=0$ for all $k$). The inertial observer in the future $\{g_k\}$ measures (Heisenberg picture again) ${}_{in}\langle 0|N_{out}(+\infty)|0\rangle_{in}$ particles, being $b_k(+\infty)=b_k$ we can then apply the well known Bogoliubov transformations and conclude the known result $${}_{in}\langle 0|N_{out}(+\infty)|0\rangle_{in}={}_{in}\langle 0|N_{out}|0\rangle_{in}=\Sigma|\beta_k|^2.$$
I hope I've been clear (enough).