---
title: Polynomial Decay of Ruin Probabilities
url: https://www.emergentmind.com/topics/polynomial-rate-of-decay-of-ruin-probabilities
type: topic
---

# Polynomial Decay of Ruin Probabilities

The polynomial rate of decay of ruin probabilities describes a fundamental asymptotic behavior in risk processes where, under various structural and distributional assumptions, the probability that the reserve of an insurer falls below zero diminishes as a negative power of the initial capital. Unlike the exponential decay observed in classical light-tailed Cramér–Lundberg models, models incorporating risky investments, heavy-tailed claims, or regime-switching dynamics exhibit power-law (polynomial) decay, characterized by specific decay exponents linked to underlying moment or renewal properties. The precise rate and leading constants are determined by spectral, moment, or implicit renewal indices associated with the model’s random coefficients and claim structures.

## 1. Discrete-Time and General IID Models: Moment Indices and Ultimate Ruin

In discrete risk models that combine both financial and insurance risks, aggregate losses over $n$ periods are modeled as $Y_n = \sum_{i=1}^n A_1\cdots A_{i-1} B_i$ with $(A_i, B_i)$ IID, where $A_i$ models the annual stochastic discount (financial risk) and $B_i$ the net insurance loss. The ruin probability is $\psi(u) = P(\sup_n Y_n > u)$, the probability that cumulative (discounted) losses ever exceed initial capital $u$. The polynomial decay is quantified via so-called *moment indices*. For a real random variable $X$, define
$$I(X) = \sup\{s \geq 0 : E[(X^+)^s] < \infty\}, \quad X^+ := \max(X,0),$$
and a Lundberg-type index $I^1(X) = \sup\{s \geq 0 : E[X^s] \leq 1\} \leq I(X)$. Under the assumption that $P(\sup_n Y_n = \infty) = 1$ (e.g., $P(A > 1) > 0$, $P(A < 1) > 0$, $P(B > 0) > 0$),
$$
\limsup_{u \to \infty} \frac{\log \psi(u)}{\log u} = -\min \{I^1(A), I(B)\},
$$
i.e., $\psi(u) \sim \text{const} \times u^{-\kappa}$ for large $u$, $\kappa = \min\{I^1(A), I(B)\}$, with the constant not explicit under general assumptions [1303.0522].

For the finite-horizon probability $\psi_n(u) = P(\max_{1 \leq k \leq n} Y_k > u)$, similar asymptotics hold with an explicitly computable decay exponent in terms of joint moment indices across $k$-step sign patterns. Under mild independence or positivity (e.g., $B_i \geq 0$), the decay rate simplifies to $\min\{I(A), I(B)\}$ for $n \geq 2$. The existence and calculation of moment indices thus completely determine the polynomial decay regime in such models.

## 2. Stochastic Investment Models: Geometric Brownian Motion and Lévy Processes

When the insurer’s surplus is invested in a risky asset, often modeled via geometric Brownian motion (gBm), the ruin probability decay is governed by a critical exponent $\beta = 2a/\sigma^2 - 1$, where $a$ is the mean return, $\sigma^2$ the variance. For a surplus process
$$
dX_t = a X_t\,dt + \sigma X_t\,dW_t + c_t\,dt - \sum_{i=1}^{N_t} \xi_i, \qquad X_0 = u \geq 0,
$$
and under $\beta > 0$ ($a > \sigma^2/2$), sharp bounds hold:
$$
C_*(\beta) u^{-\beta} \leq \Psi(u) \leq C^*(\beta) u^{-\beta}, \quad u \gg 1,
$$
with $C_*, C^*$ explicit via renewal formulas. If $c_t = c^*e^{\gamma t}$ with $\gamma \leq 0$, one has the exact asymptotic $\Psi(u) \sim C(\beta) u^{-\beta}$ as $u \to \infty$. If $\beta \leq 0$, ruin becomes certain, i.e., $\Psi(u) = 1$ for all $u \geq 0$ [1011.1329].

This polynomial regime is robust to the presence of jumps, alternative sources of randomness, and varying claim arrival intensities. Generalization to reserve evolution driven by independent Lévy processes for asset returns and premium/claims yields analogous tail asymptotics. The critical exponent $\beta$ becomes the unique positive root of $H(\beta) = 0$ for
$$
H(q) = \log E[e^{-q V_1}],
$$
$V_t$ being the log-asset process. Explicit asymptotics $\Psi(u) \sim C_\infty u^{-\beta}$ hold under non-arithmeticity and finite moments, with $C_\infty$ computable via Kesten–Goldie renewal formulas [1604.06370].

## 3. Markov-Modulated and Switching Investment Models

In models where asset returns switch according to Markovian or randomly resetting regimes, the ruin probability decay exponent $\beta$ is determined by implicit renewal theory. Under a regime-switching gBm with drift/variance parameters $(a_i, \sigma_i^2)$ and transition intensity $\lambda^{ij}$, the decay exponent $\beta$ solves an equation of the form
$$
\sigma_0^2 \sigma_1^2 \beta(\beta_0 - \beta)(\beta_1 - \beta) + 2\sigma_0^2 (\beta_0 - \beta) \lambda^{10} + 2\sigma_1^2 (\beta_1 - \beta) \lambda^{01} = 0,
$$
with $\beta_i = 2 a_i/\sigma_i^2 - 1$, and the ultimate ruin probability satisfies $\psi_i(u) \sim C_i u^{-\beta}$ [2012.05083]. If drift and volatility coefficients are time-varying and reset at each claim epoch, the ruin probability decays as $u^{-\beta}$, where $\beta$ is the unique solution to $E[e^{\beta \nu}] = 1$ with $\nu$ encoding the cumulative “dilution” of investment over a random cycle [2302.11682].

These results remain valid under only mild conditions: finite higher claim moments, nondegenerate volatility in each regime, and a non-degenerate switching process.

## 4. Ruin Under Heavy-Tailed Claims and Non-Existence of Lundberg Exponents

In portfolios subjected to heavy-tailed claims (e.g., regularly varying tails with index $\alpha > 0$), even in the absence of risky investments, ruin probabilities decay only polynomially. For i.i.d. claim sizes $X_1,\ldots,X_n$ with $P(X_1 > x) = x^{-\alpha} L(x)$ and capital allocation scaling with $u_n = a_n n^\beta$, one finds
$$
\psi(u_n) \sim n^{-\alpha \beta}, \quad \text{or } \psi(u) \sim u^{-\frac{\alpha\beta}{\beta + 1/\alpha}}
$$
for large $u$ [2512.24352]. In the case of discrete-time processes with proportional reinsurance and investment, and Pareto-claims ($\alpha > 1$), ruin probability decays as $C(i_s) u^{-\alpha}$, where $C(i_s)$ solves a matrix renewal equation reflecting both reinsurance retention and interest-rate dynamics [1306.3479].

Notably, in heavy-tailed models, the adjustment (Lundberg) coefficient does not exist, and the decay exponent is dictated by the tail index of the claim size distribution and systemically by the risk-sharing or investment mechanism.

## 5. Integro-Differential Equations and Polynomial Tails in Risky Investment/Annuity Models

In surplus models with mixed risky and riskless investment and Cramér–Lundberg-type jumps, the survival probability $\Phi(x) = 1 - \psi(x)$ satisfies a second-order integro-differential equation:
$$
\frac{1}{2}\kappa^2 \sigma^2 x^2 \Phi''(x) + \left((a-r)\kappa + r\right)x \Phi'(x) - c \Phi'(x) - \lambda \int_0^\infty [\Phi(x+y) - \Phi(x)] F(dy) = 0.
$$
On integrating and reducing to a Volterra equation, the solution’s derivative exhibits $g(x) \sim x^{-\gamma}$, leading to
$$
\psi(x) \sim C x^{-(\gamma - 1)} \quad \text{as } x \to \infty,
$$
with $\gamma = \frac{2((a-r)\kappa + r)}{\kappa^2\sigma^2}$ and $p = \gamma - 1 > 0$ [2601.01447]. The leading constant is given explicitly in terms of the integrated solution. Such results confirm that the presence of risky investments universally replaces the classical exponential decay by a polynomial law, even when the jump mechanism admits only minimal regularity.

## 6. Polynomial Approximation Rates in Scaled Classical Models

When analyzing the rate at which the ruin probability in the scaled classical Cramér–Lundberg risk process converges to its diffusion approximation, the error decays polynomially with respect to the scaling parameter $n$. Specifically, for the scaled process with claim arrival intensity $n \lambda$, claim size $Y/\sqrt{n}$, and initial capital $x$,
$$
|\psi_n(x) - \psi_D(x)| \leq C n^{-1/2}
$$
uniformly in $x \geq 0$, where $\psi_D(x)$ is the diffusion-limit solution [1902.00706]. In the exponential claim case, higher-order expansions yield
$$
\psi_n(x) = e^{-\gamma x} + a_1 n^{-1/2} e^{-\gamma x} + \ldots + a_k n^{-k/2} e^{-\gamma x} + O(n^{-(k+1)/2}).
$$
These polynomial rates describe the asymptotic convergence speed to the continuous approximation and highlight the robustness of polynomial error bounds in risk models.

## 7. Structural Mechanisms and Universality of Power-Law Decay

The universal mechanism underlying polynomial decay of ruin probabilities is the emergence of random recurrent affine equations for the surplus process, either at claim epochs or through embedded Markov chains. Renewal and implicit renewal theory (notably Kesten–Goldie-type results) provide the analytic foundation, linking decay exponents to roots of spectral equations for random multipliers $M$ (often $E[M^\beta] = 1$). The value of the polynomial exponent is sensitive to the balance between mean investment return and volatility, as well as the heaviness of claim tails. Risky investments, time-heterogeneous regime-switching, and heavy-tail phenomena each suppress exponential rates and replace them with explicit power laws, determined by moment or spectral criteria across all models considered [1011.1329, 2012.05083, 1604.06370, 2601.01447, 2512.24352, 1303.0522].

| Model class                 | Ruin probability decay                  | Exponent formula / regime                    |
|-----------------------------|-----------------------------------------|----------------------------------------------|
| Discrete risk model         | $\psi(u) \sim \text{const}\,u^{-\kappa}$ | $\kappa = \min\{I^1(A), I(B)\}$             |
| Risky investment (gBm)      | $\Psi(u) \sim C u^{-\beta}$             | $\beta = 2a/\sigma^2 - 1$                   |
| Lévy-driven asset returns   | $\Psi(u) \sim C_\infty u^{-\beta}$      | $H(\beta) = 0$, $H(q) = \log E[e^{-q V_1}]$ |
| Markov-modulated returns    | $\psi_i(u) \sim C_i u^{-\beta}$         | $\beta$ solves $f(\beta) = 1$ (see above)   |
| Heavy-tailed claims         | $\psi(u) \sim C u^{-\nu}$               | $\nu$ from tail/renewal equations           |

## References

- "Asymptotic behaviour of ruin probabilities in a general discrete risk model using moment indices" [1303.0522]
- "Ruin probability in the presence of risky investments" [1011.1329]
- "The ruin problem for Lévy-driven linear stochastic equations with applications to actuarial models with negative risk sums" [1604.06370]
- "On ruin probabilities with risky investments" [2012.05083]
- "On ruin probabilities in the presence of risky investments and random switching" [2302.11682]
- "Heavy-tailed distributions; extreme value theory; large deviations; ruin probabilities; solvency risk" [2512.24352]
- "Ruin probability of a discrete-time risk process with proportional reinsurance and investment for exponential and Pareto distributions" [1306.3479]
- "On the integro-differential equation arising in the ruin problem for annuity payment models" [2601.01447]
- "Rate of Convergence of the Probability of Ruin in the Cramér-Lundberg Model to its Diffusion Approximation" [1902.00706]

Source: https://www.emergentmind.com/topics/polynomial-rate-of-decay-of-ruin-probabilities