Ergodicity and Mixing of Conditional Linear Random Processes in the Problems of Information Signal Modelling and Analysis Mykhailo Fryz 1,2, Serhii Kharchenko 2 and Leonid Scherbak 2 1 Ternopil Ivan Puluj National Technical University, Ruska st. 56, Ternopil, 46001, Ukraine 2 General Energy Institute of NAS of Ukraine, Antonovycha st. 172, Kyiv, 03150, Ukraine Abstract Ergodicity is a fundamental property in statistical information signal modelling and processing because it simplifies the analysis and estimation of signal properties and system parameters. It enables practitioners to work with single realizations of signals and make meaningful statistical inferences, which is necessary when dealing with real-world data and signals. The continuous- time stationary conditional linear random process as a mathematical model of information signals has been analyzed in the paper using characteristic functions method. The mixing property of the process, from which its ergodicity follows, has been proven. Keywords 1 Mathematical model, information signal, conditional linear random process, ergodicity, mixing, kernel, Levy process, characteristic function, moment functions 1. Introduction Conditional linear random process (CLRP) is represented as a stochastic integral with a random kernel driven by the process with independent increments [1, 2]. It is used for the mathematical modelling, computer simulation, statistical analysis and forecasting of information signals and processes which can be represented as a sum of many random stochastically dependent impulses occurring at Poisson time moments. The CLRP as a mathematical model of the investigated signal considers the physical nature of its generation. The CLRP is applied in the area of information systems and technology for the problems of mathematical modelling and analysis of electrophysiological information signals, radar clutter, dynamic loads of mechanical systems, forecasting of energy loads and consumptions, water consumptions, etc. [2–4]. The conditional linear random process is a generalization of a well-known model of the linear random process [5–7] having a similar integral representation but with a nonrandom kernel, and as a result they can be used only for mathematical modelling of the signals or processes represented as a sum of independent impulses. Conditional linear random processes compared with their linear counterparts take into account the conditional heteroscedasticity of modelled signals which is important for information technology applications in economics, medicine, and energy. Ergodicity is always the important property of mathematical models which is used for information signal processing when the task is to estimate parameters of a signal or a system [8–10]. Ergodicity allows to use of time averages (averages over a single realization of a signal) to estimate these parameters. This is particularly important when dealing with non-stationary signals or time-varying systems [11, 12]. Ergodicity is closely related to the concept of stationarity. The assumption of ergodicity is fundamental in the modelling of communication systems and the analysis of random noise in electrical circuits, application to financial mathematics [13], compressive sensing [14, 15], etc. Proceedings ITTAP’2023: 3rd International Workshop on Information Technologies: Theoretical and Applied Problems, November 22–24, 2023, Ternopil, Ukraine, Opole, Poland EMAIL: mykh.fryz@gmail.com (M. Fryz); nanoavia@ukr.net (S. Kharchenko); prof_scherbak@ukr.net (L. Scherbak) ORCID: 0000-0002-8720-6479 (M. Fryz); 0000-0001-9808-7607 (S. Kharchenko); 0000-0002-1536-4806 (L. Scherbak) © 2023 Copyright for this paper by its authors. Use permitted under Creative Commons License Attribution 4.0 International (CC BY 4.0). CEUR Workshop Proceedings (CEUR-WS.org) CEUR ceur-ws.org Workshop ISSN 1613-0073 Proceedings Very often the ergodic property of the investigated signal is just a hypothesis or assumption. However it is a characteristic property of linear random processes [16, 17]. Moreover, the mixing property is inherent in linear random processes [17]. We did not find such properties of conditional linear random processes in the available literature. Thus, it is important to study the same features for the CLRP model. The main goal of the paper is to justify the conditions for the continuous-time stationary conditional linear random process to be ergodic as a consequence of mixing property using the known representation of its multidimensional characteristic function. In the following parts of the article, we define the continuous-time conditional linear random process and represent its multidimensional characteristic function. Then we consider the notions of ergodicity and mixing in terms of characteristic functions of the stationary random processes. Finally, we find the conditions of stationary CLRP to be mixing and ergodic. 2. Conditional linear random process and its properties This section is preliminary and covers the definition of continuous-time conditional linear random process and representation of its multidimensional characteristic function which is used in the next section for proving the mixing property and ergodicity. Here we follow mostly the results of [1] and [2], but Levy process is used with its Poisson jump spectrum given on the Levy-Khintchine form in the CLRP model definition. A continuous-time conditional linear random process ξ(ω, t ), ω∈ Ω, t ∈ (−∞, ∞) (where {Ω, F, P} is probability space) is defined as the following stochastic integral: ∞ ξ(ω, t )= ∫ ϕ(ω, τ, t )d η(ω, τ), ω∈ Ω, t ∈ , −∞ (1) where ϕ(ω, τ, t ), τ, t ∈  is a stochastic kernel; η(ω, τ), τ∈ (−∞, ∞) is a Levy process; random functions ϕ(ω, τ, t ) and η(ω, τ) are stochastically independent. Let Fϕ ⊂ F be a σ -subalgebra generated by the random kernel ϕ(ω, τ, t ) satisfying the condition ∞ ∫ ϕ(ω, τ, t ) < ∞ with probability 1. −∞ The m-dimensional characteristic function of conditional linear random process (1) is represented in the following form:  m  tm ) E exp i ∑ uk ξ(ω, t= F f ξ (u1 , u2 ,..., um ; t1 , t2 ,...,= k ) Ef ξ (ω, u1 , u2 ,..., um ; t1 , t2 ,..., tm ), ϕ  k =1    m   ; t1 , t2 ,..., tm ) E  exp i ∑ uk ξ(ω, tk )  Fϕ  is conditional (with respect to Fϕ ) F where f ξ ϕ (ω, u1 , u2 ,..., um=   k =1     characteristic function ( Fϕ -characteristic function) of CLRP (1), which is expressed as follows:   m ∞ exp ia  ∑ uk ∫ ϕ(ω, τ, tk )d τ + Fϕ f ξ (ω, u1 , u2 , ..., u= m ; t1 , t2 , ..., t m )   k =1 −∞   ∞ ∞   m  + ∫ ∫  exp ix  ∑ uk ϕ(ω, τ, tk )   − 1 − ix ∑ m k =1 ( uk ϕ(ω, τ, tk )   1 + x2 ) dG ( x ) d τ   , (2) −∞ −∞    k =1  1 + x2  x 2     uk , tk ∈ (−∞, ∞), k = 1, m , where G ( x), x ∈  is a real non-decreasing and bounded function such that G (−∞) = 0 (the Poisson jump spectrum in Levy-Khintchine form of infinitely divisible Levy process η(ω, τ) ); ∞ a ∈  , and if Eη(ω, τ) < ∞ then a = Eη(ω, τ) − ∫ xdG ( x). −∞ If random functions (fields) ϕ(ω, τ, t ) and ϕ(ω, τ + s, t + s ) are stochastically equivalent in the wide sense, that is, their finite-dimensional distributions are equal satisfying the following condition:  n m   n m  P  {ω : ϕ(ω, τi , t j )= < xij }  P  {ω : ϕ(ω, τi + s, t j + s ) < xij }  , xij ∈  ; (3) =i 1 =j 1  =i 1 =j 1  for any s ∈  , then conditional linear random process (1) is a strict sense stationary. Using the above representation of m-dimensional characteristic function of conditional linear random process the expressions for moment functions can be obtained which are important for information signal processing (including mathematical expectation and covariance function), properties of cyclostationarity [18] can be analyzed, mixing and ergodicity conditions can be proven. The results can be also extended for multivariate case. 3. Ergodicity and mixing of stationary conditional linear random proces In this section we define the general notion of continuous-time stationary ergodic random process. We also analyze some cases which are important for the problems of information signal processing. Finally, we define and prove the mixing property for conditional linear random process. Ergodicity is the consequence of mixing property [17, 19]. Let ξ(ω, t ), t ∈ (−∞, ∞) is a continuous-time strictly stationary random process with the values in a measurable space { X , B} and g ( x1 , x2 , ..., xm ) , m ≥ 1 is a Bm -measurable function satisfying the following condition: Eg ( ξ(ω, t1 ), ξ(ω, t2 ), ..., ξ(ω, tm ) ) < ∞, ∀t1 , t2 , ..., tm ∈  . (4) The continuous-time strictly stationary random process ξ(ω, t ), t ∈ (−∞, ∞) is called ergodic if for any function g ( x1 , x2 , ..., xm ) satisfying the above conditions, the following holds with probability 1: c 1 g ( ξ(ω, t1 + t ), ξ(ω, t2 + t ), ..., ξ(ω, tm + t ) ) dt = c →∞ c ∫ lim 0 (5) = Eg ( ξ(ω, t1 ), ξ(ω, t2 ), ..., ξ(ω, tm ) ) , ∀t1 , t2 , ..., tm ∈ . The above definition of ergodic random process is very general. But in the problems of information signal analysis several cases are the most important. The mathematical expectation and covariance function, one- and two-dimensional cumulative distribution functions and characteristic functions of the investigated signal are usually estimated, analyzed, and used for informative features detection [20]. Further, the ergodicity of stationary random process with respect to the most important probability characteristics are represented indicating the relation to the general condition (5). On the below expressions m is the dimension of the corresponding function g ( x1 , x2 , ..., xm ) , convergence is assumed with probability 1. Thus, we can define the following types of ergodicity of continuous-time strictly stationary random process ξ(ω, t ), t ∈ (−∞, ∞) [17]. Ergodicity with respect to the mathematical expectation Eξ(ω, t ) =µ : c 1 c →∞ c ∫ m= 1, g ( x)= x, t1= 0 ⇒ lim ξ(ω, t )dt =µ . 0 Ergodicity with respect to the covariance function R (= τ) E ( ξ(ω, t ) − µ )( ξ(ω, t + τ) − µ )  , τ ∈  : c 1 m = 2, g ( x1 , x2=) ( x1 − µ)( x2 − µ) , t1 = 0 , t2 = τ ⇒ lim ∫ ( ξ(ω, t ) − µ )( ξ(ω, t + τ) − µ ) dt= R (τ) . c →∞ c 0 Ergodicity with respect to the one-dimensional cumulative distribution function Fξ ( y )= P (ξ(ω, t ) < y ), y ∈  : c 1 1, g ( x) = m= 0 ⇒ lim ∫ U ( y − ξ(ω, t ) dt =Fξ ( y ) , U ( y − x), t1 = c →∞ c 0 0, y ≤ 0 where U ( y ) =  is a Heaviside function. 1, y > 0 Ergodicity with respect to the two-dimensional cumulative distribution function Fξ ( y1 , y2 ; τ)= P (ξ(ω, t ) < y1 , ξ(ω, t + τ) < y2 ), y1 , y2 ∈  : m = 2, g ( x1 , x2 ) = U ( y1 − x1 )U ( y2 − x2 ) , t1 = 0 , t2 = τ ⇒ c 1 c →∞ c ∫ ⇒ lim U ( y1 − ξ(ω, t )U ( y2 − ξ(ω, t= + τ)dt Fξ ( y1 , y2 ; τ) . 0 f ξ (u ) E exp ( iuξ(ω, t ) ) : Ergodicity with respect to the one-dimensional characteristic function= c 1 g ( x) exp ( iux ) , t1= 0 ⇒ lim c →∞ c ∫ =m 1,= exp(iuξ(ω, t ))dt =f ξ (u ) . 0 Ergodicity with respect to the two-dimensional characteristic function f ξ (u1 , u2 ; τ) = = E exp i ( u1ξ(ω, t ) + u2 ξ(ω, t + τ) )  : g ( x1 , x2 ) exp [i (u1 x1 + u2 x2 )] , t1 = 0 , t2 = τ ⇒ m = 2, = c 1 ⇒ lim ∫ exp i ( u1ξ(ω, t ) + u2 ξ(ω, t + τ= ) )  dt f ξ (u1 , u2 ; τ) . c →∞ c 0 It is easy to see that stationary random process can be ergodic with respect to mathematical expectation but not ergodic with respect to covariance function or characteristic function, etc. The notion of ergodicity as it has been defined above represents the convergence of corresponding sampling averages to probability characteristics of the random process. Mixing is another fundamental property of random process expressed the fact that events (related to investigated process) separated by long time intervals are approximately independent [19]. In terms of random process distributions, mixing property means that random vectors (ξ(ω, t1 + t ), ξ(ω, t2 + t ), ..., ξ(ω, tm + t )) and (ξ(ω, s1 ), ξ(ω, s2 ), ..., ξ(ω, sn )) (created by the samples of stationary random process) become approximately independent as t → ∞ . A mixing property of the stationary random process represented by its characteristic functions is defined as follows [17]: lim f ξ (u1 , u2 , ..., um , v1 , v2 , ..., vn ; t1 + t , t2 + t , ..., tm + t , s1 , s2 , ..., sn ) = t →∞ = f ξ (u1 , u2 , ..., um ; t1 , t2 , ..., tm ) f ξ (v1 , v2 , ..., vn ; s1 , s2 , ..., sn ), (6) u1 , u2 , ..., um , v1 , v2 , ..., vn ∈ , t1 , t2 , ..., tm , s1 , s2 , ..., sn ∈ . Ergodicity in the sense of (5) is a consequence of mixing. The hierarchy of mixing and ergodicity properties which are considered in the paper has been illustrated in Figure 1. Thus, to justify the ergodic properties of continuous-time strictly stationary conditional linear random process we need to prove the mixing property first. The idea is to use the characteristic function method using the representation (2) and relationship between conditional and unconditional characteristic functions. Obviously, the time- dependent properties of the model expressed as stochastic integral driven by Levy process heavily depend on the corresponding properties of the random kernel. Let ξ(ω, t ), t ∈ (−∞, ∞) be a continuous-time strictly stationary conditional linear random process driven by Levy process, and with the kernel satisfying (3). Let us denote Law(ξ1 (ω), ξ2 (ω), ..., ξ m (ω= )) Law(η1 (ω), η2 (ω), ..., ηm (ω)) if random vectors (ξ1 (ω), ξ 2 (ω), ..., ξ m (ω)) and (η1 (ω), η2 (ω), ..., ηm (ω)) have the same m -dimensional cumulative distribution functions (distribution laws). Figure 1: Hierarchy of mixing and ergodicity properties Let random vectors (ϕ(ω, τ, t1 + t ), ϕ(ω, τ, t2 + t ), ..., ϕ(ω, τ, tm + t )) and (ϕ(ω, τ, s1 ), ϕ(ω, τ, s2 ), ..., ϕ(ω, τ, sn )) are asymptotically independent as t → ∞ , ∀τ , t1 , t2 , ..., tm , s1 , s2 , ..., sn ∈  , that is, taking into account (3), the following condition holds: lim Law ( t →∞ ϕ(ω, τ + t , t1 + t ), ϕ(ω, τ + t , t2 + t ), ..., ϕ(ω, τ + t , tm + t ), ϕ(ω, τ, s1 ), ϕ(ω, τ, s2 ), ..., ϕ(ω, τ, sn ) ) = (7) = Law ( ϕ(ω, τ, t1 ), ϕ(ω, τ, t2 ), ..., ϕ(ω, τ, tm ) ) Law ( ϕ(ω, τ, s1 ), ϕ(ω, τ, s2 ), ..., ϕ(ω, τ, sn ) ) . Then ξ(ω, t ) is the strictly stationary CLRP satisfying the mixing condition (6) ∀t1 , t2 , ..., tm , s1 , s2 , ..., sn ∈  . Indeed, (m + n) -dimensional Fϕ -characteristic function of the stationary CLRP with probability 1 can be represented as: F f ξ ϕ (ω, u1 , u2 , ..., um , v1 , v2 , ..., vn ; t1 + t , t2 + t , ..., tm + t , s1 , s2 , ..., sn ) =   m n   = E  exp  i ∑ uk ξ(ω, tk + t ) + i ∑ ν k ξ(= ω, sk )  Fϕ    k 1 =k 1   =     m ∞ n ∞  = exp ia  ∑ uk ∫ ϕ(ω, τ, tk + t )d τ + ∑ vk ∫ ϕ(ω, τ, sk )d τ  + =   k 1 = −∞ k 1 −∞    ∞ ∞   m n + ∫ ∫  exp  ix  ∑ uk ϕ(ω, τ, tk + t ) + ∑ vk ϕ(ω=   ix m ∑ , τ, sk )   − 1 − k 1 = u k ϕ((ω, τ, t k + t ) + ∑ n k 1 vk ϕ(ω, τ, sk )   × ) 2 =−∞ −∞    k 1 = k 1   1 + x     1+ x 2  × 2 dG ( x)d τ  = x    m ∞ n ∞ = exp ia  ∑ uk ∫ ϕ(ω, τ, tk + t )d τ + ∑ vk ∫ ϕ(ω, τ, sk )d τ + =   k 1 = −∞ k 1 −∞ (  ) ∞ m n + ∫ ψ ∑ uk ϕ(ω, τ, tk + t ) + ∑ vk ϕ(ω, τ, sk ) d τ  , = −∞ k 1= k 1  where ∞  iux iux  1 + x 2 ∫ (u ) ψ= −∞  e − 1 −  1 + x2  x2 dG ( x) is the logarithm of characteristic function of infinitely divisible random variable η(ω,1) on the Levy- Khintchine form when a = 0 . The function ψ (u ) is a uniformly continuous on u ∈  and ψ (0) = 0. ∞ Since a kernel of CLRP satisfies the condition ∫ ϕ(ω, τ, t ) d τ < ∞ , ∀t ∈  with probability 1, then −∞ ϕ(ω, τ, t ) → 0 as τ → ∞ , ∀t ∈  with probability 1. Taking into account (3) and (7), properties of the above function ψ (u ) , and also properties of Fϕ -characteristic function [2], the following expressions with respect to conditional and unconditional characteristic functions of CLRP can be written: lim f ξ (u1 , u2 , ..., um , v1 , v2 , ..., vn ; t1 + t , t2 + t , ..., tm + t , s1 , s2 , ..., sn ) = t →∞ F = lim Ef ξ ϕ (ω, u1 , u2 , ..., um , v1 , v2 , ..., vn ; t1 + t , t2 + t , ..., tm + t , s1 , s2 , ..., sn ) = t →∞ F = E lim f ξ ϕ (ω, u1 , u2 , ..., um , v1 , v2 , ..., vn ; t1 + t , t2 + t , ..., tm + t , s1 , s2 , ..., sn ) = t →∞ ( F F ) E f ξ ϕ (ω, u1 , u2 , ..., um ; t1 , t2 , ..., tm ) f ξ ϕ (ω, v1 , v2 , ..., vn ; s1 , s2 , ..., sn ) = = F F Ef ξ ϕ (ω, u1 , u2 , ..., um ; t1 , t2 , ..., tm )Ef ξ ϕ (ω, v1 , v2 , ..., vn ; s1 , s2 , ..., sn ) = = = f ξ (u1 , u2 , ..., um ; t1 , t2 , ..., tm ) f ξ (v1 , v2 , ..., vn ; s1 , s2 , ..., sn ). Thus, the condition (6) holds. That is, the investigated stationary conditional linear random process undergoes the mixing property. The ergodicity of the process in the sense of (5) is an immediate corollary of mixing property, including ergodicity with respect to moment functions, one- and multidimensional cumulative distribution functions and characteristic functions which are used for the feature detection and estimation in the problems of information signal processing. 4. Conclusions The conditional linear random process driven by Levy process has been defined and described using conditional characteristic functions method, Levy-Khintchine form has been used to specify the Poisson jump spectrum of Levy process. The CLRP belongs to the class of infinitely divisible mixtures. The condition of CLRP to be strict sense stationary has been represented. The hierarchy of mixing and ergodicity properties has been analyzed, the importance of the corresponding concepts for the problems of applied information signal modelling and processing has been represented. Thus, the effective mathematical model should undergo ergodicity and mixing properties. The mixing property of the continuous-time strictly stationary conditional linear random process has been proven using the characteristic functions method. The ergodicity property is a consequence of mixing. That is why, if the information signal is modelled as stationary ergodic CLRP then it is priori justification of performing the statistical analysis using time averaging of its single realization. The prospective research is related to the study of mixing and ergodicity of discrete-time conditional linear random process. 5. References [1] M. Fryz, B. Mlynko, Properties of Stationarity and Cyclostationarity of Conditional Linear Random Processes, in: Proceedings of the 2020 IEEE 15th International Conference on Advanced Trends in Radioelectronics, Telecommunications and Computer Engineering (TCSET), Lviv- Slavske Ukraine, 2020, pp. 166–170. doi:10.1109/TCSET49122.2020.235415. [2] M. Fryz, L. Scherbak, M. Karpinski, B. Mlynko, Characteristic function of conditional linear random process, in: Proceedings of the 1st International Workshop on Information Technologies: Theoretical and Applied Problems 2021 (ITTAP 2021), Ternopil Ukraine, 2021, pp. 129-135. URL: http://ceur-ws.org/Vol-3039/short40.pdf. [3] P. A. Pierre, Central limit theorems for conditionally linear random processes, SIAM J. of Applied Math. 20.3 (1971) 449-461. doi:10.1137/0120048. [4] R. Iwankiewicz, Dynamic response of mechanical systems to impulse process stochastic excitations: Markov approach, Journal of Physics: Conference Series 721 (2016). doi:10.1088/1742-6596/721/1/012010. [5] V. P. Babak, S. V. Babak, A. O. Zaporozhets, M. V. Myslovych, V. M. Zvaritch, Methods and Models for Information Data Analysis, in: Diagnostic Systems For Energy Equipments, volume 281 of Studies in Systems, Decision and Control, Springer, Cham, 2020, pp. 23-70. [6] G. E. P. Box, G. M. Jenkins, G. C. Reinsel, G. M. Ljung, Time Series Analysis: Forecasting and Control, 5th ed., Wiley, 2015. [7] P. J. Brockwell, R. A. Davis, Introduction to Time Series and Forecasting, 3rd ed., Springer Cham, 2016. doi: 10.1007/978-3-319-29854-2. [8] D. G. Kelty-Stephen, M. Mangalam, Fractal and multifractal descriptors restore ergodicity broken by non-Gaussianity in time series, Chaos, Solitons & Fractals 163 (2022) 112568. [9] L. Xu, Random Signal Estimation by Ergodicity associated with Linear Canonical Transform, in: 2019 Asia-Pacific Signal and Information Processing Association Annual Summit and Conference (APSIPA ASC), 2019, pp. 1080–1083. doi: 10.1109/APSIPAASC47483.2019.9023088. [10] A. Baxevani, K. Podgórski, Signals Featuring Harmonics With Random Frequencies – Spectral, Distributional and Ergodic Properties, IEEE Trans. Signal Process. 69 (2021) 2779–2794. doi: 10.1109/TSP.2021.3078574. [11] A. Napolitano, W. A. Gardner, Fraction-of-Time Probability: Advancing Beyond the Need for Stationarity and Ergodicity Assumptions, IEEE Access 10 (2022) 34591–34612. doi: 10.1109/ACCESS.2022.3162620. [12] M. Fryz, L. Scherbak, Statistical analysis of random coefficient periodic autoregression and its application for short-term electricity consumption forecasting, Technical Electrodynamics 2 (2019) 38-47. doi:10.15407/techned2019.02.038. [13] A. Biswas, V. S. Borkar, Ergodic risk-sensitive control—A survey, Annu. Rev. Control 55 (2023) 118–141. doi: 10.1016/j.arcontrol.2023.03.001. [14] J. F. Silva, M. S. Derpich, On the Characterization of ℓp -Compressible Ergodic Sequences, IEEE Trans. Signal Process. 63.11 (2015) 2915–2928. doi: 10.1109/TSP.2015.2419183. [15] J. F. Silva, A Compressibility Result for AMS Processes, in: 2019 IEEE Global Conference on Signal and Information Processing (GlobalSIP), 2019, pp. 1–5. doi: 10.1109/GlobalSIP45357.2019.8969106. [16] R. Passeggeri, A. E. D. Veraart, Mixing Properties of Multivariate Infinitely Divisible Random Fields, J. Theor. Probab. 32.4 (2019) 1845–1879. doi: 10.1007/s10959-018-0864-7. [17] M. Fryz, Mixing property and ergodicity of linear random processes, in: 2009 IEEE International Workshop on Intelligent Data Acquisition and Advanced Computing Systems: Technology and Applications, Rende, Italy, 2009, pp. 343–346. doi: 10.1109/IDAACS.2009.5342967. [18] A. Napolitano, Cyclostationarity: New trends and applications, Signal Processing 120 (2016) 385– 408. doi: 10.1016/j.sigpro.2015.09.011. [19] A. Khaleghi, G. Lugosi, Inferring the Mixing Properties of a Stationary Ergodic Process From a Single Sample-Path, IEEE Trans. Inf. Theory, 69.6 (2023) 4014–4026. doi: 10.1109/TIT.2023.3248375. [20] P. S. R. Diniz, Signal Processing and Machine Learning Theory, Academic Press, 2023.