Menu

Decaying the abundances at the end of the calculation

Anonymous
2018-11-15
2018-11-16
  • Anonymous

    Anonymous - 2018-11-15

    Move from here:

    Dear professor,
    I have one question. I follow this post and get the same result of abundance like you. But it seems that the final abundance includes many unstable nuclei. Do we need a step to let these nuclei to decay into stable nuclei? Can this been done in the NucNet or need other subroutine?
    Thanks a lot, I learn a lot from your posts.

     
    • Bradley S. Meyer

      Thanks for your post. To follow decays at the end of a calculation, you can change the end time (tend) to a longer time to allow the decay to occur. Alternatively, you can use the nuclear decay project:

      https://sourceforge.net/p/nucnet-projects/wiki/nuclear_decay/

      Of course, in either case, your input data file will have to have all the relevant decays. For example, you will want to be sure you have all the alpha decays in the actinide region. You might have to add those decays, as in this post.

       
  • Anonymous

    Anonymous - 2018-11-15

    Thank you very much for the reply.
    Actually, I have tried a very long time (1.e18 s), but it seems doesn't change much in the final abundance(I compare the data, most of them are the same, only few changes a little). Still, there are many unstable nuclei.
    Does this mean that the JINA-Reaclib data doesn't include the beta decays and the alpha decays you mention above?

     
  • Anonymous

    Anonymous - 2018-11-16

    Dear professor,

    Maybe I have figure out the problem, I address it below and you can see if I am right or wrong.

    As in the r process calculation, I set the xpath as "[z<=90]", the network doesn't include the nuclei with z>90, so the unstable nuclei with z=90 can't decay into z=91 nuclei by beta decay. So when I check the final abundance, the nuclei with z=90 doesn't changed much with different time.

    Thank you very much.

     
    • Bradley S. Meyer

      Yes, that is correct. Decaying nuclei will pile up at Z = 90 if that is the limit. You might also try an xpath expression with something like "[a <= 293]".

      As for decays in the JINA data base, there may be some alpha decays in the actinide region not present. You can try updating from the JINA web site or by using the jina_to_webnucleo project.

       
  • Anonymous

    Anonymous - 2018-11-16

    Thank you very much for the reply.

     

Anonymous
Anonymous

Add attachments
Cancel





Want the latest updates on software, tech news, and AI?
Get latest updates about software, tech news, and AI from SourceForge directly in your inbox once a month.