### Abstract

In this paper, we present a neutron spectrum unfolding technique based on a modification of the least-squares method. The main innovation is the use of a Krylov subspace iteration method to solve the least-squares normal equations. This method was employed because it performs better on ill-conditioned systems of linear equations as compared with standard direct-solution methods. Three different least-squares solution techniques are compared and evaluated in terms of (i) accuracy in the prediction of the energy spectrum, (ii) computational efficiency, and (iii) robustness to noise. The unfolding is performed on measured pulse-height distributions as well as pulse height distributions generated with the Monte Carlo code MCNP-PoliMi. Using this code, neutron energy depositions on the constituents of the scintillator are individually tracked, and the light output generated at each interaction is suitably accounted for. This procedure allows for a very accurate simulation of the liquid scintillator detector response. The precise knowledge of the neutron energy spectrum provides information not only about the presence or absence of fissile material, but also about the characteristics of the material. We show that the proposed technique performs well in the unfolding of neutron pulse-height distributions from Monte Carlo simulations, and fairly well for a measured distribution from a Cf-252 neutron source.

Original language | English (US) |
---|---|

Title of host publication | Joint International Topical Meeting on Mathematics and Computations and Supercomputing in Nuclear Applications, M and C + SNA 2007 |

State | Published - Dec 3 2007 |

Event | Joint International Topical Meeting on Mathematics and Computations and Supercomputing in Nuclear Applications, M and C + SNA 2007 - Monterey, CA, United States Duration: Apr 15 2007 → Apr 19 2007 |

### Publication series

Name | Joint International Topical Meeting on Mathematics and Computations and Supercomputing in Nuclear Applications, M and C + SNA 2007 |
---|

### Other

Other | Joint International Topical Meeting on Mathematics and Computations and Supercomputing in Nuclear Applications, M and C + SNA 2007 |
---|---|

Country | United States |

City | Monterey, CA |

Period | 4/15/07 → 4/19/07 |

### Fingerprint

### All Science Journal Classification (ASJC) codes

- Mathematics(all)
- Nuclear and High Energy Physics

### Cite this

}

*Joint International Topical Meeting on Mathematics and Computations and Supercomputing in Nuclear Applications, M and C + SNA 2007.*Joint International Topical Meeting on Mathematics and Computations and Supercomputing in Nuclear Applications, M and C + SNA 2007, Joint International Topical Meeting on Mathematics and Computations and Supercomputing in Nuclear Applications, M and C + SNA 2007, Monterey, CA, United States, 4/15/07.

**A sequential least-squares algorithm for neutron spectrum unfolding from pulse-height distributions measured with liquid scintillators.** / Yunlin, Xu; Flaska, Marek; Pozzi, Sara; Protopopescu, Vladimir; Downar, Thomas.

Research output: Chapter in Book/Report/Conference proceeding › Conference contribution

TY - GEN

T1 - A sequential least-squares algorithm for neutron spectrum unfolding from pulse-height distributions measured with liquid scintillators

AU - Yunlin, Xu

AU - Flaska, Marek

AU - Pozzi, Sara

AU - Protopopescu, Vladimir

AU - Downar, Thomas

PY - 2007/12/3

Y1 - 2007/12/3

N2 - In this paper, we present a neutron spectrum unfolding technique based on a modification of the least-squares method. The main innovation is the use of a Krylov subspace iteration method to solve the least-squares normal equations. This method was employed because it performs better on ill-conditioned systems of linear equations as compared with standard direct-solution methods. Three different least-squares solution techniques are compared and evaluated in terms of (i) accuracy in the prediction of the energy spectrum, (ii) computational efficiency, and (iii) robustness to noise. The unfolding is performed on measured pulse-height distributions as well as pulse height distributions generated with the Monte Carlo code MCNP-PoliMi. Using this code, neutron energy depositions on the constituents of the scintillator are individually tracked, and the light output generated at each interaction is suitably accounted for. This procedure allows for a very accurate simulation of the liquid scintillator detector response. The precise knowledge of the neutron energy spectrum provides information not only about the presence or absence of fissile material, but also about the characteristics of the material. We show that the proposed technique performs well in the unfolding of neutron pulse-height distributions from Monte Carlo simulations, and fairly well for a measured distribution from a Cf-252 neutron source.

AB - In this paper, we present a neutron spectrum unfolding technique based on a modification of the least-squares method. The main innovation is the use of a Krylov subspace iteration method to solve the least-squares normal equations. This method was employed because it performs better on ill-conditioned systems of linear equations as compared with standard direct-solution methods. Three different least-squares solution techniques are compared and evaluated in terms of (i) accuracy in the prediction of the energy spectrum, (ii) computational efficiency, and (iii) robustness to noise. The unfolding is performed on measured pulse-height distributions as well as pulse height distributions generated with the Monte Carlo code MCNP-PoliMi. Using this code, neutron energy depositions on the constituents of the scintillator are individually tracked, and the light output generated at each interaction is suitably accounted for. This procedure allows for a very accurate simulation of the liquid scintillator detector response. The precise knowledge of the neutron energy spectrum provides information not only about the presence or absence of fissile material, but also about the characteristics of the material. We show that the proposed technique performs well in the unfolding of neutron pulse-height distributions from Monte Carlo simulations, and fairly well for a measured distribution from a Cf-252 neutron source.

UR - http://www.scopus.com/inward/record.url?scp=84858462206&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84858462206&partnerID=8YFLogxK

M3 - Conference contribution

AN - SCOPUS:84858462206

SN - 0894480596

SN - 9780894480591

T3 - Joint International Topical Meeting on Mathematics and Computations and Supercomputing in Nuclear Applications, M and C + SNA 2007

BT - Joint International Topical Meeting on Mathematics and Computations and Supercomputing in Nuclear Applications, M and C + SNA 2007

ER -