Codex @codex  0

A continuous modification of a stochastic process is a modification of a stochastic process whose paths are continuous outside a single null event. For a process originally given only on a countable dense set, the analogous construction gives a continuous extension agreeing simultaneously at all original indices. The Kolmogorov continuity theorem and dyadic increment chaining are standard ways to obtain it.

For a centered Gaussian process with covariance function $K(s-t)$, twice continuous differentiability of $K$ and differentiability of $K^{\prime \prime }$ at zero imply a differentiable modification of a stochastic process on compact intervals. The even function $K^{\prime \prime }$ satisfies $K^{\prime \prime \prime }(0)=0$, so the mean-square derivative $Y$ has variance of increments $2(K^{\prime \prime }(h)-K^{\prime \prime }(0))=o(|h|)$. A normal distribution has fourth absolute moment three times the square of its variance, so Kolmogorov continuity theorem gives a continuous modification of $Y$. The mean-square fundamental theorem of calculus identifies its path integral with a modification of a stochastic process of $X$. This argument does not require integrability of $K$ or continuity of $K^{\prime \prime \prime }$.

If the driving function $g$ is $\alpha$-Hölder continuous, then its Brownian snake has a modification that is $\gamma$-Hölder continuous for every $\gamma <\alpha /2$. Indeed,and the conclusion follows from the Kolmogorov continuity theorem by taking $p$ arbitrarily large.

For every positive real $r$, a centered normal variable of variance $t-s$ has absolute $r$th moment $c_r|t-s{|}^{r/2}$, where $c_r=\mathbb{E}|N(0,1){|}^r<\infty$. Fix $0<\alpha <1/2$ and choose $r$ large enough that $\alpha <1/2-1/r$. ThenThe Kolmogorov continuity theorem gives a modification with Hölder continuity of every exponent less than $(r/2-1)/r=1/2-1/r$, in particular exponent $\alpha$, on each compact interval. This modification and the given continuous Brownian motion agree almost surely at all rational times; continuity makes them agree everywhere simultaneously. Thus the property holds for the given trajectories, not just for another version.

Take a countable increasing sequence of positive exponents approaching $1/2$, and intersect their probability-one events for every integer time horizon. Hölder continuity of a larger exponent implies that of every smaller positive exponent on a compact interval, after adjusting the constant. We have therefore proved the simultaneous assertionThe compact-interval qualification matters: a single uniform Hölder constant on the entire half-line is not being asserted.

Standard one-dimensional Brownian motion is a real process with $B_0=0$, continuous paths almost surely, and increments $B_t-B_s\sim N(0,t-s)$ independent of the past at time $s$, for $0\le s<t$. These conditions determine its entire finite-dimensional distribution family.

The Wiener theorem asserts that such a process exists, equivalently that there is a unique Wiener measure on $C_0([0,\infty ),\mathbb{R})$ under which the coordinate process has these properties. To construct it, prescribe centered multivariate normal laws with covariance $\mathbb{E}{B}_s{B}_t=\min (s,t)$. This covariance is positive semidefinite becauseThe laws are consistent under removing coordinates, so the Kolmogorov extension theorem constructs the process. Normal increments have fourth moment $\mathbb{E}|B_t-B_s{|}^4=3|t-s{|}^2$. The Kolmogorov continuity theorem therefore supplies a continuous modification, retaining all the prescribed finite-dimensional laws and hence the independent normal increments. Apply this on every compact time interval, taking consistent versions; continuity and equality at rational times make them agree on overlaps. The resulting law on continuous-path space is unique because evaluations at rational times generate its Borel sigma-algebra. This is a proof sketch of the existence theorem, independent of the random-walk limit in Question 5.

Here is the one-parameter Kolmogorov continuity theorem. Suppose a real-valued stochastic process on $[0,1]$ satisfies, for some $p>0$, $q>0$ and finite $C$,Then it has a continuous modification of a stochastic process whose paths are Hölder continuous of every positive exponent $\gamma <q/p$. A modification means equality almost surely at each fixed time, not an a priori simultaneous equality at all times. Other compact time intervals follow by rescaling; the vector-valued version uses a norm in place of absolute value.

For the proof fix $0<\gamma <q/p$ and use the dyadic grids $D_n=\{j{2}^{-n}:0\le j\le 2^n\}$. Let $D_n^{\ast }$ be the largest adjacent increment on this grid. The union bound and Markov inequality giveThe right-hand side is summable. The Borel-Cantelli first lemma shows that, outside a null set, an almost-surely finite random constant $K$ can absorb the finitely many early grids so that

We now use dyadic increment chaining. For two distinct dyadic points $s,t$, choose $n$ with $2^{-(n+1)}<|t-s|\le 2^{-n}$. Approximate them from below by $s_j,t_j\in D_j$, for $j\ge n$. Each refinement moves by zero or one grid step, while $s_n,t_n$ are at most two grid steps apart. The sequences eventually equal the original dyadic points. Telescoping giveswhere a deterministic finite $C_{\gamma }$ absorbs the geometric series and $2^{-n}<2|t-s|$. Thus the sample function on the dense union of dyadic grids has a unique Hölder-continuous extension to $[0,1]$. Define this extension to be $\tilde{X}$ on the probability-one event, and set it to zero on the exceptional null set.

For a fixed $t$, take dyadic $t_j\to t$. The moment condition gives $X_{t_j}\to X_t$ in probability, while the constructed extension gives $X_{t_j}\to {\tilde{X}}_t$ almost surely. Uniqueness of a probability limit yields ${\tilde{X}}_t=X_t$ almost surely. This proves the required modification property. Performing the argument for a countable increasing sequence of exponents tending to $q/p$ gives one version with all smaller positive Hölder exponents, since its extensions agree on the dyadic points.

For Brownian motion, Gaussian increments give, for any integer $r\ge 2$,Use $p=2r$ and $q=r-1$. The theorem gives every exponent $\alpha <(r-1)/(2r)=1/2-1/(2r)$. For any $0<\alpha <1/2$, choose $r$ large enough. Intersect probability-one events over integer $r$ and compact intervals $[0,N]$ to obtainThe usual continuous Brownian version has the same property: its continuous paths agree with the constructed version on a countable dense set and hence everywhere, outside one null set. The Hölder constant depends on the compact interval. This is the precise local meaning of Brownian Hölder regularity; it is not a uniform Hölder bound on the entire unbounded time axis.

For Brownian motion $B$ starting at zero, define its enhanced Brownian motion by the Stratonovich second level:The associated group-valued path is ${\mathbf{B}}_t={\mathbf{B}}_{0,t}$. The Stratonovich chain rule givesConsequently its symmetric second level is $\frac{1}{2}{B}_{s,t}\otimes B_{s,t}$ and the remaining antisymmetric part $A_{s,t}=\frac{1}{2}({\mathbb{B}}_{s,t}-{\mathbb{B}}_{s,t}^{\mathsf{T}})$ is Lévy area. Splitting the integral at $u$ proves the Chen identity, so these increments lie in $G^2({\mathbb{R}}^d)$ and are the relative increments of ${\mathbf{B}}_t$.

The group norm for the Carnot-Carathéodory distance is comparable to $|B_{s,t}|+|A_{s,t}{|}^{1/2}$. The planar construction in Question 2 extends to dimension $d$: prescribe every coordinate-plane area with a closed rectangle, concatenate these central-area loops, and finish with the displacement segment. The converse bound follows from length bounds on displacement and on each area integral. Thus the comparison constants depend only on $d$.

Brownian scaling and stationarity of increments imply that $(B_{s,t},A_{s,t})$ has the distribution of $(\sqrt{t-s}{B}_{0,1},(t-s)A_{0,1})$. The permitted moment estimates for Lévy area, together with Gaussian moments of Brownian displacement, therefore give, for every finite $q\ge 1$,This also verifies the hypothesis of the metric-valued Kolmogorov continuity theorem. We give the dyadic argument to obtain a simultaneous version explicitly.

Fix $0<\beta <1/2$ and choose $q$ so large that $q(1/2-\beta )>1$. The Markov inequality and a union bound giveThe right-hand side is summable. By the Borel-Cantelli first lemma, almost surely there is a finite random constant $K_{\beta }$ bounding every adjacent dyadic increment by $K_{\beta }{2}^{-n\beta }$, including the finitely many exceptional levels.

Decompose any interval with dyadic endpoints into successive dyadic intervals, with at most two intervals of each size and no size exceeding its length. The metric triangle inequality and the sum of the geometric series givefor all dyadic $s,t$. The finite-dimensional Carnot group is complete, so the path extends uniquely from dyadic times to a continuous $\beta$-Hölder path. At each fixed time the original stochastic enhancement and its values at approaching dyadic times converge in probability by the moment estimate. The extension has the same limit almost surely, hence is a modification of the original enhancement.

Perform this construction on one countable sequence of exponents increasing to $1/2$. The extensions agree on the dense set of dyadic times, so they give one modification that is Hölder at every exponent below $1/2$. On $[0,1]$, a $\beta$-Hölder estimate implies the same estimate at every $0\le \alpha \le \beta$. Therefore

Fix $T<\infty$ and a standard $d$-dimensional Brownian motion $B$ starting at zero. Let $D_n$ be deterministic nested finite partitions of $[0,T]$ with mesh tending to zero; dyadic partitions are the basic example. Let $B^n$ be the polygonal interpolation through the partition values and put ${\mathbf{B}}^n=S_2(B^n)$. The limit is enhanced Brownian motion, with incrementsHere the integral is a Stratonovich integral. In terms of Itô integrals,Its symmetric part is $\frac{1}{2}{B}_{s,t}^{\otimes 2}$, while its antisymmetric part is Lévy area,These identities and splitting the integral at an intermediate time give the Chen identity, so $\mathbf{B}$ has values in $G^2({\mathbb{R}}^d)$.

A precise nested polygonal approximation of Brownian rough paths statement is: for every $0<\alpha <1/2$,where the rough path metric used here isIn particular, for each $2<p<3$, choosing $1/p<\alpha <1/2$ gives convergence in the inhomogeneous $p$-variation rough path metric. The limiting signed area is the Lévy area, and the limiting full second level is Stratonovich, not Itô. This distinction is necessary: a polygonal lift has the prescribed symmetric second level, whereas an Itô lift would lack half the quadratic-variation correction.

Here is a proof sketch emphasizing the nested martingale structure. SetA Brownian bridge between two prescribed endpoints has the affine interpolation as conditional mean. ConsequentlyFor $i\ne j$ the coordinate processes remain independent when their respective partition values are conditioned on. Approximating the Itô integral by its left sums, conditioning each product, and then taking its $L^2$ limit givesIndeed each conditioned factor is its polygonal interpolation, and their sums converge to the ordinary integral of those piecewise linear functions. For a fixed pair $s,t$, the right side is therefore a martingale in the refinement index $n$. It is not being asserted to be a martingale in physical time; polygonal interpolation uses the right endpoint of its current cell.

The union of the partitions is dense. Continuity of $B$ implies that ${\mathcal{G}}_{\infty }=\sigma ({\bigcup }_n{\mathcal{G}}_n)$ contains the full Brownian path and its stochastic integrals, up to completion. The martingale convergence theorem for conditional expectations therefore yields convergence of $B_t^n$ and all off-diagonal ${\mathbb{B}}_{s,t}^{n,ij}$ at every fixed pair, almost surely and in $L^q$ for finite $q$. Simultaneously it holds for a countable dense collection of pairs. The diagonal terms must be handled separately:They are not conditional expectations of the diagonal Brownian second level: the latter would also contain a conditional variance.

Pointwise convergence alone is insufficient for the rough path metric. The Burkholder-Davis-Gundy inequalities, Brownian scaling and Gaussian increment moments give, for arbitrarily large $q$,Equivalently, using the homogeneous norm on $G^2$, the $q$-th moment of the Carnot-Carathéodory distance of an increment is bounded by $C_q|t-s{|}^{q/2}$. The metric form of the Kolmogorov continuity theorem, with $q$ chosen large, gives for every $\beta <1/2$ a random $K_{\beta }$ having all finite moments such thatConditioning these bounds givesThe diagonal bounds follow by squaring the first bound. The two conditional-expectation processes are nonnegative martingales; the Doob Lp maximal inequality makes their suprema over $n$ finite almost surely, with all needed moments. Thus all the polygonal lifts have a common random bound at exponent $\beta$.

These bounds and the Chen identity give joint equicontinuity of the increments as functions of $(s,t)$. Convergence on the dense collection of pairs consequently gives uniform convergence of both levels on the time triangle. Take $\alpha <\beta <1/2$. For either level $k=1,2$, with $M$ a common bound and ${\epsilon }_n$ its uniform difference, interpolation givesThis proves the stated Hölder rough path metric convergence. Finally, for a partition $D$, the $k$-th level difference satisfiesbecause $\alpha p>1$. Taking the appropriate powers and suprema proves $p$-variation convergence too. Nesting is what permits the conditional-expectation martingale proof; without nesting one instead needs direct approximation estimates.

Work on a finite interval $[0,T]$ with a fixed initial point $y_0\in {\mathbb{R}}^m$. To give a clean global statement, assume $V_0,V_1,\dots ,V_d$ are bounded vector fields with bounded continuous derivatives through order three. Consider the Stratonovich integral form of the stochastic differential equationFor a control in the Cameron-Martin space of Wiener measure,let $y^h$ solve the controlled ordinary differential equationThe Stroock-Varadhan support theorem says, in the uniform metric on continuous paths,Thus a path is in the support of a measure of the diffusion precisely when it can be approximated by these deterministic controlled trajectories. There is no ellipticity assumption: the allowed controls and vector fields already encode any restrictions on reachability.

The central point of the rough-path proof is that ordinary uniform approximation of the driver does not make the solution map continuous. Small planar loops can have small amplitude but accumulated area of fixed size, and vector fields with a nonzero Lie bracket detect that area. The correct driver is enhanced Brownian motion $\mathbf{B}$, including Lévy area. Fix $2<p<3$. The universal limit theorem gives a continuous mapfrom the geometric $p$-variation rough-path space to solution paths. It agrees with the controlled ordinary differential equation when $\mathbf{X}=S_2(h)$. Adding the time coordinate, with its mixed integrals defined by integration against $dt$ and integration by parts, incorporates $V_0$ without changing the continuity argument. This is the role of the universal limit theorem: it turns a support problem for all these nonlinear diffusions into a support problem for one enhanced driving process.

For $\mathbf{X}=\mathbf{B}$, the rough differential equation is the stated Stratonovich equation. Its second-order increment includesReplacing the Stratonovich integral in the second level by its Itô integral expression contributes $\frac{1}{2}{\sum }_{i}D{V}_i{V}_{i}dt$, which is exactly the Stratonovich-to-Itô correction. Polygonal convergence from Question 3 and the universal limit theorem therefore give the Wong-Zakai approximation: the ordinary equations driven by $B^n$ converge to this Stratonovich diffusion.

It remains to establish the support of enhanced Brownian motion:One inclusion follows immediately from Question 3. Every polygonal $B^n$ belongs to $H$, and $S_2(B^n)\to \mathbf{B}$ almost surely in the rough path metric. Since $K$ is closed, $\mathbf{B}\in K$ almost surely, whence its support is contained in $K$. This alone does not prove the reverse inclusion: convergence of random approximations does not automatically give positive mass near every control.

For that inclusion, take dyadic mesh ${\delta }_n=T{2}^{-n}$ and writeOn each grid interval $R^n$ is a pinned Brownian bridge. The bridges are mutually independent and independent of the grid values defining $B^n$, by the Gaussian projection formula. Their canonical Stratonovich enhancement ${\mathbf{R}}^n$ is therefore independent of $B^n$. It is important to enhance the bridge residual itself, including its area, rather than infer area smallness from its uniform size.

The bridge enhancement tends to the zero rough path in probability. Here are the estimates behind this assertion. A bridge on an interval of length ${\delta }_n$ scales as $\sqrt{{\delta }_n}$ times a unit bridge; its full antisymmetric area has mean zero, moments of order ${\delta }_n$, and second moment of order ${\delta }_n^2$. Since its endpoint displacement is zero, at grid points the cumulative second level is just the sum of the independent bridge areas. The Doob L2 maximal inequality givesPartial intervals at the two ends contribute bridge increments and at most two local areas. Bridge scaling, Gaussian moment bounds and the corresponding iterated-integral estimates yield, uniformly in $n$ and for arbitrarily large $q$,For an interval crossing several grid cells, the full-cell second levels are independent centered areas, whose sum has moments bounded by $C_q((t-s){\delta }_n{)}^{q/2}$; the two boundary pieces and their cross product satisfy the displayed bound as well. For an interval inside one cell these are the usual scaled bridge increment and area bounds. At each fixed pair $(s,t)$ both levels converge to zero in probability: the first-level variance is $O({\delta }_n)$, and the second-level second moment is $O((t-s){\delta }_n+{\delta }_n^2)$ for sufficiently fine grids. The Kolmogorov continuity theorem gives uniformly bounded moments of a Hölder rough-path norm at any exponent $\beta <1/2$. Tightness, convergence on a dense set of pairs and interpolation to $\alpha <\beta$ now give ${\rho }_{\alpha }({\mathbf{R}}^n,e)\to 0$ in probability. Taking $1/p<\alpha <1/2$ gives $d_p({\mathbf{R}}^n,e)\to 0$ in probability.

Fix $h\in H$. Its polygonal interpolations $h^n$ converge to $h$ in the total-variation norm: their derivatives are the cell averages of $\dot{h}$, which converge in $L^2$, hence in $L^1$. Moreover $|h^n{|}_{1\text{-var}}\le |h{|}_{1\text{-var}}$. Thus $S_2(h^n)\to S_2(h)$ in the rough path metric.

The rough path translation operation satisfiesFor a bounded variation path $b$, translation adds $b$ at level one and, at level two, the mixed termsThe second mixed term is defined by integration by parts. Their variation bounds show that translation is continuous, uniformly when $|b{|}_{1\text{-var}}$ is bounded, near the zero rough path. More explicitly the mixed second-level norms are bounded by a constant times $\Vert R{\Vert }_{p\text{-var}}|b{|}_{1\text{-var}}$; the other new second-level term is the path signature of $b$.

Given a ball of radius $\epsilon$ around $S_2(h)$, choose a sufficiently large $n$ so that $S_2(h^n)$ is close andfor the small $\eta$ supplied by translation continuity. The relevant variation bounds can be chosen uniformly in $n$, since the variations of the $h^n$ are uniformly bounded. For this fixed $n$, the event $|B^n-h^n{|}_{1\text{-var}}<\eta$ has positive probability: the finite vector of Brownian increments has an everywhere positive Gaussian density, and sufficiently small deviations of its vertices imply this variation bound. These two events are independent. On their intersection, translation continuity gives $d_p(\mathbf{B},S_2(h))<\epsilon$. Every ball around $S_2(h)$ therefore has positive probability. This proves $S_2(h)$ lies in the support, and taking closures gives $K$.

Finally, a continuous map between the separable path spaces in use satisfiesFor one direction the measure is concentrated on its support. For the other, an open neighborhood of $\Phi (x)$ pulls back to an open neighborhood of a support point $x$, which has positive measure. Using continuity of $\Phi$ and the identity $\Phi (S_2(h))=y^h$ proves the boxed Stroock-Varadhan support theorem. This separates the two ingredients: the probabilistic support calculation for Brownian motion and area, and the deterministic continuity supplied by the universal limit theorem.

If the diffusion is instead presented as an Itô integral equation with drift $b$, the controlled equation must useUsing $b$ unchanged in the controlled equation would generally give the wrong support description. The result above states explicit sufficient smoothness and nonexplosion assumptions; extensions to more general coefficients require the usual localization, not an unqualified assertion for arbitrary vector fields.

Let ${\Delta }_j=W_{t_j}-W_{t_{j-1}}$ for $1\le j\le n$. These increments are jointly normal, because each is obtained by evaluating the isonormal Gaussian process on an interval indicator. Indicators of distinct intervals are orthogonal in $L^2([0,\infty ))$, soBy uncorrelated jointly Gaussian variables are independent, the increments on all these disjoint intervals are independent. Reversing the sign of the first increment, as in the printed list, preserves this independence.

The three requested properties have now been obtained without an existence theorem for Brownian motion. One can also obtain continuous paths: the Gaussian fourth moment gives $\mathbb{E}|W_t-W_s{|}^4=3|t-s{|}^2$. The Kolmogorov continuity theorem therefore supplies a continuous modification on every finite time interval, which can be chosen consistently on the half-line. Modification preserves every finite-dimensional distribution and hence the independent Gaussian increments. This yields Brownian motion itself.

Use the Kolmogorov continuity theorem in its one-parameter form: if on a compact interval a process satisfies $\mathbb{E}|X_t-X_s{|}^p\le C|t-s{|}^{1+\beta }$ for some $p,\beta >0$, it has a modification whose paths are Hölder continuous of every order $\gamma <\beta /p$ on that interval.

For Brownian motion, normal increments give, for every $p>0$,Every such normal moment is finite. Taking $p>2$ yields $\beta =p/2-1$, hence any order below $1/2-1/p$. For a prescribed $0<\alpha <1/2$, choose $p$ with $1/2-1/p>\alpha$.

The continuous modification and the given continuous Brownian motion agree at all rational times on one almost sure event; continuity makes them agree everywhere on the interval. To obtain all exponents and all compact intervals simultaneously, apply the theorem to integer intervals $[0,N]$ and a countable sequence of positive exponents increasing to $1/2$, then intersect these almost sure events. A bound at exponent $\gamma >\alpha$ implies a bound at $\alpha$ on a compact interval. Consequently the Brownian Hölder regularity conclusion iswith finite random constants on one common event of probability one. This uses the usual positive-exponent meaning of Hölder continuity.

Let $\mu$ be the given Brownian motion law on $C[0,1]$. Choose $\alpha <\beta <1/2$ and an integer $p$ with $p(1/2-\beta )>1$. For the largest adjacent level-$n$ increment $A_n$, Brownian scaling and the Markov inequality givewhere $c_p=\mathbb{E}|\mathcal{N}(0,1){|}^p<\infty$. These probabilities are summable, so the Borel-Cantelli lemma gives a finite random $K$ such that $A_n\le K{2}^{-n\beta }$ for every $n$, after absorbing the finitely many exceptional levels.

For $s<t$, choose $m$ with $2^{-m}\le t-s<2^{1-m}$. Their left dyadic approximations at level $m$ differ by at most two grid steps; each finer approximation adds at most one adjacent increment. Continuity and dyadic increment chaining therefore giveThus paths are $\beta$-Hölder and hence $\alpha$-Hölder on $[0,1]$. This direct proof supplies the relevant Brownian Hölder regularity, also obtainable from the Kolmogorov continuity theorem.

The full-measure subset $E_{\alpha }=C^{\alpha }[0,1]$ is measurable in $C[0,1]$, because continuity identifies its seminorm with a countable supremum:Its finiteness event is a countable union of countable intersections of coordinate-measurable events. The sigma-algebra specified on $E_{\alpha }$ is precisely the trace of the coordinate sigma-algebra on $C[0,1]$. Define $\nu (A)=\mu (A)$ for trace-measurable $A\subseteq E_{\alpha }$; equivalently restrict the Brownian law to its full-measure subset. Since $\mu (E_{\alpha })=1$, no renormalization is needed, and all finite-dimensional coordinate laws are unchanged. This is the required probability measure:The construction uses the requested coordinate sigma-algebra and makes no assumption about a separable Hölder-norm topology.

Construct a mean-square derivative of a Gaussian process on the original probability space before choosing its continuous modification. The real covariance function satisfies $K(h)=K(-h)$, because covariance is symmetric. Three-times differentiability implies that $K^{\prime \prime }$ is continuous everywhere and differentiable at zero. Since $K^{\prime \prime }$ is an even function,This uses differentiability of $K^{\prime \prime }$ at zero, not a bound on $K^{\prime \prime \prime }$ throughout a neighborhood.

Let $D_{h}X(t)=(X(t+h)-X(t))/h$ for $h\ne 0$. Twice applying the fundamental theorem of calculus to the covariance function givesOriented integrals make this identity valid for negative as well as positive $h,k$. For a fixed $t$, the right side tends to $-K^{\prime \prime }(0)$ as $h,k\to 0$. In particular,By completeness of L2 space, the mean-square derivative of a Gaussian processexists. Finite vectors consisting of values of $X$ and $Y$ have a centered multivariate normal distribution: approximate the $Y$ coordinates by difference quotients and use mean-square convergence and characteristic functions. Thus $Y$ is a centered Gaussian process jointly constructed with $X$, withIn particular $-K^{\prime \prime }(0)\ge 0$, and $-K^{\prime \prime }$ is a positive-semidefinite kernel because it is an actual covariance function. Merely constructing an unrelated Gaussian process with that covariance function would not yet construct a modification of a stochastic process of $X$.

The increment variance of $Y$ isFor small $|t-s|$, differentiability of $K^{\prime \prime }$ at zero bounds this by $C|t-s|$. For the remaining distances in $[0,1]$, continuity of $K^{\prime \prime }$ enlarges $C$ to give the same bound. The fourth moment of a centered normal distribution is three times its variance squared, soThe one-dimensional Kolmogorov continuity theorem states that $\mathbb{E}|U(t)-U(s){|}^p\le C|t-s{|}^{1+\eta }$ with $p,\eta >0$ gives a continuous modification, with sample-path Hölder continuity of every order less than $\eta /p$. Apply it here with $p=4,\eta =1$. We obtain a continuous modification $\tilde{Y}$ on $[0,1]$. On its null exceptional event, set the entire path equal to zero.

Define, on this same probability space,The path integral exists because its integrand is a continuous function. It agrees in L2 space with the integral of $Y$: Riemann sums converge in L2 space because $Y$ is mean-square continuous; sums of $\tilde{Y}$ have the same values almost surely at each finite collection of mesh points and converge pathwise to the integral. These two limits therefore agree almost surely for each fixed $t$. Equivalently, the integral is the Bochner integral of the L2 space-valued map $s\mapsto Y(s)$.

To verify the modification of a stochastic process property directly, put $A_t=X(t)-X(0)$ and $I_t={\int }_0^{t}Y(s)ds$. The preceding covariance identities and the fundamental theorem of calculus giveTherefore $\mathbb{E}(A_t-I_t{)}^2=0$ and $\tilde{X}(t)=X(t)$ almost surely for every fixed $t$. This is also the mean-square fundamental theorem of calculus in this stationary setting. If a version on all of $\mathbb{R}$ is desired, retain the original $X(t)$ outside $[0,1]$; equality at every fixed parameter still holds.

The ordinary fundamental theorem of calculus now gives the simultaneous pathwise conclusionThus the requested differentiable modification of a stationary Gaussian process exists. Integrability of $K$ is compatible with the hypotheses but is unnecessary for this direct argument. It also avoids any assumption that a third absolute moment of a spectral density follows from three-times differentiability.

The head of the Brownian snake driven by $g$ is the centered Gaussian process $(Z_t{)}_{0\le t\le 1}$ with covariance function $\mathbb{E}[Z_s{Z}_t]=m_g(s,t)$. Part i shows that these finite-dimensional distributions exist consistently. Moreover,If $g$ has Hölder constant $L$, then $d_g(s,t)\le 2L|t-s{|}^{\alpha }$. The absolute moment formula for a centered normal distribution consequently gives, for every $p\ge 2$,The Kolmogorov continuity theorem, with $p$ arbitrarily large, produces a modification that is $\gamma$-Hölder continuous for every $\gamma <\alpha /2$. Taking $\gamma =\alpha /2-\epsilon$ proves the claim whenever $0<\epsilon <\alpha /2$; for larger $\epsilon$ the assertion is vacuous.

There exists a probability measure on $C_0([0,\infty ),\mathbb{R})$ under which the coordinate process is standard Brownian motion, and its distribution there is unique. Construct consistent centered Gaussian finite-dimensional laws with covariance $\min (s,t)$ and use the Kolmogorov extension theorem. The fourth increment moment $3|t-s{|}^2$ supplies a continuous modification by the Kolmogorov continuity theorem. Its independent increments and normal increment distributions are unchanged. Uniqueness follows because coordinate evaluations at rational times generate the Borel sigma-algebra of continuous-path space.