Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: the possible existence of dark matter ( dm ) is one of the most intriguing aspects of modern particle physics , astrophysics and cosmology . in order to survey the nature of dm , a number of observations and experiments have been conducted or are underway around the world . the most recent cosmological results based on _ planck _ measurements of the cosmic microwave background ( cmb ) temperature and lensing - potential power spectra indicate that dm comprises about @xmath4 of the energy density of the universe which also includes about @xmath5 baryon matter and about @xmath6 dark energy  @xcite . many theories beyond the standard model of particle physics predict natural candidates for dm , e.g. , the weakly interacting massive particles ( wimps ) which are a class of hypothetical stable neutral particles with a huge range in masses from @xmath7 gev to @xmath8 tev and interaction cross sections with normal matter ( proton ) from @xmath9 to @xmath10 @xmath11  @xcite . in terrestrial laboratory , dm might be directly detected through their elastic scattering off nuclei in particle detectors  @xcite . a number of underground dm direct detection experiments have been performed , and among them an excess of events over the expected background has been observed by cogent  @xcite , dama  @xcite , cresstii  @xcite as well as in the recent results presented by the cdms - ii(si ) collaboration  @xcite . however , these results are in strong tension with the constraints set by some other experimental groups like xenon100  @xcite , lux  @xcite and supercdms(ge )  @xcite , leaving a confusing situation for the community . this has led to a number of attempts trying to explain the discrepancy by considering atomic uncertainties  @xcite or different mechanisms that deviate from standard assumptions about dm interactions or its astrophysical distributions  @xcite . isospin - violating dark matter ( ivdm ) provides a very promising mechanism to reconcile the tension among different experiments  @xcite . within the ivdm framework , dm is assumed to couple differently with protons and neutrons , and this assumption of the isospin violation has been supported by a number of theoretical works  @xcite based on the particle physics point of view . many parameters need to be specified in the standard method of analyzing dm direct detection experiments  @xcite and recently frandsen  _ et al . _ @xcite presented a systematic discussion on the possible ways to ameliorate the tension among different experiments . they found that the tension between the cdms - ii(si ) results and the xenon100 bounds is independent of the astrophysical uncertainties concerning the dm halo and any momentum- and velocity - dependence of the cross section in particle physics , but it can be largely ameliorated or even resolved within the framework of ivdm .    besides the uncertainties in astrophysics and particle physics concerning the interaction of a dm particle scattering off a single nucleon mentioned above , the uncertainties in nuclear physics describing how the struck nucleon is distributed inside the nucleus may also play an important role in interpreting dm signals . this is because in dm direct detection experiments , the nuclear form factors are generally applied to describe the dm - nucleus cross section and the bounds on dm - nucleon cross section are then obtained accordingly . in particular , the empirical helm form factor extracted from nuclear charge distributions  @xcite has been commonly adopted in current direct detection experiments . however , the dm - nucleus interaction should be in principle described by using the form factors of both proton and neutron distributions in the nuclei rather than the charge distributions since the dm particles actually interact with the protons and neutrons in the nuclei . within the framework of relativistic mean field ( rmf ) model , chen  _ et al . _ @xcite derived the nuclear form factor for the spin - independent scattering between the wimps and nucleus , and they found that the results can deviate from the empirical helm form factor by @xmath12 to @xmath13 in a large range of recoil spectrum of @xmath14 kev . a recent work by co _ et al . _ @xcite suggests that the use of different distributions for protons and neutrons instead of the commonly used empirical charge distributions for a target nucleus could be important in interpreting dm signals , especially if ivdm is considered .    for stable nuclei , the proton distribution can be well determined from the charge distribution which can be accurately measured with electron scattering  @xcite . in contrast , the neutron distribution is usually determined from hadron scattering experiments and the results are generally highly model dependent due to the unclear non - perturbative strong interaction  @xcite . recently , the lead radius experiment ( prex ) collaboration at jefferson laboratory ( jlab ) published their results on the measurement of the parity - violating cross section asymmetry in the elastic scattering of polarized electrons from @xmath15  @xcite , which provides a model - independent determination of the neutron density distributions in @xmath15 . the prex measurement leads to a value of @xmath16 for the rms radius of the neutron distributions for @xmath15 , implying a large neutron skin thickness , i.e. , @xmath17 by assuming a point - proton rms radius of @xmath18 fm  @xcite . one can see that the obtained value of @xmath19 has a large error , indicating a large uncertainty of the neutron distribution relative to the proton distribution in @xmath15 . theoretically , it has been established  @xcite that the @xmath19 is intimately related to the nuclear matter symmetry energy which characterizes the isospin dependent part of the equation of state ( eos ) of asymmetric nuclear matter ( anm )  @xcite . in particular , it has been shown recently that the @xmath19 of heavy nuclei is uniquely determined by the density slope @xmath20 of the symmetry energy at a subsaturation cross density @xmath21 @xmath22  @xcite . these features imply that the uncertainties of @xmath19 , especially the neutron distributions , predicted by various nuclear models are essentially due to our poor knowledge about the symmetry energy . the symmetry energy is of critical importance for understanding not only the structure and reaction of radioactive nuclei , but also a number of interesting issues in astrophysics , such as the structure of neutron stars and the mechanism of supernova explosions , and has become a hot topic in current research frontiers of nuclear physics and astrophysics  @xcite . the determination of the symmetry energy provides a strong motivation for studying isospin - dependent phenomena with radioactive nuclei at a number of new / planning rare isotope beam facilities around the world , such as csr / lanzhou and brif - ii / beijing in china , ribf / riken in japan , spiral2/ganil in france , fair / gsi in germany , spes / lnl in italy , raon in korea , and frib / nscl and t - rex / tamu in usa . in this work , using the proton and neutron form factors obtained from skyrme - hartree - fock calculations by varying the symmetry energy slope parameter @xmath23 within the uncertainty region set by the prex experiment , we investigate the form factor effects in the direct detection of dark matter . this article is organized as follows . we briefly introduce in sec .  [ sec : model ] the theoretical models and methods used in the present work , and then present the results and discussions in sec .  [ sec : sym ] . finally , a conclusion is given in sec . [ sec : conclusion ] . the eos of isospin asymmetric nuclear matter , defined by its binding energy per nucleon , can be well approximated by @xmath24 in terms of baryon density @xmath25 and isospin asymmetry @xmath26 with @xmath27 and @xmath28 denoting the proton and neutron densities , respectively . @xmath29 corresponds to the eos of symmetric nuclear matter , and the nuclear symmetry energy can be expressed as @xmath30 there are no odd - order @xmath31 terms in eq . ( [ eosanm ] ) due to the exchange symmetry between protons and neutrons ( isospin symmetry ) in nuclear matter . neglecting the contribution from higher - order terms in eq . ( [ eosanm ] ) leads to the well - known empirical parabolic law for eos of anm , which has been verified by all many - body theories to date , at least for densities up to moderate values  @xcite . furthermore , around a reference density @xmath32 , the symmetry energy @xmath33 can be expanded as @xmath34 with @xmath35 the slope parameter @xmath36 characterizes the density dependence of symmetry energy around @xmath32 . it has been shown  @xcite that the neutron skin thickness @xmath37 of heavy nuclei is uniquely fixed by the density slope @xmath23 at a subsaturation cross density @xmath38@xmath22 which roughly corresponds to the average density of the nuclei . for the calculations of finite nuclei , we use in the present work the standard skyrme - hartree - fock ( shf ) approach in which the nuclear effective interaction is taken to have a zero - range , density- and momentum - dependent form , i.e.  @xcite , @xmath39 } ,   \label{v12sky}\end{aligned}\]]with @xmath40 and @xmath41 . in the above expression , the relative momentum operators @xmath42 and @xmath43 act on the wave function on the right and left , respectively . the quantities @xmath44 and @xmath45 denote , respectively , the spin exchange operator and pauli spin matrices .    the skyrme interaction in eq . ( [ v12sky ] ) includes totally @xmath46 parameters , i.e. , the @xmath47 skyrme interaction parameters @xmath48 , @xmath49 , @xmath50 , and the spin - orbit coupling constant @xmath51 . this standard shf approach has been shown to be very successful in describing the structure of finite nuclei , especially the global properties such as binding energies and charge radii  @xcite . instead of using directly the @xmath47 skyrme interaction parameters , we can express them explicitly in terms of @xmath47 macroscopic quantities , i.e. , @xmath52 , @xmath53 , the incompressibility @xmath54 , the isoscalar effective mass @xmath55 , the isovector effective mass @xmath56 , @xmath57 , @xmath58 , @xmath59 , and @xmath60 . the @xmath59 and @xmath60 are respectively the gradient and symmetry - gradient coefficients in the surface interaction part of the binding energies for finite nuclei which is defined as @xmath61 ^{2}/(2{\rho ) } .\ ] ] then , by varying individually these macroscopic quantities within their known ranges , one can examine more transparently the correlation of nuclear matter properties with each individual macroscopic quantity . recently , this correlation analysis method has been successfully applied to study the neutron skin  @xcite and giant monopole resonance of finite nuclei  @xcite , the higher - order bulk characteristic parameters of anm  @xcite , and the relationship between the nuclear matter symmetry energy and the symmetry energy coefficient in nuclear mass formula  @xcite . in the present work , we use the @xmath47 macroscopic quantities instead of using directly the @xmath47 skyrme interaction parameters . especially , we study the symmetry energy effects by varying the @xmath62 value while keeping the other macroscopic quantities unchanged .      in this work , we use the standard method of analyzing dm direct detection experiments  @xcite . the spin - independent differential event rate of nuclear recoils with a recoil energy of @xmath63 , occurring in a detector due to an elastic collision between a target nucleus of mass @xmath64 and a dm particle of mass @xmath65 , can be expressed as @xmath66 where a represents the mass number of the target nucleus and @xmath67 is the number of target nucleus per unit mass in the detector , @xmath68 is the local halo dm density , @xmath69 is the local dm velocity distribution evaluated in the galactic rest frame with @xmath70 , @xmath71 is earth velocity in the galactic rest frame , and @xmath72 is the galactic escape velocity  @xcite . @xmath73 and @xmath74 denote , respectively , the spin - independent differential dm - proton scattering cross - section and the effective form factors , and they will be detailed below . the @xmath75 is the momentum transfer between the dm particle and the struck nucleus . the lower limit @xmath76 of the integration in eq . ( [ dre ] ) corresponds to the smallest dm velocity that can give a recoil energy of @xmath63 , i.e. , @xmath77 where @xmath78 is the dm - nucleus reduced mass . the predicted normalized recoil spectrum can thus be expressed as @xmath79 where @xmath80 is the natural abundance of the isotope with mass number a for the target element . the total event rate @xmath81 expected in an energy range @xmath82 $ ] in a detector consisting of compound targets with finite detector energy resolution , usually characterized by the unit of one event @xmath83 , can be written as  @xcite @xmath84 where @xmath85 and @xmath86 are the detector acceptance and the detector response function , respectively . the @xmath63 integration region is determined by the cuts in the experiments . it is convenient , and usually adequate , to describe the matter distributions of a finite nucleus by a form factor , @xmath87 , in analyzing dm direct detection experiments  @xcite . in the first - order born approximation , the form factors of spherical nuclei are the fourier transforms of the proton and neutron density distributions , i.e. , @xmath88 where the index p ( n ) denotes protons ( neutrons ) . the density distributions , @xmath89 , for spherical nuclei can be obtained by using a mean - field approach  @xcite , e.g. , the shf approach as we discussed above . we further define the so - called effective form factors in eq . ( [ dre ] ) in terms of @xmath90 as @xmath91 where z ( n ) is the proton ( neutron ) number of the target nucleus , a = z+n , and @xmath92 is the isospin - violation factor with @xmath93 and @xmath94 denoting the effective coupling of dm to neutrons and protons , respectively . the spin - independent differential dm - proton scattering cross - section can then be expressed as  @xcite @xmath95 where @xmath96 is the reduced dm - proton mass , and @xmath97 is the elastic dm - proton cross - section at zero momentum transfer ( @xmath98 ) . integrating the differential event rate in eq . ( [ dre ] ) with respect to @xmath63 , i.e. , eq . ( [ ter ] ) , one can relate the expected number of scattering events in a direct detection experiment with the dm - proton cross - section @xmath97 . furthermore , by adopting , e.g. , a profile likelihood analysis proposed by the xenon100 collaboration  @xcite or the maximum gap method introduced by yellin  @xcite _ et al . _ , one can constrain the elastic dm - proton cross - section @xmath97 in the @xmath65-@xmath99 plane with specific confidence level .    however , in experiments usually reported is the dm - nucleon cross section @xmath100 , rather than @xmath97 . @xmath100 is an effective cross section by assuming that the dm particle couples with protons and neutrons equally ( i.e. , @xmath101 ) and the protons and neutrons have identical form factor which is in practice commonly taken as the empirical charge form factor parameterized by helm  @xcite . obviously , the dm - proton cross section @xmath97 is generally different from the dm - nucleon cross section @xmath100 , but they can be related to each other by the so - called degradation factor @xmath102 defined as  @xcite @xmath103 with @xmath104 in eq . ( [ degra ] ) , @xmath105 is the helm form factor which is obtained by using a charge density with a gaussian surface distribution and has the following simple analytical expression @xmath106 where @xmath107 is the first - order spherical bessel function , @xmath108 is an effective nuclear radius and @xmath109 is a measure of the nuclear skin thickness . for @xmath108 and @xmath109 , the following standard parametrization has been used  @xcite @xmath110 in eq . ( [ degra ] ) , it has been assumed that the variations of @xmath111 and form factors do not change the shape of recoil spectra @xmath112 in eq . ( [ spec ] ) . in this case , varying @xmath111 and form factors does not change the event distributions but simply modifies the expected number of total events , @xmath81 in eq . ( [ ter ] ) , for a given cross section . in general cases with different @xmath111 and nuclear form factors , since @xmath113 , we can simply rescale the corresponding cross sections according to the degradation factor as @xmath114 to keep the likelihood function ( the maximum gap ) in a profile likelihood analysis ( the maximum gap method ) unchanged . however , if the recoil spectrum changes significantly , the degradation factor in eq . ( [ degra ] ) should be modified according to the following identity @xmath115 where @xmath116 denotes the degradation factor defined in eq . ( [ degra ] ) , @xmath117 and @xmath118 are the number of expected events with @xmath119 and @xmath120 being the upper limits of the dm - proton and dm - nucleon cross sections with certain confidence level . we estimate @xmath117 and @xmath118 by using the `` maximum gap method ''  @xcite in this work . in this way , one can thus study the isospin - violating effects and the form factor effects simultaneously in analyzing the dm signals . as mentioned above , the neutron skin thickness of heavy nuclei is uniquely determined by the symmetry energy density slope parameter @xmath20 at a subsaturation cross density @xmath121  @xcite . therefore , the @xmath20 parameter controls the relative difference between the neutron and proton distributions ( and thus the neutron and proton form factors ) in the nuclei . in the present work , we investigate the form factor effects by varying the @xmath20 parameter in the shf calculations to fit the model - independent result of @xmath122 for @xmath0pb from the recent prex experiment at jlab . xe vs. that for @xmath0pb predicted by mean - field models with different interactions . constraints set by prex measurement  @xcite are also shown.,width=321 ]    shown in fig . [ nsts ] is the @xmath19 of @xmath123 as a function of that of @xmath15 in shf calculations with the skyrme interaction msl1  @xcite by varying @xmath20 while keeping the other 8 macroscopic quantities and the spin - orbit coupling constant @xmath51 fixed at their default values in the msl1 interaction , namely , @xmath124 @xmath22 , @xmath125 mev , the incompressibility @xmath126 mev , the isoscalar effective mass @xmath127 , the isovector effective mass @xmath128 , @xmath129 mev , @xmath130 mev@xmath131@xmath132 , @xmath133 mev@xmath131@xmath132 , and @xmath134 mev@xmath131@xmath132 . the msl1 interaction has been obtained by fitting a number of experimental data of finite nuclei , including the binding energy , the charge rms radius , the neutron @xmath135 energy level splitting in @xmath0pb , isotope binding energy difference , and neutron skin data of sn isotopes . for comparison , we also include in fig . [ nsts ] the results from shf calculations with @xmath136 other skyrme interactions ( bsk1 , bsk4 , bsk5 , bsk7 , bsk10 , bsk14 , bs15 , dutta , e - fit , esigma - fit , gsigma - fit , kde , kde0 , msl0 , ratp , rsigma - fit , sgi , sgii , sk255 , sk272 , ska , ski1 , ski2 , ski5 , skm , skmp , skm * , skt1 , skt4 , skt5 , skt6 , skz0 , skz1 , skz2 , skz3 , skz4 , skz-1 , sly4 , sly5 , sly9 , z - fit , zsigma - fit , and zsigmas - fit ) and rmf calculations with @xmath137 different interaction parameter sets ( fsu , iufsu , tm1 , pk1 , nl3 , and nl1 ) . the references of these skyrme and rmf interactions can be found in refs .  @xcite . we have selected these interactions in such a way that their @xmath23 values scatter in large region and are not close to each other . one can see clearly a nice model - independent linear correlation between the neutron skin thicknesses of @xmath15 and @xmath123 within the non - relativistic and relativistic models with different interactions . in addition , although there are few interactions that predict a @xmath19 of @xmath15 larger than @xmath138 fm measured by the model - independent prex experiment , one can easily obtain a large @xmath19 of @xmath15 by increasing the @xmath20 value in the msl1 interaction . in particular , we find a value of @xmath139 mev in the msl1 interaction predicts @xmath140 fm for @xmath0pb and @xmath141 fm for @xmath123 . therefore , we use @xmath142 mev , @xmath143 mev and @xmath144 mev in the msl1 interaction , denoted as lc40 , lc92 and lc134 , respectively , to study the symmetry energy and form factor effects in the present work . we would like to point out that the large neutron skin thickness for @xmath0pb measured by the prex experiment can also be explained within the rmf model  @xcite . .,width=510 ]    to test the skyrme interaction parameter sets lc40 , lc92 , and lc134 , we calculate the binding energies and charge rms radii for a number of closed - shell or semi - closed - shell nuclei , i.e. , @xmath145o , @xmath146ca , @xmath147ca , @xmath148ni , @xmath149zr , @xmath150sn , @xmath151sn , and @xmath0pb , as well as the nuclei @xmath152si , @xmath153ge , and @xmath151xe . the elements si , ge and xe are widely used as targets in the dm direct detection experiments . for si and ge , their natural abundances are mostly dominated by their isotopes of @xmath152si and @xmath153ge , respectively . for xe , however , it has several isotopes with comparable natural abundances and here we just select one of them , i.e. @xmath151xe , as an example to show its properties . [ ebrc ] shows the relative deviation of the binding energies and charge rms radii of these nuclei from those measured in experiments  @xcite . it is seen that the interactions lc40 and lc92 can describe the experimental data very well ( the deviations are within about @xmath154 ) except for the light nucleus @xmath145o for which the deviation of charge rms radius is about @xmath155 . it is interesting to see that , the interaction lc134 , which predicts a very strong density dependence of the symmetry energy at @xmath156 and gives a very large neutron skin thickness of @xmath157 fm for @xmath15 ( the upper limit of the prex measurement ) , still can give a reasonable description of the experimental data ( the deviations are within about @xmath158 as indicated by bands in fig . [ ebrc ] ) except the nucleus @xmath151sn for which the deviation of charge rms radius ( binding energy ) reaches about @xmath159 ( @xmath160 ) and for the light nucleus @xmath145o the deviation of charge rms radius is still about @xmath155 . these results are remarkable as lc40 , lc92 and lc134 are not obtained from fitting measured binding energies and charge rms radii of finite nuclei as in usual skyrme parametrization . it should be pointed out that our main motivation for introducing the interactions lc40 , lc92 and lc134 is not to construct new skyrme interaction parameter sets to describe data , but to use them as references to study the form factor effects in the following .     for @xmath152si ( left panels ) and @xmath151xe ( right panels ) obtained from shf with lc40 , lc92 and lc134 . the momentum transfer @xmath161 , corresponding to a value of @xmath162 kev , is indicated by dash - dot - dotted lines.,width=510 ]    to investigate the symmetry energy effects on the proton and neutron form factors in finite nuclei , we plot in fig . [ fnp ] the proton ( upper panels ) and neutron ( lower panels ) form factors as functions of momentum transfer @xmath163 for @xmath152si ( left panels ) and @xmath151xe ( right panels ) from shf calculations with lc40 , lc92 and lc134 . for comparison , the results are divided by the corresponding empirical helm form factors . the proton and neutron form factors of other isotopes or elements can be obtained from a generalized helm - like empirical parametrization given in the appendix where the proton and neutron form factors are distinguished and parameterized as functions of the neutron skin thickness of @xmath0pb by fitting the results obtained from shf calculations with lc40 , lc92 and lc134 . in addition , also indicated in fig . [ fnp ] is the momentum transfer @xmath161 corresponding to a value of @xmath162 kev which represents the possible maximum cutoff among various experiments . however , we would like to emphasize that the upper limits of the recoil energy @xmath63 usually vary significantly in different dm direct detection experiments , and here we just select a relatively large recoil energy of @xmath8 kev as an example to see the possible symmetry energy effects on the proton and neutron form factors in different nuclei . in calculations of the scattering cross sections in eq . ( [ degra ] ) for each experiment , one must use the recoil energy region given by each experiment . it is seen from fig . [ fnp ] that , for the nucleus @xmath152si with equal proton and neutron numbers , the symmetry energy effect is tiny and the shf calculations with lc40 , lc92 and lc134 predict almost the same form factors for neutrons and protons , which are further in good agreement with the empirical helm form factor for @xmath164 . on the other hand , for the neutron - rich nucleus @xmath151xe , one can see a clear symmetry energy effect on the form factors , namely , a larger @xmath23 shifts the neutron form factor to lower @xmath163 values while shifts slightly the proton form factor to higher @xmath163 values with the empirical helm form factor in between , leading to an isospin splitting between the proton and neutron form factors with respect to the @xmath163 value . these features actually reflect the symmetry energy effects on the neutron skin thickness , namely , increasing the value of @xmath23 increases the rms radius of neutrons and reduces slightly the rms radius of protons , and thus leads to a larger neutron skin thickness as shown in fig . .     for si ( left panels ) and xe ( right panels ) with @xmath165 ( upper panels ) and @xmath166 ( lower panels ) obtained from shf with lc40 , lc92 and lc134 . the empirical helm form factor is also included for comparison.,width=510 ]    since the effective form factor ( eq . ( [ eff ] ) ) should be used to analyze the dm signals for general cases , it is expected that the isospin splitting between the proton and neutron form factors would cause significant effects on the effective form factor , especially for a negative @xmath111 which could cause strongly destructive interference of dm scattering with protons and neutrons . on the other hand , the presence of different isotopes with comparable abundances in some detector target elements ( e.g. , xenon ) suggests that it is better to study the averaged effective form factors defined as @xmath167 this is illustrated in fig . [ ffs ] where the averaged effective form factor is plotted as a function of @xmath63 for si ( left panels ) and xe ( right panels ) from shf calculations with lc40 , lc92 and lc134 . the results from the empirical helm form factor with @xmath168 are also included for comparison . in fig . [ ffs ] , two cases for the dm are considered : one is for the standard isospin - invariant dm with @xmath165 ( upper panels ) and the other is for the ivdm with @xmath166 ( lower panels ) . it should be pointed out that @xmath169 , firstly suggested by feng _ et al . _ @xcite , leads to nearly complete destructive interference of the scattering amplitudes for dm - proton and dm - neutron collisions in xe target and gives maximum suppression of the sensitivity for xe - based detectors if the empirical helm form factor is adopted . for such kind of ivdm , the confidence region of cdms - ii(si ) and the xenon100 exclusion contours are consistent with each other  @xcite . as will be shown later , for ivdm with @xmath169 , both the results from recent lux  @xcite and the most recent supercdms(ge )  @xcite are also consistent with the confidence region of cdms - ii(si ) . theoretically , such negative values of @xmath111 can arise , e.g. , in models with a new light neutral gauge boson @xmath170  @xcite . one can see from fig . [ ffs ] that the negative @xmath111 ( i.e. , @xmath171 ) indeed causes strongly destructive interference of dm scattering with protons and neutrons for both si and xe targets and thus the averaged effective form factor with @xmath169 is strongly suppressed compared to its value with @xmath172 . furthermore , it is interesting to see from fig . [ ffs ] that for @xmath172 , the averaged effective form factors of both si and xe targets from shf calculations with lc40 , lc92 and lc134 are in a very good agreement with the corresponding results from the empirical helm form factors , indicating that there are essentially no symmetry energy effects . however , for @xmath169 , one can see that the averaged effective form factors of xe target display very different behaviors for different distributions from shf calculations with lc40 , lc92 and lc134 as well as that from the empirical charge distributions , although the averaged effective form factors of si target are essentially the same for these different distributions . these features indicate that the symmetry energy effect is tiny for the averaged effective form factor of si target but is very strong for that of xe target in the case of @xmath169 . this is due to the fact that for a negative @xmath111 , the scattering amplitudes for dm - proton and dm - neutron collisions may interfere destructively , and the dm particle is almost completely decoupled from the xe isotopes , i.e. @xmath173 ( see panel ( d ) of fig .  [ ffs ] ) , especially for @xmath174  @xcite . in this case , a small difference between @xmath175 and @xmath176 may lead to a significant change for the effective form factor . as a result , the form factor of xe exhibits a very strong symmetry energy effect for @xmath169 , namely , a larger @xmath20 value leads to a larger form factor as shown in panel ( d ) of fig . in particular , one can see that the averaged effective form factors of xe target from shf calculations with lc40 , lc92 and lc134 can be significantly larger than that from the empirical helm form factor in a large range of recoil energy . this implies that the symmetry energy with a larger @xmath20 ( and thus a larger neutron skin thickness ) can give a more significant enhancement for the effective dm - nucleus form factors in eq . ( [ dre ] ) when the ivdm with @xmath169 is considered . therefore , for @xmath169 , one can expect a significant enhancement for the sensitivity of the xe - based detectors to the dm - proton cross sections by taking into account the form factor effects due to different symmetry energies compared with the corresponding result using the empirical helm form factor in analyzing the experimental data .       confidence level ( c.l . ) limits from xenon100 , lux and supercdms(ge ) together with @xmath177 c.l . favored regions from cdms - ii(si ) in @xmath65-@xmath97 plane for @xmath172 ( left panel ) and @xmath169 ( right panel ) with form factors from the empirical charge distributions ( helm ) and shf calculations with msl1 by varying @xmath23 from @xmath178 mev to @xmath144 mev.,width=510 ]    to see the effects of the form factor and the isospin - violation factor @xmath111 on the extraction of elastic spin - independent dm - proton cross sections @xmath97 in dm direct detection experiments , we show in fig . [ css ] the results of the @xmath177 confidence level ( c.l . ) limits from xenon100  @xcite , lux  @xcite and supercdms(ge )  @xcite along with the @xmath177 c.l . favored regions from cdms - ii(si )  @xcite in @xmath65-@xmath97 plane for @xmath172 ( left panel ) and @xmath169 ( right panel ) with form factors from the empirical charge distributions and shf calculations with msl1 by varying @xmath23 from @xmath178 mev to @xmath144 mev . for the results from supercdms(ge ) , only the mass region less than @xmath179 gev is shown since no data are available for @xmath180 gev  @xcite . in these analyses , we use the standard astrophysical parameters in the standard halo model , namely , a maxwell - boltzmann distribution for @xmath181 with @xmath182 and galactic escape velocity @xmath183 , and a dm density of @xmath184  @xcite .    in fig . [ css ] , the form factor effects due to the symmetry energy and the effects from @xmath111 are considered through dividing the experimental results by the degradation factor @xmath102 as discussed earlier . in particular , for xenon100 , we assume that the energy resolution is dominated by poisson fluctuations in the number of photoelectrons ( pe ) as in ref . we use s1 efficiency obtained from fig . 1 of ref .  @xcite and the scintillation efficiency @xmath185 from fig . 1 of ref .  @xcite , respectively , and integrate the differential rate over s1 from @xmath2 to @xmath179 pe to get the total rate . for the lux analysis which is similar with xenon100 s , the signal region is limited within @xmath186 to @xmath179 pe and the efficiencies are taken from fig . 1 of ref .  @xcite . for cdms - ii(si ) , we use the acceptance from fig . 1 of ref . @xcite , and choose an energy interval between @xmath187 and @xmath188 in the analysis . for supercdms(ge ) , with the efficiency from fig . 1 of ref .  @xcite , we use the signal region in the range @xmath189 in the calculations as in ref . @xcite . furthermore , since xenon100 reported @xmath186 candidate events in the low recoil energy region , we calculate the degradation factor according to eq . ( [ degra2 ] ) instead of eq . ( [ degra ] ) to take into account the effects due to shape changing of the recoil spectrum as we will show later . similarly , for lux , we follow the analysis in ref . @xcite and calculate the degradation factor according to eq . ( [ degra2 ] ) by using one observed event ( with @xmath190 pe ) . for cdms - ii(si ) and supercdms(ge ) that focus on the light dm , we have checked that the shape changing of their recoil spectra can be negligible and eq . ( [ degra ] ) is thus applied to evaluate the degradation factor .    for the isospin - invariant dm with @xmath172 , one can see from fig . [ css ] that the cdms - ii(si ) results are in some tension with the upper limit placed by the xenon100 experiment , and they become in even strong disagreement with the upper limits set by the recent lux  @xcite experiment and the more recent supercdms(ge )  @xcite experiment . however , for the ivdm with @xmath169 , one can see that the tension between cdms - ii(si ) and xenon100 is essentially disappeared as found by frandsen  _ et al . _ furthermore , it is remarkable to see from the right panel of fig . [ css ] that the disagreement of the cdms - ii(si ) results with lux and supercdms(ge ) can also be largely ameliorated within the framework of ivdm with @xmath169 ( we note that a similar conclusion is also obtained in ref . therefore , these results indicate that ivdm indeed provides a very promising mechanism to reconcile the tension among various experiments . it should be mentioned here that although @xmath169 leads to a maximum suppression for the sensitivity of xe - based detector for the empirical helm form factor , it does not for the more realistic shf form factors as will be shown later . furthermore , it is interesting to see that , although there are essentially no form factor effects on the extracted spin - independent dm - proton cross sections @xmath97 in different experiments in the case of isospin - invariant dm with @xmath172 , the form factor effects can significantly affect the extraction of @xmath97 for the xe - based experiments ( xenon100 and lux ) in the case of ivdm with @xmath169 . therefore , our results indicate that , for isospin - invariant dm , the widely used empirical helm form factor is well grounded and this is consistent with our previous discussions about the form factor effects . on the other hand , for ivdm with @xmath169 , using the form factors obtained from shf calculations in msl1 with varied @xmath23 values generally increases the sensitivity of the xe - based detectors ( xenon100 and lux ) compared with using the empirical helm form factor , and a larger @xmath23 value , which leads to a larger neutron skin thickness in xe isotopes , generally leads to a stronger sensitivity of the xe - based detectors . this form factor effect becomes more pronounced with the increment of the dm mass @xmath65 . while the form factor effect is about @xmath191 for the dm mass of @xmath192 gev , it can become very significant for the dm with mass above a few tens gev . for @xmath193 gev , which roughly corresponds to the dm mass upper limit of cdms - ii(si ) bounds , the sensitivities of xenon100 and lux can be enhanced by a factor of @xmath2 using the form factor with @xmath194 mev compared with that using the empirical helm form factor . in particular , the sensitivities can have a factor of more than @xmath46 improvement for @xmath195 gev where many constraints for supersymmetric wimps have been put by the data from lhc  @xcite . on the other hand , the relative variation between the results extracted from the helm form factor and the shf calculations is small ( @xmath196 ) for supercdms(ge ) and it can be negligible ( @xmath197 ) for cdms - ii(si ) for which the results from the shf calculations are not shown in fig . [ css ] . with @xmath169 for @xmath198 ( panel ( a ) ) , @xmath199 ( panel ( b ) ) and @xmath200 ( panel ( c ) ) , respectively . cases for different form factors from the empirical helm charge distribution ( black solid lines ) and shf calculations with msl1 by varying @xmath23 from @xmath178 mev to @xmath144 mev are shown in the figure . for comparison , spectra calculated with helm form factors for isospin - invariant cases ( @xmath201 ) are also shown.,width=529 ]    the dm mass dependence of the form factor effects observed in xenon100 and lux for @xmath169 can be understood as a result of the shape changing of the expected recoil spectra for varied @xmath65 due to the enhancement effects of the form factors when we take the symmetry energy effects into account as discussed in the previous subsection . shown in fig . [ spe ] is the normalized dm - xenon recoil spectrum @xmath112 for @xmath202 , @xmath203 and @xmath204 with @xmath169 and various form factors from shf calculations with lc40 , lc92 and lc134 as well as the empirical helm form factor . we note that for @xmath172 , the @xmath112 essentially does not depend on the variation of the form factors and the results are almost the same as that from the empirical helm form factor with @xmath169 . it is very interesting to see from fig . [ spe ] that for the case with @xmath169 , although the form factor effect essentially does not change or only slightly modifies the predicted normalized dm - xenon recoil spectra for light dm ( see , e.g. , figs . [ spe ] ( a ) and ( b ) ) , it can drastically change the shape of the recoil spectra for heavy dm ( see , e.g. , fig . [ spe ] ( c ) ) . for both xenon100 and lux , we have considered this spectrum shape changing effect in eq . ( [ degra2 ] ) by multiplying eq . ( [ degra ] ) by another term @xmath205 for the degradation factor , which will result in , e.g. , @xmath206 ( @xmath207 ) enhancement of the upper limit of @xmath97 from xenon100 ( lux ) relative to the results obtained from eq . ( [ degra ] ) for the lc134 case with @xmath208 . more specifically , one can see from fig . [ spe ] that for heavier dm , the expected dm - xenon recoil spectra become harder and thus the effective form factors are typically probed at higher momentum transfer , and this will significantly enhance the effects caused by different form factors ( see , e.g. , fig . [ ffs ] ( d ) ) adopted for the @xmath63 integration in the numerator of r.h.s of eq . ( [ degra ] ) , leading to the observed dm mass dependence of the form factor effects in xenon100 and lux for @xmath169 . it is also very interesting to see that for the case with @xmath169 , the bounds of xenon100 and lux overlap with each other in the mass region of @xmath209 gev although the lux gives much more stringent limit on the spin - independent elastic dm - proton scattering cross section than xenon100 in the lower mass region as shown in the right panel of fig . this can also be explained as a result of the shape changing of the expected recoil spectrum for varied @xmath65 as shown in fig . [ spe ] . particularly , fig . [ spe ] suggests that a significant fraction of the events is expected at larger recoil energies for heavier dm , and this is especially the case for the ivdm with @xmath169 and shf form factors with lc92 and lc134 where one can see the non - standard recoil spectra peaked around @xmath210 kev . since lux adopts a relatively small energy range , i.e. @xmath211 kev  @xcite , compared with xenon100 where @xmath212 kev is adopted  @xcite , it thus exhibits relatively weaker sensitivity than xenon100 for heavier dm . as a matter of fact , even for the isospin - invariant case where the form factor effects are negligible as shown in fig . [ css ]  ( a ) , the sensitivity of lux becomes relatively weaker with the increment of @xmath65 compared with that of xenon100 . it should be mentioned that the small form factor effect observed for supercdms(ge ) in fig . [ css ]  ( b ) is also partially due to the smaller threshold energy cut , i.e. , @xmath213 kev  @xcite . our results imply that the non - standard recoil spectra shown in fig . [ spe ]  ( c ) may provide an extremely important experimental evidence for non - standard dark matter interactions and they also give valuable information on the direct detection of ivdm .      in above calculations for the case of ivdm , we have fixed the isospin - violation factor at @xmath169 . as mentioned earlier , for the specific ivdm with @xmath169 , the sensitivity of xe - based detectors is maximumly suppressed when the empirical helm form factor is used and the tension among various experiments using different target elements has been shown to be largely ameliorated . for more realistic shf form factors , it is interesting to see how the form factor effects change with more general values of the isospin - violation factor @xmath111 . following the method by feng  _ et al . _ @xcite , we study in the following how the form factor effects change for various isospin - violation factor @xmath111 . shown in fig . [ xepho ] is the degradation factor @xmath102 as a function of the isospin - violation factor @xmath111 obtained by using the empirical helm form factors as well as the form factors from shf calculations with lc40 , lc92 and lc134 for @xmath214 gev ( panel ( a ) ) and @xmath215 gev ( panel ( b ) ) . since both xenon100 and lux are xe - based experiments and their results have similar dependence on the @xmath111 and the symmetry energy , we thus only consider cdms - ii(si ) , supercdms(ge ) and xenon100 in fig . [ xepho ]  ( a ) for @xmath214 gev . in addition , for @xmath215 gev , the form factor effect is negligible for cdms - ii(si ) and relatively small for supercdms(ge ) which adopts a small threshold energy cut and has been designed especially for light dm searching , we thus only show the results of xenon100 ( panel ( b ) ) .     as a function of the isospin - violation factor @xmath111 obtained by using the empirical helm form factors as well as the form factors from shf calculations with lc40 , lc92 and lc134 for @xmath214 gev ( panel ( a ) ) and @xmath215 gev ( panel ( b ) ) . for @xmath214 gev , the results of cdms - ii(si ) , supercdms(ge ) and xenon100 are shown while only the results of xenon100 are shown for @xmath215 gev.,width=510 ]    it is seen from fig . [ xepho ] that generally the sensitivity of each detector is significantly suppressed for a negative @xmath111 within the range of @xmath216 due to destructive interference between dm - neutron and dm - proton scattering amplitudes as pointed out in ref . however , the @xmath102 reaches its minimum value ( where the detector is most insensitive ) at different values of @xmath111 for different detectors due to the different compositions and distributions of neutrons and protons in each target element of the detector . in particular , one can see that the corresponding value of @xmath111 which leads to a minimum value of @xmath102 , denoted as @xmath217 , is about @xmath218 , @xmath219 and @xmath220 for cdms - ii(si ) , supercdms(ge ) and xenon100 , respectively , for the case of @xmath214 gev when the empirical helm form factor is used . for @xmath214 gev , the form factor effect is very small for cdms - ii(si ) and the largest relative variation with respective to the results from the empirical helm form factor when @xmath23 varies from @xmath221 mev to @xmath144 mev is @xmath222 at @xmath223 . the largest relative variation due to the form factor effect is @xmath224 at @xmath225 for supercdms(ge ) , and it becomes @xmath226 at @xmath227 for xenon100 . as a matter of fact , a clear shift of curves for the xe - based detectors calculated with different form factors has been already seen even for such a small dm mass in fig . [ xepho ]  ( a ) . these features indicate that the form factor effect can be already pronounced for xenon100 even for a small mass dm with @xmath214 gev . for xenon100 , one can see from fig . [ xepho](b ) that the form factor effect becomes very strong for @xmath215 gev . in particular , a larger @xmath23 value gives significantly stronger sensitivity of the detectors and the largest relative variation due to the form factor effect is @xmath228 ( almost one order ) at @xmath229 .    furthermore , one can see from fig . [ xepho ] that the value of @xmath217 also depends significantly on the form factor , and it is @xmath230 , @xmath231 and @xmath232 for the shf form factor with @xmath233 , @xmath143 and @xmath144 mev , respectively , smaller than @xmath169 for the empirical helm form factor . in addition , it is very interesting to see that , while using the shf form factors enhances the sensitivity of the xe - based detector relative to the case with helm form factor for @xmath234 , it suppresses the sensitivity for @xmath235 . these features indicate that the form factor effects can either enhance or suppress the sensitivity of the detector relative to the helm s case depending on the specific value adopted for @xmath111 . it is also very interesting to see that one can not find a @xmath111 value such that @xmath236 , leading to zero sensitivity for scattering off the elements , and this is mainly due to the fact that the elements have multiple isotopes and completely destructive interference can not be simultaneously achieved for all isotopes  @xcite . we would like to point out that even for elements with only one naturally abundant isotope , one still can not find a @xmath111 value such that @xmath236 since the neutron form factor is generally different from the proton form factor in a nucleus . the above results indicate that , the form factor effect can not only significantly change the sensitivities of the detectors , but also affect the interaction behaviors of the target element with the isospin - violating dm in various isospin - violating regions we are interested in for dm direct detection experiments . in future when different direct detection experiments using various target elements had specify the dm signals precisely , one may determine simultaneously the isospin - violation factor @xmath111 and the form factor ( and thus the neutron skin thickness as well as the symmetry energy ) , especially if the ivdm particle has a larger mass of above tens of gev . conversely , our results indicate that precise knowledge on the symmetry energy or the neutron skin thickness ( and thus the neutron and proton form factors ) of @xmath0pb is of critical importance for the direct detection of ivdm using xe - based detectors . in particular , the future experiment prex - ii  @xcite at jlab is expected to improve significantly the measurement accuracy of the neutron skin thickness for @xmath0pb and thus could make important contribution to this issue . in the present work , we have shown that isospin - violating dark matter ( ivdm ) indeed provides a possible mechanism to ameliorate the tension among recent direct detection experiments , including cdms - ii(si ) , xenon100 , lux , and supercdms(ge ) . for ivdm , we have demonstrated that the results of the dm direct detection experiments based on neutron - rich target nuclei , e.g. , xe - based detector , may strongly depend on the density slope @xmath23 of the symmetry energy at a subsaturation cross density @xmath38@xmath22 , which is presently largely unknown and uniquely determines the neutron skin thickness and thus the relative difference of neutron and proton form factors of the target nuclei .    in particular , using the proton and neutron form factors obtained from skyrme - hartree - fock calculations by varying the @xmath23 within the uncertainty region set by the latest model - independent measurement of the neutron skin thickness from prex experiment at jlab , we have found that although the form factor effects on the extracted bounds on dm - proton cross sections are negligible in the direct detection for isospin - invariant dm , they could become critically important in the detection for ivdm . especially , for ivdm with neutron - to - proton coupling ratio fixed to @xmath237 in the mass region constrained by cmds - ii(si ) , the form factor effect may enhance the sensitivity of xe - based detectors ( e.g. , xenon100 and lux ) to the dm - proton cross section by a factor of @xmath2 , compared with the results using the empirical helm nuclear form factors extracted from charge distributions . this form factor effect can even enhance the sensitivity by more than a factor of @xmath46 for such kind of ivdm with mass larger than @xmath3 gev . furthermore , we have found that the form factor effect can significantly modify the recoil spectrum of xe - based detectors for heavy ivdm with @xmath237 . we have also studied how the form factor effects change with the variation of @xmath238 and found that the @xmath238 value maximumly suppressing the sensitivity of the detector may depend on the form factor , and the form factor effects can either enhance or suppress the sensitivity of the detector relative to the helm s case depending on the specific value adopted for @xmath238 . our results imply that the precise determination of the symmetry energy or the neutron skin thickness ( and thus the neutron and proton form factors ) of @xmath0pb is extremely useful for the direct detection of ivdm based on detectors with neutron - rich targets ( e.g. , xenon ) . the authors would like to thank fei gao for useful discussions . this work was supported in part by the nnsf of china under grant nos . 11275125 and 11135011 , the shanghai rising - star program under grant no . 11qh1401100 , the `` shu guang '' project supported by shanghai municipal education commission and shanghai education development foundation , the program for professor of special appointment ( eastern scholar ) at shanghai institutions of higher learning , and the science and technology commission of shanghai municipality ( 11dz2260700 ) . for si ( upper panels ) , ge ( middle panels ) and xe ( lower panels ) with @xmath165 ( left panels ) and @xmath166 ( right panels ) obtained from shf with lc40 , lc92 and lc134 ( thick lines ) and from the empirical parametrization for proton and neutron form factors in eq . ( [ fitting ] ) ( thin lines ) with @xmath239 , @xmath240 and @xmath241 . the empirical helm form factor is also included for comparison.,width=12 ]    shown in fig . [ fnp ] are the proton and neutron form factors as functions of momentum transfer for @xmath152si and @xmath151xe from shf calculations with lc40 , lc92 and lc134 . the results for other isotopes of si and xe as well as ge isotopes can also be calculated in the same approach . thanks to the approximately linear relationship between the symmetry energy density slope parameter @xmath20 and the neutron skin thickness @xmath242 of @xmath0pb , it is possible to approximate the proton and neutron form factors obtained from shf calculations with varied @xmath20 by analytical parameterizations with explicit dependence on the @xmath242 of @xmath0pb . in particular , here we construct the analytical expressions for the proton and neutron form factors separately according to helm - like form factor with an additional term to consider isospin dependence and the symmetry energy effects ( and thus @xmath242 of @xmath0pb ) , i.e. , @xmath243 where the effective proton ( neutron ) distribution radius is parameterized as @xmath244 with @xmath245 compared with the parametrization of helm form factor ( i.e. , eq . ( [ helm ] ) ) , the parametrization here includes an additional term , namely , the last term in eq . ( [ cp ] ) and eq . ( [ cn ] ) . in addition , our parametrization here distinguishes protons and neutrons , and depends on explicitly the proton and neutron numbers of the nucleus and the @xmath242 of @xmath0pb .        by fitting a number of proton and neutron form factors calculated from shf approach with varied symmetry energies , including the 9 stable xe isotopes , 5 stable ge isotopes and 3 stable si isotopes with the symmetry energy density slope parameter @xmath20 varying from @xmath246 to @xmath247 , we obtain all parameters in the above expressions and they are listed in table  [ table ] . in order to test the quality of the above parametrization , we compare in fig . [ ffs02 ] the averaged effective form factors as functions of @xmath63 for si , ge and xe using the above parameterized proton and neutron form factors with the results obtained from the shf approach . it is very interesting to see that one can indeed simply use just one set of parameters given in table  [ table ] to well describe the proton and neutron form factors in different nuclei in a large mass region with very different values of the symmetry energy density slope parameter @xmath20 ( i.e. , @xmath242 of @xmath0pb ) . * cogent * collaboration , c. e. aalseth et al . , _ results from a search for light - mass dark matter with a p - type point contact germanium detector _ , _ phys * 106 * ( 2011 ) 131301 [ arxiv:1002.4703 ] .                    c. savage , g. gelmini , p. gondolo , and k. freese , _ xenon10/100 dark matter constraints in comparison with cogent and dama : examining the @xmath248 dependence _ , _ phys . _ * d 83 * ( 2011 ) 055002 [ arxiv:1006.0972 ] . g. duda , a. kemper , and p. gondolo , _ model independent form factors for spin independent neutralino - nucleon scattering from elastic electron scattering data _ , _ jcap _ * 0704 * ( 2007 ) 012 [ hep - ph/0608035 ] . chen , y.z . luo , l. li , h. shen , and x.q . li , _ determining nuclear form factor for detection of dark matter in relativistic mean field theory _ , _ commun . * 55 * ( 2011 ) 1059 [ arxiv:1101.3049 ] . a. ong , j. c. berengut , and v. v. flambaum , _ the effect of spin - orbit nuclear charge density corrections due to the anomalous magnetic moment on halonuclei _ , _ phys . _ * c 82 * ( 2010 ) 014320 [ arxiv:1006.5508 ] . brown , _ neutron radii in nuclei and the neutron equation of state _ , _ phys * 85 * ( 2000 ) 5296 ; s. typel and b.a . brown , _ neutron radii and the neutron equation of state in relativistic models _ , _ phys . rev . _ * c 64 * ( 2001 ) 027302 . z. zhang and l.w . chen , _ constraining the symmetry energy at subsaturation densities using isotope binding energy difference and neutron skin thickness _ , _ phys . * b 726 * ( 2012 ) 234 [ arxiv:1302.5327 ] . p. klpfel , p .- reinhard , t.j . brvenich , and j.a . maruhn , _ variations on a theme by skyrme : a systematic study of adjustments of model parameters _ , _ phys . _ * c 79 * ( 2009 ) 034310 [ arxiv:0804.3385 ] . f. le blanc , l. cabaret , e. cottereau , j.e . crawford , s. essabaa , j. genevey , r. horn , g. huber , j. lassen , j.k.p . lee , g. le scornet , j. lettry , j. obert , j. oms , a. ouchrif , j. pinard , h. ravn , b. roussiere , j. sauvage , and d. verney , _ charge - radius change and nuclear moments in the heavy tin isotopes from laser spectroscopy : charge radius of * * 132sn _ , _ phys . rev . _ * c 72 * ( 2005 ) 034305 .            c. strege et al . updated global fits of the cmssm including the latest lhc susy and higgs searches and xenon100 data _ , _ jcap _ * 1203 * ( 2012 ) 030 [ arxiv:1112.4192 ] ; a. fowlie , m. kazana , k. kowalska , s. munir , l. roszkowski , e. sessolo , s. trojanowski , and y .- l . tsai , _ the cmssm favoring new territories : the impact of new lhc limits and a 125 gev higgs _ , _ phys . rev . _ * d 86 * ( 2012 ) 075010 [ arxiv:1206.0264 ] ; o. buchmueller et al . , _ higgs and supersymmetry _ , _ eur . j. _ * c 72 * ( 2012 ) 2020 [ arxiv:1112.3564 ] .
Abstract:

Answer: isospin - violating dark matter ( ivdm ) provides a possible mechanism to ameliorate the tension among recent direct detection experiments . for ivdm , we demonstrate that the results of direct detection experiments based on neutron - rich target nuclei may depend strongly on the density dependence of the symmetry energy which is presently largely unknown and controls the neutron skin thickness that reflects the relative difference of neutron and proton form factors in the neutron - rich nuclei . in particular , using the neutron and proton form factors obtained from skyrme - hartree - fock calculations by varying the symmetry energy within the uncertainty region set by the latest model - independent measurement of the neutron skin thickness of @xmath0pb from prex experiment at jlab , we find that , for ivdm with neutron - to - proton coupling ratio fixed to @xmath1 , the form factor effect may enhance the sensitivity of xe - based detectors ( e.g. , xenon100 and lux ) to the dm - proton cross section by a factor of @xmath2 in the dm mass region constrained by cmds - ii(si ) and even by more than an order of magnitude for heavy dm with mass larger than @xmath3 gev , compared with the results using the empirical helm form factor . our results further indicate that the form factor effect can significantly modify the recoil spectrum of xe - based detectors for heavy ivdm with @xmath1 .