Vamos a empezar con el ansatz (voy a asumir su mayoría de más de la métrica de la firma)
\begin{equation}
\hat\phi(x) = \int \frac{\mathrm{d}^3 \mathbf{k}}{(2\pi)^{3/2}}\left(\hat A_\mathbf{k} e^{i k\cdot x} + \hat B_\mathbf{k}e^{-ik\cdot x}\right)
\end{equation}
donde $\hat A_\mathbf{k}$ $\hat B_\mathbf{k}$ son algunos arbitraria operadores y $k^0 = E_k$. Primero nos tenga en cuenta que $\phi$ es un verdadero campo que significa que debemos tener $B_\mathbf{k} = A_\mathbf{k}^\dagger$. También un campo cuántico debe obedecer a la igualdad de tiempo de conmutación relación
\begin{equation}
\left[ \hat \phi(\mathbf{x},t),\hat{\dot\phi} (\mathbf{y},t)\right] = i\delta(\mathbf{x-y}).
\end{equation}
Conectar nuestro ansatz en esta relación podemos obtener la condición de
\begin{equation}
2E_k \left[\hat A_\mathbf{k}, \hat A_\mathbf{k'}^\dagger \right] = \delta(\mathbf{k-k'})
\end{equation}
Esto es sólo la escalera operador de conmutación relación con un extra de factor de normalización. Para mayor comodidad podemos definir un reescalado operador $\hat a_\mathbf{k} \equiv \sqrt{2E_k}\hat A_\mathbf{k}$ a que el convencional de conmutación relación.
Cálculo detallado
Tomar el tiempo derivado de la esfera de obtener
\begin{equation}
\hat{\dot \phi}(x) = \int \frac{\mathrm{d}^3 \mathbf{k}}{(2\pi)^{3/2}}\left(-iE_k\hat A_\mathbf{k} e^{-i E_kt + i\mathbf{k\cdot x}} + iE_k\hat A_\mathbf{k}^\dagger e^{i E_kt - i\mathbf{k\cdot x}}\right).
\end{equation}
A continuación, el colector entre el campo y su tiempo derivativo (más generalmente, su conjugado) es
\begin{eqnarray}
\left[ \hat \phi(\mathbf{x},t),\hat{\dot\phi} (\mathbf{y},t)\right] = \int \frac{\mathrm{d}^3 \mathbf{k}\mathrm{d}^3 \mathbf{k'}}{(2\pi)^{3}}\left\{iE_{k'}\left[\hat A_\mathbf{k},\hat A_\mathbf{k'}^\dagger\right] e^{-i(E_k-E_{k'})t+i(\mathbf{k\cdot x - k'\cdot y)}} + iE_k\left[\hat A_\mathbf{k'},\hat A_\mathbf{k}^\dagger\right] e^{i(E_k-E_{k'})t-i(\mathbf{k\cdot x - k'\cdot y)}}\right\} = \int \frac{\mathrm{d}^3 \mathbf{k}\mathrm{d}^3 \mathbf{k'}}{(2\pi)^{3}}\left\{iE_{k'}\left[\hat A_\mathbf{k},\hat A_\mathbf{k'}^\dagger\right] e^{-i(E_k-E_{k'})t} + iE_k\left[\hat A_\mathbf{-k'},\hat A_\mathbf{-k}^\dagger\right] e^{i(E_k-E_{k'})t}\right\}e^{i(\mathbf{k\cdot x - k'\cdot y)}}
\end{eqnarray}
donde hemos utilizado el hecho de que los operadores conmutan con ellos mismos y cambiado las variables $\mathbf{k\rightarrow -k}$, $\mathbf{k'\rightarrow -k'}$ en el segundo término. Esto debe ser igual a la función delta de
\begin{equation}
i\delta(\mathbf{x-y}) = i \int \frac{\mathrm{d}^3 \mathbf{k}}{(2\pi)^3}e^{i\mathbf{k\cdot(x-y)}} = i \int \frac{\mathrm{d}^3 \mathbf{k}\mathrm{d}^3 \mathbf{k'}}{(2\pi)^3}\delta(\mathbf{k-k'})e^{i(\mathbf{k\cdot x - k'\cdot y)}}
\end{equation}
lo que significa que debemos tener
\begin{equation}
\left[\hat A_\mathbf{k},\hat A_\mathbf{k'}^\dagger\right] + \left[\hat A_\mathbf{-k'},\hat A_\mathbf{-k}^\dagger\right] = \frac{\delta(\mathbf{k-k'})}{E_k}
\end{equation}
a que $\left[\hat A_\mathbf{k},\hat A_\mathbf{k'}^\dagger\right]=\delta(\mathbf{k-k'})/2E_k$ es la solución.