Time_dont_match

Physics Processes, Models and Cross Sections
1

Hi every one,
I’m trying to run scavenger example. At energies above 150 eV, the example runs well, but at energies less than 150 eV (I considered 75 eV) , threads are failing with the following fatal exception:

-------- EEEE ------- G4Exception-START -------- EEEE -------
*** G4Exception : TIME_DONT_MATCH
      issued by : G4MoleculeCounter::RemoveAMoleculeAtTime
Time of species e_aq^-1 is 1 ps  while  global time is 1 ns 

*** Fatal Exception *** core dump ***
Backtrace:
[PID=6821, TID=1][ 0/23]> /usr/lib/x86_64-linux-gnu/libc.so.6(pthread_kill+0x12c) [0x7feba0c96a7c]
[PID=6821, TID=1][ 1/23]> /usr/lib/x86_64-linux-gnu/libc.so.6(raise+0x16) [0x7feba0c42476]
[PID=6821, TID=1][ 2/23]> /usr/lib/x86_64-linux-gnu/libc.so.6(abort+0xd3) [0x7feba0c287f3]
[PID=6821, TID=1][ 3/23]> /home/ali/Geant4/Geant4_Install/lib/libG4global.so(_Z11G4ExceptionPKcS0_19G4ExceptionSeverityS0_+0x1117) [0x7feba14d5e87]
[PID=6821, TID=1][ 4/23]> /home/ali/Geant4/Geant4_Install/lib/libG4global.so(_Z11G4ExceptionPKcS0_19G4ExceptionSeverityRNSt7__cxx1119basic_ostringstreamIcSt11char_traitsIcESaIcEEE+0xa5) [0x7feba14d6025]
[PID=6821, TID=1][ 5/23]> /home/ali/Geant4/Geant4_Install/lib/libG4processes.so(_ZN17G4MoleculeCounter21RemoveAMoleculeAtTimeEPK24G4MolecularConfigurationdPKN5CLHEP10Hep3VectorEi+0x121c) [0x7feba28963ac]
[PID=6821, TID=1][ 6/23]> /home/ali/Geant4/Geant4_Install/lib/libG4processes.so(_ZN10G4MoleculeD1Ev+0x80) [0x7feba28971d0]
[PID=6821, TID=1][ 7/23]> /home/ali/Geant4/Geant4_Install/lib/libG4processes.so(_ZTv0_n24_N10G4MoleculeD0Ev+0x17) [0x7feba28973a7]
[PID=6821, TID=1][ 8/23]> /home/ali/Geant4/Geant4_Install/lib/libG4track.so(_ZN7G4TrackD1Ev+0x44) [0x7feba0bd5d14]
[PID=6821, TID=1][ 9/23]> /home/ali/Geant4/Geant4_Install/lib/libG4processes.so(_ZN10G4FastListI7G4TrackE12DeleteObjectEPS0_+0x73) [0x7feba2882e13]
[PID=6821, TID=1][10/23]> /home/ali/Geant4/Geant4_Install/lib/libG4processes.so(_ZN15G4ITTrackHolder5ClearEv+0x89) [0x7feba2868039]
[PID=6821, TID=1][11/23]> /home/ali/Geant4/Geant4_Install/lib/libG4processes.so(_ZN11G4Scheduler9ClearListEv+0x14) [0x7feba287c204]
[PID=6821, TID=1][12/23]> /home/ali/Geant4/Geant4_Install/lib/libG4processes.so(_ZN11G4Scheduler7ProcessEv+0x2cc) [0x7feba287f90c]
[PID=6821, TID=1][13/23]> /home/ali/Geant4/Geant4_Install/lib/libG4processes.so(_ZN21G4DNAChemistryManager3RunEv+0x28e) [0x7feba27fb51e]
[PID=6821, TID=1][14/23]> /home/ali/Geant4/Geant4_Install/lib/libG4event.so(_ZN14G4StackManager12PopNextTrackEPP13G4VTrajectory+0x56) [0x7feba3d9aff6]
[PID=6821, TID=1][15/23]> /home/ali/Geant4/Geant4_Install/lib/libG4event.so(_ZN14G4EventManager12DoProcessingEP7G4Event+0x52c) [0x7feba3d52f9c]
[PID=6821, TID=1][16/23]> /home/ali/Geant4/Geant4_Install/lib/libG4run.so(_ZN18G4WorkerRunManager11DoEventLoopEiPKci+0x1cc) [0x7feba3e1ac7c]
[PID=6821, TID=1][17/23]> /home/ali/Geant4/Geant4_Install/lib/libG4run.so(_ZN12G4RunManager6BeamOnEiPKci+0x62) [0x7feba3e0c472]
[PID=6821, TID=1][18/23]> /home/ali/Geant4/Geant4_Install/lib/libG4run.so(_ZN18G4WorkerRunManager6DoWorkEv+0x36b) [0x7feba3e1e28b]
[PID=6821, TID=1][19/23]> /home/ali/Geant4/Geant4_Install/lib/libG4run.so(_ZN20G4MTRunManagerKernel11StartThreadEP14G4WorkerThread+0x19e) [0x7feba3e41f2e]
[PID=6821, TID=1][20/23]> /usr/lib/x86_64-linux-gnu/libstdc++.so.6(+0xdc2b3) [0x7feba10dc2b3]
[PID=6821, TID=1][21/23]> /usr/lib/x86_64-linux-gnu/libc.so.6(+0x94b43) [0x7feba0c94b43]
[PID=6821, TID=1][22/23]> /usr/lib/x86_64-linux-gnu/libc.so.6(+0x126a00) [0x7feba0d26a00]
G4WT0 > G4Track (0x7feb95b0c7b0) - track ID = -10, parent ID = -2
G4WT0 >  Particle type : H3O - creator process : H2O_DNAMolecularDecay, creator model : Undefined

: Aborted (Signal sent by tkill() 6821 1000)
G4WT0 >  Kinetic energy : 38.5195 meV - Momentum direction : (0.679851,0.125113,-0.722599)
G4WT0 >  Step length : 1.02418e-297 fm  - total energy deposit : 2.70738 eV 
G4WT0 >  Pre-step point : (0,0,-6.59298e-09) - Physical volume : World (G4_WATER)
Aborted (core dumped)


How can it be fixed?


_Geant4 Version: 4.11.0
_Operating System: Linux
2

Dear,

I am investigating the issue. However, please, don’t go to so low energies cause less 150 eV of electrons, we may have only the elastic scattering process. Ionisation and excitation may not be created to form the dissociation processes and chemical species.

3

Dear @ngoc
thank you for your investigation.
I checked the physics used in the PhysicsList.cc class (G4EmDNAPhysics_option2,3,4,5,6,7,8 and EmDNAChemistry). There are no limitations in energy for e-particles; only in line 58 of this class, in the following command, we have SetEnergyRange(100 * eV, 1 * GeV); yes, I changed 100 eV to 35 eV, but the problem still persists.

4

I means you should not shoot electrons so low energies cause they cannot make a lot chemical species due to less ionisation and excitation. Can you try with G4.11.2beta ?

5

thank you @ngoc I installed G4.11.beta and works well, I shoot these electrons that are energy averages of electrons that are produced from 35 Kev x-ray sources in water with the Klein–Nishina formula