Self-similar processes are crucial for modeling complex systems that exhibit scale-invariant or fractal-like behavior. These processes have been studied extensively due to their theoretical significance and practical applications in various fields, including network traffic modeling, financial time series analysis, geophysics, and hydrology. Significant contributions to this area include the works of Mandelbrot and Van Ness [14], Samorodnitsky and Taqqu [19], and Embrechts and Maejima [9]. A key feature of self-similar processes is their capacity to capture long-range dependence, where correlations persist over large time intervals.

On the other hand, processes exhibiting family-wise scaling are vital in modern stochastic modeling because they offer a more flexible and realistic framework for systems where scaling behavior isn’t simple and fixed. While traditional self-similarity describes a system that looks statistically the same at all scales, family-wise scaling accounts for a deeper, more complex relationship where the very nature of the process changes when it is scaled, all while remaining part of the same “family” of processes. Jørgensen, Martínez and Demétrio [12] introduced a framework for self-similarity and the Lamperti convergence in one-parameter families of processes. Also, Borgnat, Amblard and Flandrin [5] discussed scale invariances and Lamperti transformations as tools for understanding how time rescaling interacts with parameter changes. These perspectives justify distinguishing “family-wise scaling” from the standard notion of self-similarity of a fixed-parameter process. Processes exhibiting family-wise scaling have importance and practical benefits. From a theoretical point, studying processes with family-wise scaling expands the scope of stochastic analysis beyond the well-understood self-similar processes like fractional Brownian motion. It allows researchers to develop new mathematical tools and theories to understand a broader range of complex, nonstandard behaviors. This contributes to the development of more sophisticated models that can handle a wider array of real-world data. For instance, in financial markets, volatility and other statistical properties of asset returns can change over different time scales. A traditional self-similar model might miss these nuances. A family-wise scaling model, however, can capture this by allowing a parameter (e.g., related to volatility) to change with the time scale, providing a more accurate and robust model for risk assessment. These processes enable us to deal with non-Gaussian phenomena. For example, many natural and human processes, such as network data traffic and seismic activity, exhibit sudden jumps or shocks that cannot be modeled by Gaussian (natural) processes.

Long-memory and scaling phenomena play a central role in stochastic processes and time series. Classical models often rely on self-similar processes and their stationary increments, such as fractional Brownian motion and fractional Gaussian noise. Extensions driven by Poisson-type or Lévy-type noises provide a richer framework for modeling discontinuities and heavy-tailed behavior.

Stationary long-range dependent (LRD; this abbreviation will also stand for ‘long-range dependence’) series can be broadly classified into three main types, each with distinct mechanisms for LRD, jumps, and scaling. Lévy-based series, including fractional Lévy noise and Lévy-driven moving averages, exhibit LRD via slowly decaying kernel correlations, with jumps intrinsic to the Lévy process, and show approximate scaling: X ( c t ) ∼ c H X ( t ) for large c [6, 7]. Time-changed Gaussian processes generate LRD through memory in the stochastic time-change τ ( t ); if τ ( t ) is self-similar with index ν, then X ( c t ) = d c ν / 2 X ( t ) [8]. Counting-based Poisson-type series, such as fractional Poisson and mixed Poisson models, produce LRD through heavy-tailed interarrival times or correlated random intensity, with jumps intrinsic to the counting mechanism, but generally lack scaling except in certain asymptotic regimes [2, 13, 15].

In this paper, we introduce a family-type scaling process, where scaling the time variable by a factor a > 0 produces a scaled version of another process from the same family, rather than the same process. That is, the process parameters themselves are affected by scaling, distinguishing this construction from classical self-similar processes. Based on this process, we construct, under certain conditions, a stationary time series exhibiting long-range dependence. This framework extends classical results on a fractional Gaussian noise and Lévy or Poisson type series, providing a versatile tool for spectral analysis of one-sided long-memory processes.

The paper is organized as follows. In Section 1, we review the concepts of self-similarity and long-range dependence, followed by an overview of classical models such as fractional Gaussian noise (fGn). We then introduce the notion of pseudo-Poisson type processe, which serve as the driving process in our construction. In Section 2, we introduce our approach to family-type scaling processes, and then present the construction of such a process as the integral of a pseudo-Poisson type process. The Gamma distribution of the random intensity in the underlying Poisson process provides a natural mechanism for heavy-tailed behavior and significantly influences the statistical properties of the model. We refer to this process as fractional PSI-process. Its covariance structure, as given in Proposition 1, reflects long-range dependence with a power-law decay governed by the exponent parameter H ∈ 1 2 , 1 . This implies persistent correlations over long time intervals, making the model suitable for capturing memory effects observed in various real-world phenomena such as network traffic, financial time series, and environmental data.

The fractional PSI-process Ψ θ ( t ) , t ∈ R + , inherits several important characteristics from both the underlying fractional Brownian motion and the pseudo-Poisson time change. In Section 3, we study a time series defined on Z + , constructed from the stationary increments of a fractional PSI-process. The main contribution of this paper is to establish the long-memory behavior of the resulting stationary time series. Specifically, we show that its autocovariance function exhibits a power-law decay of the form γ X ( h ) = L ( h ) h 2 d − 1 , d ∈ ( 0 , 1 / 2 ) , where L ( h ) is a slowly varying function at infinity. Also, we prove the functional central limit theorem for LRD series X n n ∈ Z + in the Skorokhod space D [ 0 , 1 ]. Finally, in Section 4, we derive the spectral density and discuss its implications for the long-memory behavior. Although a stationary time series defined on Z + is not strictly causal, its one-sided domain allows for the application of the unilateral Fourier transform to derive its spectral density. Nevertheless, the lack of symmetry and causality introduces additional technical challenges in both the spectral representation and its interpretation. By using distributional Fourier analysis, we show that the resulting spectral density is, in general, a complex-valued and nonsymmetric function. As a result, the corresponding spectral measure is not classical but must be understood in the sense of a generalized spectral measure, specifically, as a tempered distribution.

In this section, we will briefly review key concepts and results that are essential for our construction and analysis. All concepts and definitions presented here can be found in [16]. If a different source is referenced, it will be indicated explicitly.

Typically, one selects k = 1 when examining the properties of the stationary time series Y n = X ( n + 1 ) − X ( n ).

As is known, the most famous self-similar process with stationary increments is a fractional Brownian motion (fBm). The fBm, denoted by B H = B H ( t ) t ∈ R , is Gaussian SSSI process with Hurst parameter H ∈ ( 0 , 1 ), zero mean and covariance function, for s , t ∈ R, given by E B H ( t ) B H ( s ) = E B H ( 1 ) 2 2 | t | 2 H + | s | 2 H − | t − s | 2 H .

A stochastic process X ( t ) t ∈ T is called weakly(or second-order) stationary if for all s , t ∈ T, E X ( t ) = E X ( 0 ) , and Cov ( X ( t ) , X ( s ) ) = Cov ( X ( t − s ) , X ( 0 ) ) , where T = R or Z.

The autocovariance function of a weakly stationary time series X n n ∈ Z is defined as γ X ( h ) = Cov ( X n + h , X n ) , h ∈ Z . By weak stationarity, γ X ( h ) depends only on the lag h. The autocovariance function serves as a fundamental measure of dependence in time series analysis and it has a symmetry property, i.e., γ X ( h ) = γ X ( − h ) , h ∈ Z.

The spectral density function is a fundamental tool in analyzing weakly stationary time series, providing a frequency-domain perspective of their dependence structure. It is especially important for detecting and characterizing long-range dependence, as it reveals how the variance of the series is distributed over low-frequency components. This insight helps identify persistent correlations over long time horizons, which are crucial for accurate modeling and forecasting in many fields of application.

The parameter d ∈ ( 0 , 1 / 2 ) is called a long-range dependence parameter.

Thus, although the four conditions capture different perspectives – covariance decay, spectral density, variance of partial sums, and nonsummability – they are tied together through a network of implications, providing a unified but flexible framework for analyzing long-memory behavior. This is shown in Table 1.

The fractional Gaussian noise is a zero-mean stationary process with variance E X n 2 = E [ B H ( 1 ) 2 ] , and the autocovariance function given by (6) γ ( k ) = E [ B H ( 1 ) 2 ] 2 | k + 1 | 2 H + | k − 1 | 2 H − 2 | k | 2 H (7) ∼ E [ ( B H ( 1 ) ) 2 ] · H ( 2 H − 1 ) k 2 H − 2 , as k → ∞ . According to [16, 20], the spectral density of fGn is given by f X ( λ ) = c H | 1 − e − i λ | 2 ∑ k = − ∞ ∞ | 2 π k + λ | − 1 − 2 H , for λ ∈ ( − π , π ] , where c H = E [ B H ( 1 ) 2 ] 2 π Γ ( 2 H + 1 ) sin ( π H ).

Moreover, the spectral density exhibits a power-law divergence near the origin and behaves asymptotically as f X ( λ ) ∼ c H λ 1 − 2 H , as λ → 0 , for H > 1 2 .

In his foundational monograph [11], Feller introduced the concept of a pseudo-Poisson process, through what he called Poissonian randomization of time. This approach involves subordinating the index set of a sequence of random variables ( ξ n ) n ≥ 0 using a Poisson process, typically assumed to be independent of the underlying sequence, which itself exhibits the Markov property.

In this paper, we consider a generalization of this construction in which the index time of the sequence ( ξ n ) n ≥ 0 is subordinated by a Poisson process Π λ ( ω ) with a random intensity λ ( ω ). This type of model has been studied by Rusakov, Yakubovich and Baev [17], as well as by Rusakov, Yakubovich, and Laskin [18].

Note that the standardized Poisson process with rate λ ( ω ) can be obtained by a random time-scaling of the unit-rate process: Π λ ( ω ) ( s ) = d Π 1 ( λ ( ω ) s ) .

Throughout this work, we assume that ( ξ n ) n ≥ 0 is a sequence of independent and identically distributed (i.i.d.) random variables with zero mean and finite variance E ξ 0 2 = σ 2 < ∞ and that λ ( ω ) is a positive random variable. As shown in [17, 18], the mean and covariance function of the PSI-process are given by (9) E ψ ( s ) = E ξ 0 = 0 , Cov ( ψ ( s ) , ψ ( t ) ) = E ξ 0 2 L λ ( | t − s | ) , s , t ∈ R + , where L λ ( · ) denotes the Laplace transform of the random variable λ ( ω ).

For simplicity, we denote λ ( ω ) by λ.

In this paper, we introduce a new class of stochastic models that extends the classical notion of self-similarity, referred to as family-wise scaling processes. Up to now, there has been no rigorous or unified definition for processes whose scaling property involves not only time rescaling but also a simultaneous rescaling of an additional parameter. The existing literature mainly provides precise definitions for self-similar processes (characterized by a single Hurst exponent). However, none of these frameworks captures the family-wise relation that links the dynamics of the process with the simultaneous transformation of its family parameter. Therefore, the definition we introduce here provides the formulation of such a family-wise scaling processes.

Instead of the classical self-similarity condition, here the scaling also acts on the parameter θ.

In this section, we introduce a stochastic model that combines two important features: family-wise scaling and randomness in the time index through a pseudo-Poisson-type process. This model introduced by Azzo and Rusakov in [1] is defined by (11) Ψ θ ( t ) = ∫ 0 t θ κ / 2 ξ Π θ λ ( s ) d s , t ∈ [ 0 , T ] , where θ ∈ ( 0 , ∞ ) and κ ∈ ( 0 , 1 ).

The process Ψ θ ( t ) can be interpreted as a randomly time-changed integral of a scaled discrete-time noise sequence { ξ n } n ≥ 0 , where the time distortion is governed by a Poisson process Π θ λ ( s ) with random intensity λ, scaled by the parameter θ. The parameter θ > 0 controls the degree of aggregation or dispersion in the time deformation, while κ ∈ ( 0 , 1 ) determines the scaling behavior of the integrand. Moreover, since the Poisson process and the noise sequence are independent, the process Ψ θ ( t ) is well defined and integrable under suitable assumptions on ξ 0 and λ.

This model can be applied to the analysis of stochastic systems in which random events generate workload or transmit information over time. Examples include telecommunications networks, financial systems, and queueing models. The process reflects the total load that is the accumulated effect or quantity of events occurring within a specified period – offering important insights into system efficiency, resource allocation, and risk evaluation.

Moreover, the stochastic process defined in (11), indexed by parameter θ, exhibits a family-wise scaling property with exponent parameter H ∈ ( 1 2 , 1 ), as opposed to the classical H-self-similarity. Specifically, a linear scaling of time by a factor a corresponds, in the distributional sense, to a deterministic rescaling of another process from the same family, indexed by a θ, i.e., for any a > 0, (12) Ψ θ ( a t ) = d a H Ψ a θ ( t ) , t ∈ R + . This property can be verified through a straightforward calculation. Ψ θ ( a t ) = ∫ 0 a t θ κ / 2 ξ Π θ λ ( s ) d s = d ∫ 0 a t θ κ / 2 ξ Π 1 ( θ λ s ) d s = a ∫ 0 t θ κ / 2 ξ Π 1 ( θ λ a u ) d u = a 1 − κ / 2 ∫ 0 t ( a θ ) κ / 2 ξ Π 1 ( a θ λ u ) d u = a H Ψ a θ ( t ) , where the scaling exponent H = 1 − κ 2 ∈ 1 2 , 1 .

We are interested in analyzing the covariance structure of the process defined in (11), under the assumption that the random intensity λ of the Poisson process follows a gamma distribution with unit scale and shape parameter κ ∈ ( 0 , 1 ).

Denote the covariance function of Ψ θ ( t ) by K Ψ θ ( s , t ) : = Cov Ψ θ ( s ) , Ψ θ ( t ) .

Observe that for H ∈ 1 2 , 1 , the constant factor σ 2 2 H ( 2 H − 1 ) is strictly positive. Consequently, we refer to the process Ψ θ ( t ) , t ∈ R + , defined in (11) and satisfying the assumptions of Proposition 1, as the fractional PSI-process with exponent parameter H = 1 − κ 2 ∈ 1 2 , 1 .

The following corollary follows directly from the stationary increments property of the fractional PSI-process and highlights that the increment process is weakly stationary, with constant variance and an autocovariance function depending solely on the time lag. Corollary 1.

Let Ψ θ ( t ) , t ∈ R + , be the fractional PSI-process. For any fixed h > 0, the increment process Ψ θ ( t + h ) − Ψ θ ( t ) is weakly stationary. Its autocovariance depends only on the lag | t − s | and the fixed increment length h.