Message: Li6+ n(thermal) -> alpha + triton: incorrect final states (extra gammas ?) Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question Li6+ n(thermal) -> alpha + triton: incorrect final states (extra gammas ?) 

Forum: Hadronic Processes
Date: 20 Sep, 2009
From: Yuri Shitov <Yuri Shitov>

Simulating thermal neutrons in Li-6 loaded material I have expected well known reaction:

n + Li6 -> H-3(2.75 MeV) + He-4 (2.05 MeV)

But I have found that G4 (ver. 4.9.2.p01 + G4NDL3.13, (QGSP_BIC_HP, LHEP_PRECO_HP, LBE phys lists) ) simulates strange list of secondary particles:

*********************************************************************************************************
* G4Track Information:   Particle = neutron,   Track ID = 1,   Parent ID = 0
*********************************************************************************************************

Step#    X(mm)    Y(mm)    Z(mm) KinE(MeV)  dE(MeV) StepLeng TrackLeng  NextVolume ProcName
   55      417     6.69      465  1.54e-06        0        0       141 scinstripcore.phys initStep
   56      415     11.6      463  1.54e-06        0     5.58       146 scinstriprefl.phys Transportation
   57      415     11.8      463  1.54e-06        0    0.273       147 scinstriprefl.phys Transportation
   58      415     12.1      463  1.54e-06        0    0.273       147 scinstripcore2.phys Transportation
   59      415     12.9      463         0        0        1       148 scinstripcore2.phys NeutronInelastic
MyUserSteppingAction.UserSteppingAction: Current particle: neutron: N-interaction in path='/vessel.phys_0/chamber.phys_0/scinstrip.phys[93,39]_39093/scins\
triprefl.phys_0/scinstripcore.phys_0/scinstripcore2.phys_0' (defined material trigger=4) at position (X=414.532, Y=12.9458, Z=462.992), Strip[93,39], capt\
ure time=1900.3 (global), 1900.3 (local) ns
MyUserSteppingAction.UserSteppingAction: process: NeutronInelastic: list of secondary particles:
  triton: 2.81048e+06 keV, time=1900.3 (global), 1900.3(local) ns
  alpha: 3.72855e+06 keV, time=1900.3 (global), 1900.3(local) ns
  gamma: 3.20305 keV, time=1900.3 (global), 1900.3(local) ns
  gamma: 1025.31 keV, time=1900.3 (global), 1900.3(local) ns
  gamma: 1025.31 keV, time=1900.3 (global), 1900.3(local) ns

How can I correct this? Problem is also confirmed by other user, see thread http://hypernews.slac.stanford.edu/HyperNews/geant4/get/phys-list/530.html It was proposed there to shift this discussion to this category of forum.

Cheers, Yuri.

Inline Depth:
 1 1
 All All
Outline Depth:
 1 1
 2 2
 All All
Add message: (add)

1 Sad: Re: Li6+ n(thermal) -> alpha + triton: incorrect final states (extra gammas ?)   (Yuri Shitov - 20 Sep, 2009)
2 None: Re: Li6+ n(thermal) -> alpha + triton: incorrect final states (extra gammas ?)   (Calvin Giles - 22 Sep, 2009)
 Add Message Add Message
to: "Li6+ n(thermal) -> alpha + triton: incorrect final states (extra gammas ?)"

 Subscribe Subscribe

This site runs SLAC HyperNews version 1.11-slac-98, derived from the original HyperNews


[ Geant 4 Home | Geant 4 HyperNews | Search | Request New Forum | Feedback ]