From tsushima@kn3.physast.uga.edu Mon Sep 22 17:31:23 2003
Date: Wed, 03 Sep 2003 15:43:39 -0500 (CDT)
From: Tsushima <tsushima@kn3.physast.uga.edu>
To: Arie Bodek <bodek@pas.rochester.edu>
Cc: Dr Koichi Saito <ksaito@tohoku-pharm.ac.jp>, hung@phya.yonsei.ac.kr,
     gzeller@fnal.gov, hbudd@fnal.gov, ksmcf@pas.rochester.edu,
     manly@pas.rochester.edu
Subject: code for bound/free ffs. 12C and 56Fe


Dear Prof. A. Bodek (and others (sorry to call this!)):

I send you an interpolation code and input files to get bound/free 
form factor ratios calculated in the QMC model by attachment.

If you have any questions, please inform me.

Sincerely,  Kazuo Tsushima. 

P.S Some of the parameters used are, current quark mass m_q=5 MeV and 
the bag radius of nucleon R=0.8 fm, both in free space.


On Mon, 25 Aug 2003, Arie Bodek wrote:

> Greetings
>    Thanks for the answers to my questions
> I have a couple of more requets
> 
> 1, In your paper, you give the ratio of Bound/Free ratios for
> F1, F2 and FA versus Q2.  Can you provide numbers for
> ratio of Bound/Free in Carbon and Iron for
> Gep, Gen,  Gmp,  Gmn  separately versus Q2 (or a computer program)
> 
> 2. I understand that your model is not valid for Q2 gt 1.5 where
>   one expects the ratio of bound/free form factors to be 1.0.   If one
> were to implement your ratio in our model (to predict bound rather
> than free) neutrino cross section. Therefore, can you provide
> your best estimate versus Q2 as to the ratio of bound/free that goes 
> to 1.0 at large Q2.
> 
> Best regards
>   Arie Bodek
> 
> At 3:52 PM +0900 7/7/03, Dr Koichi Saito wrote:
> >Dear Prof. A. Bodek,
> >
> >Thanks a lot for your interests to our paper.
> >We below attach a reply, and we hope this can answer to your
> >questions.
> >
> >Sincerely, 
> >K. Saito, Hungchong Kim, K. Tsushima,
> >
> >================================
> >    A/Prof. Koichi Saito
> >    Tohoku Pharmaceutical University
> >    Sendai, 981-8558, JAPAN
> >    ksaito@tohoku-pharm.ac.jp
> >    (or ksaito@nucl.phys.tohoku.ac.jp)
> >    tel: japan+22-234-4181 (ext 3201)
> >    fax: japan+22-275-2013
> >================================
> >
> >%%%%%
> >>
> >>  Greetings
> >>
> >>  I have read with interest your recent preprint below:
> >>
> >>  http://arxiv.org/abs/nucl-th/0307013
> >>
> >>  I have several questions.
> >>
> >>  1. I am wondering if you tested if your model satisfies the Coulomb sum
> >rule
> >>  from closure  (integral of nucleus elastic and inelastic form factors).
> >>
> >
> >As Kazuo mailed to you last week, we have already calculated the
> >Coulomb sum rule based on the relativistic RPA (RRPA) approach plus
> >the bound nucleon form factors (P.L. B465 (1999) 27-35).
> >We then found that the sum rule is reduced by about 20% (relative to the
> >standard relativistic Hartree result), which originates in the usual RPA
> >correlation and the bound nucleon form factors, F1 and F2.
> >This reduction is of the right magnitude to understand the recent
> >experimental analysis (P.L. B515, 269).
> >
> >>  2 In order to test this, on needs to include the Pauli suppression
> >>  factor, which
> >>  makes the quasielastic scattering zero at Q2=0.
> >>
> >
> >Since in RRPA  the Pauli exclusion effect is fully taken into account
> >in the nucleon propagator used in the intermediate states,  we believe, in
> >addition to it, it is not necessary to apply another "Pauli suppression
> >factor". Furthermore, at Q^2 (=q_0^2 - {\vec q}^2)=0 the cross section
> >is certainly zero in our neutrino result (see for example, Figs.3-6 in
> >P.R.C51 (1995) 2739 by Hungchong Kim et al.).
> >Note that at Q^2=0 the modification of the vector form factor due to
> >the nuclear medium should vanish (even in matter) because of baryon
> >number (or charge) conservation.
> >
> >>  3. Since the contributions to the
> >>  quasielastic  cross section of  of F1, F2 and FA are Q2 dependent,
> >>  and since the Q2=0 part is zero from Pauli suppression, I wonder how
> >>  including this effect would change your conclusions about the magnitude
> >>  of the reduction in the cross section.
> >>
> >
> >As we have mentioned above, our calculation is based on RRPA (for
> >Coulomb sum rule) or relativistic Fermi gas model (for the present
> >neutrino calculation), in which the Pauli exclusion effect and the Q^2
> >dependence of the form factors are fully taken into account.
> >So, we believe that our present conclusion on the neutrino reaction does
> >not change.
> >
> >Thank you very mcuh for your interests in our work.
> >
> >%%%%%
> 

    [ Part 2, "main code"  Text/PLAIN (Name: "qmc_CFe_ffratio.f")  153 ]
    [ lines. ]
    [ Unable to print this part. ]


    [ Part 3, "12C input data file"  Text/PLAIN (Name: ]
    [ "GEMAmed_qmcrat_c12.indat")  6,136 lines. ]
    [ Unable to print this part. ]


    [ Part 4, "56Fe input data file"  Text/PLAIN (Name: ]
    [ "GEMAmed_qmcrat_fe56.indat")  6,136 lines. ]
    [ Unable to print this part. ]

