in the activation measurements using the nuclear reactor , the neutron capture reaction rates are the product of the neutron flux at the sample position and the neutron capture cross - section of irradiation nuclide .
the neutron capture cross - section at particular neutron energy is the probability that a neutron at that energy will be captured by the nuclide .
the standard energy for tabulation of thermal neutron cross - section ( 0 ) is that of room temperature ( 293.6k or 20.43  c ) , corresponding to a neutron energy 0.0253  ev or to a neutron velocity of 2200  m / s .
however , in many reactors , the temperature ( t ) in the irradiation channels is not 20.43  c , but at higher temperature .
westcott developed a method for converting 0 to (t ) , which is the effective cross - section at the actual temperature ( t ) of the irradiation channel , by describing the neutron spectrum as a combination of a maxwellian distribution function , which is characterized by a temperature ( t ) , and a component of epithermal energy neutrons .
according to westcott , the reaction rate per target atom is given by:1\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ r = \phi_{0 } \mathop \sigma \limits^{\uplambda } ( t ) = nv_{0 } \sigma_{0 } \left [ { g(t ) + r\sqrt { t / t_{0 } } s_{0 } } \right ] $ $ \end{document}where 0  =  nv0 neutron flux defined as the total neutron density times the 2200  m / s neutron velocity ; \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \mathop \sigma \limits^{\uplambda } ( t ) $ $ \end{document } effective cross - section ; 0 thermal cross - section for 2200  m / s neutron ; g(t ) parameter which represents the departure of the cross - section from the 1/v law in the thermal region ( g(t )  =  1 if nuclide obeys the 1/v law in this energy region ) and which can be calculated from expression:2\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ g(t ) = { \frac{2}{{\sqrt { \pi e_{0 } } \sigma_{0 } } } } \int\limits_{0}^{\infty } { \sqrt e } \sigma ( e)\sqrt { e / e_{t } } e^ { { - e / e_{t } } } de / e_{t } $ $ \end{document}where \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ e_{t } = e_{0 } t / t_{0 } $ $ \end{document } , \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ e_{0 } = 0.0253 $ $ \end{document } ev , \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ t_{0 } = 293.6 $ $ \end{document}k , e  neutron energy , and v  neutron velocity .
c ) to 673.16  k ( 400  c ) for many nuclides were calculated by holden .
\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ r\sqrt { t / t_{0 } } $ $ \end{document}epithermal index which denotes the strength of the epithermal flux ( it is zero for a pure thermal flux ) and s0parameter which represents the ratio of the resonance integral and thermal cross - section .
according to the hogdahl s convention , a method , theoretically less rigorous , but commonly used defines the reaction rate by:3\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ r = \phi_{th } \phi_{th } + \phi_{epi } i $ $ \end{document}where th , i sub - cadmium cross - section and the epi - cadmium cross - section , respectively ; th neutron flux defined as the thermal neutron density times the 2200  m / s neutron velocity ; epi epithermal flux per unit ln(e ) . from eqs .  1 and 3 it follows that ( see ) : 4\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				
\usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \sigma_{th } = g(t)\sigma_{0 } + { \frac{{\phi_{epi } } } { { \phi_{th } } } } ( 1 - 2\sqrt { e_{0 } - e_{c } } ) g(t)\sigma_{0 } $ $ \end{document}or5\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \sigma_{th } = g(t)\sigma_{0 } + { \frac{{\phi_{epi } } } { { \phi_{th } } } } [ \updelta i ' - i(1/v ) ] $ $ \end{document}with6\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ i(1/v ) = \int\limits_{{e_{c } } } ^{\infty } { g(t)\sigma_{0 } \sqrt { e_{0 } /e } } de / e \cong 2g(t)\sigma_{0 } \sqrt { e_{0 } /e_{c } } $ $ \end{document}where ec  cd - cut - off energy ( this cut - off energy is a function of thickness of cd ; ec  =  0.55  ev if a small cylindrical cd - filter of 1  mm thickness ) , i(1/v)1/v - part of resonance integral , i and i is the epi - cadmium resonance integral excluding the 1/v part and part , shielded by cd - filter , which depends on the neutron temperature , respectively .
clearly , looking at eqs .  1 , 3 , 4 , 5 , and 6
there is one notice that the cross - sections and the resonance integral vary as a function of g(t ) , which are temperature dependence .
many authors often assume th  =  0 to calculate the resonance integral , it is not completely justified . only the nuclides with cross - section
obey the 1/v law up to 12  ev , g(t ) would be unity as mentioned above and s0  =  0 ; then th and the resonance integral are no longer depends on temperature at the irradiation position .
thus , in the activation measurements using the nuclear reactor , neutron temperature ( t ) at the irradiation position should be known to preserve the accuracy of the results . in experiment , a well - known method for the determination of the neutron temperature in the irradiation channels is co - irradiation of lu by using lu ( n ,  ) lu reaction and 1/v - monitors   with acceptable accuracy . in this work ,
we describe a method for the determination of the neutron temperature based on thermal neutron spectrum calculated by mcnp code .
as we know , neutron spectrum at an irradiation position can be expressed as a sum of a thermal equilibrium spectrum ; th(e ) ( maxwellian distribution , neutron energy from 0 to 0.55  ev ) and epithermal spectrum in the slowing down region ; epi(e ) ( neutron energy from 0.55  ev to maximal ) . in the paper [ 6 , 7 ]
, the neutron spectra were written by semi - empirical functions:7\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \varphi_{th } \left ( e \right ) = i_{th } \left ( { { \frac{e}{{e_{t}^{2 } } } } } \right)e^ { { - { \frac{e}{{e_{t } } } } } } $ $ \end{document}8\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \varphi_{epi } \left ( e \right ) = i_{epi } \left ( { \frac{e}{1ev } } \right)^ { - \alpha } .\frac{1}{e }
$ $ \end{document}where ith and iepi are scaling constants for the thermal and epithermal portion ; et is a characteristic energy of maxwellian portion of the spectrum .
the neutron temperature is given by:9\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ e_{t } ( ev ) = 8.617 \times 10^ { - 5 }
from the temperature in degrees kelvin , we can convert to temperature in degrees celsius by : t273.15 .
thus , if the neutron spectrum at the irradiation position is known , et would be determined by fitting maxwellian portion of the spectrum according to eq . 
this method was applied to determinate neutron temperature in the irradiation channels 7 - 1 , 1 - 4 , the rotary specimen rack and neutron trap of dalat reactor , vietnam .
it is a general purpose monte carlo code , which facilitates independent or coupled neutron , photon and electron transport calculations .
the code treats an arbitrary three - dimensional configuration of material and geometric cell and provides a versatile description of the source , the variance reduction techniques , a flexible tally structure and an extensive collection of cross - section data in continuous energy representation . for neutron ,
all reactions given in a particular cross - section data evaluation are accounted for and cover the energy range from 10 to 20  mev .
used mcnp version 4b for the calculation of spectral parameters in typical irradiation channels of the triga reactor , ljubljuna .
the comparison between the simulated and experimentally determined parameter f is in good agreement except for channel ic40 .
[ 6 , 7 ] also used mcnp for calculation of some characteristic parameters of neutron flux in irradiation channels of dalat reactor . the comparison of results obtained by experimental and monte carlo simulation showed good agreement for all channels .
the obtained results in that papers showed that monte carlo calculation of neutron spectra in irradiation channels of reactor provide useful information about the nature of spectra , the neutron flux , as well as reactor s spectral characteristics .
the cylindrical reactor is arranged symmetrically by several kinds of material elements . in the center of active region
the fuel elements , beryllium and graphite blocks , control rods and irradiation channels are the active region of reactor and form a hexagonal prism lattice in the cylindrical reactor with radius 20.8  cm and height 60.0  cm . outside of the active region is a graphite reflector with thickness of 32.6  cm .
.  1.fig .  1configuration of reactor core of dalat reactor configuration of reactor core of dalat reactor the detailed mcnp computation model used for the calculation of neutron spectra in irradiation channels of 7 - 1 , 1 - 4
, the rotary specimen rack and neutron trap of dalat reactor is based on a model thoroughly described by hung et al .
kcode ( criticality problem ) and ksrc ( initial spatial distribution of source points ) cards are used for the source description . in this case
, mcnp offers several options for the keff values estimate , including the combined average of the absorption , collision and track - length estimator . for all calculations of neutron spectra , the f4 tally was used .
f4 is a track - length tally card for calculating average neutron flux over a cell .
the s (  ,  ) thermal treatment was used for h2o , beryllium and graphite .
it is a general purpose monte carlo code , which facilitates independent or coupled neutron , photon and electron transport calculations .
the code treats an arbitrary three - dimensional configuration of material and geometric cell and provides a versatile description of the source , the variance reduction techniques , a flexible tally structure and an extensive collection of cross - section data in continuous energy representation . for neutron ,
all reactions given in a particular cross - section data evaluation are accounted for and cover the energy range from 10 to 20  mev .
used mcnp version 4b for the calculation of spectral parameters in typical irradiation channels of the triga reactor , ljubljuna .
the comparison between the simulated and experimentally determined parameter f is in good agreement except for channel ic40 .
[ 6 , 7 ] also used mcnp for calculation of some characteristic parameters of neutron flux in irradiation channels of dalat reactor . the comparison of results obtained by experimental and monte carlo simulation showed good agreement for all channels .
the obtained results in that papers showed that monte carlo calculation of neutron spectra in irradiation channels of reactor provide useful information about the nature of spectra , the neutron flux , as well as reactor s spectral characteristics .
the cylindrical reactor is arranged symmetrically by several kinds of material elements . in the center of active region
the fuel elements , beryllium and graphite blocks , control rods and irradiation channels are the active region of reactor and form a hexagonal prism lattice in the cylindrical reactor with radius 20.8  cm and height 60.0  cm . outside of the active region is a graphite reflector with thickness of 32.6  cm .
1.fig .  1configuration of reactor core of dalat reactor configuration of reactor core of dalat reactor the detailed mcnp computation model used for the calculation of neutron spectra in irradiation channels of 7 - 1 , 1 - 4
, the rotary specimen rack and neutron trap of dalat reactor is based on a model thoroughly described by hung et al .
kcode ( criticality problem ) and ksrc ( initial spatial distribution of source points ) cards are used for the source description . in this case , mcnp offers several options for the keff values estimate , including the combined average of the absorption , collision and track - length estimator . for all calculations of neutron spectra , the f4 tally was used .
f4 is a track - length tally card for calculating average neutron flux over a cell .
the s (  ,  ) thermal treatment was used for h2o , beryllium and graphite .
therefore , this value is quite satisfactory . the energy spectra with 115 groups of energy from 0  ev to 20  mev at 1 - 4 ,
7 - 1 channels and rotary specimen rack and neutron trap are calculated . as an example , fig . 
2 shows a calculated energy spectrum in 7 - 1 channel and the fitted et - value in the thermal neutron region .
the relative error of the calculated spectrum intensity in thermal neutron region is about than 1% and in the epithermal is about 3% for all irradiation channels.fig . 
2calculated neutron flux spectrum in 7 - 1 channel of dalat reactor using mcnp4c code ( dark 
 line ) and fitted et - value in the thermal neutron region ( light 
 line ) calculated neutron flux spectrum in 7 - 1 channel of dalat reactor using mcnp4c code ( dark 
 line ) and fitted et - value in the thermal neutron region ( light 
 line ) from the neutron spectra calculated by mcnp code , the et values are determined by fitting the neutron spectra in the thermal neutron region ( neutron energy from 0 to 0.55  ev ) according to eq .
the et - fitting values and the temperatures in degrees celsius in the irradiation channels are presented in table  1 .
the uncertainty of the neutron temperature values determined from fitting the calculated energy spectra in the thermal neutron region is about 2.5% ; it defends on the uncertainty of the calculated spectrum intensity in thermal neutron region.table  1the fitted et values and the neutron temperatures in irradiation channels of dalat reactor determined from the neutron spectra calculated by mcnp code and experimentirradiation channelsfitted et- values ( 10 ev)calculated temperatures (  c)experimental temperatures (  c)lu  au monitorneutron trap2.714    0.02841.8    1.01 - 4 channel2.727    0.02943.3    1.07 - 1 channel2.827    0.03152.0    1.053.0    2.0rotary specimen rack2.698   
0.02740.0    1.0 the fitted et values and the neutron temperatures in irradiation channels of dalat reactor determined from the neutron spectra calculated by mcnp code and experiment as a comparison to fitted value , an experiment of the t - determination in 7 - 1 channel was performed .
the choice of 7 - 1 channel is suitable for the experiment due to 7 - 1 channel coupled with a pneumatic sample transfer system .
the t - determination method was presented in paper . after measuring the irradiated monitors , the westcott factor for lutetium [ g ( t)]lu could be obtained .
thus , the neutron temperature also was determined by using the tables [ g ( t)]lu vs. tn in paper .
the comparison of the result of the neutron temperature in 7 - 1 channel shows that the value of the neutron temperature agrees well with each other .
, it also shows that the neutron temperatures in the irradiation channels of dalat reactor are from 40  c ( in the rotary specimen rack ) to 53  c ( 7 - 1 channel ) in the core of reactor . the neutron temperature in 7 - 1 channel is highest , because it is a dry channel using for neutron activation analysis ( naa ) only . note that , in general , the recommended k0 , au using in the neutron activation analysis in reports is the average k0 , au - value corresponding to an average neutron temperature of 60  c .
hence , in k0 naa standardization method , the neutron temperature at the irradiation position is importance for nuclides having ( n, ) reaction which do not exactly obeys the (v )  ~  1/v law up to 12  ev . the cross - section of ( n, ) reaction of these nuclides depends clearly on the neutron temperature .
for this reason , 7 - 1 channel is suitable in naa using k0 standardization method , because there is no need of the correction of neutron temperature .
the neutron temperature values in the channels of 7 - 1 , 1 - 4 , the rotary specimen rack and neutron trap of dalat reactor were determined from the neutron spectra in the thermal neutron region , which was calculated by mcnp code .
it was compared with the experimental value in 7 - 1 channel and showed that the values are agreement with each other .
the results from this study suggest that base on the neutron spectrum calculated by mcnp code , the neutron temperature in the irradiation channels of reactor will be determined by fitting the spectrum in the thermal neutron region according to eq . 
it is simple , effective and is a good method for the verification of experimental results .
 the neutron temperature is a characteristic parameter in irradiation channels of reactor . for nuclides which have resonances in the thermal neutron energy range , 
 their westcott g - factors are different from unity . 
 the values of g - factors and cross - sections of ( n ,  ) reaction of these nuclides are temperature dependence . 
 the standard energy for tabulation of thermal neutron cross - section ( 0 ) is that of room temperature ( 293.59  k or 20.43  c ) , corresponding to a neutron energy 0.0253  ev or to a neutron velocity of 2200  m / s . 
 however , in the irradiation channels of reactor , the temperature is not exact at 20.43  c . 
 thus , the temperature at the irradiation position must be known to convert 0 to (t ) . 
 a method for determination of the neutron temperature in irradiation channels of dalat reactor is presented by fitting the thermal neutron spectrum obtained from the calculation using mcnp code . 
Introduction
Base of method
MCNP code and configuration of Dalat reactor core
MCNP code
Configuration of Dalat reactor core
Results and discussion
Conclusion