Semiclassical CEV Option Pricing Model: an Analytical Approach

Jose A. Capit\'an; Jose Lope-Alba; Juan J. Morales-Ruiz

REVIEW 2 major objections 4 minor 21 references

Semiclassical CEV Option Pricing Model: an Analytical Approach

T0 review · 2 major / 4 minor · reviewed 2026-08-12 · deepseek-v4-flash

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

Pith's one-line read A closed-form WKB heat kernel prices CEV options, restoring a missing exponential factor.

desk verdict Useful semiclassical CEV kernel derivation undermined by a likely typo in the central formula (3.28), which needs correction and re-validation before the paper's claims can be trusted. read the letter →

arxiv 2411.18154 v5 pith:4UJX4XHY submitted 2024-11-27 q-fin.MF

Jose A. Capitán , Jose Lope-Alba , Juan J. Morales-Ruiz This is my paper · ORCID

classification q-fin.MF MSC 91G2081Q20

keywords CEVmodelsemiclassicalapproximationWKBheatkernelVanVleck-Morettedeterminantoptionpricingvolatilityskewclosed-formsolution

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 derives a closed-form semiclassical (WKB) approximation for the heat kernel of the CEV diffusion, the option-pricing model whose volatility varies as a power of the asset price. It claims that using the Van Vleck-Morette determinant, instead of the Van Vleck determinant used in earlier work, yields a much simpler explicit formula for the propagator. The calculation also uncovers an exponential prefactor that previous semiclassical CEV pricing formulas had dropped. If correct, the formula gives a fast analytical route to option prices in regimes of strong volatility skew and tail risk, where the approximation is shown to be accurate for short maturities.

What carries the argument

The load-bearing object is the Pauli-Morette semiclassical formula $K_{\mathrm{WKB}} = (2\pi J)^{-1/2}\, e^{\frac{1}{2}\int \frac{\partial^2 H}{\partial p\partial x}\,d\tau}\, e^{-S(\gamma)}$, with $J$ the Van Vleck-Morette determinant. For the CEV Hamiltonian $H = 2xp^2 + (bx-a)p$, the paper evaluates the action $S(\gamma)$, the cross-derivative integral, and $J$ in closed form, expressing all three through two integration constants $D_1$ and $D_2$ fixed by the endpoints. The determinant is computed two independent ways — from $\partial x/\partial p_T$ and from the variational equations — and both give the same expression, which is what allows the prefactor to be written explicitly.

What would settle it

Find a triple $(x_T, x, T)$ in the pricing domain where the determinant $J$ in (3.25) vanishes or changes sign; at such a point the WKB kernel diverges, so the closed-form price could not match Monte Carlo results, and the formula would need a caustic or multi-path correction.

Watch

Extended reading notes

Core claim

The central claim is that equation (3.28), with the constants $D_1$ and $D_2$ given by (3.20)-(3.21), is the closed-form semiclassical heat kernel for the CEV model, expressed directly in terms of the initial and final positions and the time to maturity. The kernel is built from the Pauli-Morette WKB formula with the Van Vleck-Morette determinant $J = \partial x/\partial p_T(T)$, computed both from the Hamiltonian flow and from the variational equations. The novelty over the earlier WKB treatment is an additional exponential factor $e^{\frac{1}{2}\int_0^T \frac{\partial^2 H}{\partial p\partial x}\,d\tau}$, which the paper shows is required when the Hamiltonian mixes position and momentum. The authors also give the European call price as the convolution of this kernel with the payoff and verify the approximation against Monte Carlo simulation, finding good agreement for short maturities and for exponents $\alpha$ in $[-1, -0.4]$.

Load-bearing premise

The entire construction rests on the Pauli-Morette formula being valid for the CEV diffusion Hamiltonian and on there being one real, non-caustic classical path between every endpoint pair and every maturity; the paper does not analyze points where the Van Vleck-Morette determinant $J$ vanishes.

Editorial extensions

If this is right

The closed-form kernel allows European call prices under CEV to be evaluated by a one-dimensional integral rather than by solving the PDE numerically.
The corrected prefactor implies that earlier WKB prices carry a systematic bias that grows with the cross-derivative term of the Hamiltonian.
Because the kernel is explicit in $x_T$ and $x$, option Greeks and parameter sensitivities can in principle be derived analytically from (3.28).
The approximation is most reliable exactly where the CEV model departs from Black-Scholes: short maturity and strong volatility skew ($\alpha$ near $-1$).
For $\alpha$ close to zero, the singular change of variables degrades accuracy, so the formula should not be used as a Black-Scholes replacement.

Reading between the lines

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

The same Van Vleck-Morette-plus-cross-derivative recipe could be applied to other one-factor diffusions whose Hamiltonian is polynomial, provided the associated variational equations are integrable.
A natural stress test is to search parameter regions where $J(T)$ changes sign or vanishes; the authors note $J(T)>0$ for small $T$ but do not rule out zeros at finite $T$, where the kernel would blow up and a caustic treatment would be needed.
The missing exponential factor is likely present, not only in CEV, but in any WKB pricing formula derived from a mixed position-momentum Hamiltonian; earlier diffusion WKB results in other models may deserve the same correction.
One could use the exact confluent-hypergeometric solution of the CEV Cauchy problem as a benchmark to quantify when the semiclassical kernel breaks down, rather than relying only on Monte Carlo averages.

Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

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

Desk Editor's Note

read the letter

The paper does something genuinely useful: it derives a semiclassical heat kernel for the CEV model using the Van Vleck–Morette determinant and a Hamiltonian formulation, rather than the Van Vleck determinant used by Araneda et al. The derivation is mostly checkable, and the two independent computations of the Van Vleck–Morette determinant (direct and via variational equations) agreeing is a real strength. The Black–Scholes sanity check is also a nice touch. If the final formula is correct, it is simpler than the earlier WKB result and can be evaluated quickly, which matters for a model used in practice.

The problem is that the central printed result, Eq. (3.28), does not follow from the preceding equations. The exponent in (3.27) contains b/(8D2)(D1^2 - d^2)(e^{-bT}-1). When you combine the logarithms and the -1/2 bT term, you get a clean factor (1 - bd/2) log(...) plus b/(8D2)(D1^2 - d^2)(e^{-bT}-1) plus 1/2 bT(db-1). There is no algebraic step that replaces d^2 with b^2. And dimensionally, D1 and d carry the dimension of the transformed variable, while b has inverse-time dimension, so D1^2 - b^2 mixes incompatible quantities. This looks like a typo, but it is a load-bearing one: Eq. (3.28) is the formula used in the pricing integral (3.31), and presumably in the Monte Carlo comparison in Appendix C. If the numerics actually used (3.28) with b^2, then the reported accuracy is for a different, dimensionally inconsistent kernel.

A second issue: the paper claims an exponential factor was missing in Araneda et al., but it never reproduces that paper's formula or shows the comparison explicitly. That claim needs to be substantiated, not just asserted. A third, softer gap is the caustic question: the paper notes J>0 for small T but does not analyze whether J can vanish for finite T in the pricing domain. That is not fatal, but it should be acknowledged.

The derivation and the two computations of J show clear thinking, and the paper deserves a serious referee. But the authors need to fix the typo, confirm which formula the numerics actually evaluated, and either show the direct comparison with [3] or soften the claim. If those are addressed, this would be a solid contribution for people working on semiclassical option pricing or CEV approximations.

Referee Report

2 major / 4 minor

Summary. This paper derives a closed-form semiclassical (WKB) estimate of the heat kernel for the CEV diffusion after the change of variables S_t -> X_t, using the Pauli-Morette formula (1.8). The authors compute the action on the classical path, the time integral of the mixed second derivative of the Hamiltonian, and the Van Vleck-Morette determinant by two independent methods (direct differentiation of the endpoint map and the variational equations), obtaining Eq. (3.25) in both cases. The resulting kernel, Eq. (3.28), is then used in a convolution integral, Eq. (3.31), to price European calls, and the prices are compared with Monte Carlo simulations in Appendix C. The paper claims that the exponential prefactor e^{1/2∫∂²H/∂p∂x dτ} is missing in the earlier WKB formula of Araneda et al.

Significance. If the printed formula is corrected, the paper offers a concise analytical approximation to the CEV pricing kernel that is considerably simpler than the previous WKB expression, with a prefactor determined by a Van Vleck-Morette determinant computed in two independent ways and with no parameters fitted to data. The Monte Carlo comparisons are an independent check of the approximation, not a calibration. However, the central printed equation (3.28) is inconsistent with the preceding derivation (3.27), so the main claim as stated is not supported until corrected.

major comments (2)

[§3.2, Eq. (3.28)] The step from Eq. (3.27) to Eq. (3.28) is invalid as printed. Combining the exponentials in (3.27) gives the factor R^{1-bd/2} e^{(1/2)bT(db-1)} e^{b/(8D2)(D1^2-d^2)(e^{-bT}-1)} with R=(2D2e^{bT}+D1-d)/(2D2+D1-d); no algebraic rearrangement can change d^2 into b^2. Since Eq. (3.28) is the kernel used in Eq. (3.31) and in the numerical tests of Appendix C, the central result and its validation depend on the wrong printed formula. Please restore d^2, and re-run or confirm the Monte Carlo comparisons with the corrected expression.
[§3.2, Eqs. (3.20)-(3.21) and (3.25)] The Pauli-Morette formula (1.8) presupposes a single real, non-caustic classical path connecting x_T and x. The endpoint formulas (3.20)-(3.21) take the positive square root without specifying the branch or the parameter conditions under which D1 and D2 produce a real path with x(τ)>0, and no proof is given that J in (3.25) is nonzero on the whole integration domain of (3.31); J(0)=0 shows that zeros are not excluded a priori. The authors should either prove J≠0 for all relevant endpoints and T, or state the restriction to the caustic-free domain and adjust the pricing integral accordingly.

minor comments (4)

[§3.2, Eqs. (3.20)-(3.21)] The variables in Section 3 are not explicitly nondimensionalized; as written, formulas such as (3.20) combine quantities that would have different units unless a scaling is understood. Please state the units or the scaling used.
[§1, Eq. (2.14)] The Black-Scholes formula for d1 in Eq. (2.14) is garbled: it should be [log(S_T/E)+(r+σ^2/2)T]/(σ√T), with a plus sign and a factor 1/2. Please fix the formula.
[§3.1 and Appendix A] The claim in Appendix A that the non-integrability of the confluent hypergeometric equation leaves 'no hope' of a closed-form solution is too categorical; special functions are standard closed forms, and the cited Martinet-Ramis conditions concern integrability by quadratures. Please rephrase.
[§3.2, after Eq. (3.16)] The statement that C1 and C2 are positive 'because they are given by suitable exponential functions' is not justified by the displayed expressions and requires either a proof or an explicit condition on the parameters.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity found: the semiclassical kernel is derived from the standard Pauli-Morette formula with no fitted inputs and Monte Carlo checks are independent; the algebraic mismatch between Eqs. (3.27) and (3.28) is a correctness issue, not circularity.

full rationale

The core derivation is self-contained rather than circular. Equation (1.8) is taken from the external semiclassical literature ([12, 11]) and applied to the CEV Hamiltonian (3.6), which is obtained from the Feller-form PDE (3.4) via the standard change of variables. The action S(gamma) in Eq. (3.23), the integral in Eq. (3.24), and the Van Vleck-Morette determinant J in Eq. (3.25) are all computed explicitly from the classical Hamiltonian flow and are cross-checked by a second computation via the variational equations in Appendix B. No parameter is fitted to option prices before the kernel is used; the Monte Carlo simulations in Appendix C serve as an independent numerical benchmark, not as a calibration input. The self-citations to [15] and [17] are not load-bearing: the paper does not rely on a uniqueness theorem to force its formula, and the closed form is exhibited by direct calculation. The claimed missing exponential factor is obtained by evaluating the Hamiltonian mixed derivative integral, Eq. (3.24), not by renaming an empirical pattern or by redefining a fitted quantity. A separate issue is that Eq. (3.28) does not follow algebraically from Eq. (3.27): the printed formula replaces d^2 by b^2 in the term b/(8D2)(D1^2 - d^2)(e^{-bT}-1), and no step in the paper justifies this change. That is an internal consistency or correctness problem, not a circularity, so it does not raise the circularity score. Similarly, the unanalyzed caustics where J vanishes would be a domain-of-validity gap rather than a circular reduction. Overall, the derivation chain does not reduce to its own inputs.

Assumptions & free parameters 0 free parameters · 3 assumptions · 0 invented entities

The CEV parameters (sigma, alpha, mu) are model inputs, not fitted. The WKB derivation itself introduces no adjustable parameters. The main external inputs are the semiclassical formula and the integrability framework, both cited from prior literature.

assumptions (3)

domain assumption The Pauli-Morette semiclassical formula (1.8) with the additional exponential prefactor e^{(1/2)∫∂²H/∂p∂x dτ} is valid for diffusion equations of the form (1.6).
This formula is cited from [12,11] and is the foundation of the whole derivation; the paper does not prove it for the CEV Hamiltonian.
standard math The classical Hamiltonian system for CEV is integrable and its variational equations are solvable in closed form (guaranteed by differential Galois theory for one degree of freedom).
Used to justify that J can be computed in closed form; supported by references [15,17].
ad hoc to paper The solution of the endpoint equations selects a real branch of the square root in Eqs (3.20)-(3.21) such that D1 and D2 define the classical path; the paper does not specify the branch or conditions for positivity.
The kernel formula depends on D1,D2 written with an ambiguous square root; without a branch rule the formula is not fully defined.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Semiclassical CEV Option Pricing Model: an Analytical Approach." pith.science (2026). https://pith.science/paper/4UJX4XHY

@misc{pith2026241118154,
  author       = {Pith},
  title        = {Pith review of: Semiclassical CEV Option Pricing Model: an Analytical Approach},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/4UJX4XHY}},
  note         = {Machine review of arXiv:2411.18154}
}

read the original abstract

This paper is devoted to obtain closed form solutions for the semiclassical (or WKB) approximation of the heat kernel propagator of the diffusion equation defined by the constant elasticity variance (CEV) option pricing model. One of the key points is that our calculations are based on the Van Vleck-Morette determinant instead of the Van Vleck determinant used by other authors. In fact, we compute this determinant in two different ways: by means of the solution of the classical Hamiltonian equations, and by solving the variational equations. Furthermore, the calculation reveals an exponential factor in the prefactor of the kernel not considered in previous works.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

21 extracted references · 21 canonical work pages

[3]

Journal of Computational and Applied Mathematics388(2021): 113244

Araneda, A.A., and Villena, M.J., Computing the CEV option pricing formula using the semiclassical approxima- tion of path integral. Journal of Computational and Applied Mathematics388(2021): 113244. 15

work page 2021
[1]

and Morales-Ruiz, J.J., Integrability of stochastic birth-death processes via differential Galois theory.Mathematical Modelling of Natural Phenomena15(2020) 70

Acosta-Hum´ anez, P.B., Capit´ an, J.A. and Morales-Ruiz, J.J., Integrability of stochastic birth-death processes via differential Galois theory.Mathematical Modelling of Natural Phenomena15(2020) 70

work page 2020
[2]

and Pantazi, C., Semiclassical quantification of some two degree of freedom potentials: A differential Galois approach

Acosta-Hum´ anez, P., L´ azaro, J.T., Morales-Ruiz, J.J. and Pantazi, C., Semiclassical quantification of some two degree of freedom potentials: A differential Galois approach. Journal of Mathematical Physics65(2024)

work page 2024
[4]

Journal of Political Economy81(1973) 637–654

Black, F., and Scholes, M., The pricing of options and corporate liabilities. Journal of Political Economy81(1973) 637–654

work page 1973
[5]

The Journal of Portfolio Management23 (1996), 15–17

Cox, J.C., The Constant Elasticity of Variance Option Pricing Model. The Journal of Portfolio Management23 (1996), 15–17

work page 1996
[6]

Journal of Financial Economics3(1976) 145–166

Cox, J.C., and Ross, S.A., The valuation of options for alternative stochastic processes. Journal of Financial Economics3(1976) 145–166

work page 1976
[7]

The semiclassical expansion

DeWitt-Morette, C. The semiclassical expansion. Annals of Physics,97(1976), 367–399

work page 1976
[8]

R.,Quantum Mechanics and Path Integrals

Feynman R.P., Hibbs A. R.,Quantum Mechanics and Path Integrals. McGraw-Hill, New York, 1965

work page 1965

Show all 21 references

[9]

Feller W., Two singular diffusion problems. Ann. Math54(1961) 173–182

1961
[10]

Glasserman, P.Monte Carlo methods in financial engineering. Vol. 53. Springer, New York, 2004

2004
[11]

Springer, 2007

Kolokoltsov, V.N.Semiclassical analysis for diffusions and stochastic processes. Springer, 2007

2007
[12]

Reidel, Dortrech, 1982

Langouche F., Rockaerts D., Tirapegui E.Functional Integration and Semiclassical Expansion. Reidel, Dortrech, 1982

1982
[13]

The constant elasticity of variance model.Encyclopedia of Quantitative Finance328–

Linetsky, V., Mendoza, R. The constant elasticity of variance model.Encyclopedia of Quantitative Finance328–
[14]

and Martinet, J., 1990

Ramis, J.P. and Martinet, J., 1990. Th´ eorie de Galois diff´ erentielle et resommation.Computer Algebra and Dif- ferential Equations, Comput. Math. Appl., Academic Press, London, pp. 117–214

1990
[15]

Journal of Mathematical Physics61(2020)

Morales-Ruiz, J.J., A differential Galois approach to path integrals. Journal of Mathematical Physics61(2020)

2020
[16]

Modern Birkh¨ auser Classics, Reprint 2013 of the 1999 Edition by Birkh¨ auser, Springer, Basel

Morales-Ruiz, J.J.,Differential Galois Theory and Non-integrability of Hamiltonian Systems. Modern Birkh¨ auser Classics, Reprint 2013 of the 1999 Edition by Birkh¨ auser, Springer, Basel

2013
[17]

Methods and Applications of Analysis8(2001) 33–96

Morales-Ruiz J.J., Ramis J.-P., Galoisian obstructions to integrability of Hamiltonian systems. Methods and Applications of Analysis8(2001) 33–96

2001
[18]

On the definition and approximation of Feynman’s path integrals

Morette, C. On the definition and approximation of Feynman’s path integrals. Physical Review81(1951) 848

1951
[19]

Pauli, W.,Selected topics in field quantization, Pauli Lecures in Physics vol. 6. Dover, 2000 (based on lectures given by Pauli at Zurich in 1950–51)

2000
[20]

The Journal of Finance44 (1989), 211–219

Schroder, M., Computing the constant elasticity of variance option pricing formula. The Journal of Finance44 (1989), 211–219. 16

1989
[334]

John Wiley & Sons, 2010

2010

Pith tools

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

[1] [3]

Journal of Computational and Applied Mathematics388(2021): 113244

Araneda, A.A., and Villena, M.J., Computing the CEV option pricing formula using the semiclassical approxima- tion of path integral. Journal of Computational and Applied Mathematics388(2021): 113244. 15

work page 2021

[2] [1]

and Morales-Ruiz, J.J., Integrability of stochastic birth-death processes via differential Galois theory.Mathematical Modelling of Natural Phenomena15(2020) 70

Acosta-Hum´ anez, P.B., Capit´ an, J.A. and Morales-Ruiz, J.J., Integrability of stochastic birth-death processes via differential Galois theory.Mathematical Modelling of Natural Phenomena15(2020) 70

work page 2020

[3] [2]

and Pantazi, C., Semiclassical quantification of some two degree of freedom potentials: A differential Galois approach

Acosta-Hum´ anez, P., L´ azaro, J.T., Morales-Ruiz, J.J. and Pantazi, C., Semiclassical quantification of some two degree of freedom potentials: A differential Galois approach. Journal of Mathematical Physics65(2024)

work page 2024

[4] [4]

Journal of Political Economy81(1973) 637–654

Black, F., and Scholes, M., The pricing of options and corporate liabilities. Journal of Political Economy81(1973) 637–654

work page 1973

[5] [5]

The Journal of Portfolio Management23 (1996), 15–17

Cox, J.C., The Constant Elasticity of Variance Option Pricing Model. The Journal of Portfolio Management23 (1996), 15–17

work page 1996

[6] [6]

Journal of Financial Economics3(1976) 145–166

Cox, J.C., and Ross, S.A., The valuation of options for alternative stochastic processes. Journal of Financial Economics3(1976) 145–166

work page 1976

[7] [7]

The semiclassical expansion

DeWitt-Morette, C. The semiclassical expansion. Annals of Physics,97(1976), 367–399

work page 1976

[8] [8]

R.,Quantum Mechanics and Path Integrals

Feynman R.P., Hibbs A. R.,Quantum Mechanics and Path Integrals. McGraw-Hill, New York, 1965

work page 1965