1 Introduction and main result
Stochastic partial differential equations (SPDEs) do not only provide key models in modern probability theory, but also become increasingly popular in applications, for instance, in neurobiology or mathematical finance. Consequently, statistical methods are required to calibrate SPDE models from given observations. However, in the statistical literature on SPDEs, see
[5] for a recent review, there are still basic questions which are not yet settled.A natural problem is parameter estimation based on discrete observations of a solution of an SPDE which was first studied in
[10] and which has very recently attracted considerable interest. Applying similar methods the three related independent works [6, 2, 4] study parabolic SPDEs including the stochastic heat equation, consider highfrequency observations in time, construct estimators using power variations of timeincrements of the solution and prove central limit theorems. As we shall see below, the marginal solution process along time at a fixed spatial point is not a (semi)martingale such that the wellestablished highfrequency theory for stochastic processes from [8] cannot be (directly) applied. In view of this difficulty, different techniques are required to prove central limit theorems. Interestingly, the proof strategies in [6, 2, 4] are quite different. Cialenco and Huang [6] consider the realised fourth power variation for the stochastic heat equation with both an unbounded spatial domain , or a bounded spatial domain . In the first setting they apply the central limit theorem by Breuer and Major [3] for stationary Gaussian sequences with sufficient decay of the correlations. For, they use Malliavin calculus instead and the fourth moment theorem from
[12]. Also in case of a bounded domain , with Dirichlet boundary conditions, Bibinger and Trabs [2] study the normalized discrete quadratic variation and establish its asymptotic normality building upon a theorem by Peligrad and Utev [14] for triangular arrays which satisfy a covariance inequality related to mixing. Finally, Chong [4] has proved (stable) central limit theorems for power variations in the case , based on a nonobvious martingale approximation in combination with the theory from [8]. The strategy of proofs by [2] and [4] do not directly rely on a purely Gaussian model and can be transferred to more general settings. While [2] considers further nonparametric inference on a timevarying deterministic volatility, [4] already provides a proof beyond the Gaussian framework including stochastic volatility.This note presents a concise analysis which transfers the asymptotic theory from [2] to an unbounded spatial domain , and from the normalized discrete quadratic variation to general power variations. Contrarily to [2], we do not start with the illustration of a solution as an infinitedimensional SDE but exploit the explicit representation of the solution with the heat kernel thanks to the continuous spectrum of the Laplace operator on the whole real line. We stick here to the simplest Gaussian setting to illustrate the main aspects and deviations from the classical theory. Our findings show that the central limit theorem under a mixing type condition used in [2] for the case with a bounded spatial domain can be used likewise for this different model with unbounded spatial domain. We moreover expect that it provides a perspective to prove central limit theorems very generally, although many approximation details, for instance, to address stochastic volatility, remain far from being obvious. We consider the stochastic heat equation in one spatial dimension
(1) 
for spacetime white noise
, and with parameters , and some initial condition which is independent of . is defined as a centred Gaussian process with covariance structure , and is in terms of a distribution the spacetime derivative of a Brownian sheet. Since the Laplace operator on the whole real line does not have a discrete spectrum and we do not have to discuss boundary problems, the asymptotic analysis actually simplifies compared to [2] and allows for more transparent proofs.A mild solution of (1) is a random field that admits the representation
(2) 
for , where the integral is welldefined as the stochastic Walsh integral and with
is the heat kernel, the fundamental solution to the heat equation. Let us refer to [9, Ch. 2.3.1] for an introduction to the heat equation and SPDEs in general. Suppose we observe this solution on a discrete grid , at equidistant observation times . We consider infill or highfrequency asymptotics where . For statistical inference on the parameters in (1), the key quantities to study are power variations
with . The normalization of with takes into account the (almost) Hölder regularity in time of , see [9, Ex. 2.3.5]. By homogeneity in space, statistics to consider for volatility estimation are spatial averages
(3) 
The main result of this note is a central limit theorem for in the double asymptotic regime where and (possibly) . An important role in our asymptotic analysis is played by the secondorder increment operator for some function , being well defined on . For brevity we assume , but the result readily extends to sufficiently regular initial conditions which are independent of .
Theorem 1.
Consider (1) with . For assume that as . Then the power variations from (3) with satisfy as and
with , and with for
jointly normally distributed with expectation 0, variances 1 and correlation
.Note the explicit formula , also referred to as for even. In particular for , that is, for the normalized discrete quadratic variation, we have and the asymptotic variance is
in analogy with Example 2.11 in [4] and with [2]. This coincides with the variance of the normalized discrete quadratic variation of a fractional Brownian motion with Hurst exponent and scale parameter , see also Theorem 6 in [1] and [11].
The above result allows for a growing time horizon and, more general than in [2], the number of spatial observations in the unbounded spatial domain can be larger than the number of observation times . The relevant condition that induces decorrelated observations in space is , tantamount to a finer observation frequency in time than in space. Based on Theorem 1, one can construct estimators and confidence statements for the parameters and , if the other one is known, see [6, 2, 4]. If no parameter is known apriori, [2, Sec. 5] show that the “viscosityadjusted volatility” can be estimated consistently, also noted in [4, Sec. 2.3].
2 Highfrequency asymptotic analysis of power variations
Our analysis builds upon the following result, whose proof is postponed to Section 3.
Proposition 2.
For with , we have that
The increments thus have nonnegligible covariances and is not a (semi)martingale. The terms will turn out to be asymptotically negligible in the variance of the power variations. Since secondorder differences of the square root decay as its second derivative, we observe that . This motivates an asymptotic theory exploiting mixing arguments. From the proposition and joint normality of the increments, we readily obtain the expectation and variance of the power variations at one spatial point .
Corollary 3.
For any , we have that
E[V_n^p(x)]=(2πϑ)^p4σ^pμ_p+O(n^1) and
Var(V_n^p(x))=1n(2πϑ)^p2σ^2p((μ_2pμ_p^2)+2∑_r=2^∞ρ_p(12 D_2(⋅,r)))+O(1n)
with , and with for jointly centred Gaussian with variances 1 and correlation .
Proof.
For , Proposition 2 yields . Since , we obtain by a Taylor expansion that
Using the joint normality of the increments , and writing , with a tight sequence , we deduce for any that
(4)  
By the above bound, the term with is negligible such that up to this negligible term. For the covariance terms, we use Proposition 2 to obtain
The first equality comes from approximating the variances by one and the second approximation is based on the Hermite expansion of absolute power functions (15) with Hermite rank 2, see also [1, (A.6)]. The last estimate follows from
∎
As we can see from the previous proof, the term in the variance would also appear for independent increments, while the additional term involving comes from the nonvanishing covariances. Proposition 2 moreover implies that the covariance of and decreases with a growing distance of the spatial observation points and . In particular, averaging over all spatial observations in (3) reduces the variance by the factor , as long as the highfrequency regime in time dominates the spatial resolution. The next corollary determines the asymptotic variance in Theorem 1.
Corollary 4.
Under the conditions of Theorem 1, we have that
Proof.
For bivariate Gaussian with correlation and variances and , we exploit the inequality
with some constant , which is based on the Hermite expansion (15) and given in Equation (4) of [7], see also Lemma 3.3 of [13].
By this inequality and Proposition 2 for , we deduce that
With the estimate
we obtain in combination with Corollary 3 that
under the condition , where we use that
(5) 
∎
We turn to the proof of the central limit theorem transferring the strategy from [2] to our model. Define the triangular array
Peligrad and Utev [14, Thm. B] established the central limit theorem , with variance , under the following conditions:

The variances satisfy and there is a constant , such that

The Lindeberg condition is fulfilled:

The following covariance inequality is satisfied. For all , there is a function satisfying , such that for all integers :
Therefore, Theorem 1 follows if the conditions (A) to (C) are verified. (C) is a mixing type condition generalizing the more restrictive condition from [16] that the triangular array is mixing with a certain decay of the mixing coefficients.
Proof of Theorem 1 (A) follows from Proposition 2. More precisely, we can verify analogously to the proofs of the Corollaries 3 and 4 that
and we obtain that
(B) is implied by the Lyapunov condition, since the normal distribution of yields with some constant that
(C) Define . For a decomposition , where is independent of , an elementary estimate with the CauchySchwarz inequality shows that
(6) 
see [2, (52)]. To determine such a suitable decomposition, we write for
(7)  
Then, we set and
where is indeed independent from .
Lemma 5.
Under the conditions of Theorem 1, holds.
3 Remaining proofs
In this section, we write for .
3.1 Proof of Proposition 2
Since , with from (7) and
(8) 
with and centred and independent for , we derive for that
(9)  
Noting that is the density of , we obtain for based on the identity for the convolution that
(10)  
We moreover obtain for and :
Based on that, we determine the terms in (9). Setting
(11) 
we obtain for by the generalization of Itô’s isometry for Walsh integrals
(12)  
Similarly, we have for that
(13)  
For , with , we obtain that
(14)  
Inserting (12), (13) and (14) in (9) yields
For we have and obtain the result in Proposition 2. Since the second derivative of is bounded by for all , we deduce for . Similarly, implies for . With for , we conclude that for :
Comments
There are no comments yet.