In the microcanonical ensemble $(E,V,N)$, the density operator is $$\hat{\rho}=\frac{\delta(\hat{H}-E\,\hat{I})}{Tr(\delta(\hat{H}-E\,\hat{I}))}$$
Where $\hat{H}$ is the Hamiltonian of the system and $\hat{I}$ is the identity operator.
In the canonical ensemble $(T,V,N)$, the density operator is $$\hat{\rho}=\frac{exp\{\hat{H}/(k_B\,T)\}}{Tr(exp\{\hat{H}/(k_B\,T)\})}$$
Where $k_B$ is Boltzmann constant.
In the macrocanonical ensemble $(T,V,\mu)$, the density operator is $$\hat{\rho}=\frac{exp\{(\hat{H}-\mu\,\hat{N})/(k_B\,T)\}}{Tr(exp\{(\hat{H}-\mu\,\hat{N})/(k_B\,T)\})}$$
Where $\hat{N}$ is the particle number operator.
Now, in the isothermal–isobaric ensemble $(T,p,N)$ the expression for the density operator should be similar to the expression used in the canonical ensemble but with the Boltzmann factor multiplied by $exp\{-p\,V/(k_B\,T)\}$.
My question is: How do you write the density operator for an isothermal–isobaric ensemble $(T,p,N)$ in terms of operators? Shall the volume $V$ of the system be written as an operator?