Pith. sign in

REVIEW 2 major objections 3 minor 36 references

Explosion in the quasi-Gaussian HJM model

T0 review · 2 major / 3 minor · reviewed 2026-08-14 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read The quasi-Gaussian HJM short rate explodes in finite time with positive probability for CEV exponents in $(1/2,1]$, and almost surely for large initial rates.

desk verdict Worth a serious referee, but Theorem 1 has a missing hitting-Γ argument that must be repaired. read the letter →

arxiv 1908.07102 v1 pith:TT2HS46G submitted 2019-08-19 q-fin.MF math.PR

Dan Pirjol , Lingjiong Zhu This is my paper · ORCID
classification q-fin.MFmath.PR MSC 60H1060J6091G30
keywords HJMmodelshortrateexplosionCEVvolatilityLyapunovfunctionmultidimensionaldiffusionzerocouponbondEurodollarfuturesfinite-time
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

The paper studies the one-factor quasi-Gaussian Heath-Jarrow-Morton interest-rate model with an $\varepsilon$-CEV volatility specification, in which the short rate's volatility behaves like $\sigma r^\gamma$ for large $r$. It establishes that for elasticity parameter $\gamma\in(1/2,1]$, including the log-normal case $\gamma=1$, the short rate explodes to infinity in finite time with positive probability whenever one of two explicit parameter conditions holds; for $\beta>0$ and a sufficiently large initial short rate, the explosion occurs with probability one. The practical consequence is that a finite-time explosion makes zero-coupon bond prices collapse to zero and Eurodollar futures prices infinite at sufficiently long maturities, so the model's useful pricing horizon is bounded by the explosion time. The paper also proves that for $\gamma\in(0,1/2]$ the model is non-explosive.

What carries the argument

The argument is carried by a Lyapunov function $V(r,y)=C_1-\frac{C_2}{(1+y)^{\delta_1}}-\frac{C_3}{(1+r)^{\delta_2}}$ with $(1+\delta_1)(1+\delta_2)=2\gamma$, together with a comparison theorem that reduces the time-dependent forward-curve model to a time-homogeneous system when $\lambda'(t)+\beta\lambda(t)\ge\beta\lambda(0)$. Applying the infinitesimal generator of the diffusion to $V$ and invoking the Lyapunov-function explosion criterion of [12] turns the explosion problem into two algebraic inequalities in $\beta,\sigma,\gamma$, giving the explicit conditions (i) and (ii). A second auxiliary function $V_0=e^{-r}+e^{-y}$ is used to verify the hitting-time condition needed for almost-sure explosion.

What would settle it

Simulate the time-homogeneous SDE (20)-(21) with parameters satisfying condition (i) or (ii), for example $\gamma=1$, $\beta=0.05$, $\sigma=0.2$, $r_0=0.1$, over a long horizon and estimate the probability of explosion from a large ensemble; the proof guarantees this probability is positive, so a carefully controlled run finding no explosion at all would contradict the claim. To test the comparison step directly, simulate the time-inhomogeneous model (8)-(9) with a forward curve that violates $\lambda'(t)+\beta\lambda(t)\ge\beta\lambda(0)$ and check whether explosions persist.

Watch

Extended reading notes

Core claim

The central discovery is that explosion is a genuine property of the quasi-Gaussian HJM model with a CEV-type volatility, not an artefact of a small-noise approximation. For the time-homogeneous model with constant initial forward rate $r_0$, Theorem 1 shows that the two-dimensional diffusion $(r_t,y_t)$ with $\varepsilon$-CEV volatility explodes with positive probability for $\gamma\in(1/2,1]$ provided that either $\sup_{R\ge\varepsilon}F(R;\beta,\sigma)>0$ or $\sup_{R\ge\varepsilon}\bigl(G(R)-(2\beta+\tfrac12\sigma^2\delta_2(\delta_2+1))\bigr)\ge0$, where $F$ and $G$ are explicit functions of the model parameters and $\delta_1,\delta_2>0$ satisfy $(1+\delta_1)(1+\delta_2)=2\gamma$. Under the extra assumptions $\beta>0$ and a large enough initial short rate satisfying inequality (40), Theorem 2 upgrades this to almost-sure explosion. The log-normal case $\gamma=1$ is included, and the authors note that this matches explosions observed numerically and previously derived in a deterministic small-noise limit.

Load-bearing premise

The most fragile load-bearing premise is the comparison step near equation (19): the explosion proved for the time-homogeneous auxiliary model is transferred to the original time-inhomogeneous model only when the forward curve satisfies $\lambda'(t)+\beta\lambda(t)\ge\beta\lambda(0)$, an inequality that can fail for a declining forward curve.

Editorial extensions

If this is right

  • Zero-coupon bonds with maturity beyond the explosion time have price zero with positive probability, because $P(T,T+\delta)$ contains the factor $\exp(-G(T,T+\delta)x_T-\tfrac12 G(T,T+\delta)^2y_T)$ and $x_T,y_T$ explode.
  • Eurodollar futures prices, which involve $\mathbb{E}^Q[P^{-1}(T,T+\delta)]$, become infinite for sufficiently large $T$ when the Theorem 1 conditions hold, and LIBOR-linked caps, swaptions and CMS products inherit the singularity.
  • Explosion is confined to $\gamma\in(1/2,1]$; for $0<\gamma\le1/2$ the coefficients satisfy a linear-growth condition and the solution is non-explosive and square-integrable.
  • With mean reversion $\beta>0$ and a sufficiently high initial short rate satisfying inequality (40), explosion occurs almost surely rather than merely with positive probability.
  • Capping the short-rate volatility at a finite level restores sub-linear growth and eliminates explosion, a practical remedy for pricing applications.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The two conditions in Theorem 1 are sufficient, not necessary; the true explosion region in the $(\beta,\sigma,\gamma)$ parameter space is likely larger, so explosions may also occur for parameter sets outside the proved region.
  • The comparison assumption $\lambda'(t)+\beta\lambda(t)\ge\beta\lambda(0)$ means the proved explosion applies most directly to flat or rising forward curves; for downward-sloping yield curves the mechanism could be delayed or suppressed, and this is testable by simulation.
  • The same Lyapunov ansatz may extend to multi-factor quasi-Gaussian models, where no necessary-and-sufficient explosion criterion is known, potentially yielding explicit sufficient conditions for explosion in higher-dimensional term-structure models.
  • The almost-sure explosion in Theorem 2 requires a large initial short rate; for small initial rates only positive-probability explosion is guaranteed, so low-rate calibrations may remain on the non-explosive side in practice.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

2 major / 3 minor

Summary. The paper analyzes the ε-CEV modified one-factor quasi-Gaussian HJM short-rate model (8)-(9). Its main results are for the time-homogeneous case λ(t)=λ0: Theorem 1 proves non-explosion for γ∈(0,1/2] and, for γ∈(1/2,1], finite-time explosion with positive probability under either of the parameter conditions (i) sup_R F(R;β,σ)>0 or (ii) sup_R (G(R)-(2β+1/2σ^2δ2(δ2+1)))≥0. Theorem 2 upgrades this to almost-sure explosion when β>0 and r0 satisfies (40). The proof uses a Lyapunov function V(r,y)=C1-C2/(1+y)^{δ1}-C3/(1+r)^{δ2} together with the Chow-Khasminskii explosion criteria. The paper closes with implications for zero-coupon bond prices and Eurodollar futures and a small numerical illustration.

Significance. The result, if fully established, is a useful rigorous addition to the interest-rate modeling literature: it identifies a wide and practically relevant parameter region in which the quasi-Gaussian HJM model with CEV-type volatility has explosive short-rate paths, and it quantifies the consequences for derivatives pricing. The construction of the Lyapunov function is explicit and the parameter conditions are checkable; the numerical section illustrates that the proved region is not vacuous. The main obstruction is that the proof of Theorem 1 as written does not connect the initial condition (r0,0) to the set Γ where the Khasminskii criterion applies, and the comparison reduction at Eq. (19) is asserted without proof. Both are localizable and appear repairable, so the manuscript merits a major revision rather than rejection.

major comments (2)
  1. [§3, proof of Theorem 1(b)] The proof verifies hypotheses (A.1)-(A.3) of Proposition 1 and then concludes P(τ<∞)>0. Proposition 1, as quoted in the paper, gives this conclusion only when the process starts at a point of Γ=[2R,∞)×[2R,∞). The initial condition of Theorem 1 is (r0,0); this point never lies in Γ because y0=0, and for the parameter values used in Figure 1 (r0=0.1, R=1/δ2=1) it lies in D. The paper gives no argument that the diffusion reaches Γ with positive probability from (r0,0). The inequality LV≥CV is only used on D^c and says nothing about the crossing from D to Γ. Please add a reachability lemma, for example via the support theorem and an explicit control path from (r0,0) to Γ, or restrict the theorem to initial conditions in Γ.
  2. [§3, Eq. (19)] The reduction from the time-inhomogeneous system (8)-(9) to the time-homogeneous system (20)-(21) is asserted in one sentence with a reference to Yamada [36], but no comparison proof is given. For the statement to hold one needs both λ'(t)+βλ(t)≥βλ(0) and the monotone dependence of y on r, and this should be shown explicitly. Since a declining forward curve can violate (19), the paper should either prove the comparison or clearly state that the rigorous explosion results are confined to λ(t)≡r0 and that Theorem 1 does not cover general forward curves.
minor comments (3)
  1. [Appendix, Eq. (58)] The definition of x should be x = r^{δ2+1}/y, not x = r^{δ1+1}/y. With the printed definition the two terms in (57) do not collapse to a single function of x; with the corrected exponent the derivation is valid.
  2. [Appendix, Eq. (69)] For the record, Eq. (69) is consistent with Eq. (68): substituting a and b from (71) into (68) reproduces (69) exactly; the apparent discrepancy in the ordering of the terms inside the minimum is only a notational matter.
  3. [§3.1, Eq. (38)] The small-δ2 expansion of G(R0(δ2)) appears to have a sign error: from G(δ)=δ^{1+δ}(1+δ)^{-(1+δ)} one obtains G(δ)=δ+δ^2(logδ-1)+O(δ^3(logδ)^2), not δ2+δ2^2(logδ2+1). The leading-order conclusion σ_max=√2 is unaffected.

Circularity Check

0 steps flagged · score 2.0 of 10

No material circularity: the central explosion theorem rests on an independent Lyapunov-function proof; the only self-citation is motivational/illustrative and is not load-bearing.

full rationale

The main derivation chain, Theorems 1 and 2, is self-contained apart from the use of external explosion criteria. The paper explicitly constructs a Lyapunov function (47), verifies the conditions (A.1)-(A.3) of Proposition 1 in the proof of Theorem 1, and reduces the parameter conditions (i) and (ii) to the algebraic inequalities (89) and (90) obtained from Lemma 2. No parameter is fitted to explosion data, and no 'prediction' is obtained by renaming a fitted input. The comparison step around Eq. (19) is asserted via an external comparison theorem of Yamada [36], and the main theorems in fact assume the time-homogeneous case λ(t) ≡ r0, so this comparison is not the load-bearing step. The only self-citation [30] appears in the introduction as motivation and in the numerical example as the small-noise explosion time and critical value β_C; it is not used to prove Theorem 1(b) or Theorem 2. Thus there is no circular reduction by construction. The score of 2 reflects the presence of a minor self-citation that is not load-bearing; it does not indicate circularity in the central proof. A possible proof gap concerning the application of Proposition 1 to initial data outside Γ is a correctness issue, not a circularity issue, and does not raise this score.

Assumptions & free parameters 4 free parameters · 6 assumptions · 0 invented entities

The central claim rests on standard probabilistic theorems (comparison, Lyapunov explosion criteria), on the model structure of the quasi-Gaussian HJM framework, and on the forward-curve inequality (19). The only ad hoc element is the ε-CEV regularization inherited from prior literature. The Lyapunov constants and auxiliary exponents δ1,δ2,R are proof devices, not empirical fits.

free parameters (4)
  • C2, C3 = not fixed numerically; existence only
    Positive constants in the Lyapunov function (47); conditions (76)/(77) in Lemma 2 ensure a non-empty range, but the theorem does not need their actual values.
  • δ1, δ2 = any positive pair satisfying (1+δ1)(1+δ2)=2γ
    Exponents in the Lyapunov function; the theorem requires existence of at least one pair for which condition (i) or (ii) holds, so they are tuneable auxiliary parameters.
  • R = any R≥ε in the sup in (24)/(25)
    Domain size for the Lyapunov argument; conditions (i) and (ii) quantify over R, so R is an auxiliary variable, not a model input.
  • ε = small positive cutoff
    Model regularization parameter in (7); the theorem requires ε>0 sufficiently small and r0>ε. It is a model modification inherited from [1], not fitted to data.
assumptions (6)
  • standard math The comparison theorem of Yamada [36] applies to the SDEs (8)-(9) and (20)-(21) under the drift ordering condition (19).
    Used in Section 2 and near Eq. (19) to bound the original short rate below by the time-homogeneous solution; the exact hypotheses are not verified in the text.
  • standard math Theorem 1 and Theorem 2 of Chow-Khasminskii [12] provide sufficient conditions for explosion with positive probability and almost surely.
    The entire proof of Theorem 1 rests on verifying conditions (A.1)-(A.3); Theorem 2 adds (A.4)-(A.5). These criteria are taken as given.
  • domain assumption The one-factor quasi-Gaussian HJM model admits the two-state Markov representation (3) under the separable volatility form (2).
    This is the modeling setup from [31,3]; the SDEs (8)-(9) describe the short rate under the risk-neutral measure Q.
  • domain assumption The forward curve satisfies λ'(t)+βλ(t) ≥ βλ(0) (Eq. 19).
    Needed for the comparison step that reduces the time-inhomogeneous model to the time-homogeneous SDE; it is an assumption on the input yield curve.
  • ad hoc to paper The ε-CEV modification (7) leaves the large-r behavior unchanged and keeps the origin unattainable.
    The modification is standard in the literature [1], but the proof of explosion is stated for the modified process; transferring to the plain CEV model relies on this behavior.
  • domain assumption The origin (0,0) is unattainable and the boundary ∂D is regular.
    Used in the proof of Theorem 1 to apply Proposition 1; argued via Eq. (46) and the positive drift βr0 in r, though for β=0 the argument relies on y>0.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Explosion in the quasi-Gaussian HJM model." pith.science (2026). https://pith.science/paper/TT2HS46G

@misc{pith2026190807102,
  author       = {Pith},
  title        = {Pith review of: Explosion in the quasi-Gaussian HJM model},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/TT2HS46G}},
  note         = {Machine review of arXiv:1908.07102}
}
read the original abstract

We study the explosion of the solutions of the SDE in the quasi-Gaussian HJM model with a CEV-type volatility. The quasi-Gaussian HJM models are a popular approach for modeling the dynamics of the yield curve. This is due to their low dimensional Markovian representation which simplifies their numerical implementation and simulation. We show rigorously that the short rate in these models explodes in finite time with positive probability, under certain assumptions for the model parameters, and that the explosion occurs in finite time with probability one under some stronger assumptions. We discuss the implications of these results for the pricing of the zero coupon bonds and Eurodollar futures under this model.

Figures

Figures reproduced from arXiv: 1908.07102 by the authors.

Figure 1
Figure 1. Sample paths for {rt}t≥0 for σ = 0.2 and λ = 0.1 in the log-normal quasi￾Gaussian HJM model. Left plot: sample paths with β = 0. The red vertical line is at texp = 47.03y, which is the small-noise explosion time following from Proposition 4 in [30]. Right plot: β = 0.05. corresponding explosion time can be found in closed form for λ(t) = λ0, and is given in Proposition 4 of [30]. The prediction is shown in [PITH_FU… view at source ↗
Figure 2
Figure 2. Regions D and Γ for the application of Proposition 1. The infinitesimal generator of this diffusion is LεV (r, y) = (σ 2 r 2 min(r 2γ−2 , ε2γ−2 (22) ) − 2βy)∂yV +(y − βr + βr0)∂rV + 1 2 σ 2 r 2 min(r 2γ−2 , ε2γ−2 )∂ 2 rV . We would like to study the explosion time of this diffusion, defined as (23) τ := sup{t > 0 : yt < ∞, rt < ∞}. We present a few preliminary results which will be used in our proof. The following t… view at source ↗
Figure 3
Figure 3. Region in the (σ, β) plane allowed by the condition (ii) in Theorem 1. For given γ ∈ ( 1 2 , 1], this condition is satisfied for β below the curves shown. (ii) For case (ii) of Theorem 1, C2, C3 > 0 satisfy the inequality δ2C3 δ1C2σ 2  R 1 + R δ2−δ1 ≤ min  R δ2(δ1+2) , κ1 R−δ2 − κ2  (33) , where κ1, κ2 are defined in (73), (74), and R = 1 δ2 . Remark 2. The constraints of Theorem 1 can be made stronger by replac… view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

36 extracted references · 35 canonical work pages

  1. [36]

    (1973) On a comparison theorem for solutions of stochastic differential equations and its applications

    Yamada, T. (1973) On a comparison theorem for solutions of stochastic differential equations and its applications. J. Math. Kyoto Univ. 13 497-512. E-mail address: dpirjol@gmail.com E-mail address: ling@cims.nyu.edu

    work page 1973

  2. [1]

    Andersen, L. and J. Andreasen. (2000). Volatility skews and extensions of the LIBOR market model. Applied Mathematical Finance 7(1), 1-32

    work page 2000

  3. [2]

    Andersen, L. and V. Piterbarg. (2007). Moment explosions in stochastic volatility models. Finance and Stochas- tics. 11, 29-50

    work page 2007

  4. [3]

    Andersen, L. and V. Piterbarg. (2010). Interest Rate Modeling. Atlantic Press, London

    work page 2010

  5. [4]

    Andreasen, J. (2010). Markovian term structure models, in Encyclopedia of Quantitative Finance , Ed. R. Cont, Wiley, New York

    work page 2010

  6. [5]

    (1993) Generalized Vasicek models of the term structure, Applied Stochastic Models and Data Analysis

    Babbs, S. (1993) Generalized Vasicek models of the term structure, Applied Stochastic Models and Data Analysis. 1, 49-62

    work page 1993

  7. [6]

    Cakici, N. and J. Zhu. (2001) Pricing Eurodollar futures options with the Heath-Jarrow-Morton model. Journal of Futures Markets. 21, 655-680

    work page 2001

  8. [7]

    Cao, L. and P. Henry-Labordere. (2016) Interest rate models enhanced with local volatility. SSRN id=2793125

    work page 2016

Show all 36 references
  1. [8]

    (1992) Markov representation of the Heath-Jarrow-Morton model

    Cheyette, O. (1992) Markov representation of the Heath-Jarrow-Morton model. Working paper, Barra

    1992

  2. [9]

    (2012) Explicit volatility specification for the linear Cheyette model

    Chibane, M. (2012) Explicit volatility specification for the linear Cheyette model. Working paper, Shinsei Bank Ltd

    2012

  3. [10]

    Chibane, M. and D. Law. (2013) A quadratic volatility Cheyette model. Risk, 60-63 (July)

    2013

  4. [11]

    Chien, H. H. (2001) On the complexity of the Ritchken Sankarasubramanian interest rate model. Ph.D. Thesis, National Taiwan University

    2001

  5. [12]

    Khasminskii

    Chow, P.-L., and R. Khasminskii. (2014) Almost sure explosion of solutions to stochastic differential equations. Stoch. Processes and their Applications 124, 639-645

    2014

  6. [13]

    (1996) Notes on option pricing I: Constant elasticity of diffusion

    Cox, J. (1996) Notes on option pricing I: Constant elasticity of diffusion. J. Portfolio Management 124, 639-645

    1996

  7. [14]

    (1952) The parabolic differential equations and the associated groups of transformations

    Feller, W. (1952) The parabolic differential equations and the associated groups of transformations. Annals of Mathematics 55, 468-519

    1952

  8. [15]

    Jablecki and D

    Gatarek, D., J. Jablecki and D. Qu. (2016) Non-parametric local volatility formula for interest rate swaptions. Risk, 1-5 (February)

    2016

  9. [16]

    Jarrow and A

    Heath, D., R. Jarrow and A. Morton. (1992) Bond pricing and the term structure of interest rates: A new methodology for contingent claim valuation. Econometrica 60, 77-105

    1992

  10. [17]

    (1993) Problems in certain two-factor term structure models

    Hogan, M. (1993) Problems in certain two-factor term structure models. Ann. Appl. Prob. 3(2), 576-581

    1993

  11. [18]

    Hogan, M. and K. Weintraub. (1998). The lognormal interest rate model and Eurodollar futures. Citibank Working Paper. QUASI-GAUSSIAN HJM MODEL 21

    1998

  12. [19]

    Hoorens, B. (2011). On the Cheyette short rate model with stochastic volatility, MSc Thesis, Delft University of Technology

    2011

  13. [20]

    Hunt, P. J. and J. E. Kennedy. (2004) Financial Derivatives in Theory and Practice . Wiley Series in Probability and Statistics, Chichester

    2004

  14. [21]

    (1991) Bond and option evaluation in the Gaussian interest rate model

    Jamshidian, F. (1991) Bond and option evaluation in the Gaussian interest rate model. Research in Finance 8, 131-170

    1991

  15. [22]

    Karatzas, I. and S. Shreve. (1991) Brownian Motion and Stochastic Calculus . Graduate Texts in Mathematics, Volume 113, Springer Science & Business Media, Second Edition, New York

    1991

  16. [23]

    Karatzas, I. and J. Ruf. (2016) Distribution of the time to explosion for one-dimensional diffusions. Prob. Theory and Related Fields. 164, 1027-1069

    2016

  17. [24]

    (2012) Stochastic Stability of Differential Equations

    Khasminskii, R. (2012) Stochastic Stability of Differential Equations . Springer, New York

    2012

  18. [25]

    (1969) Stochastic Integrals

    McKean, H.P.Jr. (1969) Stochastic Integrals. Academic Press, New York

    1969

  19. [26]

    Morton, A. J. (1989) Arbitrage and martingales. Doctoral dissertation, Cornell University

    1989

  20. [27]

    Acar and M

    Natcheva-Acar, K., S.K. Acar and M. Krekel (2015). Modeling credit spreads with the Cheyette model and its applications to credit default swaptions. Journal of Credit Risk 5, 47-71

    2015

  21. [28]

    Pinsky, R. G. (2008) Positive Harmonic Functions and Diffusion. Cambridge Studies in Advanced Mathematics v. 45, Cambridge University Press, Cambridge

    2008

  22. [29]

    (2014) Hogan-Weintraub singularity and explosive behaviour in the Black-Derman-Toy model

    Pirjol, D. (2014) Hogan-Weintraub singularity and explosive behaviour in the Black-Derman-Toy model. Quan- titative Finance 15, 1243-1257

    2014

  23. [30]

    Pirjol, D. and L. Zhu. (2017) Small-noise limit of the quasi-Gaussian lognormal HJM model. Operations Research Letters. 45, 6-11

    2017

  24. [31]

    Ritchken, P. and L. Sankarasubramanian. (1995) Volatility structures of forward rates and the dynamics of the term structure. Math. Finance 5, 55-72

    1995

  25. [32]

    Sandmann, K. and D. Sondermann. (1997) A note on the stability of lognormal interest rate models and the pricing of Eurodollar futures. Math. Finance 7(2), 119-125

    1997

  26. [33]

    Stroock, D. W. and S. R. S. Varadhan. (1979) Multidimensional Diffusion Processes. Grundlehren der mathema- tischen Wissenschaften, Volume 233. Springer-Verlag, Berlin

    1979

  27. [34]

    (2007) Some results on strong solutions of SDEs with applications to interest rate models

    Wissel, J. (2007) Some results on strong solutions of SDEs with applications to interest rate models. Stochastic Processes and their Applications 117, 720-741

    2007

  28. [35]

    Xing, J., and Y. Li. (2016) Nonlinear Lyapunov criteria for stochastic explosive solutions, Stat. Prob. Lett. 109 63-67

    2016

Pith tools

Reviewed August 14, 2026 · model on record in the stance chip above.