REVIEW 2 major objections 3 minor 45 references
Asymptotics of survival probabilities and lower tail probability problem
T0 review · 2 major / 3 minor · reviewed 2026-08-10 · deepseek-v4-flash
Pith's one-line read The paper proves that for Stieltjes-Lévy processes, as $t\to\infty$ the survival probability $P[\bar X_t<x]$ is $1-p_\infty(x)+O(t^{-1})$ when $\mu_1<0$, $p_{\infty,0}(x)t^{-1/2}+O(t^{-1})$ when $\mu_1=0$, and $O(t^{-1})$ when $\mu_1>0$…
desk verdict Explicit SL-process asymptotics are a real contribution, but the proof of the O(t^{-1}) rates has a genuine contour-deformation gap near q=0 that needs fixing. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The machinery is the contour-integral representation of the survival probability together with the Wiener-Hopf factorization identity. For a Lévy process with characteristic exponent $\psi$, the probability is written as a double contour integral; the paper deforms the $q$-contour to $\operatorname{Re} q=\sigma/t$, uses integration by parts to show the large $|\operatorname{Im} q|$ part is $O(t^{-1})$, and then studies $q$ near $0$. The small-$q$ behavior of the positive Wiener-Hopf factor $\varphi_q^+$ is captured through a modified factor $\varphi_q^{-,0}$ defined by (6.13), for which the paper proves uniformity in $q$ near the origin and an explicit representation (6.15) in terms of zeros of $\psi$ and the Stieltjes-Lévy measure $G_-$. The link between $\psi$ and the Lévy measure is Theorem 3.16, which expresses the jump density as an integral of $\operatorname{Im}\psi$ on the boundary of the strip of analyticity; this is what turns the abstract asymptotic coefficients into explicit SL-measure integrals.
What would settle it
Take a KoBoL process with $\mu_1<0$ and numerically evaluate $P[\bar X_t<x]$ for large $t$ by inverting the Fourier-Laplace representation with high precision; if the difference from $1-p_\infty(x)$ is not $O(t^{-1})$ or has a different constant than (6.16), the claim fails. A more direct check is to search the strip $S[\omega'_{-,0}]$ for zeros of $q+\psi(\xi)$ with $\operatorname{Re} q=\sigma/t$ and $|\operatorname{Im} q|\ge \epsilon$: any zero invalidates the integration-by-parts step.
Extended reading notes
Core claim
The central discovery is that the large-time survival probability for an SL-process is asymptotically determined by local behavior of the characteristic exponent at zero and by the absolutely continuous component of the Stieltjes-Lévy measure. Specifically, if $\mu_1<0$, then $P[\bar X_t<x] = 1-p_\infty(x)+O(t^{-1})$ with $p_\infty(x)$ given by (6.16); if $\mu_1=0$, then $P[\bar X_t<x] = p_{\infty,0}(x)t^{-1/2}+O(t^{-1})$ with $p_{\infty,0}(x)$ given by (6.20); and if $\mu_1>0$, then $P[\bar X_t<x] = O(t^{-1})$. The coefficients $p_\infty(x)$ and $p_{\infty,0}(x)$ are expressed in terms of the density of the absolutely continuous component of the Stieltjes-Lévy measure $G_+$, the modified Wiener-Hopf factor $\varphi_0^{-,0}$, the first two instantaneous moments $\mu_1,\mu_2$, and the zeros of the characteristic exponent on the imaginary axis. The proof works by deforming the $q$-integration contour in the Laplace-Fourier inversion formula, isolating the small-$q$ contribution, and expanding the Wiener-Hopf factors in terms of the roots of $q+\psi(\xi)=0$ and the SL-measures. The lower tail problem is solved at the same level of generality: for SINH-regular processes, $P[\bar X_t<x] = \kappa_k(t)x^{\nu_+}+O(x^{\nu_++s})$ as $x\downarrow 0$, where $\nu_+$ is the decay rate of the positive Wiener-Hopf factor and $\kappa_k(t)$ is a Laplace inversion integral.
Load-bearing premise
The argument assumes that $q+\psi(\xi)$ has no zeros on the entire deformed contour $\operatorname{Re} q=\sigma/t$ with $|\operatorname{Im} q|\ge \epsilon$ inside the strip $S[\omega'_{-,0}]$; the paper verifies this only in a small ball around $q=0$, so an unnoticed zero elsewhere would spoil the error estimate.
Editorial extensions
If this is right
- For SL-processes with $\mu_1<0$, survival to a fixed level $x$ converges to a positive limit $1-p_\infty(x)$ at rate $1/t$, giving a first-order correction to ruin probabilities in insurance and barrier option prices.
- At zero drift, the survival probability decays as $t^{-1/2}$ with an explicit prefactor, so the persistence exponent is $1/2$ and the constant can be computed from the Lévy measure.
- For positive drift, survival up to a fixed level is $O(t^{-1})$, so the probability of never exceeding $x$ is asymptotically small at least at that rate.
- The lower tail result gives $P[\bar X_t<x] \approx \kappa(t)x^{\nu_+}$ for small $x$, so the exponent $\nu_+$ and prefactor $\kappa(t)$ determine the near-barrier behavior of no-touch options for any fixed maturity.
- The explicit formulas allow direct computation of asymptotics once the Stieltjes-Lévy representation is available, without solving the Wiener-Hopf equation numerically.
Reading between the lines
- The zero-drift $t^{-1/2}$ prefactor depends on $\psi$ only through $\mu_2$ and $\varphi_0^{-,0}$, suggesting a universality class: processes with equal $\mu_2$ and matched small-$x$ structure share the same leading persistence rate, which could be checked against simulations.
- Because the paper proves zero-freeness of $q+\psi(\xi)$ only near $q=0$, one could numerically search for stray zeros on the deformed contour; if found for a particular SL-process, the $O(t^{-1})$ error rate would be replaced by a different rate, and the formulas would need modification.
- The same small-$q$ Wiener-Hopf expansion may apply to other extremal functionals, such as the joint law of the running maximum and the process value, opening a route to asymptotic joint densities.
- For mixtures of stable processes and regular SL-processes, the paper notes the survival problem is open; one could attempt a rescaling argument to reduce it to the lower tail problem, as suggested in Section 6.1.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper derives asymptotic formulas for the survival probability P[\bar X_t < x] as t→∞ for fixed x>0, for Lévy processes with exponentially decaying tails, and refines the coefficients for Stieltjes-Lévy (SL) processes in terms of SL-measures and zeros of the characteristic exponent. It also reformulates lower tail probability asymptotics as x↓0 for SINH-regular processes. The main results are (6.1)–(6.3), with explicit expressions for p∞(x) and p∞,0(x) in (6.16) and (6.20), obtained by contour deformation and residue calculus.
Significance. If the proofs are completed, the paper gives explicit, parametric formulas for the leading terms of survival probabilities, expressed in terms of objects that are often explicitly computable for popular Lévy models (KoBoL, NIG, NTS, VG). A strong point is that the coefficient formulas (6.16) and (6.20) reproduce known Brownian motion results when specialized, providing a consistency check. The paper is an addendum to the authors' prior work, but the survival-probability asymptotics for SL-processes appear to be new.
major comments (2)
- [§6.2, after (6.9)] The O(t^{-1}) estimate for the integrals over {q: Re q = σ/t, |Im q| ≥ ε} uses the bound (6.9) on the entire deformed contour. Lemma 6.2 establishes a lower bound for |q + ψ(ξ)| only for q in a small ball B(0,ε) and ξ in a small strip S(-ω',ω') around 0; it says nothing about the line Re q = σ/t when |Im q| ≥ ε, where q is not small. For an SL-process, the zero-location statement of Proposition 4.2 applies to real q > 0, so for complex q with Re q = σ/t > 0 and |Im q| ≥ ε there may be zeros of q + ψ(ξ) inside the strip S[ω'_-,0]. If such a zero is crossed during the deformation, the residue calculus in (6.10), (6.17) and (6.18) omits a contribution, and the claimed O(t^{-1}) rate could fail or pick up a t^{-1/2} term, which would alter the leading coefficient in (6.1) and (6.3). The authors should prove a global zero-free region for the deformed contour, or otherwise avoid deforming over the full line.
- [§6.2, 'Integrating by parts'] The derivation of the O(t^{-1}) rate by integration by parts is only sketched. For each of the integrals in (6.4), (6.6) and (6.8), one needs to identify the boundary terms as |Im q| → ∞ and justify that they vanish or are negligible, and one needs uniform estimates on the q-derivative of the entire integrand, including the inner ξ-integral and the factor φ_q^+(ξ). Without these estimates, the statement 'Integrating by parts and using (6.9), we obtain' is not a complete proof. Since the O(t^{-1}) accuracy is part of the theorems' claims, the details should be supplied or a reference given where these estimates appear.
minor comments (3)
- [Theorem 6.8] Theorem 6.8 states that (6.1) holds with p∞(x), but the case μ1 = 0 should correspond to (6.3) with p∞,0(x); the displayed formula is the coefficient of the t^{-1/2} term.
- [§6.2, after (6.8)] The notation 'S(ω'_- ,0]' should be 'S[ω'_-,0]' for consistency with the rest of the paper.
- [§2.2.4] The abbreviations 'VGP' and 'VG' are used interchangeably; using one consistent label would improve readability.
Circularity Check
No significant circularity: the asymptotic coefficients are derived by contour deformation from integral representations, and prior SL-process results are used as theorems, not as fitted inputs.
full rationale
The paper's central claims (6.16) and (6.20) are obtained by substituting analytic representations of the Wiener-Hopf factors into the Laplace-Fourier inversion formulas (6.4), (6.6), (6.8), then deforming contours and applying residue calculus. The coefficients p_infty(x) and p_infty,0(x) are expressed in terms of the characteristic exponent psi, its zeros, and the SL-measures G_±; the latter are themselves defined from psi via Theorem 3.16(b), so the final formulas are explicit rewritings of the same model ingredients, not independently fitted quantities. No parameter is calibrated to the survival probability being predicted, and no 'prediction' is obtained by substituting a quantity that was fit to that same target. The paper relies heavily on the authors' previous papers [11, 12, 13, 18] for definitions of RLPE, SINH-regular and (s)SL processes and for integral representations of Wiener-Hopf factors; this is normal use of prior mathematical work. The load-bearing zero-location statement for SL-processes, 'for q > 0, the equation q + psi(xi) = 0 has no zeros outside the imaginary axis iR' (attributed to [18], Section 1), is a prior theorem with stated assumptions that do not include the target asymptotic result; citing it does not make the present derivation circular. The genuine weakness identified in the proof is the claim in Section 6.2 that the integrals over |Im q| >= epsilon are O(t^{-1}) using (6.9), which requires q + psi(xi) to be zero-free on the whole deformed contour, whereas Lemma 6.2 only establishes the zero structure for q in a small ball around 0. This is a correctness gap in the proof, not a circularity, because even if the gap is filled, the derivation is still from the characteristic exponent and SL-measures rather than from the quantity being predicted. The paper also explicitly notes in Remark 4.3.1 that for sSL-processes zeros may lie outside iR and a proof is not supplied; again, this is a stated limitation, not a circular step. No step in the derivation reduces by construction to its own inputs, and no fitted input is renamed as a prediction. Score 0 is therefore appropriate.
Assumptions & free parameters
assumptions (6)
- standard math Standard complex analysis: Cauchy integral formula, analytic continuation, residue theorem.
- domain assumption The process X is an SL-process (or sSL/SINH-regular), so ψ is analytic in a strip and satisfies the asymptotic conditions of Definition 2.1.
- domain assumption The integral representation (6.4)-(6.6) holds for Lévy processes with ψ analytic in a strip (cited from [7]).
- ad hoc to paper There exists ω'_- < 0 such that q + ψ(ξ) ≠ 0 for all q on the deformed contour Re q = σ/t and all ξ in S[ω'_-,0].
- ad hoc to paper Technical integrability condition (b) in Theorems 6.5 and 6.8: ∫_{U+} |ψ(-iw-0)|^{-2} G+(dw) < ∞.
- domain assumption Zeros of ψ on i(-∞,0) are simple.
Cite this review
Pith. "Pith review of Asymptotics of survival probabilities and lower tail probability problem." pith.science (2026). https://pith.science/paper/2P3DFU4L
@misc{pith2026250104218,
author = {Pith},
title = {Pith review of: Asymptotics of survival probabilities and lower tail probability problem},
year = {2026},
howpublished = {\url{https://pith.science/paper/2P3DFU4L}},
note = {Machine review of arXiv:2501.04218}
}
abstract
The present paper is an addendum to the paper ``L\'evy models amenable to efficient calculations", where we introduced a general class of Stieltjes-L\'evy processes (SL-processes) and signed SL processes defined in terms of certain Stieltjes-L\'evy measures. We demonstrated that SL-processes enjoyed all properties that we used earlier to develop efficient methods for evaluation of expectations of functions of a L\'evy process and its extremum processes, and proved that essentially all popular classes of L\'evy processes are SL-processes; sSL-processes fail to possess one important property. In the present paper, we use the properties of (s)SL-processes to derive new formulas for the Wiener-Hopf factors $\phi^\pm_q$ for small $q$ in terms of the absolute continuous components of SL-measures and their densities, and calculate the leading terms of the survival probability also in terms of the absolute continuous components of SL-measures and their densities. The lower tail probability is calculated for more general classes of SINH-regular processes constructed earlier.
Reference graph
Works this paper leans on
- [1]
-
[2]
S. Asmussen, F. Avram, and M.R. Pistorius. Russian and American put options under exponential phase- type L´ evy models.Stochastic Processes and their Applications , 109(1):79–111, 2004
work page 2004
-
[3]
Persistence probabilities \& exponents
F. Aurzada and T. Simon. Persistence probabilities & exponents. Working paper, March 2012. https://arXiv.org/abs/1203.6554
work page Pith review arXiv 2012
-
[4]
O.E. Barndorff-Nielsen. Processes of Normal Inverse Gaussian Type. Finance and Stochastics , 2:41–68, 1998
work page 1998
-
[5]
O.E. Barndorff-Nielsen and S.Z. Levendorski ˇi. Feller Processes of Normal Inverse Gaussian type. Quantita- tive Finance, 1:318–331, 2001
work page 2001
-
[6]
M. Boyarchenko, M. de Innocentis, and S. Levendorski ˘i. Prices of barrier and first-touch digital options in L´ evy-driven models, near barrier.International Journal of Theoretical and Applied Finance , 14(7):1045– 1090, 2011. Available at SSRN: http://papers.ssrn.com/abstract=1514025
work page 2011
-
[7]
M. Boyarchenko and S. Levendorski ˘i. Prices and sensitivities of barrier and first-touch digital options in L´ evy-driven models.International Journal of Theoretical and Applied Finance , 12(8):1125–1170, December 2009
work page 2009
-
[8]
M. Boyarchenko and S. Levendorski ˘i. Ghost Calibration and Pricing Barrier Options and Credit Default Swaps in spectrally one-sided L´ evy models: The Parabolic Laplace Inversion Method.Quantitative Finance, 15(3):421–441, 2015. Available at SSRN: http://ssrn.com/abstract=2445318
work page 2015
Show all 45 references
-
[9]
Boyarchenko and S
S. Boyarchenko and S. Levendorski ˘i. Generalizations of the Black-Scholes equation for truncated L´ evy processes. Working Paper, University of Pennsylvania, April 1999
1999
-
[10]
Boyarchenko and S
S. Boyarchenko and S. Levendorski ˘i. Option pricing for truncated L´ evy processes.International Journal of Theoretical and Applied Finance, 3(3):549–552, July 2000
2000
-
[11]
Boyarchenko and S
S. Boyarchenko and S. Levendorski˘i. Barrier options and touch-and-out options under regular L´ evy processes of exponential type. Annals of Applied Probability , 12(4):1261–1298, 2002
2002
-
[12]
Boyarchenko and S
S. Boyarchenko and S. Levendorski ˘i. Non-Gaussian Merton-Black-Scholes Theory , volume 9 of Adv. Ser. Stat. Sci. Appl. Probab. World Scientific Publishing Co., River Edge, NJ, 2002
2002
-
[13]
Boyarchenko and S
S. Boyarchenko and S. Levendorski ˘i. Perpetual American options under L´ evy processes.SIAM Journal on Control and Optimization , 40(6):1663–1696, 2002
2002
-
[14]
Boyarchenko and S
S. Boyarchenko and S. Levendorski ˘i. Sinh-acceleration: Efficient evaluation of probability distri- butions, option pricing, and Monte-Carlo simulations. International Journal of Theoretical and Applied Finance , 22(3):1950–011, 2019. DOI: 10.1142/S0219024919500110. Available a...
1950 doi
-
[15]
Boyarchenko and S
S. Boyarchenko and S. Levendorski ˘i. Conformal accelerations method and efficient evaluation of stable distributions. Acta Applicandae Mathematicae , 169:711–765, 2020. Available at SSRN: https://ssrn.com/abstract=3206696 or http://dx.doi.org/10.2139/ssrn.3206696. 24 SVETLANA B...
2020 doi
-
[16]
Boyarchenko and S
S. Boyarchenko and S. Levendorski ˘i. Static and semi-static hedging as contrarian or conformist bets. Mathematical Finance , 3(30):921–960, 2020. Available at SSRN: https://ssrn.com/abstract=3329694 or http://arXiv.org/abs/1902.02854
2020 arXiv
-
[17]
Boyarchenko and S
S. Boyarchenko and S. Levendorski˘i. Efficient evaluation of expectations of functions of a stable L´ evy process and its extremum. Working paper, September 2022. Available at SSRN: http://ssrn.com/abstract=4229032 or http://arxiv.org/abs/2209.12349
2022 arXiv
-
[18]
Boyarchenko and S
S. Boyarchenko and S. Levendorski ˘i. L´ evy models amenable to efficient calculations. Working paper, June
-
[19]
Boyarchenko and S
S. Boyarchenko and S. Levendorski ˘i. Efficient evaluation of joint pdf of a L´ evy process, its ex- tremum, and hitting time of the extremum. Working paper, December 2023. Available at SSRN: http://ssrn.com/abstract=4656492 or http://arxiv.org/abs/2312.05222
2023 arXiv
-
[20]
Boyarchenko and S
S. Boyarchenko and S. Levendorski ˘i. Simulation of a L´ evy process, its extremum, and hitting time of the extremum via characteristic functions. Working paper, December 2023. Available at SSRN: http://ssrn.com/abstract=4656466 or http://arxiv.org/abs/2312.03929
2023 arXiv
-
[21]
Boyarchenko and S
S. Boyarchenko and S. Levendorski ˘i. Efficient evaluation of double barrier options. Interna- tional Journal of Theoretical and Applied Finance , 27(2):2450007, 2024. Available at SSRN: http://ssrn.com/abstract=4262396 or http://arxiv.org/abs/2211.07765
2024 arXiv
-
[22]
Boyarchenko and S
S. Boyarchenko and S. Levendorski ˘i. Efficient evaluation of expectations of functions of a L´ evy process and its extremum. Finance and Stochastics, 2024. To appear. Working paper version is available at SSRN: https://ssrn.com/abstract=4140462 or http://arXiv.org/abs/2209.02793
2024 arXiv
-
[23]
Boyarchenko and S
S. Boyarchenko and S. Levendorski˘i. Faster than fft. Proceedings icfda2024.sciencesconf.org 12th Conference on Fractional Differentiation and its Applications 9-12 July 2024, Bordeaux, France, October 2024
2024
-
[24]
Cai and S.G
N. Cai and S.G. Kou. Option pricing under a mixed-exponential jump diffusion model. Operations Research, 60(1):64–77, 2012
2012
-
[25]
P. Carr, H. Geman, D.B. Madan, and M. Yor. The fine structure of asset returns: an empirical investigation. Journal of Business , 75:305–332, 2002
2002
-
[26]
Eberlein, Y
E. Eberlein, Y. Kabanov, and T. Sshmidt. Ruin probabilities for a sparre andersen model with investment. Stochastic Processes and their Applications , 144:72–84, 2022
2022
-
[27]
Furling, P
J.R. Furling, P. Salminen, and P. Vallois. On joint hit distribution of the running maximum of Brownian motion. Stochastic Processes and Their Applications , 150:1204–1221, 2022
2022
-
[28]
S.G. Kou. A jump-diffusion model for option pricing. Management Science, 48(8):1086–1101, August 2002
2002
-
[29]
Kuznetsov
A. Kuznetsov. Wiener-Hopf factorization and distribution of extrema for a family of L´ evy processes. Ann.Appl.Prob., 20(5):1801–1830, 2010
2010
-
[30]
Kuznetsov, A.E
A. Kuznetsov, A.E. Kyprianou, and J.C. Pardo. Meromorphic L´ evy processes and their fluctuation identi- ties. Annals of Applied Probability , 22(3):1101–1135, 2012
2012
-
[31]
Levendorski ˘i
S. Levendorski ˘i. Degenerate Elliptic Equations , volume 258 of Mathematics and its Applications . Kluwer Academic Publishers Group, Dordrecht, 1993
1993
-
[32]
Levendorski ˘i
S. Levendorski ˘i. Pricing of the American put under L´ evy processes. Research Report Ma- PhySto, Aarhus, 2002. Available at http://www.maphysto.dk/publications/MPS-RR/2002/44.pdf, http://www.maphysto.dk/cgi-bin/gp.cgi?publ=441
2002
-
[33]
Levendorski ˘i
S. Levendorski ˘i. Pricing of the American put under L´ evy processes. International Journal of Theoretical and Applied Finance, 7(3):303–335, May 2004
2004
-
[34]
Levendorski ˘i
S. Levendorski ˘i. Convergence of Carr’s Randomization Approximation Near Barrier. SIAM FM , 2(1):79– 111, 2011
2011
-
[35]
Levendorski ˘i
S. Levendorski ˘i. Method of paired contours and pricing barrier options and CDS of long maturities. Inter- national Journal of Theoretical and Applied Finance , 17(5):1–58, 2014. 1450033 (58 pages)
2014
-
[36]
Levendorski ˘i and B
S. Levendorski ˘i and B. Paneyakh. Degenerate elliptic equations and boundary value problems. In M.S. Agranovich and M.A. Shubin, editors, Encyclopeadia of Mathematical Sciences, Vol.63 , pages 131–202. Springer, Berlin, 1994
1994
-
[37]
A. Lipton. Assets with jumps. Risk, pages 149–153, September 2002
2002
-
[38]
Madan and F
D.B. Madan and F. Milne. Option pricing with V.G. martingale components. Mathematical Finance , 1(4):39–55, 1991. ASYMPTOTICS OF SURVIVAL PROBABILITIES AND LOWER TAIL PROBABILITY PROBLEM 25
1991
-
[39]
Madan and M
D.B. Madan and M. Yor. Representing the CGMY and Meixner L´ evy processes as time changed Brownian motions. Journal of Computational Finance , 12(1):27–47, 2009
2009
-
[40]
R.C. Merton. Option pricing when underlying stock returns are discontinuous. Journal of Financial Eco- nomics, 3:125–144, 1976
1976
-
[41]
Pitman and M
J. Pitman and M. Yor. Infinitely divisible laws associated with hyperbolic functions. Canadian Journal of Mathematics, 55(2):292–330, 2003
2003
-
[42]
K. Sato. L´ evy processes and infinitely divisible distributions , volume 68 of Cambridge Stud. Adv. Math. Cambridge University Press, Cambridge, 1999
1999
-
[43]
Schilling, R
R. Schilling, R. Song, and Z. Vodrachˇ ek.Bernstein Functions: Theory and Applications. De Gruyter, Berlin, 2d edition, 2012
2012
-
[44]
Schoutens and J.L
W. Schoutens and J.L. Teugels. L´ evy processes, polynomials and martingaless. Communications in Statis- tics: Stochastic Models , 14:335–349, 1998
1998
-
[2022]
Available at SSRN: https://ssrn.com/abstract=4116959 or http://arXiv.org/abs/2207.02359
Reviewed August 10, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.