Wu et al. study Log S-fBM model in finance; small intermittency drives rough/multifractal link

The Log S-fBM model: Statistical analysis

The Log S-fBM model is a stochastic volatility model introduced by Wu et al. in [ 25 ] . It is characterized by a stochastic log volatility being a stationary fractional brownian motion (S-fBM) process: a stationary gaussian process whose autocovariance function has a power decay driven by the Hurst exponent H H and variance having a multiplicative coefficient namely the intermittency coefficient. One of the main properties of the Log S-fBM model is that it conciliates the rough volatility setting where the Hurst exponent is typically near 0.1 0.1 (see [ 24 ] ) and the multifractal volatility setting where the Hurst exponent is of order 0 0 as introduced in [ 23 , 29 ] as follows: the volatility measure of the Log S-fBM model tends to the one of multifractal volatility when the Hurst exponent goes to 0 0 . According to the numerical findings in [ 25 ] , this intermittency was observed to be of order 0.02 0.02 across a large range of financial assets which justifies the meaningfulness of considering a small intermittency approximation of log volatility moments as a model calibration method (so called general method of moment known as GMM). In this work, we perform a statistical analysis of the Log S-fBM model. First, we derive scaling properties related to the S-fBM process as well as the Log S-fBM integrated volatility measure. Then, we present deviation inequalities of the Log S-fBM process with a precise description of its tail distribution emphasizing on its sensitivity with respect to the Hurst exponent and the intermittency coefficient. Moreover, we develop a statistical hypothesis testing in order to test the null hypothesis of the Hurst exponent, in other words, test the rough VS multifractal hypothesis. Last but not least, we revisit the scale invariance properties of the log volatility increment process with explicit formulas leveraging the small intermittency approximation, enabling to reproduce the analoguous properties to the rough as well as the multifractal volatility settings.

Keywords: Log S-fBM model, S-fBM process, small intermittency approximation, deviation inequalities, scale invariance.

Volatility is a critical component in the field of asset pricing, influencing the valuation of financial instruments, the assessment of risk and the formulation of investment strategies. In other words, the accurate modeling of volatility is crucial for several reasons. It directly affects the estimation of risk premia, influences the pricing of derivatives and provides insights into the market’s response to macroeconomic events. Moreover, understanding the dynamics of volatility is essential for effective risk management and the development of trading strategies that can adapt to changing market conditions. Traditional asset pricing models, such as the Capital Asset Pricing Model also known as the CAPM introduced in [ 42 ] and the Black-Scholes option pricing model (see the seminal work of [ 43 ] ), often assume that volatility is either constant or deterministic. However, real-world data suggest that volatility is inherently stochastic and exhibits time-varying characteristics, such as clustering, mean reversion and asymmetry (for more details see [ 44 ] ). These observations have led to the development of advanced volatility models that better capture the complexity of financial markets. The introduction of stochastic volatility models, such as the Heston model in [ 47 ] , marked a significant advancement in pricing theory. These models incorporate an additional stochastic process defined through its own dynamic to describe volatility, allowing more accurate pricing of derivatives and a better understanding of market dynamics. Furthermore, they address key empirical phenomena like the volatility smile observed in options markets, which cannot be explained by constant volatility models [ 48 ] . In fact, many stochastic volatility models, such as the Heston model, allow to derive a smile dynamic controlled via a finite number of parameters together with a clear interpretation of every parameter and its impact on the smile. On the other hand, the Multifractal Random Walk (MRW) is a model that has gained prominence in the modeling of financial markets and turbulent systems. Introduced by Bacry, Muzy and Delour [ 23 ] (see also [ 23 , 28 , 22 ] ), the MRW is designed to capture the complex scaling properties and long-range dependence observed in empirical data. Traditional models like geometric Brownian motion, used in the Black-Scholes framework, assume homogeneity in the underlying stochastic processes. However, empirical studies have demonstrated that financial time series exhibit multifractality, meaning their statistical properties vary across different time scales [ 16 ] . The main advantage of the MRW is its ability to reproduce the heavy tails, volatility persistence and scaling laws commonly observed in financial returns. These characteristics are important for accurately modeling risk and pricing derivatives, particularly in markets that deviate from the assumptions of constant volatility and Gaussian returns. From a mathematical perspective, the MRW builds on the theory of multifractals and cascades, which were originally developed in the context of turbulence and chaos theory [ 49 ] . In the MRW, the volatility process is itself stochastic, driven by a Gaussian process with long memory, resulting in a non-linear, multifractal structure for the asset returns. In recent years, the concept of rough volatility, as proposed by Gatheral et al. [ 24 ] , has gained traction. Rough volatility models represent volatility as a process with rough paths, characterized by fractional Brownian motion with a Hurst parameter less than 0.5. These models have been shown to provide a superior fit to high-frequency data, offering new insights into the microstructure of financial markets and improving the accuracy of derivative pricing. At this stage, one may wonder about the relationship between multifractal models and rough volatility models. In fact, Wu et.al in [ 25 ] introduced the Log S-fBM model: a stochastic volatility model whose log volatility is a gaussian process, with the following autocovariance function:

Where H H and T T are respectively the Hurst exponent and the decorrelation limit. Its variance is parameterized through the parameter ν \nu . It has been proved in [ 25 ] that the volatility measure also known as the multifractal random measure of the Log S-fBM model converges to multifractal random measure of the MRW model when H H is closed to 0 0 . Consequently, the Log S-fBM model adopts a rough volatility behaviour when H H is different from 0 0 and a multifractal behaviour when H H is closed to 0 0 . The calibration of the Log S-fBM model is done in the small intermittency regime, where the squared intermittency parameter is around 0.02 0.02 which is inline with the one measure from observed market data according to the evidence in figure 11 in [ 25 ] . Leveraging the general method of moments in [ 13 ] , one can derive a small intermittency approximation of the autocovariance log volatility, then, the model parameters are determined by minimizing its discrepancy with the one measure from the observed data.

This paper aims to provide a comprehensive study of the Log S-fBM model, building upon the framework introduced in [ 25 ] . We begin by presenting the definition and main structural properties of the S-fBM process, and we establish its fundamental scaling and self-similarity features, which play a central role in understanding its multiscale behaviour. Exploiting the Gaussian structure of the S-fBM marginals, we then derive quantitative large deviation inequalities for both the S-fBM process itself and the associated Log S-fBM random measure, obtaining explicit upper bounds that characterize the concentration of logarithmic fluctuations. Then, we propose a statistical test based on Log S-fBM random measure observations in the small intermittency regime to assess the null Hurst exponent hypothesis as a way to disentangle the rough and multifractal volatility regimes. In a subsequent step, we leverage the small intermittency regime in order to derive closed-form scale invariance relations for the moments, thereby clarifying the connection between the model and its Gaussian approximation. Such scale invariance properties have previously been studied in the context of rough volatility models and multifractal random measures; here, we investigate that these properties are well reproduced in the Log S-fBM framework and show how analogous scaling structures naturally emerge in this setting. Finally, numerical experiments based on synthetically generated data are conducted to illustrate and validate the established closed form scaling formulas.

1 The Log S-fBM model

This section introduces the Log Stationary fractional Brownian motion (Log S-fBM) model proposed by Wu et al. [ 25 ] . Its construction provides a unified framework linking rough and multifractal volatility models. The presentation proceeds by first defining the Stationary fractional Brownian motion (S-fBM) relying on a geometric representation of the so called time–scale. We present its autocorrelation structure parametrized by the triplet ( H , λ 2 , T ) (H,\lambda^{2},T) respectively the Hurst exponent, intermittency coefficient and the correlation scale. Then we introduce the associated Log S-FBM random measure and its connection with multifractal random measures as its weak limit as H → 0 H\rightarrow 0 and finally embedding this construction to the Log S-fBM process.

The model is built upon a time–scale domain encoding how fluctuations at different temporal scales contribute to the value of the process at a given time. For any ( ℓ , T ) ∈ ℝ + 2 (\ell,T)\in\mathbb{R}_{+}^{2} and reference time t ∗ ∈ ℝ t^{\ast}\in\mathbb{R} , define:

This set corresponds geometrically to a truncated cone in the ( t , h ) (t,h) -plane, where t t represents time localization and h h the observation scale. The parameter T T acts as a maximal interaction scale, while ℓ \ell introduces a small-scale cutoff. In the limiting case ℓ = 0 \ell=0 , one recovers the cone C T ( . ) := C 0 , T ( . ) C_{T}(.):=C_{0,T}(.) :

which naturally appears in the construction of log-correlated Gaussian fields and multifractal cascades [ 27 , 23 ] .

Based on this geometric structure, the Stationary fractional Brownian motion (S-fBM) is defined as an integrated Gaussian noise over the time-scale domain C T ( . ) C_{T}(.) :

where d G H \mathrm{d}G_{H} denotes a non-homogeneous Gaussian white noise whose variance measure is given by

Here H ∈ ] 0 , 1 2 [ H\in\left]0,\tfrac{1}{2}\right[ is a Hurst exponent governing the local roughness of the process and λ > 0 \lambda>0 is an intermittency parameter controlling fluctuation intensity. The process is stationary since the covariance depends only on the overlap of cones translated in time. The normalization constant μ H \mu_{H} is fixed by the exponential normalization condition:

The S-fBM process is stationary and Gaussian with autocovariance function admiting the explicit expression as demonstrated in Appendix A.1 in [ 25 ] :

showing that the parameter T T induces a finite decorrelation horizon while the exponent H H determines the short-time rough behavior.

Exponentiating the Gaussian field naturally leads to the Log S-fBM random measure for H ∈ ] 0 , 1 2 [ H\in\left]0,\tfrac{1}{2}\right[ defined as follows:

which belongs to the class of Gaussian multiplicative chaos measures. The normalization ensures that the measure preserves the mean time scale, namely 𝔼 ⁡ [ M H , T ( d t ) ] = d t \mathbb{E}[M_{H,T}(\mathrm{d}t)]=\mathrm{d}t . A fundamental property established in [ 25 ] is that this construction provides a continuous interpolation between rough