REVIEW 3 major objections 3 minor 41 references
A recursive formula for the $n^\text{th}$ survival function and the $n^\text{th}$ first passage time distribution for jump and diffusion processes. Applications to the pricing of $n^\text{th}$-to-default CDS
T0 review · 3 major / 3 minor · reviewed 2026-08-05 · deepseek-v4-flash
Pith's one-line read A recursive formula gives the nth survivor's survival function for many-body Markov processes with killing boundaries.
desk verdict Two-coordinate core and examples are solid, but the many-coordinate theorem is asserted rather than proved; worth a referee, not a desk reject. 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 central object is the directed graph Γ whose nodes are binary strings labeling alive (A) and dead (D) coordinates, and the path contributions I_g(t) (Definition 1, Eq. 44). Each I_g is a nested integral over the ordered kill times τ1 < ... < τn−1, built from the conditional density P^N_n of the surviving coordinates given the full kill history, the marginal first-passage densities F_{n+1}, ..., F_N of successively reduced systems, and the transition density of the final survivor. The path sum over G^n, the set of paths from the all-alive state to states with exactly n alive coordinates, carries the recursion in Theorem 2 (Eq. 45).
What would settle it
Take two correlated Brownian motions with a killing boundary and a post-kill drift for the survivor that depends on which coordinate was killed first (e.g., the survivor's drift is +1 if coordinate 1 dies first and −1 if coordinate 2 dies first). Compute the last-survivor survival function by Monte Carlo simulation and compare with the Theorem 1 formula evaluated numerically; a statistically significant discrepancy in the regime where the two kill times are close would show the factorization fails.
Extended reading notes
Core claim
The central claim is that for N Markovian coordinates with a killing boundary, the survival function of the nth surviving coordinate satisfies the recursion S^n(t) = S^{n+1}(t) + sum_{g in G^n} I_g(t), where each I_g is a nested integral over the ordered kill times of a path through the graph of alive/dead configurations, involving the product of conditional densities of the survivor given the history and the marginal first-passage densities of successively reduced systems. The two-coordinate case (Theorem 1) is the base case, writing the last-survivor survival function as the joint survival function plus integrals over the first kill time of the survivor's subsequent independent evolution.
Load-bearing premise
The joint density of the ordered kill history is assumed to factorize into the product of marginal first-passage densities of successively reduced systems times a conditional density of the surviving coordinates, and the survivor's post-kill evolution is assumed to be Markovian with the reduced dynamics——this is the premise that carries the nested-integral formula.
Editorial extensions
If this is right
- The recursion provides an exact, if combinatorially heavy, method to compute order statistics of first passage times for any Markovian many-body system with killing, without resorting to Monte Carlo simulation.
- It yields semi-analytical pricing formulas for nth-to-default credit default swaps, with the fee and protection legs expressed directly through the survival functions S^{N-n+1}.
- In the single-file diffusion example, the formula reproduces the known reflection-principle result, confirming the method on a nontrivial interacting diffusion.
- The path-sum structure makes explicit how the kill history enters the dynamics, highlighting that the reduced survivor is Markovian only conditionally on the full ordered sequence of deaths.
- The framework naturally handles simultaneous killings, as in the multivariate Poisson model, where multiple coordinates can hit the boundary at the same time.
Reading between the lines
- The recursive path structure could suggest approximate schemes for large N by truncating or averaging over paths in Γ, potentially reaching mean-field-like descriptions of absorbing particle systems.
- If the factorization of the kill-history density fails for systems where the survivor's post-kill law depends on the identity rather than just the position of the killed coordinate, corrections to the formula would be needed; such cases may motivate generalized conditional densities.
- The paper's emphasis on path dependence connects naturally to persistence exponents and extreme-value statistics of correlated processes, where the order of extreme events matters.
- The CDS application could be extended to collateralized debt obligations or to models with default-dependent intensity term structures, as the paper itself suggests, while retaining the path-sum form.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper claims a recursive formula for the survival function of the n-th surviving coordinate of an N-coordinate Markov process with a killing boundary. Theorem 1 treats the two-coordinate case, expressing the last-survivor survival function as the joint survival function plus integrals over the first-kill time and the survivor's position. Theorem 2 and Definition 1 extend this to N coordinates via a graph of killing paths, and applications are given to a bivariate Poisson process, a two-particle single-file diffusion, a trivariate Poisson process, and n-th-to-default CDS pricing. The examples are accompanied by Monte Carlo comparisons, and the single-file result is asserted to reproduce the known Locatelli formula.
Significance. If the general formula were correct, it would provide an exact order-statistics representation for first-passage problems in many-body Markov systems, with potential applications in physics and credit risk. The paper's strengths are its explicit analytic formulas for the Poisson and single-file examples and the Monte Carlo verification of those examples. The two-coordinate single-file result being consistent with an independent known result is a useful sanity check. However, the claimed generality for N≥3 is not established, and the main theorem is not supported by the proof supplied.
major comments (3)
- [§3, Definition 1, Eq. (44)] The path contribution I_g(t) integrates only over the final survivor position y and the kill times τ1,...,τ_{n-1}. For N≥3, the law of later kill times depends on the positions of the surviving coordinates at earlier kill times. For example, for three independent Brownian particles killed at the boundary, the second kill time distribution after the first kill is the first-passage time of a Brownian started from the conditional position at τ1, not from the original initial distribution. Eq. (44) uses F_{n+1}(τ_{n-1})...F_N(τ_1) as marginals from the original system, so it does not represent the ordered kill-time density. The one-sentence induction proof of Theorem 2 does not address this. The trivariate Poisson example is special because the independent exponential structure makes the factorization exact; it does not validate the general claim. This is a load-bearing gap for Eq. (45).
- [§2, Theorem 1, Eq. (7)] Eq. (7) states P(τ2=τm)=1−S2(t). This is not an identity: the left side is the time-independent probability that coordinate 2 is the first to be killed, while S2(t) was defined as the joint survival P(t≤τm). If S2(t) is instead read as the marginal survival of X2, the equality would still be false because P(τ2=τm) does not equal P(τ2≤t). The denominator in Eq. (8) should be P(τ2=τm), so the cancellation used to obtain Eq. (4) is invalid. The two-coordinate formula may be correct, as the examples suggest, but the proof as written is not rigorous.
- [§2.2, Eq. (37)] The statement that Eq. (37) is 'identical to the result reported in [27]' is not substantiated. Eq. (37) is an infinite series whose terms contain an integral with a ratio of infinite series 1/C(τ), whereas Eq. (40) is a closed-form expression. The agreement is only shown visually in Fig. 2. A numerical comparison of the two expressions should be provided, or the claim should be softened to 'numerically agrees with'.
minor comments (3)
- [§3.1, Eq. (50)] In the displayed computation for I(t), the second integral is written with lower limit 0 and upper limit τ, but the outer variable is t and the first integral is over τ. The second upper limit should be t.
- [§3.1, Eq. (49)] In the third displayed path, 'X0=1' should presumably be 'X1=1' to match the notation of the other paths.
- [§2, Eq. (7)] The notation S2(t) is used inconsistently: at the beginning of Section 2 it is defined as P(t≤τm), but in Eq. (7) it appears to denote a marginal survival probability. This ambiguity contributes to the proof error.
Circularity Check
No significant circularity: the recursive survival formula is a path-decomposition identity validated against independent benchmarks and Monte Carlo.
full rationale
The claimed derivation is a decomposition of the event that exactly n coordinates survive. Theorem 1 (Eq. 4) conditions on the first kill time; the object P^2_1(y,s|X2∈∂R,s)F2(s) is by definition the joint density of the kill time and survivor position, so the formula does not use the target survival function as an input. Theorem 2 (Eqs. 44–45) recursively expresses S^n through S^{n+1} plus path contributions, with base case S^N taken as given; this is a downward recursion, not a circular definition. The examples are checked against Monte Carlo, and the single-file result is explicitly identical to the independent formula of Locatelli et al. [27]. The two self-citations [21,22] provide the eigenfunction machinery for the specific single-file example, but the recursion formula itself does not reduce to those papers and is not calibrated on them. The main weakness is rigor, not circularity: the proof of Theorem 2 is a one-sentence induction sketch, and Definition 1 (Eq. 44) omits the required integrations over survivor positions at intermediate kill times, so the general N≥3 statement is underived as written. The conclusion also admits that Eq. (37) is 'infinitely slower and more complicated' than [27], which is a limitation of usefulness, not circularity of derivation. No fitted parameter is renamed as a prediction and no uniqueness or ansatz claim is imported via self-citation.
Assumptions & free parameters
free parameters (3)
- Poisson intensities (bivariate: lambda1, lambda2, lambda12; trivariate: lambdai, lambdaij) =
lambda1=1, lambda2=2, lambda12=0.8, M=5 (Fig 1); lambda1=1.2, lambda2=0.5, lambda3=3.3, lambda12=1.4, lambda13=3.1, lamb
- CDS pricing inputs r, T, payment frequency =
r=0.02, T=5, semiannual fee payments
- Single-file diffusion coefficient and box geometry =
diffusion coefficient 1, box (0,1) with reflecting boundary at 0 and killing boundary at 1
assumptions (5)
- domain assumption All transition densities (P^1_i and P^2_i), conditional densities P^2_i(x,t|Xj in dR, t), and marginal first passage densities exist and are computable analytically or numerically.
- domain assumption All hitting times exist, are finite, and are inaccessible stopping times; coordinates are killed instantaneously upon hitting the boundary.
- ad hoc to paper The density of the kill history factorizes: P^2_1(y,s|X2 in dR, s) F2(s) is the joint density of (first kill at s, survivor position y), and in general the product F_(n+1)(tau_(n-1)) ... F_N(tau_1) with P^N_n(...) is the joint density of the ordered kill history.
- standard math The two-particle single-file propagator is the Bethe-ansatz eigenfunction representation of Lapolla and Godec [22] with an ordering operator (Eqs. 25-32), and the first passage density F^2(t) comes from Lapolla [21] (Eq. 36).
- domain assumption For the single-file example the initial positions are iid uniform in the box, so the ordered particles map to independent particles and only the leftmost/rightmost ordering matters.
Cite this review
Pith. "Pith review of A recursive formula for the $n^\text{th}$ survival function and the $n^\text{th}$ first passage time distribution for jump and diffusion processes. Applications to the pricing of $n^\text{th}$-to-default CDS." pith.science (2026). https://pith.science/paper/ZOYZXMW4
@misc{pith2026250902347,
author = {Pith},
title = {Pith review of: A recursive formula for the $n^\textth$ survival function and the $n^\textth$ first passage time distribution for jump and diffusion processes. Applications to the pricing of $n^\textth$-to-default CDS},
year = {2026},
howpublished = {\url{https://pith.science/paper/ZOYZXMW4}},
note = {Machine review of arXiv:2509.02347}
}
abstract
We derive some rather general, but complicated, formulae to compute the survival function and the first passage time distribution of the $n^\text{th}$ coordinate of a many-body stochastic process in the presence of a killing barrier. First we will study the case of two coordinates and then we will generalize the results to three or more coordinates. Even if the results are difficult to implement, we will provide examples of their use applying them to a physical system, the single file diffusion, and to the financial problem of pricing a $n^\text{th}$-to-default credit default swap ($n^\text{th}$-CDS)
Figures
Figures from the paper (1 more)
Reference graph
Works this paper leans on
-
[27]
Active Brownian particles escaping a channel in single file
Emanuele Locatelli, Fulvio Baldovin, Enzo Orlandini, and Matteo Pierno. Active Brownian particles escaping a channel in single file. Phys. Rev. E , 91(2):022109, Feb- ruary 2015
work page 2015
-
[1]
Persistence Probabilities and Exponents
Frank Aurzada and Thomas Simon. Persistence Probabilities and Exponents. In L´ evy Matters V , volume 2149, pages 183–224. Springer International Publishing, Cham,
-
[2]
N. Balakrishnan. Discrete Multivariate Distributions. In Samuel Kotz, Campbell B. Read, Narayanaswamy Balakrishnan, and Brani Vidakovic, editors, Encyclopedia of Statistical Sciences. Wiley, 2 edition, December 2005
work page 2005
-
[3]
Alan J. Bray, Satya N. Majumdar, and Gr´ egory Schehr. Persistence and first-passage properties in nonequilibrium systems. Advances in Physics, 62(3):225–361, June 2013
work page 2013
-
[4]
Interest rate models: theory and practice: with smile, inflation, and credit
Damiano Brigo and Fabio Mercurio. Interest rate models: theory and practice: with smile, inflation, and credit . Springer finance. Springer, Berlin ; New York, 2nd ed edition, 2006. nth SUR VIV AL FUNCTION ANDnth FIRST PASSAGE TIME DISTRIBUTION 17
work page 2006
-
[5]
First Passage Time Distribution of a Two- Dimensional Wiener Process with Drift
Marco Domin´ e and Volkmar Pieper. First Passage Time Distribution of a Two- Dimensional Wiener Process with Drift. Prob. Eng. Inf. Sci. , 7(4):545–555, October 1993
work page 1993
-
[6]
A Simple Exponential Model for Dependent Defaults
Kay Giesecke. A Simple Exponential Model for Dependent Defaults. JFI, 13(3):74–83, December 2003
work page 2003
-
[7]
Monte Carlo methods in financial engineering
Paul Glasserman. Monte Carlo methods in financial engineering . Number 53 in Ap- plications of mathematics : stochastic modelling and applied probability. Springer, New York, NY, softcover version of original hardcover ed. 2003 edition, 2010
work page 2003
Show all 41 references
-
[8]
Gradshteyn, I. S. and Ryzhik, I. M. Table of integrals, series, and products . Else- vier/Academic Press, Amsterdam, seventh edition, 2007
2007
-
[9]
Dependent Default Modeling through Multi- variate Generalized Cox Processes, 2025
Djibril Gueye and Alejandra Quintos. Dependent Default Modeling through Multi- variate Generalized Cox Processes, 2025. Version Number: 1
2025
-
[10]
Mean-field limit of particle systems with absorp- tion, 2023
Gaoyue Guo and Milica Tomasevic. Mean-field limit of particle systems with absorp- tion, 2023. Version Number: 2
2023
-
[11]
collisions
T. E. Harris. Diffusion with “collisions” between particles. Journal of Applied Prob- ability, 2(2):323–338, December 1965
1965
-
[12]
P. Holgate. Estimation for the Bivariate Poisson Distribution. Biometrika, 51(1/2):241, June 1964
1964
-
[13]
Valuation of a CDO and an n -th to Default CDS Without Monte Carlo Simulation
John C Hull and Alan D White. Valuation of a CDO and an n -th to Default CDS Without Monte Carlo Simulation. JOD, 12(2):8–23, November 2004
2004
-
[14]
Hummer, J
G. Hummer, J. C. Rasaiah, and J. P. Noworyta. Water conduction through the hy- drophobic channel of a carbon nanotube. Nature, 414(6860):188–190, 2001
2001
-
[15]
D. W. Jepsen. Dynamics of a Simple Many-Body System of Hard Rods. J. Math. Phys., 6(3):405–413, 1965
1965
-
[16]
Analysis of sports data by using bivariate Poisson models
Dimitris Karlis and Ioannis Ntzoufras. Analysis of sports data by using bivariate Poisson models. J Royal Statistical Soc D , 52(3):381–393, October 2003
2003
-
[17]
The structure of trivariate Poisson distribution
Kazutomo Kawamura. The structure of trivariate Poisson distribution. Kodai Math. J., 28(1), January 1976
1976
-
[18]
The structure of multivariate Poisson distribution
Kazutomo Kawamura. The structure of multivariate Poisson distribution. Kodai Math. J. , 2(3), January 1979
1979
-
[19]
Kloeden and Eckhard Platen
Peter E. Kloeden and Eckhard Platen. Numerical Solution of Stochastic Differential Equations. Springer Berlin Heidelberg, Berlin, Heidelberg, 1992
1992
-
[20]
First-passage times of two-dimensional Brownian motion
Steven Kou and Haowen Zhong. First-passage times of two-dimensional Brownian motion. Adv. Appl. Probab., 48(4):1045–1060, December 2016
2016
-
[21]
The First Exit Time Statistics and the Entropic Forces in Single File Diffusion, 2022
Alessio Lapolla. The First Exit Time Statistics and the Entropic Forces in Single File Diffusion, 2022. Version Number: 1
2022
-
[22]
BetheSF: Efficient computation of the exact tagged- particle propagator in single-file systems via the Bethe eigenspectrum
Alessio Lapolla and Aljaˇ z Godec. BetheSF: Efficient computation of the exact tagged- particle propagator in single-file systems via the Bethe eigenspectrum. Comput. Phys. Commun, page 107569, August 2020
2020
-
[23]
Basket default swaps, CDOs and factor copulas
Jean-Paul Laurent and Jon Gregory. Basket default swaps, CDOs and factor copulas. JOR, 7(4):1–20, June 2005
2005
-
[24]
Sean D. Lawley. Distribution of extreme first passage times of diffusion. J. Math. Biol., 80(7):2301–2325, June 2020
2020
-
[25]
David X. Li. On Default Correlation: A Copula Function Approach. JFI, 9(4):43–54, March 2000
2000
-
[26]
Berg, and Johan Elf
Gene-Wei Li, Otto G. Berg, and Johan Elf. Effects of macromolecular crowding and DNA looping on gene regulation kinetics. Nature Phys, 5(4):294–297, April 2009
2009
-
[28]
Single-File Escape of Colloidal Particles from Microfluidic Channels
Emanuele Locatelli, Matteo Pierno, Fulvio Baldovin, Enzo Orlandini, Yizhou Tan, and Stefano Pagliara. Single-File Escape of Colloidal Particles from Microfluidic Channels. Phys. Rev. Lett. , 117(3):038001, July 2016
2016
-
[29]
Credit and basket default swaps
Dilip Madan, Michael Konikov, and Mircea Marinescu. Credit and basket default swaps. JCR, 2(1):67–87, 2006. 18 ALESSIO LAPOLLA
2006
-
[30]
Majumdar, Arnab Pal, and Gr´ egory Schehr
Satya N. Majumdar, Arnab Pal, and Gr´ egory Schehr. Extreme value statistics of cor- related random variables: A pedagogical review. Physics Reports, 840:1–32, January 2020
2020
-
[31]
Marshall and Ingram Olkin
Albert W. Marshall and Ingram Olkin. A Multivariate Exponential Distribution. Journal of the American Statistical Association , 62(317):30–44, March 1967
1967
-
[32]
Semi-Analytical Valuation of Basket Credit Derivatives in Intensity- Based Models
Allan Mortensen. Semi-Analytical Valuation of Basket Credit Derivatives in Intensity- Based Models. JOD, 13(4):8–26, May 2006
2006
-
[33]
Navarro and Ian G
Daniel J. Navarro and Ian G. Fuss. Fast and accurate calculations for first-passage times in Wiener diffusion models. Journal of Mathematical Psychology, 53(4):222–230, August 2009
2009
-
[34]
Stopping times occurring simultaneously
Philip Protter and Alejandra Quintos. Stopping times occurring simultaneously. ESAIM: PS , 28:110–131, 2024
2024
-
[35]
A guide to first-passage processes
Sidney Redner. A guide to first-passage processes. Cambridge Univ. Press, Cambridge, digitally printed version (with corrections) 2007 edition, 2007. OCLC: 830642837
2007
-
[36]
Single-file diffusion in an interval: First passage properties
Artem Ryabov. Single-file diffusion in an interval: First passage properties. J. Chem. Phys., 138(15):154104, April 2013
2013
-
[37]
Stochastic Integrate and Fire Models: A Review on Mathematical Methods and Their Applications
Laura Sacerdote and Maria Teresa Giraudo. Stochastic Integrate and Fire Models: A Review on Mathematical Methods and Their Applications. In Mostafa Bachar, Jerry Batzel, and Susanne Ditlevsen, editors, Stochastic Biomathematical Models , volume 2058, pages 99–148. Springer Ber...
-
[38]
First passage times of two-dimensional correlated processes: Analytical results for the Wiener process and a numerical method for diffusion processes
Laura Sacerdote, Massimiliano Tamborrino, and Cristina Zucca. First passage times of two-dimensional correlated processes: Analytical results for the Wiener process and a numerical method for diffusion processes. Journal of Computational and Applied Mathematics, 296:275–292, A...
2016
-
[39]
Herbert S. Wilf. What is an Answer? The American Mathematical Monthly , 89(5):289–292, May 1982
1982
-
[40]
An Analysis of Default Correlations and Multiple Defaults
Chunsheng Zhou. An Analysis of Default Correlations and Multiple Defaults. Rev. Financ. Stud. , 14(2):555–576, April 2001. Email address : alessiolapolla@gmail.com
2001
-
[2015]
Series Title: Lecture Notes in Mathematics
Reviewed August 5, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.