DD reaction channels have different coefficients

[shimwell]

Just looking at the reac_DD function and noticed that is handles both reaction channels dd->p+t and dd->n+3He with the same equation. I notice bosch and hale paper has slightly different coefficients for each reaction

I think they are very very close to each, the other channel does not emit a neutron so perhaps this is not important

However at 20,000eV I get these two different reactivity values when using the different coefficients

D+D=p+T reactivity 2.398954986744111e-24 cm^3/s
D+D=n+He3 reactivity 2.602651467423054e-24 cm^3/s <- nest is using this one

[aidancrilly]

Correct that I only look at the neutron producing channel currently. Would be straightforward to include the p+T channel.

I have charged particle extensions to NeSST in the works (will have to rethink the name!) but if you need this sooner I can add this early?

[shimwell]

No hurry from my side.