Uncertainty analysis and order-by-order optimization of chiral nuclear interactions
| Autor: | Carlsson, B. D., Ekstr��m, A., Forss��n, C., Str��mberg, D. Fahlin, Jansen, G. R., Lilja, O., Lindby, M., Mattsson, B. A., Wendt, K. A. |
|---|---|
| Jazyk: | angličtina |
| Rok vydání: | 2015 |
| Předmět: |
Nuclear Theory
Computer science Automatic differentiation QC1-999 Monte Carlo method FOS: Physical sciences General Physics and Astronomy 01 natural sciences Nuclear Theory (nucl-th) Subatomic Physics 0103 physical sciences Convergence (routing) Effective field theory Statistical physics Nuclear Experiment (nucl-ex) Uncertainty quantification Probability Theory and Statistics 010306 general physics Nuclear Experiment Uncertainty analysis 010308 nuclear & particles physics Physics Observable Regression analysis Computational Physics (physics.comp-ph) Physics - Computational Physics |
| Zdroj: | Physical Review X, Vol 6, Iss 1, p 011019 (2016) Physical Review X Physical Review X (2160-3308) vol.6(2016) |
| Popis: | Chiral effective field theory (chiEFT) provides a systematic approach to describe low-energy nuclear forces. Moreover, chiEFT is able to provide well-founded estimates of statistical and systematic uncertainties -- although this unique advantage has not yet been fully exploited. We fill this gap by performing an optimization and statistical analysis of all the low-energy constants (LECs) up to next-to-next-to-leading order. Our optimization protocol corresponds to a simultaneous fit to scattering and bound-state observables in the pion-nucleon, nucleon-nucleon, and few-nucleon sectors, thereby utilizing the full model capabilities of chiEFT. We study the effect on other observables by demonstrating error-propagation methods that can easily be adopted by future works. We employ mathematical optimization and implement automatic differentiation to attain efficient and machine-precise first- and second-order derivatives of the objective function with respect to the LECs. We use power-counting arguments to estimate the systematic uncertainty that is inherent to chiEFT and we construct chiral interactions at different orders with quantified uncertainties. Statistical error propagation is compared with Monte Carlo sampling showing that statistical errors are in general small compared to systematic ones. In conclusion, we find that a simultaneous fit to different sets of data is critical to (i) identify the optimal set of LECs, (ii) capture all relevant correlations, (iii) reduce the statistical uncertainty, and (iv) attain order-by-order convergence in chiEFT. Furthermore, certain systematic uncertainties in the few-nucleon sector are shown to get substantially magnified in the many-body sector; in particlar when varying the cutoff in the chiral potentials. The methodology and results presented in this Paper open a new frontier for uncertainty quantification in ab initio nuclear theory. Accepted for publication in Physical Review X. This version includes Supplemental Material |
| Databáze: | OpenAIRE |
| Externí odkaz: |
|