Wigner distributions for an electron

We study the Wigner distributions for a physical electron, which reveal the multidimensional images of the electron. The physical electron is considered as a composite system of a bare electron and photon. The Wigner distributions for unpolarized, longitudinally polarized and transversely polarized electron are presented in transverse momentum plane as well as in impact parameter plane. The spin-spin correlations between the bare electron and the physical electron are discussed. We also evaluate all the leading twist generalized transverse momentum distributions (GTMDs) for electron.


Introduction
What is the shape of an electron? Most probably our answer is spherical but electron is a point particle having no underlying substructure. In quantum electrodynamics (QED) electron is considered as an elementary field. When we talk about the electron "structure", what we are actually doing is probing the fluctuations in the quantum theory. According to quantum theory, one can consider the fluctuations of electron into electron-photon pair i.e., e → eγ → e with same quantum number. The virtual photon further broke up into virtual electrons and positrons with all possible combinations. Thus, as a result bare electron is no more an isolated particle but it is surrounded by virtual cloud of electrons, positrons and photons. Therefore, bare electron becomes a dressed electron and one can consider electrons, positrons and photons as composite particles in the original electron.
The structure of the electron can be revealed when virtual cloud interacts with a probe, the parton content of the electron is resolved. Recently Wigner distributions have been used as powerful tools to understand the multi-dimensional images of the hadrons. Ji [1] first introduced the Wigner distributions in quantum chromodynamics (QCD), commonly known as phase-space distributions. These distributions reduce to generalized parton distributions (GPDs) [2] and transverse momentum distributions (TMDs) after certain phase-space reductions. GPDs and TMDs are measurable in high energy experiments (For review on these distributions and experiments to measure them see [3][4][5][6]). GPDs provide us the method for spatial imaging of nucleon [7][8][9][10][11] by using the definition of impact-parameter dependent parton distribution functions (ipdpdfs) which resolve the correlation between the parton distribution in impact-parameter plane and longitudinal momentum of different parton and target polarizations. On the other side, TMDs contains three dimensional information regarding the spin-spin and spin-orbit correlations in momentum plane [12][13][14].
Therefore, complete understanding of internal structure can be gained by combining the position and momentum distribution called as Wigner distribution. Wigner distributions contain the one-body information of partons. Wigner distributions have been used in many areas of physics for example in quantum molecular dynamics, quantum information, heavy ion collision and signal analysis [15][16][17] and measurable in some experiments [18][19][20]. Wigner distributions have been studied in different models. e.g., light-cone chiral quark soliton model [21], light-front dressed quark model [22], light-cone spectator model [23], AdS/QCD inspired quark-diquark model [24,25] etc. Wigner distributions are also related with generalized parton correlation functions (GPCFs) [26] and integrating over the light-cone energy of the quark provide us the generlaized transverse momentum distributions (GTMDs) [27][28][29]which have been widely studied in light-front quark model and AdS/ QCD inspired quark model [21,24,25] models. GTMDs are also called as mother distributions as they give rise to GPDs and TMDs. Therefore, Wigner distributions are important tools to understand the internal structure of nucleons. A nice attempt has been made to understand the structure of electron in QED in Ref. [30] where several properties such as form factors, generalized parton distributions, transverse densities, Wigner distributions and the angular momentum content have been evaluated for the electron-photon component of the electron wave function. The transverse shape of electron has also been discussed in Ref. [31]. From last many years, it was believed that the quantum mechanical version of the Wigner distributions cannot be measured in experiments. But recent series of papers put this conception to an end. Recently, it was shown that GTMDs of the gluons can be accessed by diffractive di-jet production in deep-inelastic lepton-nucleus scattering [32][33][34] and in virtual photon-nucleus quasi-elastic scattering [35]. A saturation model carrying numerical calculations on the gluons at small x is also presented in Ref. [35,36]. But for the first time, a physical process which can give access to quark GTMDs by entering the exclusive pion-nucleus double Drell-Yan process is presented in Ref. [37]. An upcoming experimental facility Electron Ion Collider [38] in USA and JLab 12 GeV [39] upgrade is gearing up to measure these distributions . In addition to this, azimuthal Wigner distributions for twisted photons is also measured [40].
In the present work, we investigate the Wigner distribution of an electron, which provide the multidimensional images on the distribution of QED partons in the dressed electron. We consider the self-consistent model of spin-1 2 composite system based on the quantum fluctuations of electron in QED [41]. We represent the spin-1 2 composite system with a mass M consist of spin-1 2 fermion (electron) and spin-1 vector boson (photon) with respective masses m and λ. This model is also act as a guideline to nucleon strcture because electron-photon component of the physical electron can be used as a component of the nucleon wave function consists of a quark and a vector diquark and therefore it probes the structure of nucleon and widely used for the calculations of gravitational form factors and spin and orbital angualr momentum of a composite relativistic system [41], GPDs [42,43], ipdpdfs [44], charge and magnetization densities [45]. Recently this approach has been used to understand the electron TMDs in momentum plane [46]. Using the overlap of light-front wave functions (LFWFs) we evaluate the Wigner distributions for electron dressed with another electron having different polarizations (unpolarized, longitudinal polarized and transversely polarized). The QED-based LFWFs can be useful since they are frame-independent and satisfy properties such as J z conservation. The spin-spin correlations for longitudinally and transversely polarized spin-1 2 composite system and fermion constituent, transversely polarized spin-1 2 composite system with longitudinally polarized fermion constituent and longitudinally polarized spin-1 2 composite system with transversely polarized fermion constituent are also investigated. Further, all the leading twist GTMDs for a physical electron are evaluated from the Wigner distributions. Presently we do not have any experimental data over the distribution of QED partons in both transverse position and monetum space and therefore a direct comparison of our calcuclations is not possbile. Experimental facilities like ACME collaboration [47,48] is putting efforts to understand the shape of electron by means of atomic studies i.e., measuring the charge distributions and electric dipole moment (EDM) [49]. They improved the upper limits for electric dipole moment of the electron and focused on the non-spherical shape of electron. New limits obtained there suggests the deviation of the shape of electron from spherical symmetry. Within the standard model, a non-spherical shape of the electron arises naturally via the quantum fluctuation in QED. The manuscript is organized as follows. In section 2, we discuss the QED model. In section 3, we discuss the Wigner distributions and present the analytic expressions for the different Wigner distribution functions. Results and discussion are given in section 4. In section 5, we discuss the spin-spin correlations and in section 6, we calculate all the leading twist GTMDs. Transverse shape of the electron is given in section 7. Conclusions are drawn in section 8.

QED model
We evaluate the results for the Wigner distribution of the physical electron by considering it as a two particle state (electron and photon). The two particle Fock state for an electron with J z = + 1 2 has four possible combinations [41]: The wavefunctions can be represented as (2) Similarly, two particle Fock's state for an electron with J z = − 1 2 also has four possible combinations: where ψ ↓ This one loop model is self consistent since it has correct correlation of different Fock's components of the state as given by light-front eigen value equation.

Wigner Distributions
Wigner distribution in the light-front framework is defined as [21,26,50,51] where with Γ, for example γ + , γ + γ 5 , iσ j+ γ 5 and S is the spin of the composite system. We have defined the P = (P + , P − , ∆ ⊥ 2 ) and P = (P + , P − , − ∆ ⊥ 2 ) are the initial and final momentum of the composite system. From Eq. 7, we can write the quark-quark correlator by assuming different operators (γ + , γ + γ 5 , iσ 1+ γ 5 , iσ 2+ γ 5 ). Wigner distribution is a fivedimensional(b x , b y , p x , p y , x) quantity. To obtain the phase-space distribution, we will integrate it over x so that the distribution is in impact-parameter and transverse momentum plane and therefore called as transverse Wigner distribution. In this work, we present the Wigner distribution in impact-parameter plane with fixed transverse momentum and in transverse momentum plane with fixed value of impact-parameter. One can also obtain the TMDs and ipdpdfs from the Wigner distributions. Wigner distributions can be presented in mixed plane i.e., by integrating over a transverse momentum and an impact-parameter along two transverse directioñ These mixing distributions are real distributions and the remaining variables are not violated by the uncertainty principle. They describe the correlation of impact-parameter and transverse momentum in perpendicular directions. Following, we have presented the results for the Wigner distribution by combining the various polarization configurations of unpolarized (U), longitudinal polarized (L) and transversely polarized (T) physical electron and internal electron. Among 16 independent twist-two Wigner distributions, we define the unpolarized Wigner distribution the unpolarized-longitudinally polarized Wigner distribution the unpolarized-transversely polarized Wigner distribution the longitudinal-unpolarized Wigner distribution the longitudinal polarized Wigner distribution the longitudinal-transversely polarized Wigner distribution the transverse-unpolarized Wigner distribution the transverse-longitudinal Wigner distribution the transverse Wigner distribution and the pretzelous Wigner distribution Here, first subscript in ρ represent the polarization state of physical electron and second subscript represents the polarization state of internal electron. Using Eqs. 1 and 5 in Eq. 7, the correlator function W [Γ] (∆ ⊥ , p ⊥ , x; S) can be expressed in terms of LFWFs as where σ i are the Pauli spin matrices. GTMDs in unpolarized, longitudinally polarized and transversely polarized cases can be represent in the form of LFWFs by Eqs. 19, 20 and 21. Using the overlap of LFWFs, we can explicitly represent the expression for the different Wigner distributions as for i = j = 1 and i = j = 2, for i = 1, j = 2 and for i = 2, j = 1, where and

Numerical results and discussion
In this section, we discuss the numerical results for Wigner distributions of electron with different polarizations configurations. The Wigner distributions in the pure transverse plane are obtained by integrating them over x. To get the complete contribution correctly at x = 1, one needs to include the contribution from single particle Fock state expansion of the physical electron state i.e. |e − = |e − + |e − γ . The form of the single particle contribution to Wigner distribution function is [41,52]. This corresponds that the single particle carries all the longitudinal momentum at b ⊥ = 0 and the transverse momentum p ⊥ is also zero. In the present work, we fix the non-zero value of b ⊥ in p ⊥ plane and vice-versa, which makes ρ single−particle (b ⊥ , p ⊥ , x) = 0, thus, the Wigner distributions do not get the contribution from the single particle component. The LFWFs include the three mass parameters M , m and λ where M is the mass of physical electron, m is the bare electron mass and λ is the photon mass. The value of masses are taken as M = 0.51M eV , m = 0.5M eV and λ = 0.02M eV . The non-zero mass of photon is used to regulate the infrared divergence [30]. In order to obtain the results in impactparameter plane (b x − b y ), one needs to fix the value of transverse momentum and integrate over x. Here, we perform the numerical integration by setting the higher limit of x = 0.95. Similarly, we need to fix the value of transverse position in order to obtain the results in 9 M eV and p y = 0.8 M eV to show convergence using Levin's integration method [53,54]. momentum plane (p x − p y ). In the case of mixed plane representation one can integrate over p x − b y or p y − b x plane to get the distributions. In both cases i.e. (in momentum and mixed plane) we again impose the same higher cut-off on the integration over x.
Further, the integration in Eqs. (22)(23)(24)(25)(26)(27)(28)(29) is set over ∆ ⊥ having range from -∞ to ∞. We perform the numerical integration by applying a suitable cut off on the ∆ ⊥ (∆ max ). We employ the Levin's integration technique [52][53][54] to handle the oscillatory behaviour of integrands. The ∆ max dependence of the Wigner distributions are shown in Fig. 1 upto 1000 M eV . It can be noticed that ρ U U and ρ LL is constant after the ∆ max = 6 M eV and hence this represents the cut-off on ∆ ⊥ integration. Similarly for the other Wigner distributions, the cut-off on ∆ ⊥ integration is ∆ max = 20 M eV . The Wigner distributions ρ T U and ρ T L vanish in this QED model. Recently, the Wigner distributions of a quark dressed with gluon have been studied in light-front dressed quark model [55], where the authors have also used the same cut-off procedure over the ∆ ⊥ integration to handle the oscillatory behaviour of the integrands.
We show the first Mellin moment of the Wigner distributions ρ U U (x, b ⊥ , p ⊥ ) and ρ LU (x, b ⊥ , p ⊥ ) for a composite system of bare electron and photon in Fig.3 which represent the transverse phase-space distribution of the unpolarized and longitudinally polarized physical electron when the constituent bare electron is unpolarized. Fig.3(a) and Fig.3(d) show the distributions ρ U U (b ⊥ , p ⊥ ) and ρ LU (b ⊥ , p ⊥ ) in impact-parameter plane respectively with fixed transverse momentum p ⊥ alongx for p x = 0.8 M eV whereas the variation of the distributions in the transverse momentum plane are shown in Fig.3(b) and Fig.3 with fixed impact-parameter b ⊥ alongx for b x = 0.8 M eV −1 . The mixing distributions of ρ U U and ρ LU are shown in Fig.3(c) and Fig.3(f), respectively. For the mixed plane representation, we present the results in (p x , b y ) plane for both ρ U U and ρ LU . In the present model, we obtain ρ U L = ρ LU , therefore results are not presented for ρ U L . One can notice that the distributions ρ U U in transverse momentum plane as well as in impact-parameter plane are circularly symmetric but in momentum plane, the peak of the distribution is in negative direction. The distributions ρ LU in both the impact-parameter and the momentum planes exhibit dipolar patterns due to privilege direction but the peaks of the dipolar pattern are opposite to each other. We observe the quadrupole structure for ρ LU in mixed plane representation. Such patterns are also observed in recent study of quark Wigner distributions in light-front dressed quark model where a quark is dressed with gluon [55]. This model is analogous to the present QED model. In the GPDs limit, ρ U U reduces to the GPD H in impact-parameter plane which by integrating over x provides charge density of the unpolarized physical electron. In the TMDs limit, ρ U U reduces to the unpolarized TMD f 1 which gives the electron density in momentum plane. The Wigner distribution ρ LU is connected with the Fourier transform of the GTMD F 1,4 [24] and the multipole structure is due to the presence of factor ij ⊥ p i ⊥ ∂ ∂b j ⊥ which usually break the right-left symmetry.
Futher they are also used to study the orbital angular momentum issues [21,24,51]. For longitudinal-unpolarized Wigner distribution, no twist-2 TMDs or IPDs are related and therefore they will vanish at the TMD/IPD limit. Concentrating on the physical interpretation, these distributions reflects the correlation between the electron spin and orbital angular momentum if C z > 0, electron spin and OAM is parallel to each other and if C z < 0 then electron spin and OAM are anti-parallel to each other. Further, OAM can be calculated from the F 1,4 GTMD as Wigner distributions do not have a probabilistic interpretation due to uncertainty principle. We integrate over p y and b x dependence which in result giving us the probability densities and correlating p x and b y , and this correlation is not violated by the uncertainty principle. As Wigner distributions have only quasi-probabilistic interpretation, a original probabilistic interpretation can be obtained by integrating over b ⊥ and p ⊥ which reduce them to TMDs and GPDs respectively. In Fig. 4, we present the results for the Wigner distribution ρ U T (b ⊥ , p ⊥ ). They describe the distributions when the bare electron is transversely polarized in an unpolarized physical electron. In impact-parameter plane, ρ U T shows the dipole behaviour whereas in momentum plane it shows the quadruple structure. In the IPD limit, ρ U T reduces toH T together with some other distributions and in the TMD limit it reduces to h ⊥ 1 (Boer-Mulder's function) whereas former corresponds to T-even part and latter corresponds to T-odd part. In the mixed plane representation, ρ U T gives the dipole structure. The physical interpretation of ρ U T is connected with the spin structure. ρ U T only corresponds the correlation between the the transverse polarization of the bare electron with its perpendicular transverse coordinate. It is clear from Eq.(24) that this distribution vanishes when electron intrinsic transverse coordinate is parallel to polarization i.e. electron transverse spin has no correlation with its parallel transverse coordinate. In Fig. 5, we show the longitudinal and transverse Wigner distributions i.e. ρ LL and ρ T T which describe the distributions when both the bare and physical electrons are longitudinally and transversely polarized, respectively. In impact-parameter plane, we fix the transverse momentum p ⊥ = p ⊥êx with p ⊥ = 0.8M eV and ρ LL and ρ T T shows a circularly symmetric behavior as similar as ρ U U . For the distributions in momentum plane, we fix impact-parameter b ⊥ = b ⊥êx with b ⊥ = 0.8M eV −1 and observe that both ρ LL and ρ T T again exhibit a similar nature as ρ U U and the peak is in negative direction in momentum plane. In Fig. 5(c) and (f), we show the distributions in the mixed plane for ρ LL and ρ T T respectively. In the TMD limit, transverse Wigner distribution reduces to transversity distribution h 1 when the polarization of the bare electron is parallel to the physical electron. Here, we choose the polarizations alongx-axis. Further, the bare electron and the physical electron can be polarized perpendicular to each other which corresponds to pretzelous Wigner distribution. The pretzelous Wigner distribution vanishes in this QED model. However in QCD spectator model [23] it is nonzero but it also vanishes in the dressed quark model [55]. In Fig. 6, we present the longitudinal-transverse ρ LT Wigner distributions. ρ LT describes the correlation between transverse spin of bare electron and longitudinal spin of physical electron. This distribution vanishes when the transverse momentum of the bare electron is perpendicular to its polarization. This illustrates that the transverse spin of bare electron has strong correlation only with the transverse momentum parallel to its polarization. In impact-parameter plane, we observe a peak at b x = b y = 0 whereas in momentum plane distribution is dipolar in nature. In mixed-plane, it again shows a dipolar structure. An analogous result has also been observed in a QCD spectator model [23] and in dressed quark model [55].

Spin spin correlations
The Wigner distribution with the composite system helicity Λ and fermion constituent helicity λ is defined for Γ = γ + 1+λγ 5 2 and S = ΛŜ z as [21] which can be decomposed as with Λ =↑, ↓ and λ =↑, ↓ where ↑ and ↓ are corresponding to +1 and −1 for longitudinal polarizations respectively. In this model ρ U L = ρ LU , thus, for Λ =↑ and λ =↑ and for Λ =↑ and λ =↓ Similarly, one can also consider the case when composite system has transverse polarization Λ T =⇑, ⇓ and the fermion constituent has transverse polarization λ T =⇑, ⇓. The expression of Wigner distribution for this situation can be written as [25] which can be decomposed as We can also define the Wigner distribution for longitudinally polarized fermion constituent in transversely polarized composite system, ρ Λ T λ (b ⊥ , p ⊥ , x) and transversely polarized fermion constituent in a longitudinally polarized composite system, ρ Λλ (b ⊥ , p ⊥ , x) provides information about the correlations between the spin of composite system and the spin of bare electron in the longitudinal direction. The Wigner distributions, ρ Λλ (b ⊥ , p ⊥ , x) with the polarization of composite system Λ =↑ and bare electron polarization λ =↑ are shown in Fig.7 which implies that the physical electron is align parallel to the bare electron. The lower panel represents top view of the same distributions. We observe that in impact-parameter plane, the peak of the distribution is slightly shifted sideways toward the +b x direction which is due to ρ LU , as it is dipole in nature and therefore it distorts the distribution significantly. No significant shifting of peak is observed in momentum plane but a distortion is observed in the mixed plane. The results for the spin-spin correlation for the case of Λ =↑ and λ =↓, (which implies that the physical electron is polarized opposite to bare electron) are shown in Fig. 8. In the present case, the effective spin-spin correlation is given by Eq.(37), i.e., distribution is not getting any contribution from ρ LU , thus there is no shifting of peaks. In the mixed plane representation (Fig. 8(b) and (e)), the peak is symmetric and there is no shifting from central axes. It should be mentioned here that in the light-front quark diquark model ρ LU = −ρ U L [25] which in result does not contribute towards the distortion in ρ ↑↑ distribution. In Fig. 9, we plot the transverse Wigner distribution ρ ⇑⇑ using Eq. (39). It represents the spin-spin correlation of a transversely polarized fermion constituent in transversely polarized composite system. This correlation consist of contributions from ρ U U , ρ U T , ρ T U and ρ T T . Out of these four contributions, ρ T U vanishes in this model. One finds that in impact-parameter plane, the distribution is symmetric about the center of the plane as both ρ U U and ρ T T are symmetric whereas ρ U T is also contributing which is dipolar in nature but with very small magnitude and therefore does not contribute to provide any distortion. In momentum plane the distribution has a sharp peak in negative direction but no shifting of peak has been obtained. Effective contributions of ρ U U , ρ U T and ρ T T towards the mixed distributions is shown in Fig. 9(c). The Wigner distribution, ρ ⇑⇓ which corresponds the spin-spin correlation of a transversely polarized fermion constituent with λ T =⇓ in transversely polarized composite system with Λ T =⇑ has been shown in Fig.  10. In this case, again ρ T U vanishes. The qualitative nature of ρ U U and ρ T T are the same but due to the subtraction they cancel each other and therefore ρ U T plays vital role in providing the dipole nature of the distribution in impact-parameter plane. Howevere, ρ ⇑⇓ exhibit the quadrupole structure in momentum plane. In mixed plane, we obtained dipole structure of the distribution. Results for the ρ Λ T λ which describes the spin-spin correlation between a longitudinally polarized fermion constituent with polarization λ and transversely polarized composite system with polarization Λ T are presented in Fig. 11 and Fig. 12 for Λ T =⇑, λ =↑ and Λ T =⇑, λ =↓ respectively. Since, ρ T U and ρ T L vanish in present QED model therefore effective contributions are only from ρ U U and ρ U L . In impact-parameter plane, we observe that ρ ⇑↑ shifts towards +b y direction whereas ρ ⇑↓ gets shifted in opposite direction. It is due to the fact that for Λ T =⇑ λ =↑, the contribution from ρ LU distorts the symmetric nature of ρ U U . Similarly when Λ T =⇑ λ =↓ one obtains the same distortion in opposite direction. However in momentum plane representation, no distortion has been observed for both cases.
In Fig. 13 and Fig. 14, we show the distribution ρ Λλ T where the spin-spin correlation is considered of a transversely polarized fermion constituent with polarization λ T in a longitu- dinally polarized composite system with the polarization Λ. Fig. 13 is for Λ =↑, λ T =⇑ and Fig. 14 describes the case when Λ =↑, λ T =⇓. In both the cases, unlike ρ Λ T λ , in impactparameter plane ρ Λλ T gets distorted along the −b x direction. ρ Λλ T is getting contribution from ρ LU , ρ U T and ρ LT which break the symmetric nature of ρ U U . Spin-spin correlation distributions immediately reduce to spin densities in impactparameter plane [10,11,56] when one integrates them over p ⊥ . In general, spin density ρ(x, b ⊥ , λ, Λ) and ρ(x, b ⊥ , λ T , Λ T ) gives the probability of finding a parton with momentum fraction x and transverse position b ⊥ with light-cone helicity λ in a composite system with longitudinal polarization Λ or with a transverse spin λ T in the composite system with transverse spin Λ T . We find that the longitudinal spin density ρ(x, b ⊥ , λ, Λ) is related to the Wigner distribution ρ Λλ (b ⊥ , p ⊥ , x) defined in Eq. (35) and the transverse spin density is connected with the transverse spin-spin correlation ρ Λ T λ T (b ⊥ , p ⊥ , x) defined in Eq. (39). ρ(x, b ⊥ , λ, Λ) involves the GPDs H(x, b 2 ) andH(x, b 2 ), thus comparing with Eq.(35), one can easily find out that ρ U U and ρ LL reduce H(x, b 2 ) andH(x, b 2 ) respectively and ρ U L and ρ LU give zero when we integrate them over p ⊥ . Similarly, in the case of transverse polarization On the other hand, integrating over b ⊥ the correlation distributions for electron of definite longitudinal or transverse polarizations, one can obtain spin densities in the momentum plane [11,46] which are again expressed in terms of TMDs. Hence with these connections, the correlations between the electron distributions in impact-parameter plane and longitudinal momentum as well as the three-dimensional information about the strength of different spin-spin and spin-orbit correlations in the momentum plane for different polarizations of bare electron and physical electron can be obtained from the Wigner distributions.

GTMDs for electron
Generalized transverse momentum distributions [12,13,26] which are known as the mother distributions of GPDs and TMDs can be extracted from different Wigner distributions. For the leading twist, the Wigner correlator, Eq. 7, can be parametrized in terms of (a) (e) (f) GTMDs as [26] W In Fig. 15, we show the GTMDs F 1,1 , F 1,3 and F 1,4 with fixed value of p ⊥ = 0.3M eV but with different values of ∆ ⊥ (left panel) and fixed value of ∆ ⊥ = 0.3M eV with different values of p ⊥ (right panel) respectively. It can be noticed that F 1,1 and F 1,3 diverges at x = 1. This is due to the fact that we consider only the quantum fluctuation into the |eγ Fock's sector which encodes the information on the structure of the physical electron, but if one considers the state of a physical electron as |e phy → |e + |eγ , the single particle(|e ) also contributes to the F 1,1 at exactly x = 1 which, with proper normalization, cancels the divergence. Since the GPD H is related to the GTMD F 1,1 , a similar behavior of F 1,1 has also been observed in H for physical electron [57]. We find F 1,2 = 0 in this model. For fixed p ⊥ with increasing values of ∆ ⊥ the peaks in F 1,4 shift towards higher values of x and the height of the peaks increase. But for fixed ∆ ⊥ with increasing values of p ⊥ , peaks move towards lower values of x however magnitude of GTMDs decreases. In Fig. 16, we present the results for G 1,1 , G 1,2 and G 1,4 GTMDs with same parameters used in Fig. 15. We observe that G 1,1 = −F 1,4 and G 1,3 = 0. We show the GTMDs H 1,2 , H 1,3 , H 1,6 and H 1,7 in Fig. 17

Transverse shape of electron
The shapes of an electron can be obtained by using the following relation which was first introduced in [58], where φ is the angle between p ⊥ and S ⊥ and φ n is the angle between n and S ⊥ . n is the unit vector which describes the arbitrary spin of the particle in a fixed direction. S ⊥ is the physical electron polarization in the transverse direction. Further, f 1 , h 1 and h ⊥ 1T are the unpolarized electron distribution, transversity, and pretzelous distributions respectively and f (p 2 ⊥ ) = dxf (x, p 2 ⊥ ). The definition of the TMDs (f 1 , h 1 , h ⊥ 1T ) can be found in Ref. [46] and their expressions in the present QED model are given by where Since the pretzelous distribution h ⊥ 1T is zero in this model, thus the shape of an electron explicitly depends upon f 1 and h 1 . In Fig. 18 we present the shape of electron by considering different angle φ n . We show that how the shape of an electron emerges when we take the different angles between n and S ⊥ . We consider the value of φ n ranging from 0 to 2π (φ n = 0 or 2π reflects that n is parallel to S ⊥ and φ n = π reflects that n is anti-parallel to S ⊥ ) and increase the value of p ⊥ . We observe that with increasing the value of p ⊥ , the shape gets distorted (moving from black color→ brown color). Since we are excluding the p ⊥ = 0 point, we present the results starting from p ⊥ = 0.01M eV which comes out to be spherical and thereafter by increasing the value of p ⊥ , we observe that distortion increases in the shape of electron.

Conclusions
We discussed all the leading twist Wigner distributions of a physical electron which provide the multi-dimensional images of electron. Here, we used the perturbative QED model where one can consider physical electron as an effective composite system of electron and photon. Using the overlap of LFWFs the Wigner distributions have been evaluated for the different polarization configurations e.g. unpolarized, longitudinally polarized and transversely polarized composite system and fermion constituent. The spin-spin correlations have been investigated for the system. We have also evaluated the GTMDs of electron. In this model we have found 10 nonzero leading twist GTMDs. Results provide rich and interesting information on the distribution of QED partons in impact-parameter, momentum plane and mixed planes representation. We have also discussed the shape of electron.