Message: Re: neutron beam hit Zinc target Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question Re: neutron beam hit Zinc target 

Forum: Hadronic Processes
Re: Question neutron beam hit Zinc target (Tuoyu Cao)
Re: None Re: neutron beam hit Zinc target (Michael H. Kelsey)
Date: 01 Mar, 2013
From: Tuoyu Cao <Tuoyu Cao>

Hi Michael,

Thanks for reply.

Here is the code I use to define Zinc.

G4Isotope* iso_Zn64 = new G4Isotope("Zn64", iz=30, ia=64, density= 64*g/mole);
  G4Isotope* iso_Zn66 = new G4Isotope("Zn66", iz=30, ia=66, density= 66*g/mole);
  G4Isotope* iso_Zn67 = new G4Isotope("Zn67", iz=30, ia=67, density= 67*g/mole);
  G4Isotope* iso_Zn68 = new G4Isotope("Zn68", iz=30, ia=68, density= 68*g/mole);
  G4Isotope* iso_Zn70 = new G4Isotope("Zn70", iz=30, ia=70, density= 70*g/mole);

  G4Element* ele_Zn = new G4Element("ele Zn", symbol="eZn",ncomponents=5);

  ele_Zn->AddIsotope(iso_Zn64,abundance=48.60*perCent);
  ele_Zn->AddIsotope(iso_Zn66,abundance=27.90*perCent);
  ele_Zn->AddIsotope(iso_Zn67,abundance=04.10*perCent);
  ele_Zn->AddIsotope(iso_Zn68,abundance=18.80*perCent);
  ele_Zn->AddIsotope(iso_Zn70,abundance=00.60*perCent);

  G4Material* mat_ele_Zn = new G4Material("mix_Zn", density= 7.14*g/cm3, ncomponents=1);
  mat_ele_Zn->AddElement(ele_Zn, fractionmass =1);

But G4 uses nat_zinc instead. What should I do if I want to simulate one isotope individually?

The printout:

*********************************************************************************************************
* G4Track Information:   Particle = neutron,   Track ID = 1,   Parent ID = 0
*********************************************************************************************************
Step#      X         Y         Z        KineE    dEStep   StepLeng  TrakLeng    Volume     Process
    0     -5 mm   2.72 mm  -4.68 mm    2.5 MeV     0 eV      0 fm      0 fm   Absorber   initStep
    1   3.57 mm   2.72 mm  -4.68 mm   2.48 MeV     0 eV   8.57 mm   8.57 mm   Absorber   hadElastic
    :----- List of secondaries ----------------
       Zn65[0.0]:  energy =  16.3 keV  time = 392.5 ps 
    :------------------------------------------
    2   4.01 mm   2.58 mm     -5 mm   2.48 MeV     0 eV    565 um   9.13 mm   Absorber   Transportation
    3   5.42 mm   2.14 mm     -6 mm   2.48 MeV     0 eV   1.78 mm   1.09 cm      World   OutOfWorld

shows there is only one secondary, which is Zn65, and it seems like it comes from hadElastic.... Which is very confusing. Is this an error in G4. Could you please explain this.

Thanks a lot.

Tuoyu

On Thu, 28 Feb 2013 06:00:27 GMT, Michael H. Kelsey wrote:

> On Thu, 28 Feb 2013 03:56:02 GMT, Tuoyu Cao wrote:
> 
> > I used a 2.5 MeV neutron beam to hit a zinc target and analyse the
> > gammas. I defined zinc material based on zinc isotopes.
> 
> And what was your material definition? Natural zinc (material G4_Zn from
> the so-called "NIST" database) is 48.63% Zn-64, 27.9% Zn-66, 18.75%
> Zn-68, and a little bit of other isotopes.
> 
> > I saw the biggest peak in the
> > gamma spectrum is on 1115.5 keV, which comes from Zn65->Cu65+gamma.
> 
> From reading your tracking dump:
> 
> > * G4Track Information: Particle = Zn65[0.0], Track ID = 2, Parent ID =
> 1
> 
> The Zn65 came from the initial neutron, either by capture on Zn-64 or by
> spallation (n,2n) on Zn-66. Did that event contain a second neutron?
> 
> In other events, you should have also seen some Zn-67 and Zn-69
> production in the same way.
> 
>   -- Michael Kelsey
> 

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

1 None: Re: neutron beam hit Zinc target   (Mike Kelsey - 01 Mar, 2013)
(_ Question: Re: neutron beam hit Zinc target   (Tuoyu Cao - 01 Mar, 2013)
1 None: Re: neutron beam hit Zinc target   (Mike Kelsey - 01 Mar, 2013)
2 None: Re: neutron beam hit Zinc target   (Koi, Tatsumi - 01 Mar, 2013)
 Add Message Add Message
to: "Re: neutron beam hit Zinc target"

 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 ]