Photonuclear Neutron Production in OpenMC: Verification Against MCNPX, FLUKA, and a First-Collision Analytical Solution

Introduction to the show: ident: Physics. The day in general physics: instruments, methods and results across the field.

Kai: We're looking at the paper: "Photonuclear Neutron Production in OpenMC".

Mira: Here's the gist of the paper. The modeling of photonuclear reactions is increasingly important for applications involving highenergy photon fields, including accelerator-driven neutron sources, radiation shielding, medical physics, and fusion technologies <ref:2607.26045#pg2>.

Kai: This work presents a systematic verification of the OpenMC photonuclear implementation using six single-collision broomstick benchmarks based on 2H, 9Be, and 238U targets irradiated by monoenergetic 5 MeV and 15 MeV photons and by a continuous 1–20 MeV linear accelerator (LINAC)-representative spectrum <ref:2607.26045#pg4>.

Mira: Verification Scope The study systematically verifies the OpenMC photonuclear implementation against MCNPX, FLUKA, and a first-collision analytical solution <ref:2607.26045#pg2>.

Kai: This verification is conducted under idealized, single-collision conditions to isolate the treatment of individual particle interactions from secondary transport and macroscopic geometrical effects <ref:2607.26045#pg9>.

Title and authors: Kai: Now, looking at what they actually wrote in the abstract, "Photonuclear Neutron Production in OpenMC: Verification Against MCNPX, FLUKA, and a First-Collision Analytical Solution" <ref:2607.26045#pg3>. They say that modeling these photonuclear reactions is increasingly important for applications involving highenergy photon fields like accelerator-driven neutron sources, radiation shielding, medical physics and fusion technologies <ref:2607.26045#pg3>.

Mira: That sets the stage; it’s not just academic curiosity, this is for things like accelerator-driven neutron sources or radiation shielding <ref:2607.26045#pg3>.

Lev: So what they are doing is taking OpenMC's photonuclear physics, which they say is currently only in an unofficial development branch and needs verification before broader scientific use or possible integration into the official code distribution <ref:2607.26045#pg3>.

Kai: They present a systematic verification using six single-collision broomstick benchmarks based on 2H, 9Be and 238U targets irradiated by monoenergetic five MeV and fifteen MeV photons and by a continuous one to twenty MeV linear accelerator LINAC representative spectrum <ref:2607.26045#pg4> <ref:2607.26045#pg3,using six single-collision broomstick benchmarks based on 2H, 9Be and 238U>.

Mira: That's the test setup; they are hitting these specific targets with specific photon energies to see how the code reacts <ref:2607.26045#pg4>.

Lev: They compared OpenMC with MCNPX using common ENDF7u photonuclear data, FLUKA using its own native models, and a first-collision analytical solution for the integrated neutron yield <ref:2607.26045#pg4>.

Kai: And they also performed calculations using the IAEA/PD-two thousand nineteen library to quantify nuclear-data sensitivity <ref:2607.26045#pg4> <ref:2607.26045#pg1>.

Mira: Basically, they are proving that this implementation reproduces the integrated neutron-production response of an established reference code like MCNPX when using common data <ref:2607.26045#pg3>.

Lev: It’s a direct comparison to see if OpenMC is doing what it should be doing when the physics matches existing standards <ref:2607.26045#pg3>.

Kai: And they also used FLUKA's own models to see how it compares, noting that FLUKA differs from the analytical solution by approximately six to nine percent for monoenergetic cases and by no more than approximately four percent for continuous-source cases <ref:2607.26045#pg4>.

Mira: Those are some specific numbers showing where the discrepancies lie between different modeling approaches <ref:2607.26045#pg4>.

Lev: It’s a good baseline to start with, but we need to see how OpenMC stacks up when the data itself changes, which is what they test next <ref:2607.26045#pg4>.

The paper's summary: Kai: So the paper points out that one of the key things they did was quantifying nuclear-data sensitivity <ref:2607.26045#pg4>. They performed calculations using the IAEA/PD-two thousand nineteen library to see how much difference there is when you replace ENDF7u with IAEA/PD-two thousand nineteen <ref:2607.26045#pg4> <ref:2607.26045#pg1>.

Mira: They found that replacing ENDF7u with IAEA/PD-two thousand nineteen produces differences that are generally larger than the OpenMC MCNPX code-to-code differences obtained under common data <ref:2607.26045#pg4> <ref:2607.26045#pg1>.

Lev: That’s a big deal for users; it means if you're using the IAEA/PD-two thousand nineteen library, you have to expect bigger shifts in your predictions than just switching from ENDF7u to whatever ENDF7u was originally calibrated with <ref:2607.26045#pg4> <ref:2607.26045#pg1>.

Kai: For the monoenergetic cases, the results using IAEA/PD-two thousand nineteen differ from the analytical ENDF7u solution by approximately plus zero point zero three percent for 9Be, plus eleven point three percent for 2H and minus four point six percent for 238U <ref:2607.26045#pg6> <ref:2607.26045#pg3>.

Mira: Those specific deviations show that the sensitivity depends strongly on both the nuclide and the source spectrum <ref:2607.26045#pg6>.

Lev: That tells us you can't just pick one data set and assume it’s fine; you have to check what happens across different isotopes and energy ranges <ref:2607.26045#pg6>.

Kai: And they also demonstrated that agreement in the integrated yield doesn't necessarily imply agreement in the emitted neutron spectrum <ref:2607.26045#pg2>.

Mira: That’s a crucial point for anyone doing shielding or dosimetry; you can get the same total number of neutrons even if the way those neutrons are distributed across energy is totally different <ref:2607.26045#pg9>.

Lev: So, while they agree on the total yield, the shape of that spectrum is where things get messy and needs more work <ref:2607.26045#pg9>.

The paper's improvements: Kai: To wrap up, the paper "Photonuclear Neutron Production in OpenMC: Verification Against MCNPX, FLUKA, and a First-Collision Analytical Solution" <ref:2607.26045#pg1> it shows that the OpenMC photonuclear implementation reproduces the integrated neutron-production response of an established reference code when using common ENDF7u photonuclear data <ref:2607.26045#pg2>.

Mira: So in simple terms, they’ve confirmed that the new part of OpenMC works as expected against MCNPX, which is a key piece of evidence for its use <ref:2607.26045#pg2>.

Lev: It supports the use of this unofficial OpenMC branch for integrated photonuclear neutron production under the verified configurations <ref:2607.26045#pg3>.

Kai: But they also flag secondary-particle kinematics as an area requiring additional assessment <ref:2607.26045#pg9>.

Mira: So they aren't saying it’s perfect; they are highlighting that we still need to look at how those secondary particles behave <ref:2607.26045#pg9>.

Lev: Further work should extend the verification to additional nuclides, reaction thresholds, angular distributions and geometries in which secondary transport is significant <ref:2607.26045#pg9>.

Kai: So the overall idea here is that while this verifies the integrated yield against MCNPX under common data, we need more tests on things like angular distributions and where secondary transport matters <ref:2607.26045#pg9>.

Conclusion: Mira: We've talked about how this paper, "Photonuclear Neutron Production in OpenMC: Verification Against MCNPX, FLUKA, and a First-Collision Analytical Solution" <ref:2607.26045#pg1> confirms the integrated yield agreement with MCNPX under common ENDF7u photonuclear data <ref:2607.26045#pg2>.

Kai: And they also showed that using IAEA/PD-two thousand nineteen data can change things by up to eleven point three percent in some cases like 2H or minus four point six percent for 238U <ref:2607.26045#pg4> <ref:2607.26045#pg3>.

Lev: That means the sensitivity depends strongly on both the nuclide and the source spectrum <ref:2607.26045#pg6>.

Mira: It’s a strong indication that we have to be careful about which data set we use because it impacts things differently across different isotopes and energies <ref:2607.26045#pg6>.

Kai: And they also highlighted that agreement in total yield doesn't necessarily imply agreement in the emitted neutron spectrum <ref:2607.26045#pg2>.

Lev: That means for applications like shielding, you can get the same total number of neutrons even if the energy distribution is completely different <ref:2607.26045#pg9>.

Mira: So while this verification provides confidence in the integrated yield against MCNPX, we still need to look at those spectral details and secondary transport effects <ref:2607.26045#pg9>.

Kai: The paper "Photonuclear Neutron Production in OpenMC: Verification Against MCNPX, FLUKA, and a First-Collision Analytical Solution" <ref:2607.26045#pg1> is a solid piece of evidence that the new implementation reproduces the integrated neutron-production response of an established reference code <ref:2607.26045#pg2>.

Lev: So we need to keep checking those angular distributions and geometries where secondary transport is significant in future work <ref:2607.26045#pg9>.