The macroscopic Maxwell's equations can be put in terms of differential forms as
$$\mathrm{d}\mathrm{F}=0,\quad\delta \mathrm{D}=j\implies \delta j=0,\quad \mathrm{D}=\mathrm{F}+\mathrm{P}.$$
$\mathrm{F}$, $\mathrm{D}$, and $\mathrm{P}$ are 2-forms; $j$ is a 1-form. $\mathrm{F}$ is the $electromagnetic\ field$, which has 3-vector components $E_i,B_i$, usually called the $electric\ field$ and the $magnetic\ field$; $\mathrm{D}$ is a generalized electromagnetic field, generalized in the sense that $\mathrm{dD}\not=0$, which has 3-vector components $D_i,H_i$, called, amongst other names, the $electric\ displacement\ field$ and the $magnetizing\ field$; $j$ is a conserved 4-current; and $\mathrm{P}$ is the electromagnetic polarization/magnetization, which has 3-vector components $P_i,M_i$, usually called, up to a factor $\pm 1$, the $polarization$ and the $magnetization$, which in the macroscopic form of Maxwell's equations model the internal degrees of freedom, whereas in the microscopic form of Maxwell's equations we have $\mathrm{P}=0$. The Maxwell's equation Wiki page, together with its many links, does a reasonable job of summarizing all this.
In QED, there are internal degrees of freedom associated with the Dirac field, so should we be working with $\mathrm{D}$ when we quantize? Classically, given a differentiable bivector field such as $\mathrm{D}$, there is inevitably a conserved 4-current $\delta \mathrm{D}$, unless the equations of motion are such that this quantity is identically zero. [Edit, to reflect Cristi Stoica's comment] Given an arbitrary 2-form such as $\mathrm{D}$, if Hodge decomposition could be applied, we could decompose $\mathrm{D}$ uniquely as $$\mathrm{D}=\mathrm{dA}+(\mathrm{D}-\mathrm{dA}).$$
For Minkowski space, such decompositions are possible but are not unique because there is no positive-definite bilinear form on the function space; for the forward light-cone component $\mathrm{D}^+$ the decomposition can be restricted by minimizing the positive semi-definite integral $$\int [\widetilde{\delta(\mathrm{D}-dA)}]_\mu(k)[\widetilde{\delta(\mathrm{D}-dA)}]^\mu(k)\theta(k^2)\theta(k_0)\frac{\mathrm{d}^4k}{(2\pi)^4},$$
and similarly for the backward light-cone, but no such restriction is available for space-like components.
Looked at this way, the electromagnetic field $\mathrm{F}$ is just a component of the real object of interest, $\mathrm{D}$, that can be expressed as $\mathrm{F}=\mathrm{dA}$, for which, inevitably, $\mathrm{dF}=0$. The nonunique component $\mathrm{D}-\mathrm{dA}=\mathrm{P}$ of $\mathrm{D}$ cannot be expressed in terms of a single potential. [End of edit — constructed somewhat in the moment so it may have to be modified.]
If $\mathrm{D}$ is taken to be the real object of interest, with $\delta\mathrm{D}=j$ an independent degree of freedom, then does that mean that we should look for second-order differential equations for $\mathrm{D}$ in the first instance, probably in association with additional degrees of freedom, in view of the additional degrees of freedom that are present in the Standard Model?
I should point out that, for the purposes of this question, I take the Dirac wave function not to be observable in the classical Maxwell-Dirac theory, with $\overline{\psi(x)}\psi(x)$, $(\delta\mathrm{D})^\mu=j^\mu=\overline{\psi(x)}\gamma^\mu\psi(x)$, $\overline{\psi(x)}\gamma^{[\mu}\gamma^{\nu]}\psi(x)$, $\overline{\psi(x)}\gamma^{[\mu}\gamma^\nu\gamma^{\rho]}\psi(x)$, and $\overline{\psi(x)}\gamma^5\psi(x)$ being (non-independent) 0-, 1-, 2-, 3-, and 4-form observables, respectively.
This is almost not a question, for which apologies, however I am interested in any critique of details or of whatever minimally elaborated plan there might seem to be to be here, or in any reference that seems pertinent. The ideas lie somewhat in the same territory as an earlier question of mine, Gauge invariance for electromagnetic potential observables in test function form, in that I prefer to disavow unobservable fields unless there are absolutely clear reasons why they are unavoidable.
Best Answer
My recent article might be relevant: http://link.springer.com/content/pdf/10.1140%2Fepjc%2Fs10052-013-2371-4.pdf (published in the European Physical Journal C, open access). An excerpt from the abstract: "after introduction of a complex four-potential of electromagnetic field, which generates the same electromagnetic fields as the initial real four-potential, the spinor field is algebraically eliminated from the equations of spinor electrodynamics. It is proven that the resulting equations for electromagnetic field describe independent evolution of the latter"
The article may be relevant, although I do not use the formalism of differential forms, as you "take the Dirac wave function not to be observable in the classical Maxwell-[Dirac] theory" and "prefer to disavow unobservable fields unless there are absolutely clear reasons why they are unavoidable."