Problem with photonuclear physics phys:p in MCNP

  • Thread starter Thread starter Nour_wahban
  • Start date Start date
  • Tags Tags
    Mcnp
Click For Summary
SUMMARY

The discussion centers on troubleshooting photonuclear physics settings in MCNP while attempting to activate cobalt-59 using neutrons from californium-252. Users reported warnings regarding the necessity of photonuclear physics (phys:p) and the non-utilization of material 3. The recommended solution involves adjusting the PHYS:p card to PHYS:p 100 0 0 -1 0 J 0 for unbiased photonuclear treatment, which may slow down computations but is deemed safe. Additionally, the user confirmed that material 3 is correctly defined and utilized in the input file.

PREREQUISITES
  • Understanding of MCNP (version 6.3) and its input file structure
  • Knowledge of photonuclear reactions and their thresholds, particularly regarding cobalt-59
  • Familiarity with neutron activation processes and gamma emissions from californium-252
  • Experience with interpreting warnings and errors in MCNP simulations
NEXT STEPS
  • Research the implications of using the PHYS:p card in MCNP simulations
  • Study the photoneutron reaction mechanisms and their thresholds in nuclear physics
  • Explore the MCNP manual for version-specific changes related to photonuclear physics
  • Investigate the activation spectrum of cobalt-60 and its decay processes
USEFUL FOR

Nuclear physicists, radiation safety professionals, and researchers working with MCNP simulations for neutron activation studies.

Nour_wahban
Messages
5
Reaction score
0
TL;DR
I'm trying to activate cobalt 59 to cobalt 60 and obtain its peaks, but the phys:p is not working.
I'm working on a very basic setup to activate cobalt 59 with neutrons emitted from californium 252. MCNP is giving me the following warnings:
warning. photonuclear physics may be needed (phys:p).
warning. material 3 is not used in the problem.
could anyone help me understand where the problem is. I've attached the input file.
 

Attachments

Last edited:
Engineering news on Phys.org
Nour_wahban said:
TL;DR Summary: I'm trying to activate cobalt 59 to cobalt 60 and obtain its peaks, but the phys:p is not working.

warning. photonuclear physics may be needed (phys:p).
I suspect that one may have to consider photoneutron reactions from the gammas emitted by 252Cf. Most photoneutron thresholds are above 7 MeV.
See - https://journals.aps.org/prc/abstract/10.1103/PhysRevC.10.795

This is a useful reference: https://www.osti.gov/servlets/purl/1323135
MatMCNP: A Code for Producing Material Cards for MCNP

About photoneutrons
https://www.sciencedirect.com/topics/physics-and-astronomy/photoneutrons
A gamma photon with energy sufficiently large to overcome the neutron binding energy (about 7 MeV in most nuclides) may cause a (γ, n) reaction.
 
  • Like
Likes   Reactions: Alex A
You will see a lot of warnings and you should chase them. Warnings often show up a bigger problem. That said, this one is usually a comment, not a warning. It can usually be ignored, unless you know it's a common event in the problem, but again I don't think I've seen it as an actual warning, so I don't know. Might be related to the use of a distribution function for your neutron source energy.

phys:p is a more complicated card that doesn't just turn on and off photonuclear physics, but a lot of other stuff. It also tends to change between versions so you have to check the right manual to make sure the switch is the one you want.
The default for 6.3, according to the manual is PHYS:p 100 0 0 0 0 J 0
I am guessing the safest way to turn on photonuclear physics would be,
PHYS:p 100 0 0 -1 0 J 0
That according to the manual is an unbiased 'sample every collision' photonuclear treatment.
This might slow things down and probably won't affect the results much.

I do not see an issue with material 3. In the def there are only 3 spaces for the extra lines, but & is used to continue the previous line, so that is absolutely fine. Material 3 is used in cell 3. You should not be getting that error with that input file as best as I can tell.

Hope that helps!

Edit, Just to check, you are after the fast activation spectrum (n,gamma) right? You are not going to see Co-60 decay, because the program doesn't handle time like that.
 
Last edited:
  • Like
Likes   Reactions: Astronuc

Similar threads

  • · Replies 8 ·
Replies
8
Views
2K
  • · Replies 4 ·
Replies
4
Views
826
  • · Replies 4 ·
Replies
4
Views
2K
  • · Replies 15 ·
Replies
15
Views
3K
  • · Replies 5 ·
Replies
5
Views
5K
  • · Replies 1 ·
Replies
1
Views
4K
  • · Replies 2 ·
Replies
2
Views
4K
  • · Replies 20 ·
Replies
20
Views
19K
Replies
1
Views
2K
  • · Replies 4 ·
Replies
4
Views
1K