J. Miller, “Analytical Inverse Model for Post Event Attribution of Plutonium”, Proceedings of the 2009 Annual Meeting of the Institute of Nuclear Materials Management, Tucson, AZ, July 12-15, 2009.
An integral part of deterring nuclear terrorism is the swiftattributionof anyeventto a particular state or organization. By quickly being able to identify the responsible party after a nuclearevent, appropriate people may be held accountable for their actions. Currently, there is a system in place to determine the origin of nuclear devices and materials frompost-eventdata; however, the system requires significant time to produce an answer within acceptable error margins. Described here is a deterministic approach derived from first principles to solve theinverseproblem. The derivation starts with the basic change rate equation and ends in relationships for important nuclear concentrations and device yield. This results in a computationally efficient and timely method for producing an estimate of the material attributes. This estimate can then be used as a starting point for other more detailed methods and reduce the overall computation time of thepost-eventforensics. This work focused on a specific type of nuclearevent: aplutoniumimprovised nuclear device (IND) explosion. Frompost-eventisotopic ratios, this method determines the device’s pre-eventisotopic concentrations of special nuclear material. From the original isotopic concentrations, the field of possible origins for the nuclear material is narrowed. In this scenario, knowing where the nuclear material did not originate may be as important as knowing where it did. The derived methodology was tested using several cases of interest including simplified and realistic cases. For the simplistic cases, only two isotopes comprised the material being fissioned. In the realistic cases, both Weapons Grade and Reactor Gradeplutoniumwere used to cover the spectrum of possible fissile material. The methodology performed very well over the desired energy range. Errors were under two percent from the expected values for all yields under 50 kT. In the realistic cases, competing reactions caused an increase in error; however, these stayed under five percent. As expected, with an increased yield, the error continued to rise, but these errors increased linearly. A sensitivity analysis was performed on the methodology to determine the impact of uncertainty in various physical constants. The result was that theinversemethodology is not overly sensitive to perturbations in these constants.