The polymer representations, which are partially motivated by loop quantum gravity, have been suggested as alternative schemes to quantize the matter fields. Here we apply a version of the polymer representations to the free electromagnetic field, in a reduced phase space setting, and derive the corresponding effective (i.e., semiclassical) Hamiltonian. We study the propagation of an electromagnetic pulse and we confront our theoretical results with gamma ray burst observations. This comparison reveals that the dimensionless polymer scale must be smaller than 4×10−35 , casting doubts on the possibility that the matter fields are quantized with the polymer representation we employed.Loop quantum gravity (LQG) [1][2][3], which is a prominent quantum gravity candidate, has inspired alternative matter quantization methods, known as polymer representations [4][5][6][7][8]. These alternative methods resemble LQG in that they are nonperturbative and unitarily inequivalent to the Schrödinger representation. Also, the formal way the states and the fundamental operators are expressed in the polymer representations, mimics the cylindrical functions and the holonomy-flux algebra of LQG, respectively. Moreover, the polymer representations have been considered, by themselves, as interesting alternatives to the Schrödinger quantization [9][10][11][12].Notably, most works on the polymer representations of matter fields use scalar fields or do not make contact with experimental data [13][14][15][16][17][18][19]. In contrast, our goal is to study the empirical consequences of applying such a quantization scheme to the free electromagnetic field in the framework of Ref. 7. To that end, we polymer quantize the Maxwell theory and then use well-known methods to extract the corresponding effective dynamics.As it is well known, the electromagnetic field A ν (x), ν being a spacetime index, has a U (1) gauge symmetry and, to quantize it, we utilize a reduced phase space quantization (see, for example, Ref. 20). Furthermore, we work in the Minkowski spacetime with a global Cartesian coordinate frame where t represents the time index and i, j are spatial indices. We fix the gauge by taking A t = 0 = ∂ i A i , which can be consistently imposed when there are no sources [21, chapter 6.3]. In this case the action takes the formIn this work we use a metric with signature +2 and adopt natural units, i.e., Lorentz-Heaviside units with the additional conditions c = 1 = . we use the spacetime foliation associated with constant t hypersurfaces and denote the canonically conjugated momenta by E i , resulting inThe fact that no constraints arise reflects that there is no remaining gauge freedom.To properly implement the polymer quantization we turn to Fourier space. Notice, however, that a priori we cannot assume Lorentz invariance, and thus, we do not use the standard four-dimensional Fourier transform. Instead we only perform such a transformation on the spatial coordinates. Furthermore, to have a countable number of modes, we consider the system t...