Please fill out the following information to help in answering your question, and also see tips for posting code snippets. If you don’t provide this information it will take more time to help with your problem!
_Geant4 Version:_11.3.0
_Operating System:_ubuntu
Compiler/Version:
_CMake Version:_3.28.3
Hello, I am currently working on a GEANT4 simulation for the inverse kinematic reaction p(⁷Li, n)⁷Be and have encountered some challenges that I would greatly appreciate your guidance on.
I have tested different combinations of reference physics lists, including G4HadronicPhysicsQGSP_BIC_AllHP, G4IonPhysicsPHP, and G4IonINCLXXPhysics, among others. However, the simulation yields counts only for ⁷Li energies up to 22 MeV, and beyond this energy, I do not observe any counts—even when simulating 10⁹ ⁷Li ions.
Comparing my results with prior published data for the direct kinematics reaction ⁷Li(p, n)⁷Be, I found that the cross-section data is available for proton energies up to 10 MeV. Simulating this reaction in GEANT4 produces results consistent with the expected cross-section. Since the center-of-mass (COM) cross-sections should be the same for both direct and inverse kinematics, I would expect to observe counts beyond 22 MeV in the inverse kinematics simulation, but this is not the case. Additionally, the simulated cross-section does not match the direct kinematics cross-section trend.
I am trying to simulate the p(⁷Li, n)⁷Be reaction at ⁷Li = 28 MeV. Given the limited experimental data for inverse kinematics, I am wondering if GEANT4 is not correctly handling the data for this reaction. Is there a possibility that GEANT4 does not correctly convert direct kinematics cross-sections into the COM frame for inverse kinematics? Could this issue be related to the physics list selection?
I would greatly appreciate any insights or suggestions on how to refine my simulation to obtain accurate results.
Looking forward to your response.