Problem 2162 - fusion reactions
Summary: fusion reactions
Status: RESOLVED FIXED
Alias: None
Product: Geant4
Classification: Unclassified
Component: physics_lists (show other problems)
Version: 10.5
Hardware: All All
: P4 minor
Assignee: Vladimir.Ivantchenko
URL:
Depends on:
Blocks:
 
Reported: 2019-04-24 11:51 CEST by michel.maire
Modified: 2020-11-06 10:32 CET (History)
3 users (show)

See Also:

Attachments
Add an attachment (proposed patch, testcase, etc.)

Note You need to log in before you can comment on or make changes to this problem.
Description michel.maire 2019-04-24 11:51:24 CEST 
Fusion reaction : deuteron + tritium can be simulated with G4IonPhysics but no more with G4IonPhysicsXS.
Same remark for dd and td reactions.

Here a macro for Hadr03 :

 /control/verbose 1
 /run/verbose 1
 #
 /testhadr/det/setIsotopeMat tritium 1 3 0.1 g/cm3
 /testhadr/det/setSize 20 m
 #
 /run/initialize
 #
 /gun/particle deuteron
 /gun/energy 80 keV
 #
 /process/inactivate hadElastic
 #
 /run/beamOn 10000
Comment 1 michel.maire 2019-09-02 15:25:00 CEST 
the problem is still here in ref-08
Comment 2 Vladimir.Ivantchenko 2019-11-14 19:42:05 CET 
Hi Michel,

 80 keV seems to be below Coulomb barrier, when I set 80 MeV all reactions are seen.

Vladimir
Comment 3 michel.maire 2019-11-15 16:49:54 CET 
Did you try with G4IonPhysics or G4IonPhysicsPHP ?
Both produce the well known neutron of 14.1 MeV (for example in the sun ...)

   Michel
Comment 4 Vladimir.Ivantchenko 2019-11-16 19:43:29 CET 
Hi Michel,

in all 3 cases similar neutrons are produced with 1.5 MeV primary energy. WHat is slightly worrying me that 3 cross sections very different.

Vladimir
Comment 5 Vladimir.Ivantchenko 2020-11-06 08:54:54 CET 
Hi Michel,

with the new macro it seems the problem is fixed.

Vladimir
Comment 6 michel.maire 2020-11-06 10:30:15 CET 
Yes, it seems ok.
I scanned few energies. Cross sections agreed with the plot send by Baptiste Guegan.

 Thanks,  Michel