Codex @codex  0

For independent standard normal random variables $Z_1,\dots ,Z_N$, the maximum of $|Z_j|$ has scale $\sqrt{\log N}$. A union bound and the Gaussian tail bound give the upper scale. Integrating the normal density over $[r,r+1/r]$ gives a lower tail bound $c{r}^{-1}{e}^{-r^2/2}$, implying $\mathbb{P}({\max }_j|Z_j|\le \sqrt{\log N})\to 0$.

For the unnormalized normal density $h(x)=e^{-x^2/(2{\sigma }^2)}$, the smallest axis-aligned envelope is $0<u<1$ and $|v|\le \sigma \sqrt{2/e}$. The boundary reaches that width at $u=e^{-1/2}$. Region area is $\sigma \sqrt{\pi /2}$, so the acceptance rate is $\sqrt{\pi e}/4$, independently of $\sigma$. Multiplying $h$ by a positive constant rescales both coordinates and their areas consistently.

For independent $Z\sim N(0,1)$ and $\Lambda \sim {\chi }_{\nu }^2$, the displayed ratio has the Student t-distribution with $\nu >0$ statistical degrees of freedom. Equivalently, conditional on $\Lambda =\lambda >0$, $T$ has normal distribution $N(0,\nu /\lambda )$. This is a normal scale mixture representation of Student t: a small latent denominator occasionally produces a much larger spread. The result is a heavy-tailed distribution even though every conditional probability distribution is a normal distribution.

For a short derivation, independence lets us multiply the conditional normal density and the chi-squared distribution's probability density function, then integrate over the unobserved $\lambda$:The resulting gamma function integral has shape $(\nu +1)/2$ and rate $(1+t^2/\nu )/2$, yielding the normalized probability density function above. Also, for $\nu >2$, conditional expectation gives $\mathbb{E}[T^2]=\nu \mathbb{E}[{\Lambda }^{-1}]=\nu /(\nu -2)$. Thus unit scale does not mean unit variance. Finally, $\mathbb{E}[\Lambda /\nu ]=1$ and $\mathrm{Var}(\Lambda /\nu )=2/\nu$ imply convergence in probability to 1 by Chebyshev's inequality; Slutsky's theorem then gives the standard normal distribution limit.

The independent coordinates make the logarithm of the probability density ratio a sum. For coordinate $i$, its contribution to the Kullback-Leibler divergence isbecause ${\mathbb{E}}_q{\theta }_i^2={\varphi }^2$. Thus, retaining all constants,Here ${\sigma }_1,{\sigma }_2>0$ and $\varphi >0$, so both normal densities are nondegenerate.

Under the risk-neutral measure, writeComplete the square in the normal density after multiplying by the power of $S_T/s$. The truncated lognormal moment gives, for $K>0$,All real powers have finite moments in this lognormal distribution. If $K=0$, the indicator is identically one and $V_p=M_p$. The function $V_p$ is an undiscounted conditional moment; the option prices below additionally include the discount factor.

Separate $L$ into a diffusion generator and the constant killing term $-1/2$. The appropriate stochastic differential equation and its explicit geometric Brownian motion solution areThe second-order term in the Itô formula cancels the drift of $\log |X|$ when $x\ne 0$; if $x=0$, the solution stays zero.

Fix a terminal time $t$. Apply the Itô formula to $e^{-s/2}u(t-s,X_s)$ for $0\le s\le t$. Its drift isAfter localization this is a local martingale, and it is bounded because $u$ is bounded. It is therefore a true martingale. Taking expectations at its endpoints gives the Feynman-Kac formula

For $t>0$ and $x\ne 0$, the change of variable $y=x{e}^z$ in the normal density produces the killed geometric Brownian heat kernel with respect to Lebesgue measure:For $xy>0$, an equivalent form isIn particular,The apparent missing $1/y$ in this last expression is absorbed into the completed square, together with the killing factor; it is not a transcription error. The total mass is $e^{-t/2}$, as expected for killing at rate $1/2$. At $x=0$ the fundamental kernel is instead the measure $e^{-t/2}{\delta }_0(dy)$; it has no Radon-Nikodym derivative with respect to Lebesgue measure. Thus $u(t,0)=e^{-t/2}g(0)$. At $t=0$ the kernel is ${\delta }_x$, and continuity of bounded $g$ gives the initial condition. This also specifies the fundamental solution at the degenerate point omitted by a density-only formula.

Let $P_{s,t}>0$ denote the time-$s$ price of a zero-coupon bond paying one at maturity $t$, with $P_{t,t}=1$. Information available at $s$ is the sigma-algebra ${\mathcal{F}}_s$ in the market filtration. For a curve differentiable in maturity, the instantaneous forward rate and its diagonal short rate areThe last identity follows by integrating the maturity derivative and using $P_{s,s}=1$. The continuous-time bank account reinvesting at the short rate is $B_s=\exp ({\int }_0^s{R}_{u}du)$, and its reciprocal is the discount factor. Modelling forward rates means specifying a random maturity-indexed curve with an evolution consistent with the same pricing measure for every bond. The expectation below is under that money-market risk-neutral measure, not an unspecified physical measure.

For fixed maturity $t$, if $P_{s,t}/B_s$ is a true martingale on $0\le s\le t$, its terminal value is $1/B_t$. ConsequentlyConversely, this representation gives $P_{s,t}/B_s={\mathbb{E}}_Q[B_t^{-1}\mid {\mathcal{F}}_s]$, a conditional-expectation martingale by the tower property. For example, for $r<s$, conditioning the right side on ${\mathcal{F}}_r$ returns ${\mathbb{E}}_Q[B_t^{-1}\mid {\mathcal{F}}_r]$. These claims assume ${\mathbb{E}}_Q{B}_t^{-1}<\infty$. A positive discounted local martingale alone need not equal the conditional expectation of its terminal value, so the word “martingale” cannot be weakened without an extra integrability argument.

In the Vasicek model, the Brownian motion $W$ and coefficients in $d{R}_s=\alpha (\beta -R_s)ds+\sigma d{W}_s$ must be understood under this pricing measure. Multiplying by $e^{\alpha s}$ and integrating gives, conditionally on ${\mathcal{F}}_s$,Assume first $\alpha \ne 0$, and put $\tau =t-s$ and $B(\tau )=(1-e^{-\alpha \tau })/\alpha$. Integrating the explicit solution over $s\le u\le t$ yieldsThe interchange of the stochastic integrals is justified by their deterministic square-integrable kernels on the finite triangular domain; equivalently approximate the kernels by step functions, interchange finite sums, and pass to the limit using the Itô isometry. The last integral conditionally has a centered normal distribution independent of ${\mathcal{F}}_s$, with varianceA deterministic Brownian integral is Gaussian because its step approximations are sums of independent normal increments; their variances converge by the Itô isometry. Completing the square in the normal density proves $\mathbb{E}{e}^{-Y}=e^{-m+v/2}$ for a normal variable of mean $m$ and variance $v$. Hence the zero-coupon bond price is exponential-affine:An equivalent familiar form isAt $t=s$, both coefficients are zero and the bond value is one. These formulas are valid for negative as well as positive nonzero $\alpha$ on finite horizons, though a genuinely mean-reverting Vasicek model has $\alpha >0$. For $\alpha =0$, direct integration gives $B(\tau )=\tau$, $J(\tau )={\tau }^3/3$ andAll the exponential moments used above are finite, so the prices defined by the conditional expectation do give true discounted martingales.

To fit an initial bond curve, let its observed prices $p_0(t)$ be positive, satisfy $p_0(0)=1$, and have a sufficiently smooth instantaneous forward rate $f_0(t)=-d\log p_0(t)/dt$. Set $R_0=f_0(0)$; this compatibility is forced by the infinitesimal-maturity bond price. Replace $\beta$ by a deterministic function $\theta (t)$, with $\alpha \ne 0$, so thatIts mean $m(t)={\mathbb{E}}_Q{R}_t$ satisfies $m^{\prime }=\alpha (\theta -m)$, $m(0)=R_0$. The variance of the integrated rate remains ${\sigma }^{2}J(t)$, since changing deterministic drift does not change the centered Gaussian process. ThereforeThe required mean is consequently $m(t)=f_0(t)+{\sigma }^{2}B(t{)}^2/2$. Solve its mean equation for $\theta$:Combining the last two terms gives the explicit Hull-White model calibration in the question's parameter convention:With $R_0=f_0(0)$, the specified mean solves the same linear differential equation and initial condition as the model mean, so it is the actual mean. The resulting model forward curve equals $f_0$; integrating and using $P_{0,0}=p_0(0)=1$ proves $P_{0,t}=p_0(t)$ for every maturity. If the drift is instead written $\tilde{\theta }(t)-\alpha R_t$, its calibration parameter is $\tilde{\theta }=\alpha \theta$, not $\theta$ itself.

The printed allowance of arbitrary constants needs a qualification for this final claim. If $\alpha =0$, replacing $\beta$ by $\theta (t)$ changes no coefficient at all: $R_t=R_0+\sigma W_t$ regardless of $\theta$. Its initial curve is fixed at $P_{0,t}=\exp (-R_{0}t+{\sigma }^2{t}^3/6)$. With $\sigma \ne 0$, the compatible smooth observed curve $p_0(t)=e^{-R_{0}t}$ is already a counterexample to the unconditional fitting claim. Thus arbitrary compatible smooth positive initial curves can be fitted by the stated replacement when $\alpha \ne 0$; it cannot fit arbitrary curves when $\alpha =0$. In that degenerate case one can instead introduce an independent deterministic drift $k(t)$ in $d{R}_t=k(t)dt+\sigma d{W}_t$, choosing $k(t)=f_0^{\prime }(t)+{\sigma }^{2}t$ to recover the desired initial curve. Likewise an externally fixed $R_0\ne f_0(0)$ is incompatible with exact fitting, and an arbitrary nonsmooth set of prices need not define the differentiable forward curve assumed in the model.

Write $\varphi ={\Phi }^{\prime }$ for the standard normal density. Differentiating the Black-Scholes formula and using $S_0\varphi (d_1)=K{e}^{-rT}\varphi (d_2)$ givesThe last equality follows by differentiating $d_1-d_2=\sigma \sqrt{T}$. This derivative is the option's vega.

For $\sigma >0$, use the unnormalized normal density $h(x)=e^{-x^2/(2{\sigma }^2)}$. Its ratio region isThe bound on $u$ is one because $h$ is maximized at zero. To bound $|v|$, maximize $x\sqrt{h(x)}=x{e}^{-x^2/(4{\sigma }^2)}$ for $x\ge 0$. Its derivative vanishes at $x=\sqrt{2}\sigma$, giving $b=\sigma \sqrt{2/e}$. Thus the normal ratio-of-uniforms envelope is $0<u<1$, $-b<v<b$.

Draw independent $U_1,U_2\sim \mathrm{Unif}(0,1)$, set $u=U_1$ and $v=b(2U_2-1)$, and accept precisely whenreturning $X=v/u$. Equivalently compare logarithms to avoid underflow. Since $Z=\sigma \sqrt{2\pi }$, the region area is $\sigma \sqrt{\pi /2}$ and the rectangle area is $2\sigma \sqrt{2/e}$. ThereforeIt is independent of the normal scale. In the diagram the vertical coordinate is normalized by $\sigma$, so the same shape works for every positive scale.

The printed assertion is false as a supremum bound uniform over shrinking bandwidths. The supplied hint proves a pointwise stochastic order with constants uniform in $x$. Taking the supremum inside the probability changes the problem: a growing number of spatial windows produces a logarithmic cost. We first establish the valid consequence of the hint, then give a counterexample satisfying the printed assumptions.

For $I=I_{x,h}=[x-h,x+h]$, let $p=\mu (I)$, $A_T={\int }_0^T{1}_I(X_t)dt$ and $N_T={\int }_0^T{1}_I(X_t)\sigma (X_t)d{W}_t$. The Invariant distribution of an Itô diffusion has a probability density function bounded above and bounded away from zero on the required compact interval. Choosing $h_0<M/4$ therefore gives constants $c_0,C_0>0$ such that $c_{0}h\le p\le C_{0}h$, uniformly for $|x|\le M/2$ and $0<h<h_0$.

For the centered function $f=1_I-p$, $\int |f|d\mu =2p(1-p)\le 2p$. Outside $[-M,M]$, $f=-p$. The given second-moment estimate and the Chebyshev inequality yieldOn the event $A_T\ge Tp/2$, the drift coefficient contributes an average whose error is at most $R{h}^{\alpha }$, by Hölder continuity. The Itô isometry and stationarity give $\mathbb{E}{N}_T^2\le \Vert \sigma {\Vert }_{\infty }^{2}Tp$. Another use of the Chebyshev inequality provesThis also accounts for the zero-denominator convention through the first exceptional event. The bound proves pointwise $R{h}^{\alpha }+O_{\mathbb{P}}((Th{)}^{-1/2})$, uniformly in the location's tail probabilities. It supplies no bound for the probability of a supremum over all locations.

For an explicit counterexample, take $\sigma \equiv 1$ and the bounded globally Lipschitz functionIt has Lipschitz constant $R=1$, exponent $\alpha =1$, and satisfies the inward-drift conditions with $M=2$, $\gamma =2$. The stochastic differential equation has a unique strong solution of a stochastic differential equation. Its Invariant distribution of an Itô diffusion has densityThe exponent is zero on $[-1,1]$, equals $-(|x|-1{)}^2$ when $1<|x|<2$, and equals $3-2|x|$ outside $[-2,2]$. This verifies integrability, positivity and a constant density ${\rho }_0=Z^{-1}$ on the central interval. Start the Itô diffusion in this Invariant distribution of an Itô diffusion to obtain the required stationary process.

Set $h=T^{-1/3}$ and choose $N=\lfloor (4h{)}^{-1}\rfloor$ windows with centers $x_j=-1/2+(4j-2)h$, $1\le j\le N$. Their intervals $I_j=[x_j-h,x_j+h]$ are disjoint and lie in $[-1,1]$, where the drift is zero. WriteThen ${\hat{b}}_T(x_j,h)=M_j(T)/A_j(T)$ whenever $A_j(T)>0$. The same occupation bound as above gives $\mathbb{E}(A_j(T)-q{)}^2\le C_{1}T{h}^2$ for $T\ge 1$. With $\delta =T^{-1/6}$, a union bound and the Chebyshev inequality implyThe quadratic variations are $⟨M_j⟩=A_j$, and the cross quadratic covariations vanish since the windows are disjoint. Each clock tends to infinity almost surely: for any fixed window, the occupation estimate at times $2^k$, followed by the Borel-Cantelli lemmas, gives $A_j(2^k)/(2^k\mu (I_j))\to 1$. The Knight theorem for orthogonal martingales therefore represents $M_j(t)=B_j(A_j(t))$ using independent standard Brownian motions $B_1,\dots ,B_N$. This theorem is applied to the finite collection for each $T$; no growing-dimensional central limit theorem is being assumed.

On the preceding clock event, compare each time-changed Brownian motion with its value at deterministic time $q$. The Brownian reflection principle and a union bound show that, for every fixed $\epsilon >0$,This bound does not require the Brownian motions to be independent of their clocks. In particular, ${\max }_j|M_j(T)-B_j(q)|/\sqrt{q}\to 0$ in probability. The variables $B_j(q)/\sqrt{q}$ are independent $N(0,1)$ random variables. Their Gaussian maximum exceeds $\sqrt{\log N}$ with probability tending to one: for $r\ge 1$, integration of the normal density over $[r,r+1/r]$ gives $\mathbb{P}(|Z|>r)\ge c{r}^{-1}{e}^{-r^2/2}$, and henceSince $A_j(T)\le (1+\delta )q$ on the clock event, it follows that, for some $c_3>0$,But $\sqrt{Th}Rh=1$ for this bandwidth. Thus the error after subtraction of $Rh$ and multiplication by $\sqrt{Th}$ is not bounded in probability, contradicting the printed stochastic order. The valid pointwise bound above is provable; the uniform shrinking-bandwidth assertion is not. This counterexample shows why a uniform rate needs at least a $\sqrt{\log (1/h)}$ enlargement of the noise scale. If $h$ is instead held fixed and constants may depend on $h$, this counterexample does not contradict that different interpretation.

For the standard normal density, $P^{\prime }(x)=-xP(x)$. Using integration by parts,The boundary term vanishes because Gaussian decay dominates every power. Since the given mean is zero and variance is one, $\mathbb{E}[X^2]=1$, yielding the Gaussian fourth moment

The standard normal density is even, so $x^{3}P(x)$ is odd. Since Gaussian decay dominates every power, the integrand is absolutely integrable. Therefore

The displayed identity requires a correction to the process construction in the source. In a piecewise-deterministic Markov process, the instantaneous hazard function is evaluated at the current flowing state $X_{u-}$, and the destination is sampled from $K_u(X_{u-},dy)$. Keeping the jump times and destinations of the original Markov jump process $Y$ instead evaluates both quantities at its last reset state $Y_{u-}$, which generally differs from $X_{u-}$. The source PDF confirms that these original jump times and destinations were retained. Consequently its definition does not in general produce the asserted transition kernel, or even a Markov process in the state variable $X$ alone.

For an explicit counterexample, take $a(t,x)=1$, $\alpha (x)=2+\tanh x$, $y_0=0$, and let every $K_t(x,dy)$ be the standard normal distribution, independently of $t,x$. These are smooth coefficients with bounded positive jump rate, so nonexplosion holds. In the literal construction, the first holding time has the exponential distribution of rate $2$. Thus the probability distribution of $X_t$ has an atom of mass $e^{-2t}$ at $t$. Paths with at least one jump contribute no atom there, because the reset has a normal density. The asserted identity instead assigns to this no-jump atom the masswhich differs for every $t>0$.

Here is the full derivation for the intended piecewise-deterministic Markov process. Assume nonexplosion; smoothness alone does not guarantee it. Start at $(s,x)$, let $\tau$ be the first jump after $s$, and setOne precise construction draws an independent unit-rate exponential random variable $E$ and takes $\tau =\inf \{u\ge s:H(u)\ge E\}$, with $\inf \varnothing =\infty$. This implements a hazard function along the flow of a vector field, rather than a fixed exponential holding time at the last reset. Its survival probability and jump-time probability density function areEquivalently, $S^{\prime }=-\alpha (z)S$ with $S(s)=1$. On $\{\tau >t\}$, the terminal state is $z(t)$. Conditional on $\tau =u$, draw the new state $y$ from $K_u(z(u),dy)$ and restart the same construction from $(u,y)$. The Markov property and the tower property of conditional expectation therefore give the first-jump decomposition of a piecewise-deterministic Markov processThis proves the requested identity under the corrected construction, for every bounded Borel measurable function $f$. The two contributions partition according to whether the first jump occurs before $t$; later jumps are already contained in $P_{u,t}$.

For comparison, the literal construction does have a different renewal identity when started at a reset with $X_s=Y_s=x$. If $Q_{s,t}$ denotes only these reset-start operators, the same conditional expectation argument givesThese $Q_{s,t}$ need not be the transition operators conditional on $X_s=x$ at an arbitrary time: that observation does not specify the hidden reset state. This distinction explains exactly why replacing only the between-jump paths is insufficient.