This does not answer the original question. A reply to the comment field.
About Gauge fixing
I write some things about gauge fixing within my knowledge.
Not in electromagnetic but in elastic problems, incompressible materials such as rubber is a problem type that some authors pay great attention to. The incompressible property is described by the equation $\text{div}\vec{u}=0$, where $\vec{u}$ is the displacement vector, and Poisson's ratio is 0.5. This feature is not easy to implement within the standard FE method for elastic materials and many advanced authors are trying to solve it. As restraining incompressible elastic problem is a difficult task, I am inspecting that we encounter similar difficulty in electromagnetic gauge fixing.
As for the FEM of low-frequency electromagnetism, the reinvention of the Whitney element (vector finite element, edge element) seems to be successful. If the element shapes are limited to tetrahedron, rectangular box, rectangular prism, the Coulomb gauge $\text{div}\vec{A}=0$ is automatic inside the element. I would like to use this method.
Scalar/longitudinal approximate cancellation of photon wave.
First, remember the elastic wave problem (not the photon problem) again.
When solving the elastic wave problem, 3-component displacement vectors $u_x,u_y$ and $u_z$ are usually chosen as fundamental variables. This number 3 corresponds to S-wave (2 components) and P-wave (1 component). Since 3=2+1, it is quite OK, no problem.
However, for high frequency electromagnetic wave problems, if you choose 4-component vector $A_0,A_x,A_y$ and $A_z$ as the fundamental variable, there should be 4 wave modes: S-wave (2-component), P-wave (1-component) and scarar wave (1-component). But we know that electromagnetic waves are predominantly 2-component transverse waves, so we should hide the remaining 2-components.
As a real-world example of such hiding, the infinitesimal linear wire antenna has an analytic solution; see what I wrote in here; it is taken from Balanis' book. As you can see from equation 4-10a, the longitudinal electric field $E_r$ is not zero near a metallic antenna. In deriving 4-10a there is some cancellation of the longitudinal component and the scalar component (I have a note of my own, but it was simply an assertion of Balanis' book description, so I will not write it here). It is an example of the cancellation of the scalar photon and the longitudinal photon. Only the transverse wave survives as it moves away from the metal antenna.