26.02.2013 Views

Program - Brookhaven National Laboratory

Program - Brookhaven National Laboratory

Program - Brookhaven National Laboratory

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

[1] B. Paxton, L. Bildsten, A. Dotter et al., Astrophys. J. Suppl. Ser., 192:3 (2011). http://mesa.<br />

sourceforge.net/. [2] F. Herwig, M. E. Bennet, S. Diehl et al., “Nucleosynthesis simulations for a wide<br />

range of nuclear production sites from NuGrid” in Proc. of the 10th Symp. on Nuclei in the Cosmos (NIC<br />

X), 23 (2008). http://www.astro.keele.ac.uk/nugrid. [3] P. Möller, J. R. Nix, W. D. Myers, and W.<br />

J. Swiatecki, At. Data Nucl. Data Tables 59 (1995) 185. [4] P. Möller, J. R. Nix, and K.-L. Kratz, At.<br />

Data Nucl. Data Tables 66 (1997) 131.<br />

Session LC Covariances<br />

Wednesday March 6, 2013<br />

Room: Empire East at 3:30 PM<br />

LC 1 3:30 PM<br />

On the Use of Zero Variance Penalty Conditions for the Generation of Covariance Matrix<br />

Between Model Parameters<br />

G. Noguere, P. Archier, C. De Saint Jean, B. Habert<br />

CEA, DEN, DER, SPRC, Cadarache, F-13108 Saint-Paul-lez-Durance, France<br />

Analytic and Monte-Carlo algorithms have been developed and implemented in the nuclear data code<br />

CONRAD to generate covariances between model parameters. However, in the resolved resonance range<br />

of the neutron cross sections, thousands model parameters and experimental data points are needed to<br />

describe the main neutron reactions of interest for the nuclear applications. In practice, all of the model<br />

parameter covariances cannot be simultaneously estimated from well-defined sets of experimental data.<br />

The situation becomes more and more complex when the number of open reaction channels increases. This<br />

optimization problem is nonconvex in general, making it difficult to find the global optima. A solution to<br />

simplify the problem consists of dividing the model parameter sequence into blocks of variables, the uncertainties<br />

of which can be propagated through fixed-order sequential data assimilation procedures. According<br />

to definition and nomenclature used in statistics, we can identify the ”observable”, ”latent” and ”nuisance”<br />

variables. Observable variables can be directly determined from the experimental data. Latent variables<br />

(as opposed to observable variables) may define redundant parameters or hidden variables that cannot<br />

be observed directly. This term reflects the fact that such variables are really there, but they cannot be<br />

observed or measured for practical reasons. Nuisance variables (i.e. experimental corrections) correspond<br />

to aspect of physical realities whose properties are not of particular interest as such but are fundamental<br />

for assessing reliable model parameters. The methodology consists in calculating cross-covariance matrix<br />

between Gaussian distributed parameters by introducing ”variance penalty” terms. This approach was<br />

proposed by D.W. Muir for determining the contribution of individual correlated parameters to the uncertainty<br />

of integral quantities. It provides a useful understanding of how the model parameters are related<br />

to each other and ensures the positive-definiteness of the covariance matrix. Performances of the analytic<br />

and Monte-Carlo models will be illustrated with covariances calculated between new Am-241 resonance<br />

parameters recently established in the frame of the IRMM/CEA collaboration.<br />

LC 2 4:00 PM<br />

Evaluation of Uncertainties in βeff by Means of Deterministic and Monte Carlo Methods<br />

I. Kodeli<br />

Joˇzef Stefan Institute (JSI), Jamova 39, 1000 Ljubljana, Slovenia<br />

173

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!