Consider the bosonic Fock space $F$ and let $a^\dagger$ and $a$ denote the "usual" creation operators. As far as I know, these are defined on the dense domain $\mathcal D(N^{1/2}):=\{ \psi \in F| \sum\limits_{n=0}^\infty n\, ||\psi_n||^2 < \infty\}$, where $N$ denotes the number operator, cf. Ref. 1.
Question: What is the domain (such that all expressions are well-defined, they are the adjoint of each others and the CCR hold) of a string of length $K\in \mathbb N$ of creation and annihilation operators: $$a^\#(f_1) \, a^\#(f_2)\,\cdots\, a^\#(f_K)\quad , \tag{1}$$
where $a^\#$ denotes (in each case) either a creation or annihilation operator. I am mostly interested in the case where $K$ is even and the string takes the form
$$ a^\dagger(f_1) \,a^\dagger(f_2)\, \cdots \,a^\dagger(f_{K/2})\,a(f_{K/2+1})\,a(f_{K/2+2})\cdots a(f_{K}) \quad . \tag{2}$$
My guess here is that the domain is $\mathcal D(N^{K/2}):=\{ \psi \in F| \sum\limits_{n=0}^\infty n^{K}\, ||\psi_n||^2< \infty\}$. If so, how to prove it? Are these domains dense?
Closely related to this is the following:
What is the (largest possible) domain such that for $\psi\in F$ it holds that
$$ \langle \psi, a^\#(f_1) \,a^\#(f_2)\,\cdots\, a^\#(f_K)\,\psi\rangle_F\tag{3}$$
is well-defined?
Background:
In Ref. 1 it is stated that , although the creation and annihilation operators are unbounded, domain questions are "unproblematic", since one can take the domain of a sufficiently large power of the number operator.
In Ref. 2 a quantity is defined as $$\langle \psi,a^\dagger(f_1)\, a(f_2)\, \psi\rangle_F $$ for $\psi \in F$ with $\langle \psi, N\psi\rangle_F <\infty$. For me this condition means that we need that $\psi \in \mathcal D(N)$ such that $N\psi \in F$ and further that $\psi\in \mathcal D(N^{1/2})$ (which I think automatically follows). However, following this math stack exchange post, the form domain of an operator (which I think here is meant with the finite particle expectation value) can be larger than the operator domain. So the question remains: What is the (largest possible) domain $\mathcal D$ such that the above expression makes sense for $\psi \in \mathcal D$.
Maybe in our particular case it follows that the form domain is indeed $\mathcal D(N^{1/2})$ - or more generally for strings of length $K$ that the form domain is $\mathcal D(N^{K/2})$?
References:
Benedikter, Niels, Marcello Porta, and Benjamin Schlein. Effective evolution equations from quantum dynamics. Cham: Springer International Publishing, 2016. Chapter 3.
Solovej, J. P. Many body quantum mechanics. Lecture Notes 2007. chapter 8.