uq
Differences
This shows you the differences between two versions of the page.
Both sides previous revisionPrevious revisionNext revision | Previous revision | ||
uq [2018/11/26 13:13] – [Uncertainty Quantification] slongbottom | uq [2019/10/08 19:10] (current) – removed slongbottom | ||
---|---|---|---|
Line 1: | Line 1: | ||
- | ====== Uncertainty Quantification ====== | ||
- | |||
- | ---- | ||
- | |||
- | //potfit// uses the potential ensemble method ((Frederiksen, | ||
- | |||
- | To enable this feature, compile //potfit// with the '' | ||
- | |||
- | < | ||
- | ===== The Ensemble Method ===== | ||
- | |||
- | An ensemble of potentials are generated by taking a series on Markov chain Monte Carlo steps starting from the best fit potential parameters. The step size in each parameter direction is scaled dependant on the curvature in each parameter. This is encoded using information about the eigenvalues of the hessian at the best fit potential minimum, for potential parameters $\Theta=\{\theta_1, | ||
- | |||
- | \begin{equation} | ||
- | \Delta\theta_{i}=\sum_{j=1}^{\rm N_{\rm P}}\sqrt{\frac{\rm R}{{\rm{max}} (1, | ||
- | \end{equation} | ||
- | where $\lambda_j$ are the hessian eigenvalues, | ||
- | |||
- | The MCMC algorithm samples potentials from the distribution at a temperature, | ||
- | |||
- | ===== Parameters ===== | ||
- | |||
- | |<100% 33% 33% 33%>| | ||
- | ^ parameter name | parameter type | default value | | ||
- | | short explanation ||| | ||
- | |||
- | ==== Required parameters | ||
- | |||
- | |<100% 33% 33% 33%>| | ||
- | ^ **acceptance_rescaling*** | float | (none) | | ||
- | | R value to tune MCMC acceptance rate ||| | ||
- | |||
- | |<100% 33% 33% 33%>| | ||
- | ^ **acc_moves*** | integer | (none) | | ||
- | | Number of accepted MCMC moves required ||| | ||
- | |||
- | ==== Optional parameters | ||
- | |||
- | |<100% 33% 33% 33%>| | ||
- | ^ **sloppyfile** | string | '' | ||
- | | Potential ensemble output filename, uses '' | ||
- | |||
- | |<100% 33% 33% 33%>| | ||
- | ^ **uq_temp** | float | 1.0 | | ||
- | | Temperature scaling parameter $\alpha$ ||| | ||
- | |||
- | |<100% 33% 33% 33%>| | ||
- | ^ **use_svd** | boolean | 0 | | ||
- | | Use singular value decomposition to find Hessian eigenvalues (default is eigenvalue decomposition) ||| | ||
- | |||
- | |<100% 33% 33% 33%>| | ||
- | ^ **hess_pert** | float | (none) | | ||
- | | Percentage parameter perturbation in Hessian finite difference calculation ||| | ||
- | |||
- | |<100% 33% 33% 33%>| | ||
- | ^ **eig_max** | float | 1.0 | | ||
- | | Alternative MCMC step perturbation maximum value in max('' | ||
- | |||
- | |||