EM wave function & photon wavefunction According to this review

Photon wave function. Iwo Bialynicki-Birula. Progress in Optics 36 V (1996), pp. 245-294. arXiv:quant-ph/0508202,

a classical EM plane wavefunction is a wavefunction (in Hilbert space) of a single photon with definite momentum (c.f section 1.4), although a naive probabilistic interpretation is not applicable. However, what I've learned in some other sources (e.g. Sakurai's Advanced QM, chap. 2) is that, the classical EM field is obtained by taking the expectation value of the field operator. Then according to Sakurai, the classical $E$ or $B$ field of a single photon state with definite momentum p is given by
$\langle p|\hat{E}(or \hat{B})|p\rangle$, which is $0$ in the whole space. This seems to contradict the first view, but both views make equally good sense to me by their own reasonings, so how do I reconcile them?
 A: The expectation values
$$ \langle p | \vec E(\vec x) | p\rangle $$
and similarly for $\vec B(\vec x)$ vanish for a simple reason: the state $|p\rangle$ is by definition translational symmetric (translation only changes the phase of the state, the overall normalization) so the expectation values of any field in this state has to be translationally symmetric, too (the phase cancels between the ket and the bra). 
So if you expect to see classical waves in expectation values in such momentum eigenstates, you are unsurprisingly disappointed. Incidentally, the same thing holds for any other field including the Dirac field (in contrast with the OP's assertion). If you compute the expectation value of the Dirac field $\Psi(\vec x)$ in a one-particle momentum eigenstate with one electron, this expectation value also vanishes. In this Dirac case, it's much easier to prove so because the expectation values of all fermionic operators (to the first or another odd power) vanish because of the Grassmann grading.
The vanishing of the expectation values of fields (those that can have both signs, namely the linear functions of the "basic" fields connected with the given particle) would be true for any momentum eigenstates, even multiparticle states which are momentum eigenstates simply because the argument above holds universally. You may think that this vanishing is because the one-particle momentum eigenstate is some mixture of infinitesimal electromagnetic waves that are allowed to be in any "phase" and these phases therefore cancel.
However, the formal relationship between the classical fields and the one-particle states still holds if one is more careful. In particular, one may construct "coherent states" which are multiparticle states with an uncertain number of particles which are the closest approximations of a classical configuration. You may think of coherent states as the ground states of a harmonic oscillator (and a quantum field is an infinite-dimensional harmonic oscillator) which are shifted in the position directions and/or momentum directions, i.e. states
$$ |a\rangle = C_\alpha \cdot \exp(\alpha\cdot a^\dagger) |0\rangle  $$
This expression may be Taylor-expanded to see the components with individual numbers of excitations, $N=0,1,2,3,\dots$ The $C_\alpha$ coefficient is just a normalization factor that doesn't affect physics of a single coherent state.
With a good choice of $\alpha$ for each value of the classical field (there are many independent $a^\dagger(k,\lambda)$ operators for a quantum field and each of them has its $\alpha(k,\lambda)$), such a coherent state may be constructed for any classical configuration. The expectation values of the classical fields $\vec B,\vec E$ in these coherent states will be what you want.
Now, with the coherent state toolkit, you may get a more detailed understanding of why the momentum eigenstates which are also eigenstates of the number of particles have vanishing eigenvalues. The coherent state is something like the wave function
$$ \exp(-(x-x_S)^2/2) $$
which is the Gaussian shifted to $x_S$ so $x_S$ is the expectation value of $x$ in it. Such a coherent state may be obtained by an exponential operator acting on the vacuum. The initial term in the Taylor-expansion is the vacuum itself; the next term is a one-particle state that knows about the structure of the coherent state – because the remaining terms in the Taylor expansions are just gotten from the same linear piece that acts many times, recall the $Y^k/k!$ form of the terms in the Taylor expansion of $\exp(Y)$: here, $Y$ is the only thing you need to know.
On the other hand, the expectation value of $x$ in the one-particle state is of course zero. It's because the wave function of a one-particle state is an odd function such as
$$ x\cdot \exp(-x^2/2) $$
whose probability density is symmetric (even) in $x$ so of course that the expectation value has to be zero. If you look at the structure of the coherent state and you imagine that the $\alpha$ coefficients are very small so that multiparticle states may be neglected for the sake of simplicity, you will realize that the nonzero expectation value of $x$ in the shifted state (the coherent state) boils down to some interference between the vacuum state and the one-particle state; it is not a property of the one-particle state itself! More generally, the nonzero expectation values of fields at particular points of the spacetime prove some interference between components of the state that have different numbers of the particle excitations in them.
The latter statement should be unsurprising from another viewpoint. If you consider something like the matrix element
$$ \langle n | a^\dagger | m \rangle $$
where the bra and ket vectors are eigenstates of a harmonic oscillator with some number of excitations, it's clear that it's nonzero only if $m=n\pm 1$. In particular, $m$ and $n$ cannot be equal. If you consider the expectation values of $a^\dagger$ in a particle-number eigenstate $|n\rangle$, it's obvious that the expectation value vanishes because $a$ and $a^\dagger$, and they're just a different way of writing linear combinations of $\vec B(\vec x)$ or $\vec E(\vec x)$, are operators that change the number of particle excitations by one or minus one (the same for all other fields including the Dirac fields).
So if you want to mimic a classical field or classical wave with nonzero expectation values of the fields, of course that you need to consider superpositions of states with different numbers of particle excitations! But it's still true that all these expectation values are already encoded in the one-particle states. Let me summarize it: the right states that mimic the classical configurations are $\exp(Y)|0\rangle$ where $Y$ is a linear combination of creation operators (you may add the annihilation ones but they won't make a difference, except for the overall normalization, because annihilation operators annihilate the vacuum). Such coherent exponential-shapes states have nonzero vevs of any classically allowed form that you may want. At the same moment, the exponential may be Taylor-expanded to $(1+Y+\dots)$ and the linear term $Y$ produces a one-particle state that is the ultimate "building block" of the classical configuration. But if you actually want to calculate the vevs of the fields, you can't drop the term $1$ or others, either: you need to include the contributions of the matrix elements between states with different numbers of the particle excitations.
A: As explained by Iwo Bialynicki-Birula in the paper quoted, the Maxwell equations are relativistic equations for a single photon, fully analogous to the Dirac equations for a single electron. By restricting to the positive energy solutions, one gets in both cases an irreducible unitary representation of the full Poincare group, and hence the space of modes of a photon or electron in quantum electrodynamics.
Classical fields are expectation values of quantum fields; but the classically relevant states are the coherent states. Indeed, for a photon, one can associate to each mode a coherent state, and in this state, the expectation value of the e/m field results in the value of the field given by the mode.
For more details, see my lectures
http://arnold-neumaier.at/ms/lightslides.pdf 
http://arnold-neumaier.at/ms/optslides.pdf 
and Chapter B2: Photons and Electrons of my 
theoretical physics FAQ.
