This paper describes results concerning the construction of low-regularity solutions of nonlinear partial differential equations that depend on a random parameter. The motivations for this study are very varied. However, in the end, the results obtained and the methods used are conceptually very similar.

## 1 Multiple Fourier series and Sobolev spaces on the torus

For $x=(x_{1},…,x_{d})∈R_{d}$, we set

Let $T_{d}=(R/(2πZ))_{d}$ be the torus of dimension $d$. If $f:T_{d}→C$ is a function of class $C_{∞}$, then for all $x∈T_{d}$,

where $f (n)$ are the Fourier coefficients of $f$. For $s∈R$, the Sobolev norm of $f$ is defined by

For an integer $s≥0$, we have the equivalence of norms

In (2), $∂_{α}$ represents a partial derivative of order $≤s$. For $s=0$, we recover the norm of the Lebesgue space $L_{2}(T_{d})$.

The Sobolev space $H_{s}(T_{d})$ is defined as the completion of $C_{∞}(T_{d})$ with respect to the norm (1). In contrast to the case $s≥0$, for $s<0$ the elements of $H_{s}(T_{d})$ are not classical functions on the torus, but can be interpreted as Schwartz distributions. Note that the Sobolev spaces are nested: the larger $s$ is, the more regular the elements of $H_{s}(T_{d})$ are; the intersection of all $H_{s}(T_{d})$ is $C_{∞}(T_{d})$. On the other hand, the smaller $s$ is, the larger $H_{s}(T_{d})$ is; the union of all the spaces $H_{s}(T_{d})$ is the Schwartz space of $(2πZ)_{d}$-periodic distributions on $R_{d}$.

## 2 Probabilistic effects in fine questions of analysis

### 2.1 An almost sure improvement of the Sobolev embedding

Let $(Ω,F,p)$ be a probability space. Recall that a random variable $g:Ω→R$ belongs to $N(0,σ_{2})$, with $σ>0$, if the image of the measure $p$ under $g$ is

where $dx$ is the Lebesgue measure on $R$. The random variable $g$ then follows the centered normal distribution with variance $σ$. Similarly, a random variable $g:Ω→C$ belongs to $N_{C}(0,σ_{2})$, if $g=h+il$ with $h∈N(0,σ_{2})$ and $l∈N(0,σ_{2})$ independent.

Let $u∈L_{2}(T_{2})$ be a deterministic function. There exists a sequence $(c_{n})_{n∈Z}∈l_{2}(Z)$ (which consists of the Fourier coefficients of $u$) such that

Now consider a randomized version of $u$ given by the expression

where $(g_{n}(ω))_{n∈Z}$ is a sequence of independent variables of $N_{C}(0,1)$. The randomization has no effect on the Sobolev regularity of $u_{ω}$ (see, e.g., [5 N. Burq and N. Tzvetkov, Random data Cauchy theory for supercritical wave equations. I. Local theory. Invent. Math. 173, 449–475 (2008) ]). On the other hand, randomization has an important effect on regularity in Lebesgue spaces $L_{p}(T_{d})$. Using the rotation invariance of $N_{C}(0,1)$, we obtain that $g_{n}(ω)e_{inx}∈N_{C}(0,1)$, and then the independence of $g_{n}$ ensures that for a fixed $x$

Since Gaussian variables have finite moments of any order, we get

which implies that $u_{ω}(x)∈L_{p}(T)$ almost surely, a remarkable improvement in the $L_{p}$ regularity of $u_{ω}$ compared to that of $u$. Note that the Sobolev embedding requires that a deterministic function be $H_{1/2}(T)$-regular in order to conclude that it lies in $L_{p}(T)$ for any $p<+∞$ (and this restriction on regularity is optimal). In descriptive terms, randomization saves half a derivative compared with the Sobolev embedding. Like the Sobolev embedding, this effect has been known since the beginning of the 20th century, and it may seem surprising that the interaction between these two phenomena has not been studied more in the past.

Finally, thanks to the Khinchin inequality, in the preceding discussion one is allowed to replace Gaussian variables by a more general family of random variables (e.g., Bernoulli variables).

### 2.2 Almost sure products in Sobolev spaces of negative index

Let

be a random series, with $g_{n}$ as in the preceding section. It is easy to verify that almost surely $u_{ω}∈H_{σ}(T)$ for $σ<α−21 $, but almost surely $u_{ω}∈/H_{α−21}(T)$. In the following, we fix a number $σ$ such that $σ<α−21 $ (it should be assumed that this number is very close to $α−21 $). The series $u_{ω}$ is therefore in a Sobolev space of negative index and it is difficult to define an object like $∣u_{ω}∣_{2}$. After renormalization, it is nevertheless possible to give a meaning to $∣u_{ω}∣_{2}$, and even to determine its regularity in Sobolev spaces. Let us start by considering the partial sums

which are $C_{∞}$ functions. Now expand $∣u_{ω,N}(x)∣_{2}$ as

The first term of this expansion (the zero-order Fourier coefficient) contains all the singularity, while the second term has a limit almost surely in $H_{2σ}(T)$. We then set

and we define the renormalized partial sum as

The independence of the random variables $g_{n}$ ensures that the zero-order Fourier coefficient is well-defined. More precisely, we obtain

which has a limit as $N→+∞$ if $α>1/4$.

Similarly, the independence implies that the expectation

is bounded from above by a term of the order of

The latter sum is convergent under the condition $−4σ+4α>2$, which is equivalent to our assumption $σ<α−21 $. Consequently, the sequence

has a limit in $L_{2}(Ω;H_{2σ}(T))$. This limit is by definition the renormalization of $∣u_{ω}∣_{2}$. We can also establish (using more sophisticated arguments) the almost sure convergence in the Sobolev space $H_{2σ}(T)$ of the sequence (3). Note that since $σ<0$, the norm in $H_{2σ}(T)$ is weaker than the norm in $H_{σ}(T)$ (where the series $u_{ω}(x)$ is defined).

To sum up in a very informal way, the squared modulus of an element of $H_{σ}$ is in $H_{2σ}$, after renormalization. This is a remarkable probabilistic effect which lies at the heart of the study of evolutionary PDE, in the presence of randomness, in Sobolev spaces of negative index. We will further develop this topic in the remainder of this text.

## 3 Solving the nonlinear wave equation with low-regularity initial data

The wave equation is a typical example of a dispersive PDE. Solving dispersive PDE with low-regularity initial data has a long history, dating back to the seminal works of Ginibre and Velo and of Kato. Kenig, Ponce and Vega, Klainerman and Machedon, and most notably Bourgain have developed tools from harmonic analysis allowing to obtain solutions of very low regularity. The question of the optimality of these results then arose. It was Lebeau’s work that launched a series of results on the construction of counterexamples showing the optimality of the assumption of regularity in the previous results. It was in this context that the idea of proving a kind of probabilistic well-posedness for regularities where counterexamples were constructed was introduced in [23 N. Tzvetkov, Invariant measures for the defocusing nonlinear Schrödinger equation. Ann. Inst. Fourier (Grenoble) 58, 2543–2604 (2008) ], and then implemented in [5 N. Burq and N. Tzvetkov, Random data Cauchy theory for supercritical wave equations. I. Local theory. Invent. Math. 173, 449–475 (2008) , 6 N. Burq and N. Tzvetkov, Probabilistic well-posedness for the cubic wave equation. J. Eur. Math. Soc. (JEMS) 16, 1–30 (2014) ].

### 3.1 Solving the linear wave equation with periodic distributions as initial data

Consider the linear wave equation

where $t∈R$, $x∈T_{3}$, $u:R×T_{3}→R$ and $Δ$ is the Laplace operator. It is readily verified that for

the solution of (4) is given by the map $S(t)$ defined by

where $∣n∣=(n_{1}+n_{2}+n_{3})_{1/2}$. For $n=0$, the expression $∣n∣sin(t∣n∣) $ is naturally understood as its limit $t$.

Since $∣cos(t∣n∣)∣≤1$ and $∣sin(t∣n∣)∣≤1$, it follows from the above definition that

Since the map $S(t)$ is linear, we can define a unique extension of $S(t)$ for

and solve (4) with initial data in $H_{s}(T_{3})×H_{s−1}(T_{3})$ for arbitrary $s∈R$.

### 3.2 The nonlinear problem. Resolution by deterministic methods

The previous discussion makes it easy to solve (4) with singular initial data (in Sobolev spaces of arbitrary negative index). The argument is based on the a priori estimate (5) and the linear nature of the map $S(t)$ (or of equation (4)). The situation changes radically if we consider a nonlinear perturbation of (4). In this text, we restrict our attention to the case of a cubic nonlinear interaction. More precisely, we consider the problem

For this problem, the crucial information (5) and the linear nature of the equation are lost. Nevertheless, equation (6) is of Hamiltonian type. Therefore, formally, the solutions of (6) satisfy the algebraic relation

This relation implies that the Sobolev space $H_{1}(T_{3})$ is one of the natural settings for the study of problem (6). The starting point of this study is given by the following classical result.

**Theorem 3.1**

**.**

*For any pair $(u_{0},u_{1})∈H_{1}(T_{3})×L_{2}(T_{3})$
(real valued) there exists a unique global in time solution of (6) in the class*

*If, in addition, $(u_{0},u_{1})∈H_{s}(T_{3})×H_{s−1}(T_{3})$, for a given $s≥1$, then*

*Finally, the dependence on the initial data is continuous.*

Using compactness methods (going back to the work of Leray) we can exploit (7) and obtain a much weaker version of Theorem 3.1, without uniqueness and without the propagation of regularity (8). In Theorem 3.1, uniqueness results from the Sobolev embedding $H_{1}(T_{3})↪L_{6}(T_{3})$. The $L_{6}$-norm appears here naturally when we study the $L_{2}$-norm of the nonlinear term $u_{3}$. As for the propagation of the regularity, it derives from the estimate

The details of the proof of (9) can be found in [1 S. Alinhac and P. Gérard, Opérateurs pseudo-différentiels et théorème de Nash–Moser. EDP Sciences, InterEditions, Paris; et Editions du CNRS, Meudon (1991) ], where estimates of the type (9) are called tame. The key point in estimate (9) is that the $s$ derivatives acting on the expression $u_{3}$ are redistributed in such a way that, at the end, the right-hand side of (9) depends only linearly on the strong norm (the one that contains derivatives).

In view of the discussion of the linear problem (4), it is now natural to ask whether Theorem 3.1 generalizes to initial data in $H_{s}×H_{s−1}$ for a given $s<1$. As we shall see below, such a generalization is possible for some $s$, but not for all. By using Strichartz estimates instead of the Sobolev embedding $H_{1}(T_{3})↪L_{6}(T_{3})$, part of Theorem 3.1 generalizes to

More precisely, the local well-posedness of (6) can be established under the assumption (10). A more detailed description of Strichartz’s estimates would go beyond our objectives in this text. We merely point out that Strichartz estimates can be seen as improvements almost everywhere in time of Sobolev embeddings, when, instead of considering an arbitrary function, we consider a function that satisfies a dispersive PDE. We refer to [24 N. Tzvetkov, Random data wave equations. In Singular random dynamics, Lecture Notes in Math. 2253, Springer, Cham, 221–313 (2019) ] for more details on Strichartz estimates and the generalization of Theorem 3.1 under assumption (10). It can be conjectured that the global in time part of Theorem 3.1 remains true under assumption (10). The most advanced results towards the resolution of this conjecture are in [8 B. Dodson, Global well-posedness and scattering for the radial, defocusing, cubic wave equation with initial data in a critical Besov space. Anal. PDE 12, 1023–1048 (2019) , 21 T. Roy, Adapted linear-nonlinear decomposition and global well-posedness for solutions to the defocusing cubic wave equation on ℝ3. Discrete Contin. Dyn. Syst. 24, 1307–1323 (2009) ].

### 3.3 The limitations of deterministic methods

The restriction (10) is optimal with respect to the well-posedness in the sense Hadamard of the problem (6) with initial data in $H_{s}(T_{3})×H_{s−1}(T_{3})$. More precisely, we have the following result.

**Theorem 3.2**

**.**

*Let $s∈(0,1/2)$ and $(u_{0},u_{1})∈H_{s}(T_{3})×H_{s−1}(T_{3})$. There exists a sequence*

*such that*

*with*

*but for all $T>0$*

Well-posedness in the sense of Hadamard requires existence, uniqueness and continuous dependence on the initial data. Theorem 3.2 shows that continuous dependence on initial data fails.

The proof of Theorem 3.2 is based on an idea of Lebeau (see for example [15 G. Lebeau, Perte de régularité pour les équations d’ondes sur-critiques. Bull. Soc. Math. France 133, 145–157 (2005) ]): if the initial data are localized at high frequency, then for small times a good approximation of the solution of (6) is given by the solution of

which is obtained from (6) by neglecting the effect of the Laplacian. In other words, under the hypothesis of Theorem 3.2, nonlinear effects dominate in the regime described above. The solution of (11) manifests the phenomenon of amplification described by Theorem 3.2 and this property propagates to the solutions of (6) by a perturbative, highly non-trivial argument. A detailed proof of Theorem 3.2 can be found in [24 N. Tzvetkov, Random data wave equations. In Singular random dynamics, Lecture Notes in Math. 2253, Springer, Cham, 221–313 (2019) ].

### 3.4 Resolution by probabilistic methods beyond the limitations of deterministic theory

Despite the result of Theorem 3.2, we can ask whether a form of the well-posedness of (6) remains true for initial data in

The answer to this question is positive if we endow the space (12) with a non-degenerate probability measure such that we have existence, uniqueness and (a form of) continuous dependence almost surely with respect to this measure.

We will choose the initial data for (6) from the realizations of the following random series:

Here $(g_{n})_{n∈Z_{3}}$ and $(h_{n})_{n∈Z_{3}}$ are two families of independent random variables conditioned by $g_{n}=g_{−n} $ and $h_{n}=h_{−n} $, so that $u_{0}$ and $u_{1}$ are real valued. Furthermore, it is assumed that for $n=0$, $g_{n}$ and $h_{n}$ are complex Gaussians with distribution $N_{C}(0,1)$, while $g_{0}$ and $h_{0}$ are real Gaussians with distribution $N(0,1)$.

The partial sums associated with (13) are Cauchy sequences in $L_{2}(Ω;H_{s}(T_{3})×H_{s−1}(T_{3}))$ for all $s<α−23 $. Therefore, the initial data (13) belong almost surely to $H_{s}(T_{3})×H_{s−1}(T_{3})$ for $s<α−23 $. Furthermore, the probability of the event

is zero. It follows that for $α>5/2$ we can apply Theorem 3.1 to the data $(u_{0},u_{1})$ described by (13). For $α>2$, we can apply refined deterministic results (based on Strichartz estimates). Finally, for $α∈(3/2,2)$, Theorem 3.2 applies and we get:

**Theorem 3.3**(pathological approximations)

**.**

*Let $α∈(3/2,2)$ and $0<s<α−3/2$. For almost every $ω$,
there exists a sequence*

*such that*

*with*

*but for all $T>0$*

However, the following result also holds.

**Theorem 3.4**(probabilistic well-posedness)

**.**

*Let $α∈(3/2,2)$ and $0<s<α−3/2$. Using Theorem 3.1, define the sequence
$(u_{N})_{N≥1}$ of solutions of problem (6) with regular initial conditions given by*

*Then there exists a set $Σ⊂Ω$ of probability $1$ such that for every $ω∈Σ$ the sequence $(u_{N})_{N≥1}$ converges when $N→+∞$ in
$C_{0}(R;H_{s}(T_{3}))$ to a (unique) limit that satisfies (6) in the sense of distributions.*

Theorems 3.3 and 3.4 show that the type of approximation of the initial data is crucial when establishing the probabilistic well-posedness.

Using compactness methods (à la Leray), we can hope to obtain convergence of a subsequence of $(u_{N})_{N≥1}$. The convergence of the whole sequence $(u_{N})_{N≥1}$ is beyond the reach of weak-solutions techniques. The fact that the whole sequence converges already contains a form of uniqueness. In [6 N. Burq and N. Tzvetkov, Probabilistic well-posedness for the cubic wave equation. J. Eur. Math. Soc. (JEMS) 16, 1–30 (2014) ], one can find a form of uniqueness that can be formulated in a suitable functional framework.

In [6 N. Burq and N. Tzvetkov, Probabilistic well-posedness for the cubic wave equation. J. Eur. Math. Soc. (JEMS) 16, 1–30 (2014) ], we also obtain a probabilistic continuous dependence on the initial data, the proof of which makes use of conditioned large deviation properties which seem to be of independent interest.

We can prove the result of Theorem 3.4 for more general randomizations than (13). For example, Gaussian variables can be replaced by Bernoulli variables and the deterministic coefficients $⟨n⟩_{−α}$ by other coefficients with “similar” behaviour for $∣n∣≫1$ (see [6 N. Burq and N. Tzvetkov, Probabilistic well-posedness for the cubic wave equation. J. Eur. Math. Soc. (JEMS) 16, 1–30 (2014) ] for more details).

Theorem 3.4 provides a nice dense set $Σ$ of initial data such that for good approximations we get nice global solutions (but for bad approximations we get divergent sequences, as shown by Theorem 3.3!). On the other hand, due to [22 C. Sun and N. Tzvetkov, Concerning the pathological set in the context of probabilistic well-posedness. C. R. Math. Acad. Sci. Paris 358, 989–999 (2020) ], we also have a dense set of bad initial data, even for the natural approximations by Fourier truncation (or convolution).

**Theorem 3.5**(pathological initial data)

**.**

*Let $0<s<21 $. Then there is
a dense set $S⊂H_{s}(T_{3})×H_{s−1}(T_{3})$ such that for every $(f,g)∈S$, the sequence $(u_{N})_{N≥1}$ of (smooth) solutions of*

*with data*

*does not converge. More precisely, for every $T>0$,*

One may even prove that the pathological set $S$ contains a dense $G_{δ}$ set and consequently the good-data set is not generic in the sense of Baire.

### 3.5 Going even further

For $α<3/2$, $u_{0}$ is no longer a classical function. In this case, it can be interpreted as a distribution belonging to a Sobolev space with negative index. We cannot expect a result like that of Theorem 3.4 for $α<3/2$. A renormalization is necessary, as shown by the following result established in [19 T. Oh, O. Pocovnicu and N. Tzvetkov, Probabilistic local Cauchy theory of the cubic nonlinear wave equation in negative Sobolev spaces. Ann. Inst. Fourier (Grenoble) 72, 771–830 (2022) ].

**Theorem 3.6**

**.**

*Let $α∈(45 ,23 )$ and $s<α−3/2$. There
exists positive constants $γ,c,C,T_{0}$ and a divergent
sequence $(c_{N})_{N≥1}$ such that for any $T∈(0,T_{0})$,
there exists a set $Ω_{T}$ such that the probability of its complement is smaller than
$Cexp(−c/T_{γ})$ and such that if we denote by $(u_{N})_{N≥1}$ the solution of*

*with initial data given by (14), then for any $ω∈Ω_{T}$, the sequence $(u_{N})_{N≥1}$ converges
for $N→+∞$ in $C_{0}([−T,T];H_{s}(T_{3}))$. In particular, for almost every $ω$
there exists $T_{ω}>0$ such that $(u_{N})_{N≥1}$ converges in $C_{0}([−T_{ω},T_{ω}];H_{s}(T_{3}))$.*

Theorem 3.6 was a first step in the study of problem (6) in Sobolev spaces of negative index. In the remarkable recent work [4 B. Bringmann, Invariant Gibbs measures for the three-dimensional wave equation with a Hartree nonlinearity II: Dynamics. J. Eur. Math. Soc. (JEMS) online first (2023) DOI 10.4171/jems/1317 ] the result of Theorem 3.6 was extended to the range $α>1$.

## 4 Invariant measures for the nonlinear Schrödinger equation

Let us now consider the nonlinear Schrödinger equation, posed on the torus of dimension two:

Here the solution $u$ is complex valued, but the equation is of first order in time. We have the following analogue of Theorem 3.1 in the context of (15).

**Theorem 4.1**

**.**

*For any $u_{0}∈H_{1}(T_{2})$ there exists a unique global solution of (15) in the class
$C_{0}(R;H_{1}(T_{2}))$. If moreover $u_{0}∈H_{s}(T_{2})$ for some $s≥1$, then
$u∈C_{0}(R;H_{s}(T_{2}))$. The dependence on the initial data is also continuous.*

Equation (15) is again Hamiltonian in nature. This implies that the functional

is preserved by (15). The Gibbs measure associated with (16) is a “renormalization” of the completely formal object

This renormalization is a classic procedure in quantum field theory, which would be impossible to present in this short text. Let us just say that the measure obtained by this renormalization is absolutely continuous with respect to the Gaussian measure induced by the random series

where $(g_{n})_{n∈Z_{2}}$ is a family of independent complex Gaussians with distribution $N_{C}(0,1)$. Once the measure (17) has been rigorously defined, the natural question is whether we can define a dynamics related to (15) that leaves this measure invariant. The answer to this question is given by the work of Bourgain [3 J. Bourgain, Periodic nonlinear Schrödinger equation and invariant measures. Comm. Math. Phys. 166, 1–26 (1994) ]. The difficulty lies in the fact that (18) does not define a classical function. The object defined by (18) almost surely belongs to the Sobolev space $H_{s}(T_{2})$ for all $s<0$. Such a regularity implies that Theorem 4.1 cannot be applied in the context of an initial datum given by (18). This regularity is also beyond the reach of the most sophisticated deterministic techniques. Nevertheless, the following statement can be deduced from [3 J. Bourgain, Periodic nonlinear Schrödinger equation and invariant measures. Comm. Math. Phys. 166, 1–26 (1994) ].

**Theorem 4.2**

**.**

*Using Theorem 4.1, define the sequence $(u_{N})_{N≥1}$ of solutions of (15) for the
initial conditions of class $C_{∞}$ given by*

*Then for any $s<0$ the sequence*

*converges almost surely in $C_{0}(R;H_{s}(T_{2}))$
to a limit which satisfies (in the sense of distributions) a renormalized version of problem (15).*

There is a similarity between Theorems 3.4 and 4.2. One notable difference is the need to renormalize the sequence $(u_{N})_{N≥1}$ of Theorem 4.2 in order to obtain a limit. This renormalization is linked to the construction of the measure from the formal object (17) mentioned above.

The result of Theorem 4.2 can be formulated in the spirit of Theorem 3.6. More precisely, we can establish the convergence of the solutions of the problem

with initial data (19), where $(c_{N}(ω))_{N≥1}$ is a sequence of real numbers almost surely divergent to $+∞$.

## 5 Singular stochastic PDEs

The issues considered in the previous sections are very close to the analysis of PDE in the presence of a singular random source (noise). This topic has received a lot of attention in recent years (see, for example, [7 G. Da Prato and A. Debussche, Strong solutions to the stochastic quantization equations. Ann. Probab. 31, 1900–1916 (2003) , 9 M. Gubinelli, A panorama of singular SPDEs. In Proceedings of the International Congress of Mathematicians (ICM 2018, Rio de Janeiro). Vol. III. Invited lectures, World Scientific Publishing, Hackensack, NJ, 2311–2338 (2018) , 10 M. Gubinelli, P. Imkeller and N. Perkowski, Paracontrolled distributions and singular PDEs. Forum Math. Pi 3, article no. e6 (2015) , 12 M. Hairer, A theory of regularity structures. Invent. Math. 198, 269–504 (2014) , 11 M. Hairer, Singular stochastic PDEs. In Proceedings of the International Congress of Mathematicians (ICM 2014, Seoul). Vol. IV, Kyung Moon Sa, Seoul, 49–73 (2014) , 14 A. Kupiainen, Renormalization group and stochastic PDEs. Ann. Henri Poincaré 17, 497–535 (2016) ]). The closest equation to those in the previous sections is the nonlinear heat equation. More precisely, we consider the problem

In this equation, $ξ$ is the space-time white noise on $[0,+∞[×T_{3}$. The unknown $u$ is a real valued function. There are many physical motivations for considering a PDE perturbed by white noise. A serious discussion of these motivations is beyond the scope of this paper.

The source term $ξ$ represents the singular randomness in (20), while in (6) and (15), the initial datum is the source of the singular randomness. A little experience with the analysis of evolutionary PDE is enough to know that the two situations are very similar and even that, in some cases, for reasons of convenience, we can easily transform the problem with initial data into a problem with a source term and zero as initial data.

A representation in the spirit of (13) and (18) of white noise on $[0,+∞[×T_{3}$ is given by the formula

where $β_{n}$ are independent Brownian motions, conditioned by $β_{n}=β_{−n} $ ($β_{0}$ is real and for $n=0$, $β_{n}$ is complex valued). The derivative with respect to $t$ of $β_{n}$ in (21) is in the sense of distributions.

If $ξ∈C_{∞}([0,∞[×T_{3})$, equation (20) can be solved by deterministic methods. This is the analogue of Theorems 3.1 or 4.1 in the context of (20). For $N≫1$, an approximation of $ξ$ given by (21