There is no neutron production in simulation for DD reaction in Hadr03

Physics Processes, Models and Cross Sections
1

hi dear experts,

I’ve been trying to simulate DD and DT fusion with Hadr03 Example.
despite the approximately same cross-section between DD and DT fusion at 900keV, there is no neutron production in simulation for DD reaction.

this is the output for DT fusion

The run is 100000 deuteron of 900 keV through 100 km of tritium (density: 5 mg/cm3)

Process calls frequency:
dInelastic= 100000

MeanFreePath: 37.887 m ± 38.049 m massic: 18.943 g/cm2
CrossSection: 0.00026395 cm^-1 massic: 5.2789 mm2/g
crossSection per atom: 264.41 mbarn

Verification: crossSections from G4HadronicProcessStore:
dInelastic= 5.292 mm2/g 265.07 mbarn
total= 5.292 mm2/g 265.07 mbarn

List of nuclear reactions:

                    deuteron + tritium --> neutron + alpha:  100000   Q =   17.59 MeV

List of generated particles:
alpha: 100000 Emean = 3.9554 MeV ( 1.9056 MeV → 6.01 MeV)
neutron: 100000 Emean = 14.534 MeV ( 12.48 MeV → 16.584 MeV)

and this output for DD fusion

The run is 100000 deuteron of 900 keV through 100 km of deuteron (density: 5 mg/cm3)

Process calls frequency:
dInelastic= 100000

MeanFreePath: 40.176 m ± 40.348 m massic: 20.088 g/cm2
CrossSection: 0.00024891 cm^-1 massic: 4.9782 mm2/g
crossSection per atom: 166.49 mbarn

Verification: crossSections from G4HadronicProcessStore:
dInelastic= 4.9905 mm2/g 166.91 mbarn
total= 4.9905 mm2/g 166.91 mbarn

List of nuclear reactions:

                        deuteron + deuteron --> 2 deuteron:   94661   Q = 0.00025005 meV
             deuteron + deuteron --> N gamma or e- + alpha:    5339   Q =  23.847 MeV

                      number of gamma or e- (ic): N = 1 --> 2

List of generated particles:
alpha: 5339 Emean = 519.56 keV ( 151.2 keV → 900.77 keV)
deuteron: 189322 Emean = 450 keV ( 20.53 eV → 899.98 keV)
gamma: 6386 Emean = 20.255 MeV ( 4.0223 MeV → 24.596 MeV)

can Geant4 simulate DD fusion?
am I should use another physicsList for DD fusion?

this is cross-section of these fusions

Q
Q558×560 23.2 KB

Ongena, J., 2015. Fusion: A true challenge for an enormous reward. In EPJ Web of Conferences (Vol. 98, p. 05004). EDP Sciences.

I’m using Geant4 11.0

2

I run the same simulation but with QGSP_BIC_AllHP as phisycList for DT reaction and here is the output :

The run is 10000000 deuteron of 900 keV through 1 km of H3 (density: 5 mg/cm3 )

Process calls frequency:
dInelastic= 4592968 hadElastic= 5407032

MeanFreePath: 1.1866 cm ± 1.1861 cm massic: 8.6737 g/cm2
CrossSection: 0.84278 cm^-1 massic: 11.529 mm2/g
crossSection per atom: 577.47 mbarn

Verification: crossSections from G4HadronicProcessStore:
dInelastic= 5.292 mm2/g 265.07 mbarn
hadElastic= 6.2369 mm2/g 312.4 mbarn
total= 11.529 mm2/g 577.47 mbarn

List of nuclear reactions:

                       deuteron + H3 --> deuteron + triton: 5407032   Q = 0.00031846 meV
                         deuteron + H3 --> neutron + alpha: 4592968   Q =   17.59 MeV

List of generated particles:
alpha: 4592968 Emean = 3.9581 MeV ( 1.9055 MeV → 6.0101 MeV)
neutron: 4592968 Emean = 14.531 MeV ( 12.48 MeV → 16.584 MeV)
triton: 5407032 Emean = 422.27 keV ( 447.98 meV → 864.28 keV)

Momentum balance: Pmean = 0.027052 meV ( 0 eV → 34.528 meV)

Thibault

3

Hi Nilson,

I happened to notice this thread has been left unanswered. I am also trying to simulate DD and DT fusion using the Hadr03 example. Did you ever figure out why there was no neutron production from the DD reaction? Is GEANT4 capable of simulating DD fusion? Every thread concerning this topic has no response.

Please let me know if you have figured it out!

4

Nilson,

Just wanted to follow up to inform you that I believe I found the solution. You need to download the TENDL datafiles and set/create the G4PARTICLEHP environment variable to the dataset. The physics list you need to use is G4HadronPhysicsQGSP_BIC_AllHP and I also needed to use G4IonPhysicsPHP in order to get this to work.

List of nuclear reactions: 
                deuteron + Deuterium Isotope --> neutron + He3:    3850   Q =   3.271 MeV
              deuteron + Deuterium Isotope --> proton + triton:    3545   Q =  4.0359 MeV
                deuteron + Tritium Isotope --> neutron + alpha:  992605   Q =  17.631 MeV

List of generated particles:
            He3:    3850  Emean =  871.83 keV	( 500.31 keV --> 1.2552 MeV)
          alpha:  992605  Emean =   3.619 MeV	( 2.898 MeV --> 4.3413 MeV)
        neutron:  996455  Emean =  14.081 MeV	( 2.1298 MeV --> 14.847 MeV)
         proton:    3545  Emean =  3.0796 MeV	( 2.6627 MeV --> 3.4997 MeV)
         triton:    3545  Emean =  1.0703 MeV	( 650.22 keV --> 1.4871 MeV)

Hope this helps,

  • Brendan