Let \(Z(t)= \exp \left( \sqrt{ 2} B_H(t)- \left|t \right|^{2H}\right) , t\in \mathbb {R}\) with \(B_H(t),t\in \mathbb {R}\) a standard fractional Brownian motion (fBm) with Hurst parameter \(H \in (0,1]\) and define for x non-negative the Berman function \(\begin{aligned} \mathcal {B}_{Z}(x)= \mathbb {E} \left\{ \frac{ \mathbb {I} \{ \epsilon _0(RZ) > x\}}{ \epsilon _0(RZ)}\right\} \in (0,\infty ), \end{aligned}\) where the random variable R independent of Z has survival function \(1/x,x\geqslant 1\) and \(\begin{aligned} \epsilon _0(RZ) = \int _{\mathbb {R}} \mathbb {I}{\left\{ RZ(t)> 1\right\} }{dt} . \end{aligned}\) In this paper we consider a general random field (rf) Z that is a spectral rf of some stationary max-stable rf X and derive the properties of the corresponding Berman functions. In particular, we show that Berman functions can be approximated by the corresponding discrete ones and derive interesting representations of those functions which are of interest for Monte Carlo simulations presented in this article.
