Abstract: The secondary neutron double differential cross sections and angular distributions from heavy ion induced reactions are crucial for heavy ion beam driven inertial confinement fusion and high energy density physics. Given the wide variety of heavy ion species and the insufficiency of relevant experimental data, reliable physical models are required for accurate calculations. The neutron double differential cross sections and angular distributions of 400 MeV/u ⁸⁴Kr and ¹³²Xe projectiles interacting with Al, Cu, and Pb targets at different emission angles from 5° to 80° were calculated using GEANT4 code coupled with the BIC, G4QMD and LIQMD models, which include the SLy4, SIII and SkM* nuclear interaction potential parameters. Comparison of the calculated results with the experimental data reveals that, for the double differential cross section, the BIC model exhibited significant discrepancies with experimental data in the high energy part of the spectrum at large angles, while the results of the G4QMD model for Al, Cu, and Pb targets are in relatively good agreement with the experimental data. The calculation results of the LIQMD model with different nuclear interaction potential parameters are essentially consistent. For the angular distribution, the BIC model underestimates the experimental data at 80° for Al and Cu targets but overestimates the data for Pb targets. Both the G4QMD and LIQMD models provide a good reproduction of the experimental data.