%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% This is a sample input file for your contribution to a
% multi-author book ("proceedings-like) to be published
% by Springer Verlag.
%
% Please use it as a template for your own input, and please
% follow the instructions for the formal editing of your
% manuscript as described in the file "1readme".
%
% Please send the Tex and figure files of your manuscript
% together with any additional style files as well as the
% PS file to the editor of your book.
%
% He or she will collect all contributions for the planned
% book, possibly compile them all in one go and pass the
% complete set of manuscripts on to Springer.
%
% Vers. 07/99
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%



%%RECOMMENDED%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
\documentclass[runningheads]{svmult}
%
\usepackage{makeidx}   % allows index generation
\usepackage{graphicx}  % standard LaTeX graphics tool
%                       % for including eps-figure files
\usepackage{subeqnar}  % subnumbers individual equations
%                       % within an array
\usepackage{multicol}  % used for the two-column index
\usepackage{cropmark} % cropmarks for pages without
%                       % pagenumbers - only needed when manuscript
%                       % is printed from paper and not from data
\usepackage{physprbb}  % centered layout of diverse elements, etc.
%\makeindex             % used for the subject index
%                       % please use the style sprmidx.sty with
%                       % your makeindex program
%
%%%upright Greek letters (example below: upright "mu")
\newcommand{\greeksym}[1]{{\usefont{U}{psy}{m}{n}#1}}
\newcommand{\umu}{\mbox{\greeksym{m}}}
\newcommand{\udelta}{\mbox{\greeksym{d}}}
\newcommand{\uDelta}{\mbox{\greeksym{D}}}
\newcommand{\uPi}{\mbox{\greeksym{P}}}
\newcommand{\newc}{\newcommand} 
%\newc{\ra}{\rightarrow} 
%\newc{\lra}{\leftrightarrow} 
\newc{\beq}{\begin{equation}} 
\newc{\eeq}{\end{equation}} 
\newc{\barr}{\begin{eqnarray}} 
\newc{\earr}{\end{eqnarray}} 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


%OPTIONAL%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
%\usepackage{amstex}   % useful for coding complex math
%\mathindent\parindent % needed in case "Amstex" is used
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%AUTHOR_STYLES_AND_DEFINITIONS%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%
%Please reduce your own definitions and macros to an absolute
%minimum since otherwise the editor will find it rather
%strenuous to compile all individual contributions to a
%single book file
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

\begin{document}
%
\title* {Theoretical Directional and Modulated Rates for Direct  SUSY Dark
 Matter Detection}
%
%
\toctitle{ Theoretical Directional and Modulated Rates
\protect\newline for Direct  SUSY Dark Matter Detection}
% allows explicit linebreak for the table of content
%
%
\titlerunning{Direct SUSY Dark Matter Detection}
% allows abbreviation of title, if the full title is too long
% to fit in the running head
%
\author{J.D. Vergados{1,2}}
%
\authorrunning{J.D. Vergados}
% if there are more than two authors,
% please abbreviate author list for running head
%
%
\institute{Physics Department, UNISA, Pretoria, South Africa.
\and Permanent address: University of Ioannina, Ioannina, GR 45110, Greece.
\\E-mail:Vergados@cc.uoi.gr} 
\maketitle              % typesets the title of the contribution

\begin{abstract}
Exotic dark matter together with the vacuum energy (cosmological constant)
seem to dominate in the flat Universe.  Thus direct dark matter detection 
is central to particle physics and cosmology. Supersymmetry
provides a natural dark matter candidate, the lightest
supersymmetric particle (LSP). Furthermore from the knowledge of the density
 and velocity distribution of the LSP, 
the quark substructure of the nucleon and the nuclear structure (form factor 
and/or spin response function), one is able to evaluate
the event rate for LSP-nucleus elastic scattering. The thus obtained event rates
are, however, very low.
 So it is imperative to exploit the two signatures of the reaction, namely
the modulation effect, i.e. the dependence of
the event rate on  the Earth's motion, and the directional asymmetry, i.e.
the dependence of the rate on the the relative angle between the direction
 of the recoiling
nucleus and the sun's velocity. These two signatures are studied in this paper
employing various velocity distributions and a supersymmetric model
with universal boundary conditions at large $tan\beta$.
\end{abstract}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\date{\today}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{Introduction}
The combined MAXIMA-1 BOOMERANG and COBE/DMR Cosmic Microwave Background (CMB)
observations imply that the Universe is flat \cite{flat01},
 $\Omega=1.11\pm0.07$, while the baryonic component is very small
$\Omega_b h^2=0.032^{+0.009}_{~0.008}$.
Furthermore exotic dark matter has become necessary
in order to close the Universe. In fact about a decade ago 
 the COBE data ~\cite{COBE} suggested that CDM (Cold Dark Matter)
dominates the Universe, $\Omega_{CDM}$ being at least $60\%$ ~\cite {GAW}.
 Subsequent evidence from two different teams,
the High-z Supernova Search Team \cite {HSST} and the
Supernova Cosmology Project  ~\cite {SPF} $^,$~\cite {SCP} 
changed this view suggesting that the Universe may be dominated by 
the  cosmological constant $\Lambda$ or dark energy. In other words one
roughly finds a baryonic component $\Omega_B=0.1$ along with the exotic 
components $\Omega _{CDM}= 0.3$ and $\Omega _{\Lambda}= 0.6$.
In  a more detailed analysis by Einasto \cite{Eina01} one finds
$h=0.65\pm0.07$, $\Omega_b=0.05\pm0.02$, $\Omega_{CDM}=0.36\pm0.05$, 
$\Omega_{HDM}\le 0.05$ and $\Omega_{\Lambda}=0.65\pm0.05$.
Since the non exotic component cannot exceed $40\%$ of the CDM 
~\cite {Benne}, there is room for the exotic WIMP's 
(Weakly  Interacting Massive Particles).
  In fact the DAMA experiment ~\cite {BERNA2} 
has claimed the observation of one signal in direct detection of a WIMP, which
with better statistics has subsequently been interpreted as a modulation signal
~\cite{BERNA1}.

 In the most favored scenario of supersymmetry the
LSP can be simply described as a Majorana fermion, a linear 
combination of the neutral components of the gauginos and Higgsinos
\cite{ref1,Gomez,ref2}. 


\section{An Overview of Direct Detection - The Allowed SUSY Parameter Space.}

 Since this particle is expected to be very massive, $m_{\chi} \geq 30 GeV$, and
extremely non relativistic with average kinetic energy $T \leq 100 KeV$,
it can be directly detected ~\cite{JDV96,KVprd} mainly via the recoiling
of a nucleus (A,Z) in elastic scattering.
In order to compute the event rate one needs
the following ingredients:

1) An effective Lagrangian at the elementary particle 
(quark) level obtained in the framework of supersymmetry as described 
, e.g., in Refs.~\cite{ref2,JDV96}.

2) A procedure in going from the quark to the nucleon level, i.e. a quark 
model for the nucleon. The results depend crucially on the content of the
nucleon in quarks other than u and d. This is particularly true for the scalar
couplings as well as the isoscalar axial coupling ~\cite{drees}$-$\cite{Chen}.

3) Compute the relevant nuclear matrix elements \cite{Ress,DIVA00}
using as reliable as possible many body nuclear wave functions.
The situation is a bit simpler in the case of the scalar coupling, in which
case one only needs the nuclear form factor.

Since the obtained rates are very low, one would like to be able to exploit the
modulation of the event rates due to the earth's
revolution around the sun \cite{Verg98,Verg99}$-$\cite{Verg01}. To this end one
adopts a folding procedure
assuming some distribution~\cite{Verg99,Verg01} of velocities for the LSP.
One also would like to know the directional rates, by observing the 
nucleus in a certain direction, which correlate with the motion of the
sun around the center of the galaxy and the motion of the Earth 
\cite {ref1,UKDMC}.

 The calculation of this cross section  has become pretty standard.
 One starts with   
representative input in the restricted SUSY parameter space as described in
the literature~\cite{Gomez,ref2} (see also Arnowitt and Dutta \cite{ARNDU},
which will appear in these proceedings). 
We will adopt a phenomenological procedure taking  universal soft 
SUSY breaking terms at $M_{GUT}$, i.e., a 
common mass for all scalar fields $m_0$, a common gaugino mass 
$M_{1/2}$ and a common trilinear scalar coupling $A_0$, which 
we put equal to zero (we will discuss later the influence of 
non-zero $A_0$'s). Our effective theory below $M_{GUT}$ then 
%depends on the parameters ($\mu_0=\mu(M_{GUT})$)
depends on the parameters \cite{Gomez}:
\[
m_0,\ M_{1/2},\ \mu_0,\ \alpha_G,\ M_{GUT},\ h_{t},\ ,\ h_{b},\ ,\ h_{\tau},\ 
\tan\beta~,  
\]
where $\alpha_G=g_G^2/4\pi$ ($g_G$ being the GUT gauge coupling 
constant) and $h_t, h_b, h_\tau $ are respectively the top, bottom and 
tau Yukawa coupling constants at $M_{GUT}$. The values of $\alpha_G$ and 
$M_{GUT}$ are obtained as described in Ref.\cite{Gomez}.
For a specified value of $\tan\beta$ at $M_S$, we determine $h_{t}$ at 
$M_{GUT}$ by fixing the top quark mass at the center of its 
experimental range, $m_t(m_t)= 166 \rm{GeV}$. The value
of  $h_{\tau}$ at $M_{GUT}$ is  fixed by using the running tau lepton
mass at $m_Z$, $m_\tau(m_Z)= 1.746 \rm{GeV}$. 
The value of $h_{b}$ at $M_{GUT}$ used is such that:
\[
m_b(m_Z)_{SM}^{\overline{DR}}=2.90\pm 0.14~{\rm GeV}.
\]
after including the SUSY threshold correction. 
%$m_b(m_Z)$ is at the center experimental upper bound \cite{mb}:
%\[
%m_b(m_Z)=2.67 \pm 0.50~{\rm GeV}.
%\]
The SUSY parameter space is subject to the
 following constraints:\\
1.) The LSP relic abundance will satisfy the cosmological constrain:
\begin{equation}
0.09 \le \Omega_{LSP} h^2 \le 0.22
\label{eq:in2}
\end{equation}
 2.) The Higgs bound obtained from recent CDF \cite{VALLS}
 and LEP2 \cite{DORMAN},
%\cite{VALDOR},
 i.e. $m_h~>~113~GeV$.\\
3.) We will limit ourselves to LSP-nucleon cross sections for the scalar
coupling, which gives detectable rates
\begin{equation}
4\times 10^{-7}~pb~ \le \sigma^{nucleon}_{scalar} 
\le 2 \times 10^{-5}~pb~
\label{eq:in3}
\end{equation}
 We should remember that the event rate does not depend only
 on the nucleon cross section, but on other parameters also, mainly
 on the LSP mass and the nucleus used in target. 
The condition on the nucleon cross section imposes severe constraints on the
acceptable parameter space. In particular in our model it restricts 
$tan \beta$ to values $tan \beta \simeq 50$. We will not elaborate further
on this point, since related work has already appeared elsewhere
 \cite{Gomez,JDV02}. 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\bigskip
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%
\section{Expressions for the Differential Cross Section .} 
\bigskip

 The effective Lagrangian describing the LSP-nucleus cross section can
be cast in the form \cite {JDV96}
\beq
{\it L}_{eff} = - \frac {G_F}{\sqrt 2} \{({\bar \chi}_1 \gamma^{\lambda}
\gamma_5 \chi_1) J_{\lambda} + ({\bar \chi}_1 
 \chi_1) J\}
 \label{eq:eg 41}
\eeq
where
\beq
  J_{\lambda} =  {\bar N} \gamma_{\lambda} (f^0_V +f^1_V \tau_3
+ f^0_A\gamma_5 + f^1_A\gamma_5 \tau_3)N~~,~~
J = {\bar N} (f^0_s +f^1_s \tau_3) N
 \label{eq:eg.42}
\eeq
%and
%\beq
%J = {\bar N} (f^0_s +f^1_s \tau_3) N
% \label{eq:eg.45}
%\eeq
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

We have neglected the uninteresting pseudoscalar and tensor
currents. Note that, due to the Majorana nature of the LSP, 
${\bar \chi_1} \gamma^{\lambda} \chi_1 =0$ (identically).

 With the above ingredients the differential cross section can be cast in the 
form \cite{ref1,Verg98,Verg99}
\begin{equation}
d\sigma (u,\upsilon)= \frac{du}{2 (\mu _r b\upsilon )^2} [(\bar{\Sigma} _{S} 
                   +\bar{\Sigma} _{V}~ \frac{\upsilon^2}{c^2})~F^2(u)
%\frac{u}{\rangle~u~\rangle})\frac{\upsilon^2}{c^2}\frac{u}{\rangle~u~\rangle}) 
                       +\bar{\Sigma} _{spin} F_{11}(u)]
\label{2.9}
\end{equation}
\begin{equation}
\bar{\Sigma} _{S} = \sigma_0 (\frac{\mu_r(A)}{\mu _r(N)})^2  \,
 \{ A^2 \, [ (f^0_S - f^1_S \frac{A-2 Z}{A})^2 \, ] \simeq \sigma^S_{p,\chi^0}
        A^2 (\frac{\mu_r(A)}{\mu _r(N)})^2 
\label{2.10}
\end{equation}
\begin{equation}
\bar{\Sigma} _{spin}  =  \sigma^{spin}_{p,\chi^0}~\zeta_{spin}~~,~~
\zeta_{spin}= \frac{(\mu_r(A)/\mu _r(N))^2}{3(1+\frac{f^0_A}{f^1_A})^2}S(u)
\label{2.10a}
\end{equation}
%\begin{equation}
%\zeta_{spin}= \frac{(\mu_r(A)/\mu _r(N))^2}{3(1+\frac{f^0_A}{f^1_A})^2}S(u)
%\label{2.10b1}
%\end{equation}
\begin{equation}
S(u)=[(\frac{f^0_A}{f^1_A} \Omega_0(0))^2 \frac{F_{00}(u)}{F_{11}(u)}
  +  2\frac{f^0_A}{ f^1_A} \Omega_0(0) \Omega_1(0)
\frac{F_{01}(u)}{F_{11}(u)}+  \Omega_1(0))^2  \, ] 
\label{2.10b2}
\end{equation}
\begin{equation}
\bar{\Sigma} _{V}  =  \sigma^V_{p,\chi^0}~\zeta_V 
\label{2.10c}
\end{equation}
\begin{equation}
\zeta_V =  \frac{(\mu_r(A)/\mu _r(N))^2}{(1+\frac{f^1_V}{f^0_V})^2} A^2 \, 
(1-\frac{f^1_V}{f^0_V}~\frac{A-2 Z}{A})^2 [ (\frac{\upsilon_0} {c})^2  
[ 1  -\frac{1}{(2 \mu _r b)^2} \frac{2\eta +1}{(1+\eta)^2} 
\frac{\langle~2u~ \rangle}{\langle~\upsilon ^2~\rangle}] 
\label{2.10d}
\end{equation}
\\
$\sigma^i_{p,\chi^0}=$ proton cross-section,$i=S,spin,V$ given by:\\
$\sigma^S_{p,\chi^0}= \sigma_0 ~(f^0_S)^2~(\frac{\mu _r(N)}{m_N})^2$ 
 (scalar) , 
(the isovector scalar is negligible, i.e. $\sigma_p^S=\sigma_n^S)$\\
$\sigma^{spin}_{p,\chi^0}=\sigma_0~~3~(f^0_A+f^1_A)^2~(\frac{\mu _r(N)}{m_N})^2$ 
  (spin) ,
$\sigma^{V}_{p,\chi^0}= \sigma_0~(f^0_V+f^1_V)^2~(\frac{\mu _r(N)}{m_N})^2$ 
(vector)   \\
where $m_N$ is the nucleon mass,
 $\eta = m_x/m_N A$, and
 $\mu_r(A)$ is the LSP-nucleus reduced mass,  
 $\mu_r(N)$ is the LSP-nucleon reduced mass and  
\begin{equation}
\sigma_0 = \frac{1}{2\pi} (G_F m_N)^2 \simeq 0.77 \times 10^{-38}cm^2 
\label{2.7} 
\end{equation}
\begin{equation}
%u = q^2b^2/2~~,~~
Q=Q_{0}u~~, \qquad Q_{0} = \frac{1}{A m_{N} b^2}=4.1\times 10^4A^{-4/3}~KeV 
\label{2.15} 
\end{equation}
where
%b is (the harmonic oscillator) size parameter, 
%q is the momentum transfer to the nucleus, and
Q is the energy transfer to the nucleus and
$F(u)$ is the nuclear form factor. 

In the present paper we will concentrate on the coherent mode. For a discussion
of the spin contribution, expected to be important in the case of the light
nuclei, has been reviewed elsewhere \cite{DIVA00,JDV02}.

\section{Expressions for the Rates.} 
 The non-directional event rate is given by:
\begin{equation}
R=R_{non-dir} =\frac{dN}{dt} =\frac{\rho (0)}{m_{\chi}} \frac{m}{A m_N} 
\sigma (u,\upsilon) | {\boldmath \upsilon}|
\label{2.17} 
\end{equation}
 Where
 $\rho (0) = 0.3 GeV/cm^3$ is the LSP density in our vicinity and 
 m is the detector mass 
The differential non-directional  rate can be written as
\begin{equation}
dR=dR_{non-dir} = \frac{\rho (0)}{m_{\chi}} \frac{m}{A m_N} 
d\sigma (u,\upsilon) | {\boldmath \upsilon}|
\label{2.18}  
\end{equation}
where $d\sigma(u,\upsilon )$ was given above.

 The directional differential rate \cite{ref1},\cite{Verg01} in the
direction $\hat{e}$ is given by :
\begin{equation}
dR_{dir} = \frac{\rho (0)}{m_{\chi}} \frac{m}{A m_N} 
{\boldmath \upsilon}.\hat{e} H({\boldmath \upsilon}.\hat{e})
 ~\frac{1}{2 \pi}~  
d\sigma (u,\upsilon)
\label{2.20}  
\end{equation}
where H the Heaviside step function. The factor of $1/2 \pi$ is 
introduced, since  the differential cross section of the last equation
is the same with that entering the non-directional rate, i.e. after
an integration
over the azimuthal angle around the nuclear momentum has been performed.
In other words, crudely speaking, $1/(2 \pi)$ is the suppression factor we
 expect in the directional rate compared to the usual one. The precise 
suppression factor depends, of course, on the direction of observation.
The mean value of the non-directional event rate of Eq. (\ref {2.18}), 
is obtained by convoluting the above expressions with the LSP velocity
distribution $f({\bf \upsilon}, {\boldmath \upsilon}_E)$ 
with respect to the Earth, i.e. is given by:
\beq
\Big<\frac{dR}{du}\Big> =\frac{\rho (0)}{m_{\chi}} 
\frac{m}{A m_N}  
\int f({\bf \upsilon}, {\boldmath \upsilon}_E) 
          | {\boldmath \upsilon}|
                      \frac{d\sigma (u,\upsilon )}{du} d^3 {\boldmath \upsilon} 
\label{3.10} 
\eeq
 The above expression can be more conveniently written as
\beq
\Big<\frac{dR}{du}\Big> =\frac{\rho (0)}{m_{\chi}} \frac{m}{Am_N} \sqrt{\langle
\upsilon^2\rangle } {\langle \frac{d\Sigma}{du}\rangle }~,~ 
\langle \frac{d\Sigma}{du}\rangle =\int
           \frac{   |{\boldmath \upsilon}|}
{\sqrt{ \langle \upsilon^2 \rangle}} f({\boldmath \upsilon}, 
         {\boldmath \upsilon}_E)
                       \frac{d\sigma (u,\upsilon )}{du} d^3 {\boldmath \upsilon}
\label{3.11}  
\eeq
%where
%\beq
%\langle \frac{d\Sigma}{du}\rangle =\int
%           \frac{   |{\boldmath \upsilon}|}
%{\sqrt{ \langle \upsilon^2 \rangle}} f({\boldmath \upsilon}, 
%         {\boldmath \upsilon}_E)
%                       \frac{d\sigma (u,\upsilon )}{du} d^3 {\boldmath \upsilon}
%\label{3.12}  
%\eeq


 After performing the needed integrations over the velocity distribution,
to first order in the Earth's velocity, and over the energy transfer u  the
 last expression takes the form
\beq
R =  \bar{R}~t~
          [1 + h(a,Q_{min})cos{\alpha})] 
\label{3.55a}  
\eeq
where $\alpha$ is the phase of the Earth ($\alpha=0$ around June 2nd)
and  $Q_{min}$ is the energy transfer cutoff imposed by the detector.
In the above expressions $\bar{R}$ is the rate obtained in the conventional 
approach \cite {JDV96} by neglecting the folding with the LSP velocity and the
momentum transfer dependence of the differential cross section, i.e. by
\beq
\bar{R} =\frac{\rho (0)}{m_{\chi}} \frac{m}{Am_N} \sqrt{\langle
v^2\rangle } [\bar{\Sigma}_{S}+ \bar{\Sigma} _{spin} + 
\frac{\langle \upsilon ^2 \rangle}{c^2} \bar{\Sigma} _{V}]
\label{3.39b}  
\eeq
where $\bar{\Sigma} _{i}, i=S,V,spin$ 
%have been defined above, see Eqs (\ref {2.10}) - (\ref {2.10c}). They 
 contain all the parameters of the
SUSY models. 
 The modulation is described by the parameter $h$ .

The total  directional event rates  can be obtained in a similar fashion by
suitably modifying Eq. (\ref{3.10})
\beq
\Big<\frac{dR}{du}\Big>_{dir} =\frac{\rho (0)}{m_{\chi}} 
\frac{m}{A m_N}  
\int f({\bf \upsilon}, {\boldmath \upsilon}_E) 
\frac{{\boldmath \upsilon}.\hat{e} H({\boldmath \upsilon}.\hat{e})}{2 \pi}~  
                      \frac{d\sigma (u,\upsilon )}{du} d^3 {\boldmath \upsilon} 
\label{3.10b} 
\eeq
The integration of the above equation is difficult. So we find it convenient to
go to a coordinate system in which the polar axis is in the direction of 
of observation $\hat{e}$, which in the above coordinate system is specified
by the polar angle $\Theta$ and the azymouthal angle $\Phi$. In this new
 coordinate system polar angle specifying the velocity vector 
is simply restricted to be
 $0\le \theta \le \pi$, while the azymouthal angle $\phi$ is unrestricted.
Thus the unit vectors along the new coordinate axes, $\hat{X},\hat{Y},\hat{Z}$,
are expressed in terms of the old ones as follows:
\beq
\hat{Z}=\sin \Theta \cos \Phi \hat{x}+ \sin \Theta \sin \Phi \hat{y}+
\cos \Theta \hat{z}.
\label{3.20} 
\eeq
\beq
\hat{X}=\cos \Theta \cos \Phi \hat{x}+ \cos \Theta \sin \Phi \hat{y}-
\sin \Theta \hat{x}.
\label{3.21} 
\eeq
\beq
\hat{Y}=-\sin \Phi \hat{x}+ \cos \Phi \hat{y}.
\label{3.22} 
\eeq
Thus the LSP velocity is expressed in the new coordinate system as:
\beq
\upsilon_x=\sin \Theta \cos \Phi \upsilon_X+ \sin \Theta \sin \Phi \upsilon_Y+
\cos \Theta \upsilon_Z,
\label{3.23} 
\eeq
\beq
\upsilon_y=\cos \Theta \cos \Phi \upsilon_X+ \cos \Theta \sin \Phi \upsilon_Y-
\sin \Theta \upsilon_Z,
\label{3.24} 
\eeq
\beq
\upsilon_z=-\sin \Phi \upsilon_X+ \cos \Phi \upsilon_Y,
\label{3.25} 
\eeq
with $\upsilon_X=\upsilon \sin \theta \cos \phi$,
$\upsilon_Y=\upsilon \sin \theta \sin \phi$,
$\upsilon_Z=\upsilon \cos \theta $.
It is thus straightforward to go to polar coordinates in velocity space and
get:
\beq
\Big<\frac{dR}{du}\Big>_{dir} =\frac{\rho (0)}{m_{\chi}} 
\frac{m}{A m_N}  
\int^{\upsilon_m}_0 \upsilon^3 d \upsilon \int^1_0 d \xi \int^{2 \pi}_0 d \phi
\frac{\tilde{f}({\Theta,\Phi, \upsilon},\xi, \phi , {\boldmath \upsilon}_E)} 
           {2 \pi}
                      \frac{d\sigma (u,\upsilon )}{du}  
\label{3.26} 
\eeq
Now the orientation parameters $\Theta$ and $\Phi$ appear explicitly and not 
implicitly via the limits of integration. The function $\tilde{f}$ can be
obtained from the velocity distribution, but it will not be explicitly shown
here. Thus we obtain:
\beq
R_{dir}  =   \bar{R} [(t_{dir}/2 \pi) \, 
            [1 + (h_1-h_2)cos{\alpha}) + h_3 sin{\alpha}]
\label{4.56}  
\eeq
where the quantity $t_{dir}$ provides the non modulated amplitude,
 while $h_1,h_2$
and $h_3$ describe the modulation.
 They are functions of  $\Theta$ and $\Phi$ as well as the
parameters $a$ and $Q_{min}$.
 The effect of folding
with LSP velocity on the total rate is taken into account via the quantity
$t_{dir}$, which depends on the LSP mass. All other SUSY parameters have been
 absorbed in $\bar{R}$. 
 In the special case previously studied, i.e 
along the coordinate axes, we find that: a) in the direction of the sun's
 motion
 $h_2=h_3=0$, b)
along the radial direction (y axes) $h_3=0$ and c) in the vertical to the
 galaxy $h_2=0$.
 Instead of $t_{dir}$ itself it is more convenient
to present the reduction factor of the non modulated directional rate compared
to the usual non-directional one, i.e.
\beq
f_{red}=\frac{R_{dir}}{R}=t_{dir}/(2 \pi~t)= \kappa/(2 \pi) 
\label{kappa.eq}
\eeq
%\barr
%R_{dir}& = &  \bar{R} [\rho^{'}(0)/(\rho(0)] \, (t^0/4 \pi) \, 
%R_{dir}& = &  \bar{R} [(t^0/4 \pi) \, 
%            |(1 + h_1(a,Q_{min})cos{\alpha}) {\bf e}~_z.{\bf e}
%\nonumber\\  &-& h_2(a,Q_{min})\, 
%cos{\alpha} {\bf e}~_y.{\bf e}
%                      + h_3(a,Q_{min})\, 
%sin{\alpha} {\bf e}~_x.{\bf e}|
%\label{4.55}  
%\earr
  It turns out that the parameter $\kappa$, being the 
ratio of two rates, is less dependent on these parameters.
Another quantity, which may be of experimental interest is the asymmetry
$|R_{dir}(-)-R_{dir}(+)|/(R_{dir}(-)+R_{dir}(+))$. This is completely
 independent of all other parameters except the LSP mass and the velocity
distribution.
 The directional rates exhibit interesting pattern of modulation. Given the
functions $h_l(a,Q_{min})$, $l=1,2,3$, one can plot the the
expression in
Eqs (\ref {3.55a}) and \ref{4.56} as a function of the phase of the earth
 $\alpha$. 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

\section{Results and Discussion}
The three basic ingredients of our calculation were 
the input SUSY parameters, a quark model for the nucleon
and the velocity distribution combined with the structure of
 the nuclei involved. We will focus our our discussion on the
coherent scattering and present results for the popular target $^{127}I$.
We have utilized two nucleon models indicated by B and C, for their description
see our previous work \cite{JDV02}, which take into
account the presence of heavy quarks in the nucleon. We also considered
the effects on the rates of the energy cut off imposed by the detector,  
 by considering two  typical cases $Q_{min}=0,~10$ KeV. The thus obtained
results for the non modulated non directional event rates, $\bar{R}t$,
 in the case
 of the symmetric isothermal model for a typical SUSY parameter choice
\cite{Gomez} are shown in Fig. \ref{rate}. 
\begin{figure}
\includegraphics[height=.3\textheight]{rates.eps}
%\begin{figure}
%\epsfig{figure=rates.eps,height=3.90in,angle=0}
\caption{The Total detection rate per $(kg-target)yr$ vs the LSP mass
in GeV for a typical solution in our parameter space in the case of 
$^{127}I$ corresponding to  model B (thick line) and Model C (fine 
line). For the definitions see text.
\label{rate}
}
\end{figure}
  The two relative 
parameters, i.e. the quantities $t$ and $h$, are shown in Figs \ref{fig.t}
and \ref{fig.h} respectively in the case of isothermal
 models.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\begin{figure}
%\vspace*{-25.0 cm}
\hspace*{-0.0 cm}
\includegraphics[height=.2\textheight]{nt0.ps}
%\newline \hspace*{4cm} {\bf ($\lambda=0$)}
%\hspace{-0.5 cm}
\includegraphics[height=.2\textheight]{nt1.ps}
%\newline \hspace*{4cm} {\bf ($\lambda=1$)}
%\vspace*{-25.0 cm}
\caption{
The dependence of the quantity $t$ on the LSP mass
 for the symmetric case ($\lambda=0$) on the left as well as for the
 maximum axial asymmetry ($\lambda=1$) on the right in the case of the
 target $^{127}I$.
For orientation purposes  two  detection cutoff energies are exhibited,
$Q_{min}=0$ (thick solid line) and $Q_{min}=10~keV$ (thin solid line). 
 As expected $t$ decreases as the cutoff energy and/or the LSP 
mass increase. We see that the asymmetry parameter $\lambda$ has little
effect on the non modulated rate. 
\label{fig.t}
}
\end{figure}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\begin{figure}
\hspace*{-0.0 cm}
\includegraphics[height=.2\textheight]{nh0.ps}
\includegraphics[height=.2\textheight]{nh1.ps}
\caption{ The same as in Fig. \ref{fig.t}
for the modulated amplitude. We see that the
modulation is small and decreases with the LSP mass.
 It even changes sign
for large LSP mass. The
 introduction of a cutoff $Q_{min}$ increases the modulation (at the expense
of the total number of counts).
 It also increases
 with the asymmetry parameter $\lambda$.
\label{fig.h}
}
\end{figure}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%\begin{figure}
%\includegraphics[height=.3\textheight]{nh1.ps}
%\caption{ The same as in Fig. \ref{fig.h0} for $\lambda=1$.
%We see that the modulation increases
% with the asymmetry parameter $\lambda$.
%\label{fig.h1}
%}
%\end{figure}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 The case of non isothermal models, e.g. caustic rings, is more
complicated \cite{Verg01} and it will only briefly be discussed here.
 
In the case of the directional rates we calculated the reduction factors
and the asymmetry parameters as well as the modulation amplitudes
 as functions of the direction of observation, focusing our attention along
 the three axes, i.e along $+z,-z,+y,-y,+x$ and $-x$
\cite{Verg00}. Since $f_{red}$ is the ratio of two parameters, its dependence
on $Q_{min}$ and the LSP mass is mild. So we present results for
 $Q_{min}=0$
and give an  average as a function of the LSP mass (see Table 
\ref{table.dir}).
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

 As we have seen the modulation can be described in terms of the parameters
$h_i,~i=1,2,3$ (see Eq.  (\ref{4.56})). If the observation is done in the
direction opposite to the sun's direction of motion, the modulation amplitude
$h_1$ behaves in the same way as the non directional one, namely $h$. It is
more instructive to consider directions of observation in the plane
perpendicular to the sun's direction of motion ($\Theta=\pi/2$) even though
the non modulated rate is quite reduced
in this direction. Along the $-y$ direction ($\Phi=(3/2)\pi$)
the modulation amplitude $h_1-h_2$ is constant, $-0.20$ and $-0.30$ for
 $\lambda=0,1$ respectively. In other words it is large and leads to a maximum
rate in December. Along the $+y$ direction the modulation is exhibited in Fig.
\ref{fig.h12}.
% and \ref{fig.h13}.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\begin{figure}
\includegraphics[height=.2\textheight]{nhph12.ps}
\includegraphics[height=.2\textheight]{nhph13.ps}
\caption{ On the left figure one sees the quantity $h_1-h_2$ in the
 direction $+y$ for
 $\lambda=0$ (thick line) and $\lambda=1$ thin line. In the $-y$ direction
this quantity is constant and negative, $-0.20$ and $-0.30$ for $\lambda=0$
 and $1$ respectively. As a result the modulation effect is opposite 
(minimum in June the 3nd). On the right figure we show
%the same as in Fig. \ref{fig.h12} for
the modulation amplitude in the direction $+x$, which is essentially
$h_3$, since$|h_1|<<|h_3|$.  Thus in this case
the maximum occurs around September the 3nd and the minimum 6 months later.
The opposite is true in the $-x$ direction.
\label{fig.h12}
}
\end{figure}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%\begin{figure}
%\includegraphics[height=.3\textheight]{nhph13.ps}
%\caption{ The same as in Fig. \ref{fig.h12} for
%the modulation amplitude in the direction $+x$, which is essentially
%$h_3$, since$|h_1|<<|h_3|$.  Thus in this case
%the maximum occurs around September the 3nd and the minimum 6 months later.
%The opposite is true in the $-x$ direction.
%\label{fig.h13}
%}
%\end{figure}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\begin{table}[t]  
\caption{
The ratio $\kappa$ of the non modulated directional rate along the three
 directions to the non-directional one: z is in the direction of the sun's
 motion, x is in the radial direction  and x is perpendicular to the axis 
of the galaxy. The asymmetry is also given. $Q_{min}=0$ was assumed.
}
\label{table.dir}
%\begin{center}
\begin{tabular}{lrrrrrrr}
%\hline\hline
& & & & & & &    \\
&  & \multicolumn{3}{c|}{isothermal} & \multicolumn{3}{c|}{caustic rings}  \\ 
 &   &  &  & & & &     \\
  $\lambda$ & dir. & +  & - &asym&+&-&asym\\
\hline
    &   &  &  &      \\
0    & z   & 0.02  & 0.50  &  0.92&0.75&0.25&0.50\\
0    & y   & 0.16  & 0.16  &  0&0.22&0&1.00\\
0    & x  &  0.16  & 0.16  &  0&0.37&0.24&0.21\\
1    & z   & 0.04  & 0.58  &  0.90&-&-&- \\
1    & y   & 0.12  & 0.12  &  0&-&-&-\\
1    & x  &  0.17  & 0.17  &  0&-&-&-\\
%\hline
%\hline
\end{tabular}
%\end{center}
\end{table}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{Conclusions}
In the present paper we have discussed the parameters, which describe
the event rates for direct detection of SUSY dark matter.
Only in a small segment of the allowed parameter space the rates are above
 the present experimental goals.
 We thus looked for
characteristic experimental signatures for background reduction, i.e.
a) Correlation of the event rates
with the motion of the Earth (modulation effect) and b)
 the directional rates (their correlation  both with the velocity of the sun
 and that of the Earth.)

A typical graph for the total non modulated rate is shown Fig. \ref{rate}. 
The relative parameters $t$ and $h$ in the case of non directional experiments
are exhibited in Fig. \ref{fig.t} and  Fig. \ref{fig.h}.
We must emphasize that the two graphs of Fig. \ref{fig.t} do not contain
the entire dependence on
the LSP mass. This is due to the fact that there is the extra factor
$m^{-1}_{\chi}$ in Eq. (\ref{3.39b}) and a factor of $\mu_r^2$ arising from
 $\Sigma_i$,
$i=S,~spin,~V$, see Eqs (\ref{2.10}), (\ref{2.10a}), and  (\ref{2.10d}).
 All these factors combined lead
 to a constant. There remains, however, an LSP mass dependence, which is  due
to the fact that the nucleon cross section itself depends on the LSP mass.

 Figs \ref{fig.t} and  \ref{fig.h}  were obtained for
 the scalar interaction, but similar behavior is expected for the spin
 contribution.
 The scale of the total spin contribution,
however, is going to be very different.
 We should also mention that in the non directional
experiments the modulation $2h_1$ is small, .i.e. for $\lambda=0$ less than
 $4\%$ for
$Q{min}=0$ and $12\%$ for $Q{min}=10~KeV$ (at the expense of the total
number of counts). For $\lambda=1$ there in no change for
$Q_{min}=0$, but it can go as high as $24\%$ for $Q_{min}=10~KeV$. 

For the directional rates It is instructive to examine the reduction factors
 along the three axes, i.e 
along $+z,-z,+y,-y,+x$ and $-x$. These depend on the 
nuclear parameters, the reduced mass, the energy cutoff $Q_{min}$ and 
 $\lambda$
\cite{Verg00}. Since $f_{red}$ is the ratio of two parameters, its dependence
on $Q_{min}$ and the LSP mass is mild. So we present results for $Q_{min}=0$
and give their average as a function of the LSP mass (see Table 
\ref{table.dir}). As expected the maximum rate is  along the sun's direction of
 motion, i.e opposite to its velocity $(-z)$ in the Gaussian distribution
 and $+z$ in the case of caustic rings. In fact we find
that $\kappa(-z)$ is around 0.5 (no asymmetry) and around 0.6 (maximum
asymmetry, $\lambda=1.0$). It is not very different from the 
naively expected $f_{red}=1/( 2 \pi)=\kappa=1$.
The asymmetry along the sun's direction of motion,
$asym=|R_{dir}(-)-R_{dir}(+)|/(R_{dir}(-)+R_{dir}(+))$ is quite characteristic
and, for a given direction, it depends on the velocity distribution.
It is large in the direction of the sun's motion in the
isothermal models and a bit smaller in caustic rings. The rate in the other
directions is quite
a bit smaller (see Table \ref{table.dir}). In this case for the isothermal model
 the asymmetry
 is zero, while in the case of caustic rings it is sizable. So the two models
can be easily distinguished. The disadvantage of smaller rates in the plane
perpendicular to the sun's velocity may be compensated by the very large
and characteristic modulation.
It is interesting to note that
 even the functional dependence on the phase of the Earth $\alpha$ changes
 substantially with the direction of observation. 
 
 In conclusion: in the case of directional non modulated rates we expect
 unambiguous correlation with the motion on the sun, which can be explored by
 the experimentalists.
 The reduction factor in the direction of the
 motion
of the sun is approximately only $1/(4\pi)$ relative to the non directional
experiments.  In the plane perpendicular to the motion of the sun we expect
interesting modulation signals, but the reduction factor becomes worse.
 These difficulties may be reduced, if the TPC counters can make observations
in many directions simultaneously.
%A more complete discussion will be given elsewhere \cite{JDVb02}.


\par
This work was supported by the European Union under the contracts 
RTN No HPRN-CT-2000-00148 and TMR 
No. ERBFMRX--CT96--0090.

\def\ijmp#1#2#3{{ Int. Jour. Mod. Phys. }{\bf #1~}(#2)~#3}
\def\pl#1#2#3{{ Phys. Lett. }{\bf B#1~}(#2)~#3}
\def\zp#1#2#3{{ Z. Phys. }{\bf C#1~}(#2)~#3}
\def\prl#1#2#3{{ Phys. Rev. Lett. }{\bf #1~}(#2)~#3}
\def\rmp#1#2#3{{ Rev. Mod. Phys. }{\bf #1~}(#2)~#3}
\def\prep#1#2#3{{ Phys. Rep. }{\bf #1~}(#2)~#3}
\def\pr#1#2#3{{ Phys. Rev. }{\bf D#1~}(#2)~#3}
\def\np#1#2#3{{ Nucl. Phys. }{\bf B#1~}(#2)~#3}
\def\npps#1#2#3{{ Nucl. Phys. (Proc. Sup.) }{\bf B#1~}(#2)~#3}
\def\mpl#1#2#3{{ Mod. Phys. Lett. }{\bf #1~}(#2)~#3}
\def\arnps#1#2#3{{ Annu. Rev. Nucl. Part. Sci. }{\bf
#1~}(#2)~#3}
\def\sjnp#1#2#3{{ Sov. J. Nucl. Phys. }{\bf #1~}(#2)~#3}
\def\jetp#1#2#3{{ JETP Lett. }{\bf #1~}(#2)~#3}
\def\app#1#2#3{{ Acta Phys. Polon. }{\bf #1~}(#2)~#3}
\def\rnc#1#2#3{{ Riv. Nuovo Cim. }{\bf #1~}(#2)~#3}
\def\ap#1#2#3{{ Ann. Phys. }{\bf #1~}(#2)~#3}
\def\ptp#1#2#3{{ Prog. Theor. Phys. }{\bf #1~}(#2)~#3}
\def\plb#1#2#3{{ Phys. Lett. }{\bf#1B~}(#2)~#3}
\def\apjl#1#2#3{{ Astrophys. J. Lett. }{\bf #1~}(#2)~#3}
\def\n#1#2#3{{ Nature }{\bf #1~}(#2)~#3}
\def\apj#1#2#3{{ Astrophys. Journal }{\bf #1~}(#2)~#3}
\def\anj#1#2#3{{ Astron. J. }{\bf #1~}(#2)~#3}
\def\mnras#1#2#3{{ MNRAS }{\bf #1~}(#2)~#3}
\def\grg#1#2#3{{ Gen. Rel. Grav. }{\bf #1~}(#2)~#3}
\def\s#1#2#3{{ Science }{\bf #1~}(19#2)~#3}
\def\baas#1#2#3{{ Bull. Am. Astron. Soc. }{\bf #1~}(#2)~#3}
\def\ibid#1#2#3{{ ibid. }{\bf #1~}(19#2)~#3}
\def\cpc#1#2#3{{ Comput. Phys. Commun. }{\bf #1~}(#2)~#3}
\def\astp#1#2#3{{ Astropart. Phys. }{\bf #1~}(#2)~#3}
\def\epj#1#2#3{{ Eur. Phys. J. }{\bf C#1~}(#2)~#3}
%\doingARLO[\bibliographystyle{aipproc}]
%          {\ifthenelse{\equal{\AIPcitestyleselect}{num}}
%             {\bibliographystyle{arlonum}}
%             {\bibliographystyle{arlobib}} 
%\bibitem{GAW}E. Gawiser and J. Silk,{\it {Science}} {\bf 280}, 1405 (1988);
%@Article{Jungm:1996,
%{G. Jungman  et al.},{ Phys. Rep.} {267}, {195} , {1996}
%}
%@Article{COBE:1992,
%{Smoot, G.F. et al}, (COBE data), { Astrophys. J.}, {396 L1}, {1992}
%}
%\begin{thebibliography}{99}
\begin{thebibliography}{8.}
\addcontentsline{toc}{section}{References}
\bibitem{flat01}
Jaffe, A.H. {\it et al.},{\it  Phys. Rev. Lett.} {\bf 86}, 3475  (2001).
\bibitem{COBE} Smoot, G.F. et al., (COBE data), {\it Astrophys. J.} {\bf 396} 
(1992) L1.
\bibitem{GAW}Gawiser, E. and Silk, J.,{\it {Science}} {\bf 280}, 1405 (1988);
 Gross, M.A.K. Somerville, R.S.,  Primack, J.R.,  Holtzman , J. and  Klypin,
 A.A., {\it Mon. Not. R. Astron. Soc.} {\bf 301}, 81 (1998).
\bibitem{HSST} Riess, A.G.  {\it et al}, {\it Astron. J.}
{\bf 116} (1998), 1009.
\bibitem{SPF} Somerville, R.S., J.R. Primack and S.M. Faber, ;
{\it Mon. Not. R. Astron. Soc.} (in press).
\bibitem{SCP}Perlmutter, S. {\it et al} (1999) {\it Astrophys. J.}
{\bf 517},565; (1997) {\bf 483},565 ({\it astro-ph}/9812133).\\
Perlmutter, S., Turner, M.S.  and White, M., {\it Phys. Rev. Let.} {\bf 83},
670 (1999).
%C. Krauss, The cosmological constant and dark matter,2nd Int.
%Workshop on the Identification of Dark Matter, Buxton, UK, 7-11 September
%1998 (to appear in the proceedings).
\bibitem{Eina01} Einasto, Jaan, in Dark Matter inj Astro- and Particle Physics,
p.3, Ed. H.V. Klapdor-Kleingrothaus, Springer-Verlag Berlin Heidelberg 2001.
\bibitem{Benne} Bennett, D.P. {\it et al.}, (MACHO collaboration), A binary
lensing event toward the LMC: Observations and Dark Matter Implications, 
Proc. 5th Annual Maryland Conference, edited by S. Holt (1995);\\
Alcock, C. {\it et al.}, (MACHO collaboration), {\it Phys. Rev. Lett.} {\bf 74}
, 2967 (1995). 
\bibitem{BERNA2} Bernabei, R. et al., INFN/AE-98/34, (1998);
 Bernabei, R. et al., {it Phys. Lett.} {\bf B 389}, 757 (1996).
\bibitem{BERNA1} Bernabei, R. et al., {\it Phys. Lett.} {\bf B 424}, 195 (1998);
{\bf B 450}, 448 (1999).
\bibitem{ref1}For more references see e.g. our previous report:\\
Vergados, J.D., Supersymmetric Dark Matter Detection-
The Directional Rate and the Modulation Effect, ;
\bibitem{Gomez}  G\'{o}mez, M.E.,  Vergados, J.D., Phys. Lett. {\bf B 512}
, 252 (2001); .\\
 G\'{o}mez, M.E.,  Lazarides, G. and Pallis, C.,\pr{61}{2000}{123512} and
 {\it Phys. Lett.} {\bf B 487}, 313 (2000).
\bibitem{gtalk}  G\'{o}mez, M.E. and Vergados, J.D., .
%M.W. Goodman and E. Witten,{\it  Phys. Rev.} {\bf D 31}, 3059 (1985);\\ 
%K. Griest, {\it Phys. Rev. Lett} {\bf  61}, 666 (1988); 
%Phys. Rev. {\bf D 38}, 2357 (1988) ; {\bf D 39}, 3802 (1989);\\
%J. Ellis, and R.A. Flores, {\it Phys. Lett.} {\bf B 263}, 259 (1991);
%{\it Phys. Lett} {\bf B 300}, 175 (1993); {\it Nucl. Phys.} {\bf B 400}, 25 (1993);\\
%J. Ellis and L. Roszkowski, {\it Phys. Lett.} {\bf B 283}, 252 (1992).
\bibitem{ref2}
Bottino, A. {\it et al.}, {\it Phys. Lett} {\bf B 402}, 113 (1997).\\
Arnowitt, R. and Nath, P., {\it Phys. Rev. Lett.}  {\bf 74}, 4952 (1995);
 {\it Phys. Rev.} {\bf D 54}, 2394 (1996); ;\\
 Bednyakov, V.A.,  Klapdor-Kleingrothaus, H. V. and  Kovalenko, S.G.,
{\it Phys. Lett.}  {\bf B 329}, 5 (1994).
\bibitem{ARNDU} Arnowitt, R. and Dutta, B., Supersymmetry anmd Dark Matter,
.
\bibitem{JDV96} Vergados, J.D., {\it J. of Phys.} {\bf G 22}, 253 (1996).
\bibitem{KVprd}Kosmas, T.S. and  Vergados, J.D., {\it Phys. Rev.} {\bf D 55}, 
1752 (1997).
%\bibitem{ref2}
%A. Bottino {\it et al.}, {\it Phys. Lett} {\bf B 402}, 113 (1997);
%R. Arnowitt and P. Nath, {\it Phys. Rev. Lett.}  {\bf 74}, 4952 (1995);
% {\it Phys. Rev.} {\bf D 54}, 2394 (1996); ;\\
%V.A. Bednyakov, H.V. Klapdor-Kleingrothaus and S.G. Kovalenko,
%{\it Phys. Lett.}  {\bf B 329}, 5 (1994).
%\bibitem{als} B. Ananthanarayan, G. Lazarides and Q. Shafi,
%\pr{44}{1991}{1613}. For a more recent update see U. Sarid, 
%.
%\bibitem{pana} G. Lazarides and C. Panagiotakopoulos, 
%\pl{337}{1994}{90}.
%\bibitem{copw} M. Carena, M. Olechowski,
%S. Pokorski and C. E. M. Wagner, \np{426}{1994}{269}; 
%M. Carena and C. E. M. Wagner, .
%\bibitem{hall} L. Hall, R. Rattazzi and U. Sarid,
%\pr{50}{1994}{7048};
%R. Hempfling, \pr{49}{1994}{6168}.
%\bibitem{cdm} M. G\'{o}mez, G. Lazarides and C. Pallis, 
%%pr{61}{2000}{123512}.
%\bibitem{bsg} S. Bertolini, F. Borzumati, A. Masiero and G. Ridolfi, 
%\np{353}{1991}{591}; 
%R. Barbieri and G. F. Giudice, \pl{309}{1993}{86}; 
%R. Garisto and J. N. Ng, \pl{315}{1993}{372}; 
%A. J. Buras, M. Misiak, M. M\"{u}nz and S. Pokorski, 
%\np{424}{1994}{374}; 
%B. de Carlos and J. A. Casas, \pl{349}{1995}{300}; 
%S. Bertolini and F. Vissani, \zp{67}{1995}{513}; 
%V. Barger, M. S. Berger, P. Ohmann and R. J. N. Phillips, 
%\pr{51}{1995}{2438}.
%\bibitem{borzumati} F. M. Borzumati, M. Olechowski
%and S. Pokorski, \pl{349}{1995}{311}.
%\bibitem{bsg2} H. Baer, M. Brhlik, D. Casta\~{n}o and X. Tata, 
%\pr{58}{1998}{015007}.
%\bibitem{lahanas} A. B. Lahanas, D. V. Nanopoulos and 
%V. C. Spanos, \pl{464}{1999}{213}; .  
%\bibitem{coan} K. Griest and D. Seckel, \pr{43}{1991}{3191}.
\bibitem{drees} Drees, M. and Nojiri, M.M.,
\pr{47}{1993}{376}; 
%S. Mizuta and M. Yamaguchi, \pl{298}{1993}{120}; 
%P. Gondolo and J. Edsj\"o, \pr{56}{1997}{1879};
%M. Drees, ;  
%\bibitem{Hab-Ka}H.E. Haber and G.L. Kane, {\it Phys. Rep.} {\bf  117} 75
 (1985).
\bibitem{Dree}Drees, M. and  Nojiri, M.M., {\it Phys. Rev.} {\bf D 48}, 
3843 (1993); {\it Phys. Rev.} {\bf D 47}, 4226 (1993). 
\bibitem{Dree00}Djouadi, A. and Drees, M., Phys.\ Lett.\ B {\bf 484} (2000)
 183;
%QCD Corrections to the neutralino - no-nucleon Scattering, TUM-HEP-370-00;PM/00-14;
Dawson, S., {\it Nucl. Phys.} {\bf B359},283 (1991);
Spira, M. {it et al}, {\it Nucl. Phys.} {\bf B453},17 (1995).
\bibitem{Chen} Cheng, T.P., {\it Phys. Rev.} {\bf D 38} 2869 (1988);
 Cheng, H-Y., Phys. Lett. {\bf B 219} 347 (1989).
\bibitem{Ress} Ressell, M.T., {\it et al.}, {\it Phys. Rev.} {\bf D 48},
 5519 (1993);
%M.T. Ressell and D.J. Dean, hep-ph$/9702290$.
%\bibitem{Dimit}V.I. Dimitrov, J. Engel and S. Pittel, 
%{\it Phys. Rev.} {\bf D 51}, R291 (1995).
%\bibitem{Engel}J. Engel, {\it Phys. Lett.} {\bf B 264}, 114 (1991).
%\bibitem{Nikol}M.A. Nikolaev and H.V. Klapdor-Kleingrothaus, {\it Z. Phys. }
%{\bf A 345}, 373 (1993); {\it Phys. Lett} {\bf 329 B}, 5 (1993);
%{\it  Phys. Rev.} {\bf D 50}, 7128 (1995).
\bibitem{KVdubna} Vergados, J.D. and  Kosmas, T.S. {\it Physics of Atomic
 nuclei}, Vol. {\bf 61}, No 7, 1066 (1998) 
(from {\it Yadernaya Fisika}, Vol. 61, No 7, 1166 (1998).
\bibitem{DIVA00} Divari, P.C., Kosmas, T.S., Vergados, J.D. and Skouras, L.D.,
{\it  Phys. Rev.} {\bf C 61} (2000), 044612-1.
\bibitem{Verg98}Vergados, J.D., {\it Phys. Rev.} {\bf D 58}, 103001-1 (1998).
\bibitem{Verg99}Vergados, J.D., {\it Phys. Rev. Lett} {\bf  83}, 3597
(1999).
\bibitem{Verg00}Vergados, J.D., {\it Phys. Rev.} {\bf D 62}, 023519 (2000).
\bibitem{Verg01} Vergados, {J.D., \it Phys. Rev. } {\bf D 63}, 06351 (2001).
\bibitem{UKDMC}Buckland, K.N., Lehner, M.J. and  Masek, G.E., in 
 Proc. {\it 3nd Int. Conf. on Dark Matter
in Astro- and part. Phys.} (Dark2000), Ed. . Klapdor-Kleingrothaus, H.V.,
Springer Verlag (2000).
%\bibitem{mb} S. Mart\'{\i} i Grac\'{\i}a, J. Fuster and 
%S. Cabrera, \npps{64}{1998}{376}.
\bibitem{VALLS} {\it CDF Collaboration}, FERMILAB-Conf-99/263-E CDF;\\
http://fnalpubs.fnal.gov/archive/1999/conf/Conf-99-263-E.html.
%\bibitem{VALDOR} CDF Collaboration, FERMILAB-Conf-99/263-E CDF;
%http://fnalpubs.fnal.gov/archive/1999/conf/Conf-99-263-E.html.
\bibitem{DORMAN}
 Dorman, P.J., {\it ALEPH Collaboration} March 2000,
http://alephwww.cern.ch/ALPUB/seminar/lepc$_{}$mar200/lepc2000.pdf.
\bibitem{JDV02}
Vergados, J.D., SUSY Dark Matter in Universe- Theoretical Direct
Detection Rates, 
 Proc. {\it NANP-01, International Conference on Non Accelerator New Physics},
Dubna, Russia, June 19-23, 2001, Editors V. Bednyakov and S. Kovalenko,
.
%\bibitem{JDVb02} Vergados, J.D., to be published.
\end{thebibliography}
\end{document}

