Electromagnetic heavy-lepton pair production in relativistic heavy-ion collisions

We calculate the cross sections of electromagnetic productions of muon- and tauon-pair productions from the ultra-relativistic heavy ion collisions. Since the Compton wavelengths of muon and tauon are comparable to the radius of the colliding ions, nuclear form factors play important roles for calculating the cross sections. Recent measurement (Abrahamyan et al., Phys Rev Lett 108:112502, 2012) indicates that the neutrons are differently distributed from the protons; therefore this affects the cross section of the heavy-lepton pair production. In order to see the effects of the neutron distributions in the nucleus, we used analytical expression of the Fourier transforms of the Wood–Saxon distribution. Cross section calculations show that the Wood–Saxon distribution function is more sensitive to the parameter R compared to the parameter a.


Introduction
Colliding beams of heavy ions at energies per nucleon in the range 10 2 -10 4 GeV, corresponding to the RHIC experiments at Brookhaven, and LHC at CERN, suggested that intense electromagnetic fields of these ions can produce relatively large fluxes of muons and tauons. Since the mass of the electron (0.511 MeV) is much smaller than the muon (105.66 MeV) and the tauon (1784 MeV) masses, and the Compton wavelength of the electron (386 fm) is much larger than the muon (1.86 fm) and the tauon (0.11 fm) masses, the cross sections of producing electron pairs are much larger than the heavy leptons. In addition to this, since the Compton wavelengths of muons and tauons are smaller than the radius of the colliding heavy ions (Au, Pb), the nucleon [1] and nucleus form factors are not negligible. Therefore realistic charge form factors play important roles for calculating the cross sections of the heavy-lepton pair productions. a e-mail: melekaurora@yahoo.com In this work, we have used the Wood-Saxon nucleus form factor which is also widely used in the literature [2][3][4][5]. Instead of the monopole approximation, we have an analytical expression of the Wood-Saxon form factor in our equations. With the help of this analytical equation, we are able to account for the neutron skin and its effects to our calculations. The determination of the sizes and shapes of atomic nuclei is one of the important problems in nuclear physics. The rms radius of the charge distribution in nuclei is well known, with a high degree of accuracy achieved by the elastic electron-nucleus and muon-nucleus scattering experiments. The experimental uncertainties for the nuclear charge radius are about 1 % for several nuclei [6]. However, the neutron distribution and its rms radius in nuclei is not precisely determined theoretically and experimentally. Generally, the neutron skin thickness is defined as the neutron-proton rms radius difference in the atomic nucleus, where the measured charge radius is r 2 1/2 p = 5.43 fm for the gold nucleus and r 2 1/2 p = 5.50 fm for the lead nucleus [2]. By using the parity-violating electron scattering (PVES), the recent experiment [7] at JLab gives us the value R n as 5.78 +0. 16 −0. 18 . The difference between the radii of the neutron and proton distributions is the neutron skin and it is equal to R skin = R n − R p = 0.33 +0. 16 −0.18 fm for the lead nucleus with large total error. Therefore the difference between the distributions of protons and neutrons in the nucleus affects the cross section of the heavy-lepton pair production.
To explore a surface thickness of the neutron density in a heavy nucleus, we model the neutron density with a Wood-Saxon form: where R 0 and a n are the parameters that determine the shape of the neutron distribution. Although the same parameters are precisely determined for the proton distribution in the heavy nucleus, there are large uncertainties for the neutron distribution. In this work, we study the effect of these parameters on the heavy-lepton pair production cross sections.
The main motivation at this period was the RHIC. Dilepton production from the central collisions of heavy ions is very important to understand the hadronic interactions. However, electromagnetic production of dileptons from the peripheral collisions of heavy ions can shield those dileptons produced from the hadronic interactions. Therefore it is highly important to understand various properties of the electromagnetically produced dileptons [17].
In our previous work [11,[18][19][20], we have calculated electron-positron pair production cross sections by using the second order Feynman diagrams. In this work, we have modified the equations by including the electromagnetic form factors of nuclei in the momentum space. We have used Fermi type charge form factor. We also consider the coherent production of heavy-lepton pairs in two-photon limit for the symmetric collision of two gold nuclei at RHIC and two lead nuclei at LHC energies. In Fig. 1, the wiggly lines are the virtual photons that make the electromagnetic field of the colliding ions. Here we start with the lowest non-vanishing Feynman diagram which is called a two-photon diagram. We take the classical limit of the motion of the ions, which yields the external field model. Since the parameter Z α ≈ 0.6 is not small, Coulomb and unitarity corrections should be included to obtain exact cross sections. This calculation has been done in [21][22][23], and we will also mention it briefly in the results section.
In Sect. 2, we briefly explain the two-photon pair production mechanism. In Sect. 3, we study the effects of the form Z 1 Z 2 b>R 1 +R 2 Fig. 2 Collision of two heavy ions at ultra-relativistic velocities. Lorentz-contracted strong EM fields are produced around the ions factor that we use in our calculations and we also calculate a variety of cross sections for the single-pair process. Finally on Sect. 4, we summarize all the work that we have done in this manuscript.

Calculations and results
In this work we will investigate the two-photon mechanism for heavy-lepton pair production. In Fig. 2, the electromagnetic interaction of the colliding ions is shown. Since all the protons in the nucleus act coherently, there is a very strong electromagnetic field appearing around the heavy ions [10]. Lepton pairs are produced from this field, which is a cloud of virtual photons. Although we consider peripheral collisions of the heavy-ions where the impact parameter b > R 1 + R 2 is greater than the sum of the radii of the colliding nuclei, our cross section expressions include the whole impact parameter range from zero to infinity. We consider the electromagnetic production of lepton pairs from the collisions of heavy ions with charge Z , where l + l − can be electrons, muons or tauons. The maximum magnitude of the electric field of the ions is where γ is the Lorentz factor, b is the impact parameter, and Ze is the charge of the fully stripped heavy ion. For the impact parameters around b ∼ 7.5 fm, which is very close to the nucleus but still concerns peripheral collisions, the maximum strength of the electric field at the surface is E max ∼ 2 × 10 23 V/m for the RHIC and E max ∼ 7.14 × 10 24 V/m for LHC collisions. We can estimate the minimum electric field that can create a lepton pair of mass 2mc 2 : where the distance that we used in the above equation is twice the Compton wavelength of the particular lepton [11].
To create e + e − , μ + μ − , and τ + τ − pairs, the minimum electric fields needed are 1.32 × 10 18 , 5.66 × 10 22 and 1.62 × 10 25 V/m, respectively. We see that, even at LHC energies, the maximum electric field produced from the lead ions is below the critical electric field to create tauon pairs. Even for the impact parameters around the b ∼ 7.5 fm-30 fm ranges, the maximum electric field is smaller than the electric field that produces tauon pairs at RHIC and LHC energies. However, in QED there is still a probability to produce heavy leptons even at these energies. The source currents for producing lepton pairs are the Lorentz-boosted charges of the colliding heavy ions. Electromagnetic fields in these collisions are very strong; they can pull many lepton pairs from the vacuum. Relativistic heavy-ion collider (RHIC) at Brookhaven can accelerate the fully stripped gold nuclei at 100 GeV/nucleon energies and the large hadron collider (corresponding to run-2 at LHC) can accelerate the protons up to energies of E p = 6.5 TeV. The fully stripped lead nuclei can be accelerated at LHC run-2 up to E p Z /A ≈ 2.5 TeV/nucleon energies, which yields √ s ≈ 5.0 TeV/nucleon. In the literature, most of the calculations about electromagnetic lepton pair production have been done for the lead ions at energies 3.4 TeV/nucleon, which yields √ s ≈ 6.8 TeV/nucleon. In order to compare the results with our previous work, and also with other authors' work, we will use the energy value at LHC as 3.4 TeV/nucleon as input in our calculations. Therefore, electromagnetic fields at these collisions last a very short time, t = b/(γ cβ), and the Fourier frequencies could be up to ω max ∼ γ cβ/b [11,24]. We can write the Lagrangian density as the sum of the noninteracting fermion Lagrangian L 0 and the Lagrangian density for the coupling of the classical electromagnetic potential to the lepton fields L int [25,26], is the four-vector potential of the nucleus 1 and the scalar component of the potential can be written as and the non-zero component of the vector potential A 3 1 is in the longitudinal direction and we obtain where the transverse components of the vector potentials are A 1 1 = A 2 1 = 0. The potentials from nucleus 2 can be written by the substitutions, b → −b and β → −β. In the above equation f Z (q 2 ) is the form factor of a nucleus which gives the momentum distribution of a proton in the nucleus, and G E (q 2 ) is the form factor of a proton which represents the electric distribution of the proton [11,27]. As we mentioned above, we have used two parameter Fermi (2pF) function (or Wood-Saxon function) for the charge distribution of the protons: where R is the radius of the nucleus and a is the skin depth or diffuseness parameter. These parameters are obtained by fits to the electron scattering data [28] and ρ 0 is written by the normalization condition. For symmetric nuclei, the nuclear density for a nucleus that has mass number A, a distance r from its center is modeled in the literature with a Woods-Saxon distribution as in Eq. (9). When we normalize the Woods-Saxon charge density to unity we can obtain the parameters: ρ 0 = 0.1694 A fm −3 for the Au nucleus and ρ 0 = 0.1604 A fm −3 for the Pb nucleus [3]. The radii of the gold and lead nucleus are equal to R Au = 6.38 fm and R Pb = 6.62 fm, respectively. However, recent measurement [7] indicates that the neutrons are differently distributed from the protons, so that the parameters R and a must be modified. Therefore this affects the cross section of the heavy lepton pair production [29]. In order to see the effects of the neutron distributions in the nucleus, we need to have an analytical expression of the Fourier transforms of the Wood-Saxon distribution [30]: The parameters R and a are explicitly shown in this analytical equation and it is plotted in Fig. 3 together with the numerical result from Fourier transformation of a Woods-Saxon distribution. The agreement is excellent. Since the contribution of the last term in the above equation is very small compared to the first term, it is safely neglected. By changing the parameters R and a we can estimate the sensitivity of the nucleon distributions in the nucleus.
On the other hand, for the electric distribution function of the proton, we can write the dipole form factor for the proton where the value of 2 is 0.71 GeV 2 [27]. We have observed that the effect of this distribution is negligible. Finally, the total electromagnetic potential is the sum of the nucleus 1 and 2, To calculate the total cross section for the heavy-lepton pair production, we first represent the direct term thus: Thus, by combining both integrals we obtain for the direct free-free pair production amplitude the explicit expression: p is the energy of the intermediate state. In these equations k represents the momentum of the electron and q represents the momentum of the positron. The vector p describes the momentum of the intermediate (electron and positron) states in the field of the ion and it can be decomposed into its transverse and parallel part, p = p ⊥ + p z , relative to the motion of the ions. The momentum that is parallel to the heavy ions and the frequency in Eq. (13) are fixed by momentum conservation; For the fixed momentum and the spin states, the transition matrix element can be written as ω a and ω b are the frequencies associated with the fields of ions a and b, respectively, the function F is the scalar part of the field from each heavy ion, for the frequency ω a , and for the frequency ω b , respectively. The function T explicitly depends on the velocity of the heavy ions β, on the transverse momentum p ⊥ , and on the states k, q; Having the amplitudes for the direct and crossed diagram, we are now prepared to write down the cross section for the generation of a free-free pair in collisions of two heavy ions where S = S ab + S ba denotes the sum of the direct and crossed cases. Making use of all the simplifications from above, these cross sections can be expressed as with and these being some proper products of the transition amplitudes and scalar parts of the fields as associated with ions a and b. Table 1 displays the cross sections of electromagnetic productions of electron, muon, and tauon pairs from the ultrarelativistic gold-gold collisions at RHIC energies and leadlead collisions at LHC energies with no form factor and with a Wood-Saxon form factor distribution. In this calculation, we simply used the analytical expression of the Fourier transformation of this form factor. We have also included the electric distribution function of the proton (Eq. (11)) in our calculation, and we have seen that it has a negligible effect for the production of any lepton pairs. The cross sections of electromagnetic productions of electron, muon, and tauon pairs are also calculated in [31] by using various form factors such as uniform, Gaussian, and Woods-Saxon types, and we compared the results between them. In this present work, we are Table 1 The cross sections of electromagnetic productions of electron-, muon-, and tauon-pair productions from the ultra-relativistic gold-gold collisions at RHIC energies and lead-lead collisions at LHC energies. f NFF means the form factor is not included and f W S means that the form factor is included in the calculations

Calculations and tables
using the Woods-Saxon type form factor since it is more realistic and widely used in the literature. When we compare electron pair production cross sections with and without form factor, it is clearly seen that the differences between the results are negligibly small. Since the Compton wavelength of the electron is much larger than the radius of the colliding heavy ions, the form factors of the proton and nucleus almost have no effect for the production of electron-positron pairs.
However, when we compare muon-and tauon-pair production cross sections with and without form factors, it is clearly seen that the differences are large. The cross section results for muon-and tauon-pair production with form factors are much smaller than without form factors. Since the Compton wavelengths of muons and tauons are smaller than the radius of the colliding heavy ions, the effects of the form factors become dominant for pair production. At RHIC and LHC energies, muon-pair production is reduced by factors of about 3 and 2, respectively. On the other hand, at RHIC and LHC energies, tauon-pair production is reduced by factors of about 100 and 5, respectively. Figures 4, 5, and 6 display the differential cross sections as a function of the energy p 0 , the longitudinal momentum p z , the transverse momentum p ⊥ , and the rapidity y of the produced electron, muon, and tauon for the cases with and without form factors at RHIC and LHC energies. When we compare these figures, in Fig. 4 for the electron pairs, it is clearly seen that the form factor has almost no effect on the cross section. With and without form factors, differential cross sections show the same behavior and almost identical values. The differential cross section as a function of transverse momentum decreases rapidly compared to the longitudinal momentum and energy. The differential cross sections as a function longitudinal momentum and energy have almost identical values. This also shows that the produced leptons On the other hand, In Figs. 5 and 6 the same differential cross sections of the produced muon and tauon pairs at RHIC and LHC energies with and without form factors are shown. When we compare these figures, it is clear that the longitudinal momentum of the produced muon and tauon are much higher than the transverse momentum. Therefore, the energy E = p 2 t + p 2 z + 1 is carried by the longitudinal momentum of the produced heavy leptons. It is also clear that the differential cross sections as a function of energy and longitudinal momentum are very close to each other at RHIC and LHC energies. The rapidity distribution of the muonpair production and form factor reduce the cross sections by about ten times at the RHIC energies, and they reduce the cross section by about five times at the LHC energies. For the tauon-pair production, the reduction of the cross section is about two orders of magnitude for the RHIC energies, and it is about one order of magnitude for the LHC energies. All these calculations show that distributions of  Figure 7 displays muon and tauon total cross sections for Au + Au and Pb + Pb collisions as a function of the Lorentz factor γ of the beams for the Wood-Saxon form factors. We also compare this results with the case of no form factor. For the muon-pair production, the total cross sections are reduced by factors of about 10 for the low energies, about 3 for the RHIC energies, and about 2 for the LHC energies. For the tauon-pair production, the total cross sections are reduced by factors of about 10 9 -10 6 for the low energies, about 100 for the RHIC energies, and finally about 5 for the LHC energies. Figure 8 shows the total cross section of heavy-lepton pair production as a function of the charge of the colliding nuclei Z . Obviously the nuclear form factor reduces the total cross sections. At RHIC energies, the form factor reduces the tauon-pair production cross section by about one or two orders of magnitude, however, for the muon-pair production the reduction is about a factor of 2 or 3. On the other hand, at LHC energies, the form factors reduce the muon-pair production cross section by factors of about 1 or 2, for the tauon production cross section this reduction is slightly less than one order of magnitude. Figure 9 displays the impact parameter dependence cross section for muon-and tauon-pair production at RHIC and LHC energies with and without Wood-Saxon distribution. In these calculations, the range of the impact parameter b is between 0 and infinity. Since the Compton wavelengths of the muon are 1.86 fm and for the tauon 0.11 fm, which are smaller than the radius of the nucleus, the heavy leptons are produced mainly near the radius of the heavy ions. Maximum production of muon pairs occurs at RHIC ener-gies about impact parameters of the 20 fm range, and at LHC energies in the impact parameter of the 40-80 fm range. At a range of about 200 fm, the impact parameter dependence cross section decreases by factors of about 10 at RHIC energies, and 2 at LHC energies. On the other hand, maximum production of tauon pairs occurs for impact parameters of the 7 fm range at RHIC energies, and about 10 fm at LHC energies. At about 200 fm, the impact parameter dependence of the cross section decreases by a factor of about 100 at RHIC energies and 50 at LHC energies.  Table 2 represents the percentage of cross section contributions to the total cross section between the impact parameter ranges b < R 1 + R 2 and b > R 1 + R 2 where R 1 + R 2 ∼ 15 fm for the heavy ions. The cross section calculations are taken from Table 1 with form factor. For the electron-positron pair production cross section, pairs are overwhelmingly produced for the range of b > R 1 + R 2 . For the muon-and tauon-pair production, the range b < R 1 + R 2 has an important contribution to the total cross sections. When the energy of the colliding ions increases, the contributions from the b < R 1 + R 2 range decrease. However, when the mass of the lepton increases, the contribution from the b < R 1 + R 2 range increases. We have repeated the cal-culations for the impact parameters b for the ranges 20, 25, and 30 fm. These calculations clearly indicate that the details of the form factors, such as the distributions of the nucleons in the nucleus, are important for the calculations of the cross sections.
In all the above calculations, we have also used the proton form factor G E (q 2 ). Our results show that the contribution of this proton form factor is negligible.
The description of the matter distribution in nuclei is still under investigation in nuclear physics and precise knowledge of this distribution can help us to understand the structure of the nucleus. The charge distribution has been measured with high accuracy from the electron scattering data, so that the charge radii are usually known with uncertainties lower than 1 % [1,7]. However, the neutron distribution in the nucleus is not accurately measured as for protons. We have also investigated the neutron distribution of the nucleus. The Wood-Saxon equation gives us the proton charge distribution and neutron density distribution of the nucleus. In Eq. (8), the parameter R is the radius of the nucleus and a is the skin depth and these parameters are obtained from the electron scattering data [28]. Recent experiments show that the neutrons are distributed differently from the protons and this affects the cross sections of the heavylepton pair productions. Neutrons are electrically neutral, so they do not contribute to the electric form factor. Neutrons can contribute to the hadronic radius. In order to see the effects of the neutron distribution, we have changed the parameters R and a in the Wood-Saxon form factor. First we keep the R constant and change the parameter a between 0.3 fm ≤ a ≤ 0.8 fm. Although the parameter a has a value of about 0.55 fm, by changing it for the extreme ranges, we can see how it affects the cross sections. Next we keep the parameter a = 0.55 fm constant and the change the radius R between 5.8 fm ≤ R ≤ 7.2 fm. In Figs. 10 and 11, we have plotted these cross sections for muon and tauon productions at RHIC (gold nucleus) and at LHC (lead nucleus). When the parameter a is changed from the value of 0.4 fm to 0.6 fm , the cross section ( σ/σ ) of the produced muon pairs decreases by about 4 % at RHIC and 2 % percent at LHC. Similarly, tauon production decreases by about 4 % at LHC and 14 %    Table 3. From this calculation, we can conclude that the Wood-Saxon distribution function is more sensitive to the parameter R compared to the parameter a. It also shows that the cross sections are very sensitive to the parameter R at RHIC energies compared to the LHC energies. The proton distribution in the nucleus is well known. If the neutron surface is shifted outwards, for the proton surface this means that the skin depths of neutrons and protons are nearly equal a n ∼ a p and it is called bulk or skin type of neutron skin.
On the other hand, if the diffuseness of the neutron surface is greater than the proton, then R n ∼ R p and it is called sur f ace or halo type of the neutron skin [32].
We have compared our results with the work in Ref. [33]. The two cross section results of the RHIC for μ − μ + production nearly agree; however, the τ − τ + cross section in our work approximately is 3 times larger than their results. Our LHC cross section results of μ − μ + and τ − τ + productions are slightly higher than their results. The authors in Ref. [33] used a Gaussian distribution form factor and the equivalent photon approximation to calculate the heavy-lepton pair production cross section.
When we compare our cross section results with Ref. [27], both calculations for μ − μ + production at RHIC ener- gies are almost similar with and without the form factor. At RHIC energies, the τ − τ + cross section with no form factor in our work approximately is 5 times lower than in Ref. [27], however, the results with the form factor are close to each other.
Our cross section results of μ − μ + rates at RHIC and LHC are about slightly lower than in Ref. [21], that is, using Born calculations at RHIC and LHC energies. Since the parameter Z α ≈ 0.6 is not small, Coulomb and unitary corrections σ total = σ two−photon + σ Coulomb + σ unitary (24) should be included in the calculations to obtain exact cross sections. The Coulomb correction corresponds to multiphoton exchange between the produced lepton pairs and nuclei [22], whereas the unitary correction corresponds to the exchange of light-by-light blocks between nuclei [11]. In Refs. [11,21,22] Fig. 11 The cross section of produced muon a as a function of skin thickness a at LHC energy when R = 6.62 fm b as a function of R at LHC energy when skin thickness a = 0.55 fm. The cross section of produced tauon c as a function of the skin thickness a at LHC energy when R = 6.62 fm d as a function of R at LHC energy when the skin thickness a = 0.55 fm Table 3 The ratio of the cross sections ( σ/σ ) of the produced muon and tauon pairs. In the first column, the parameter R is kept constant and the parameter a is changed between 0. the nuclei) are small, since the form factor eliminates all the exchange photon momenta for the range of k ⊥ 2 > 1/R 2 . Therefore the cross section is suppressed by (Z α) 2 /(R M) 2 where M is the mass of the heavy lepton and R is the radius of the nuclei. The calculation shows that the contribution of the Coulomb correction is less than 1 % of the total cross sections [21].
For the unitary corrections (corresponding to the exchange of light-by-light blocks between nuclei), we can consider both the exclusive (the heavy-lepton pair production without any electron-positron pair production) and the inclusive (the heavy-lepton pair production with several electron-positron pair production) heavy-lepton pair production. The authors find that for exclusive muon-pair production for Pb+Pb collisions at LHC there are about −50 % unitary corrections to the Born cross sections; however, for inclusive muon-pair production, the unitary correction is very small. In this work, we have repeated the calculations and our results for the heavylepton pair productions coincide with the above work.
All these calculations are done for all ranges of momenta and energies of the heavy leptons. However, for a given detector, electron and muon pairs can be measured with some restricted momenta. For example, the ALICE detector can measure only forward muons with pseudorapidity 4 < y < 5 and p t > 2 GeV. Similarly, the CMS detector in LHC and the STAR and PHENIX detectors in RHIC have some rapidity cuts and restricted transverse momenta. In Ref. [3] the authors made some predictions for electromagnetic exclusive muon-pair productions for the above detectors with the relevant restrictions.
The exclusive production of mesons, pairs of quark and antiquark particles, pairs of leptons in ultraperipheral heavyion collisions (UPC) have been investigated by several authors [34][35][36]. Some of the authors calculated the neutron emission from the excited nuclear system within the Hauser-Feshbach formalism and compared the results with the recent experimental data of PHENIX, PHOBOS, BRAHMS, and ALICE collaborations. Those methods can also be applied to the nuclear excitations with lepton-pair production. In our future work, we are planning to do calculations about nuclear excitations with mainly giant dipole resonance with heavylepton pair productions with the kinematic restrictions of the above detectors and compare with the experimental results.

Conclusions
In this work, we investigate the contribution of the Wood-Saxon nuclear form factors on the heavy-lepton pair production. We study the pair productions for the RHIC and LHC energies for the Au and Pb nucleus. Even though electronpair production is not strongly affected by the form factors, the calculations show that nuclear form factors play important role for the muon-and tauon-pair production. At RHIC and LHC energies, muon-pair production is reduced by factors of about 3 and 2 by the Au and Pb form factors, respectively. On the other hand, at RHIC and LHC energies, tauon pair production is reduced by factors of about 100 and 5 by the Au and Pb form factors, respectively. These calculations show that tauon-pair production is very sensitive to the form factors at the RHIC energies. We can conclude that muon-pair production in peripheral heavy-ion collisions is coherent on distances greater than the dimensions of the nucleus. However, although for the LHC energies the tauon-pair production is coherent over distances of the dimension of the nucleus, for the RHIC energies this coherence breaks down. Since the mass of the tauon is very large when compared to the muon, incoherent pair production of tauon is dominant at the RHIC energies. When the energy is increased, the coherence pair production is restored. We expect that this result is also valid for the other heavy particles.
We have also investigated the effect of the neutron distributions in the nucleus to the cross section of the heavy-lepton pair production. The neutron skin can be produced in two ways: The first type is as a result of the shift of the neutron surface in relation to the proton one (the bulk or "skin" type a n ∼ a p ), and the second type is due to the differences in the surface thickness of proton and neutron matter in a nucleus (the surface or "halo" type R n ∼ R p ). In general, both types of mechanisms contribute to the creation of the neutron skin. We have changed the parameters a and R in the Wood-Saxon form factor, and our calculations show that the Wood-Saxon distribution function is more sensitive to the parameter R compared to the parameter a.