Message: Serious problem with isotope definition (G4Isotope) and hadronic processes Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None Serious problem with isotope definition (G4Isotope) and hadronic processes 

Forum: Hadronic Processes
Date: 08 May, 2006
From: Benoit Pirard <benoit.pirard@cesr.fr>

I am trying to simulate the activation of Germanium caused by 50 MeV protons with Geant4.7.1. In order to study the sensitivity to the Ge isotopes (70Ge, 72Ge, 73Ge, 74Ge and 76Ge), I firstly defined a material with only the 70Ge isotope. Here is my Ge material definition:

 G4Isotope* Ge70 = new G4Isotope(name="Ge70",iz=32,n=70,a=69.92425*g/mole);                            
 G4Element* elGe = new G4Element(name="Germanium",symbol="Ge",ncomponents=1);                                      
 elGe->AddIsotope(Ge70,abundance=100.00*perCent);

 density = 5.323*g/cm3;
 temperature =95.0*kelvin;
 G4Material* MGe = new G4Material("Germanium",density, ncomponents=1, kStateSolid, temperature);
 MGe -> AddElement( elGe, nAtoms=1 );
 cout << *(G4Material::GetMaterialTable()) << endl;

which prints the following table as the run starts:

 ***** Table : Nb of materials = 1 *****
  Material: Germanium     density:  5.323 g/cm3   temperature:  95.00 K  pressure:   1.00 atm  RadLength:   2.216 cm
    --->  Element: Germanium (Ge)   Z = 32.0   N =  70.0   A =  69.92 g/mole
    --->  Isotope:  Ge70   Z = 32   N =  70   A =  69.92 g/mole   abundance: 100.00 %  fractionMass: 100.00 %  Abundance 100.00 %

Surprisingly I saw in the output files that isotopes with atomic mass (nucleon count) A = 77 could be created through 70Ge + p inelastic (or even elastic !!!) collisions. By looking at the target nuclei (with G4IsoParticleChange class and GetMotherNucleus() method), I saw that all the 5 Ge stable isotopes were taken with their natural abundances. Thus the definition of the Ge element I gave was not respected in the hadronic processes (eventhough the printed isotope table was correct).

Note that this problem is not specific to the element (Ge), the particle (proton) or the physics list I used. I tried simulations with neutrons and protons on enriched Si and U, but again the program did not take in account my isotopic abundance. From the quick look I had in the source files, it seems that the class G4HadronicProcess (method ChooseAandZ) is taking only element-wise information (I mean the atomic charge number Z) given in the user DetectorConstruction class and not isotope-wise information (relative abundances, etc.). Information on isotopes is taken from the class G4StableIsotopes which contains natural abundances ...

My question is: What is the reliable way to force an isotopic abundance other than the (default) natural abundance defined in the source files ?

Some messages related to the same problem were posted on the forum in the past, in particular posting #388 in Hadronic Processes forum and posting #117 in Physics List forum. No answer was however given to them.

Could you please help me to solve this problem ? Thank you in advance,

Benoît

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

1 More: Re: Serious problem with isotope definition (G4Isotope) and hadronic processes   (Marcus Mendenhall - 15 May, 2006)
3 None: Re: Serious problem with isotope definition (G4Isotope) and hadronic processes   (Dennis H. Wright - 15 May, 2006)
1 None: Re: Serious problem with isotope definition (G4Isotope) and hadronic processes   (Koi, Tatsumi - 15 May, 2006)
(_ None: Re: Serious problem with isotope definition (G4Isotope) and hadronic processes   (Vladimir IVANTCHENKO - 15 May, 2006)
 Add Message Add Message
to: "Serious problem with isotope definition (G4Isotope) and hadronic processes"

 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 ]