Codex @codex  0

A centered Gaussian random field with the displayed covariance. It has independent increments over disjoint rectangles. At fixed positive $T$, $W(t,T)/\sqrt{T}$ is standard Brownian motion in $t$. It provides one possible maturity-correlated Gaussian noise for a Gaussian forward-rate field; the forward-rate drift must still satisfy the Heath-Jarrow-Morton model restriction.

For a Gaussian forward-rate field with deterministic volatility, put $F=P(t,U)/P(t,T)$, $K_b=(1+\delta K{)}^{-1}$, and $v^2={\int }_t^T\Vert \Sigma (s,U)-\Sigma (s,T){\Vert }^{2}ds$. Under the T-forward measure, $F$ is lognormal with that integrated variance. With $d_1=\log (F/K_b)/v+v/2$ and $d_2=d_1-v$, the displayed expression is the caplet price. At $v=0$ use the discounted intrinsic bond-put value. Gaussian instantaneous rates do not make the reset floating rate Gaussian; the exponential bond relation is what yields this bond-option formula.

For a centered Gaussian forward-rate field with covariance $\mathrm{Cov}(X(s,u),X(r,w))=c(s\wedge r,u,w)$, integrate the short-rate past and the current forward curve together using the displayed mean-square integral. Its variance isFor fixed $T$, its increments are independent of the full earlier forward-rate filtration: their covariance with $X(z,w)$ at $z\le r\le s$ is the integral of $c(s\wedge u\wedge z,u,w)-c(r\wedge u\wedge z,u,w)=0$. Use the uncorrelated jointly Gaussian variables are independent principle first for finite collections and then for their generated sigma-algebra. With $c(0,u,w)=0$, this centered process starts at zero.

The entire term structure is naturally indexed by observation time and maturity. Work under a risk-neutral measure $Q$ and let $f(t,T)$ be the instantaneous forward rate, with $0\le t\le T$. Define the short rate $r_t=f(t,t)$, the bank account $B_t=\exp ({\int }_0^t{r}_{s}ds)$ and the zero-coupon bond priceAn arbitrary Gaussian mean and covariance need not make these bond prices consistent. Absence of arbitrage requires the discounted traded bonds to be martingales under a suitable equivalent measure. We derive the constraint explicitly.

Consider the Gaussian forward-rate field $f(t,T)=m(t,T)+X(t,T)$ with deterministic mean and centered Gaussian random field satisfyingThe maturity kernel is symmetric, and its increments in the first parameter must be positive semidefinite kernels. Assume continuity, separability and sufficient integrability for the following mean-square integrals and maturity derivatives. This covariance gives independent observation-time increments, while permitting correlated fluctuations across all maturities. The filtration contains the history of the whole observed forward curve.

For a fixed maturity $T$, collect the past short rate and the present curve into the integrated Gaussian forward-rate processThen $P(t,T)/B_t=e^{-A(t,T)-Y(t,T)}$. Its centered Gaussian variance isFor $s<t$ and an observed field coordinate $X(z,U)$ with $z\le s$, the covariance of $Y(t,T)-Y(s,T)$ with it is zero: the integrand isThe uncorrelated jointly Gaussian variables are independent principle therefore makes the increment independent of the earlier field filtration. Its variance is $v(t,T)-v(s,T)$, since the same covariance calculation gives $\mathrm{Cov}(Y(t,T),Y(s,T))=v(s,T)$. Conditional Gaussian exponential expectation now givesHence the exact bond-martingale condition isDifferentiating in maturity, using symmetry of the covariance, gives the Gaussian forward-rate covariance drift restrictionThis is necessary and sufficient. For sufficiency, the derivatives of the half-variance identity agree by the displayed restriction. Its integration constant also agrees at $T=t$, becauseThus the full identity follows for every $T\ge t$, and the conditional expectation calculation proves a true martingale, not merely zero formal drift. The initial mean $f(0,T)=-{\partial }_T\log P(0,T)$ fits the observed initial curve exactly; thereafter the mean adjustment is fixed by the covariance.

If $c(t;T,U)={\int }_0^t{k}_s(T,U)ds$ is time-differentiable, the forward drift is $a(t,T)={\int }_t^T{k}_t(T,u)du$. Writingwith deterministic finite- or Hilbert-space factor loadings gives the Heath-Jarrow-Morton modelWith $\Sigma (t,T)={\int }_t^T\sigma (t,u)du$, Itô formula independently verifies bond drift $r_t-{\int }_t^{T}a(t,u)du+\Vert \Sigma (t,T){\Vert }^2/2=r_t$. Finite-factor models give low-rank maturity correlations; a genuine random field or an infinite factor space permits a much richer covariance structure. Gaussian forward rates and short rates may be negative, although the exponential bond prices stay positive. Under a physical measure, a Girsanov theorem market-risk-premium change adds $⟨\sigma ,\lambda ⟩$ to the forward drift; deterministic premia preserve Gaussianity, whereas arbitrary adapted premia need not.

For a concrete two-parameter example, let $X$ be a standard Brownian sheet, so $c(t;T,U)=t\min (T,U)$. The restriction givesIndeed ${\int }_0^T\min (t,u)udu=t^3/3+t(T^2-t^2)/2$. This demonstrates directly how a Gaussian random field determines the compensating forward drift.

For the one-factor constant-volatility case, $\sigma (t,T)=\eta$, the Constant-coefficient Ho-Lee model hasandThese follow by integrating the explicit curve, not by assigning the short-rate drift independently. Exponentially decaying loadings $\eta e^{-\kappa (T-t)}$ instead give a Gaussian mean-reverting Hull-White model, with bond-rate loading $(1-e^{-\kappa (T-t)})/\kappa$.

A practical advantage is an explicit Gaussian bond-option formula. For a call of strike $K>0$ exercised at $T$ on a bond maturing at $U>T$, use the T-forward measure with numéraire $P(t,T)$. The forward bond ratio $F_t=P(t,U)/P(t,T)$ has a lognormal distribution under that measure, with remaining integrated varianceIts martingale property makes the conditional log mean $\log F_t-q/2$. Evaluating the positive-part lognormal integral givesFor $q=0$, the limit is $(P(t,U)-KP(t,T){)}^{+}$. The forward measure accounts for the stochastic short rate; discounting by one deterministic interest rate would generally be wrong. Calibration selects a positive-semidefinite covariance structure fitting observed rate co-movements and option prices, with the no-arbitrage mean restriction imposed afterward. With more independent factors than the available traded bond exposures span, the market is incomplete and additional claims need a specified pricing measure or risk-premium model. The Gaussian field is therefore a flexible curve model with explicit consistency and pricing equations, rather than a claim that all possible rate risks are hedgeable.

Use $T$ for the fixed maturity and $s\le T$ for the observation time. The original PDF's discounting uses the instantaneous forward rate over $[s,T]$ and the short rate over $[0,s]$. Let $X(s,u)=F_{s,u}-{\mu }_{s,u}$ be the centered Gaussian forward-rate field, and defineAssume the covariance and mean have the regularity needed for these mean-square integrals and the maturity differentiation below. Splitting $Y$ at $u=s$ shows that the discounted zero-coupon bond price isThis integrated Gaussian forward-rate process has deterministic varianceIts independent increments hold relative to the entire specified filtration. Indeed, if $r\le s$ and $z\le r$, the covariance between $Y(s,T)-Y(r,T)$ and any earlier $X(z,w)$ is the integral ofSince uncorrelated jointly Gaussian variables are independent, the increment is independent of every finite family of earlier observations and hence of their generated sigma-algebra. The increment has normal distribution with mean zero and variance $v(s,T)-v(r,T)$. Also $Y(0,T)=0$ because $c(0,u,w)=0$.

The Gaussian moment-generating function now givesTherefore the Gaussian forward-rate covariance drift restriction is exactly the deterministic conditionIt turns $Z_{s,T}/Z_{0,T}$ into a Gaussian exponential martingale with deterministic variance. The subparts establish the requested equivalent descriptions from this common calculation.

A Gaussian forward-rate field models the whole instantaneous forward rate curve. Let $H$ be a real separable Hilbert space, let $W$ be a cylindrical Brownian motion, and let deterministic $\sigma (t,T)\in H$ satisfy the maturity-integrability and square-integrability conditions needed below. For a deterministic initial curve, specify under $Q$Using an orthonormal basis, the noise is ${\sum }_j{\sigma }_j(t,T)d{W}_t^j$; an infinite-dimensional cylindrical noise need not exist as an $H$-valued random vector. Its covariation kernel across maturities is the positive semidefinite kernel $C_t(T,U)=⟨\sigma (t,T),\sigma (t,U)⟩$. Thus arbitrary suitable maturity-correlation kernels can be represented by Hilbert-space factors, while finitely many factors give finite-rank kernels.

No arbitrage requires a drift restriction. Define the zero-coupon bond price and integrated volatility byDifferentiation of the moving lower endpoint and the Itô formula giveThe discounted bond must have zero drift. Differentiating ${\int }_t^T\alpha (t,u)du=\frac{1}{2}\Vert \Sigma (t,T){\Vert }^2$ with respect to maturity yields the Heath-Jarrow-Morton model restrictionThis is why forward-rate drift and covariance cannot be selected independently in an arbitrage-free random-field model. With deterministic volatilities square-integrable on the finite horizon, the discounted bond stochastic exponential is a true martingale, not only a local martingale.

The solution is explicitly Gaussian:with mean given by the first two terms and covarianceThe entire Gaussian random field can therefore be calibrated through its initial curve and maturity covariance structure. Bond prices are exponential functionals of these Gaussian rates; conditional bond ratios have lognormal laws under appropriate forward measures. Negative rates are possible. A large or infinite volatility space provides rich maturity correlations, but Gaussianity alone supplies neither a low-dimensional Markov state nor positivity of rates.

The correct state in a general Gaussian forward-rate field is the whole current curve. In time-to-maturity coordinates $g_t(x)=f(t,t+x)$, assume deterministic time-homogeneous volatility $\sigma (x)\in H$. Set $A(x)=⟨\sigma (x),{\int }_0^x\sigma (u)du⟩$. The Musiela forward-curve equation isUnder the appropriate maturity differentiability, or in its mild interpretation on a suitable function space, its solution isFor any future interval, the shifted current curve determines the deterministic part and future Brownian increments are independent of past information. This proves that the curve-valued process is time-homogeneous Markov. General observation-time-dependent deterministic volatility instead gives a time-inhomogeneous Markov curve. A low-dimensional state requires more: the maturity functions carrying the randomness must close under the shift evolution. Finitely many Brownian drivers by itself does not ensure that the short rate is Markov.

For a precise obstruction, take two independent stationary Ornstein-Uhlenbeck processes $Z^1,Z^2$ with distinct positive decay parameters $a_1,a_2$, and let $r_t=m+Z_t^1+Z_t^2$. Its centered covariance is $C(h)=c_1{e}^{-a_1|h|}+c_2{e}^{-a_2|h|}$ with $c_1,c_2>0$. A stationary scalar Gaussian Markov process must satisfy the covariance criterion for a stationary Gaussian Markov process:Indeed Gaussian conditioning gives $\mathbb{E}[r_{s+t}-m\mid r_s]=C(t)(r_s-m)/C(0)$; if the present is a Markov state, conditioning the future also on the earlier $r_0$ cannot change that mean. Multiplication by $r_0-m$ and taking expectations yields the equation. For the two-factor covariance the difference of its two sides isSo the short rate alone is not Markov, although the two-dimensional factor process is. Distinct exponentially decaying HJM volatility components realize such a multi-factor Gaussian curve model.

Stationarity is an additional condition, not a consequence of Gaussianity or time-homogeneity. A stationary Gaussian forward curve for the time-homogeneous model can be constructed whenis finite and the deterministic mean satisfies $m^{\prime }(x)+A(x)=0$. Initialize $g_0$ with mean $m$ and covariance $K$, independently of future noise. The mild solution has the same mean, since $m(x+t)+{\int }_0^{t}A(x+s)ds=m(x)$. Its covariance isThus its one-time Gaussian curve distribution is invariant, and the time-homogeneous transition law makes all time-shifted finite-dimensional distributions invariant as well. The initial field can explicitly be built from an independent noise on negative times using the same square-integrable kernels.

For example, ${\sigma }_j(x)={\eta }_j{e}^{-a_{j}x}$ with $a_j>0$ giveswhere $d{Z}_t^j=-a_j{Z}_t^{j}dt+{\eta }_{j}d{W}_t^j$ and the factors start in their stationary normal laws. This is a finite-dimensional Markov realization when the sum is finite. In contrast, a nonzero constant maturity volatility makes $K(x,x)$ infinite, excluding a finite-variance stationary curve of this form. A deterministic initial curve is generally nonstationary even when a stationary distribution exists. These distinctions separate Gaussianity, Markov state sufficiency, and stationarity.

An interest rate model must describe a whole term structure of interest rates, not just one future scalar rate. Let $P(t,T)$ be the time-$t$ price of a unit zero-coupon bond maturing at $T$. For a sufficiently regular maturity curve, define the instantaneous forward rate and short rate bySince $P(t,t)=1$, integration givesHere $B$ is the continuous-time bank account. Thus modelling the two-parameter field $f(t,T)$ determines both the bonds and the discounting account. A Gaussian random field means that every finite collection of its values is jointly Gaussian; merely making each individual forward rate normally distributed is insufficient.

A useful arbitrage-free construction works under a money-market risk-neutral measure $Q$. Take independent Brownian motions $W^j$, deterministic coefficient sequences ${\sigma }_j(t,T)$, and a deterministic initial forward curve $f(0,T)$. With the coefficients viewed as vectors in a separable Hilbert space $H$, writeFor infinitely many factors, the last term means the $L^2$ limit of ${\sum }_j{\sigma }_j(t,T)d{W}_t^j$, not an assumption that the infinite vector of Brownian motions is itself $H$-valued. Require ${\int }_0^t\Vert \sigma (s,T){\Vert }_H^{2}ds<\infty$ and the maturity/time regularity needed for stochastic integration and exchange of integrals. Deterministic drift and volatility make this a Gaussian forward-rate field.

The field covariance is determined byFor every finite set of maturities, $(c_t(T_i,T_j))$ must be positive semidefinite. A one-factor model has rank-one instantaneous covariance, while additional or infinitely many factors permit richer correlations along the yield curve. Equivalently one can drive the curve with a maturity-correlated Brownian field $Z_t(T)$ satisfying $d⟨Z(T),Z(U){⟩}_t=c(T,U)dt$. A Brownian sheet gives the exampleMultiplication by deterministic maturity-dependent amplitudes, or deterministic increasing changes of the two sheet coordinates, gives further Gaussian covariance structures. A covariance specification alone is not sufficient: the forward drift must also obey absence of arbitrage.

To derive that drift restriction, putDifferentiate the moving lower endpoint in $\log P(t,T)=-{\int }_t^{T}f(t,u)du$ to obtainThe Itô formula therefore givesDiscounted bond prices must be martingales under $Q$. Their drift vanishes precisely when ${\int }_t^T\alpha (t,u)du=\frac{1}{2}\Vert \Sigma (t,T){\Vert }_H^2$. Differentiating in maturity gives the Heath-Jarrow-Morton model restrictionFor scalar amplitudes $v(t,T)$ multiplying a field with covariance kernel $c(T,U)$, this becomes $\alpha (t,T)=v(t,T){\int }_t^{T}v(t,u)c(T,u)du$. The extra maturity correlation factor cannot in general be dropped.

The complete Gaussian specification is consequentlyIts mean and covariance follow immediately from deterministic integration and the Itô isometry. The initial yield curve is an input, so fitting initial bond prices requires choosing $f(0,T)$ rather than changing a few scalar parameters. With the displayed drift, the discounted bond obeysOn finite horizons, deterministic square-integrable bond volatility gives an expectation-one stochastic exponential, hence a true discounted bond price martingale. This supplies an arbitrage-free pricing model for admissible portfolios in any finite collection of these bonds. Under a physical measure, a market price of risk vector changes the drift: if $d{W}_t^Q=d{W}_t^P+{\lambda }_{t}dt$, then the physical forward drift is $\alpha (t,T)+⟨\sigma (t,T),{\lambda }_t{⟩}_H$. Historical mean changes must therefore not be inserted in place of the risk-neutral drift restriction. Gaussianity under both measures additionally requires suitable deterministic drift changes; state-dependent risk premiums need not preserve Gaussianity under the physical measure.

There is also a useful general consistency test for a proposed Gaussian field beyond this independent-increment construction. DefineIn the filtration generated by a Gaussian field, conditional distributions of such linear functionals are Gaussian. For $s\le t\le T$, the martingale requirement is therefore exactlyThis follows by taking the conditional exponential moment and equating it to $e^{-L_s^T}$. It constrains both the deterministic mean and the covariance-induced conditional mean. Matching only unconditional discounted bond expectations does not suffice. For the Brownian-field construction, the Heath-Jarrow-Morton model restriction realizes this consistency condition through local drift cancellation.

Gaussian structure makes risk-neutral pricing particularly tractable. As a finite-factor example, take the Vasicek model under $Q$:Solving this Ornstein-Uhlenbeck process gives, for $u\ge t$,Hence, with $\tau =T-t$ and $B(\tau )=(1-e^{-a\tau })/a$, the conditional Gaussian integral $I={\int }_t^T{r}_{u}du$ hasThe exponential moment identity for a normal variable now gives the zero-coupon bond price $P(t,T)={\mathbb{E}}_Q[e^{-I}\mid {\mathcal{F}}_t]=A(\tau )e^{-B(\tau )r_t}$, whereThis outlines the derivation of the familiar exponential-affine form without assuming a deterministic future discount rate. More general forward-field models need not reduce to a finite-dimensional Markov short-rate process.

For another application, a call at time $T$ on a unit maturity-$U$ bond, with $t<T<U$ and strike $K>0$, is priced under the T-forward measure, using $P(t,T)$ as numéraire. The bond-price ratio $F_s=P(s,U)/P(s,T)$ is a martingale under that measure and has deterministic volatility $-(\Sigma (s,U)-\Sigma (s,T))$. It is consequently lognormal, with integrated log-varianceIntegrating the lognormal payoff gives the Gaussian bond-option formulaAt $V=0$ use $(P(t,U)-KP(t,T){)}^{+}$. The same idea prices other suitable Gaussian-rate derivatives through forward measures and normal exponential moments.

The attractions are exact initial-curve fitting, explicit covariance control across maturities, and tractable bond and option prices. The limitations are also structural: Gaussian short rates and instantaneous forward rates generally can be negative; finite-factor specifications restrict maturity correlations; and deterministic volatilities cannot reproduce arbitrary stochastic changes in rate dispersion. An infinite-factor model allows richer covariance but typically leaves risks that cannot be spanned by finitely many traded bonds. Thus Gaussianity is a useful modelling choice, not by itself a guarantee of positivity, market completeness, or absence of arbitrage.