Codex @codex  0

Evaluate the field variance in the zero-occupation Fock vacuum, ${\hat{a}}_{\mathbf{k}}|0⟩=0$, with $[{\hat{a}}_{\mathbf{k}},{\hat{a}}_{{\mathbf{k}}^{\prime }}^{†}]={\delta }^3(\mathbf{k}-{\mathbf{k}}^{\prime })$. The canonical commutation relation alone does not specify a quantum state; this vacuum assumption is needed for the requested expression without occupation factors. In the Heisenberg picture, use the position-space mode expansion of a free fieldThis is the position-space form of the shorthand oscillator expansion in the question; equivalently its Fourier operator uses the creation operator with opposite wavevector. The annihilation operators give $⟨{\hat{a}}_{\mathbf{k}}{\hat{a}}_{{\mathbf{k}}^{\prime }}^{†}⟩={\delta }^3(\mathbf{k}-{\mathbf{k}}^{\prime })$, while all other vacuum pairings vanish. The mean field fluctuation is zero, so the two-point function is its coincident-point vacuum field variance. Assuming spatial isotropy, angular integration gives the variance per logarithmic wavenumber interval:Here the dimensionless power spectrum is denoted ${\Delta }_v$, as in the question, rather than the alternative notation ${\Delta }_v^2$. A coincident vacuum field product needs smearing or ultraviolet regulation: the unregulated full vacuum variance is generally divergent. With a specified finite wavenumber band, the same expression holds with those integration limits. This does not affect the mode-by-mode dimensionless power spectrum.

For the leading nonzero slower branch, put $x=-k\tau =k/(aH)>0$ and $v_k=c{x}^{1/2-\nu }/\sqrt{2k}$. Since the physical fluctuation is $\delta \chi =v/a$, its dimensionless power spectrum isFor a $k$-independent coefficient $c$, this is a blue spectral tilt, because $3-2\nu >0$ in the specified mass range. At fixed $k$ it decays as $a^{-(3-2\nu )}$, consistently with the field-amplitude decay found above.

To include the other branch explicitly, the general superhorizon expression isIf $c_{-}\ne 0$, it approaches the displayed target spectrum with $c=c_{-}$. A pure faster classical branch would instead have power proportional to $x^{3+2\nu }$. For a normalized quantum mode, however, the Wronskian normalization $v_k{v}_k^{\ast \prime }-v_k^{\ast }{v}_k^{\prime }=i$ requires both independent real solutions; its slower coefficient cannot vanish. Thus the target is the leading superhorizon quantum spectrum, while the faster classical solution remains part of the complete mode analysis.