Recently, I study quantum optics and deal with quantization of EM field in a cavity. We know we can express/quantize vector potential in terms of $\hat{a},\hat{a}^{\dagger}$ to get a quantized EM field in a cavity.
$$ \vec{A}(\vec{r},t)=\sum_{n,\sigma}\sqrt{\frac{\hbar}{2\epsilon_0\omega_n V}}\vec{e}_{n,\sigma}\Big[\hat{a}_{n,\sigma}e^{i(\vec{k}_n\cdot\vec{r}-\omega_nt)}+\hat{a}_{n,\sigma}^{\dagger}e^{-i(\vec{k}_n\cdot\vec{r}-\omega_nt)}\Big] $$
The quantized Hamiltonian is:
$$ \hat{H}=\sum_{k}\hbar\omega_k(\hat{n}_k+\frac{1}{2}) $$
The eigenstate of quantized Hamiltonian is: $\left| n_1,n_2,n_3,... \right>=\left|n_1\right>\otimes\left|n_2\right>\otimes\left|n_3\right>...$ The state means there are $n_1$ photons in the first mode and $n_2$ photons in second mode and so on...
So every mode has it own number of photons and photons in the different modes are not at the same frequencies. But why do we take photons as indistinguishable particles?