We develop an adaptive Hessian-based non-stationary Gaussian process (GP) response surface method for approximating a probability density function (pdf) that exploits its structure, particularly the Hessian of its negative logarithm. Of particular interest to us are pdfs that arise from the Bayesian solution of large-scale inverse problems, which imply very expensive-to-evaluate pdfs. The metho...

The exceptional region of the parameter space (ERPS) two Higgs doublet model (2HDM) is defined to be regime where scalar potential takes on a very special form. In standard parametrization 2HDM with squared mass parameters ${m}_{11}^{2}$, ${m}_{22}^{2}$, ${m}_{12}^{2}$, and dimensionless couplings, ${\ensuremath{\lambda}}_{1},{\ensuremath{\lambda}}_{2},\dots{},{\ensuremath{\lambda}}_{7}$, ERPS ...