Credit Default Swaps (CDS) and Credit Derivatives › Default Probabilities and Hazard Rates
Quantitative FinanceAll chapters Credit Default Swaps (CDS) and Credit Derivatives › Default Probabilities and Hazard RatesOne section · Sections in this chapter5 / 16 ↓
Reading tools Definitions & inline links
Choose your expertise level to adjust how many terms are explained. Beginners see more tooltips, experts see fewer to maintain reading flow. Hover over underlined terms for instant definitions.
Article links
Make inline references clickable
To price CDS contracts rigorously, we need a framework for modeling default. The standard approach uses hazard rates (also called default intensities) to characterize the instantaneous probability of default. This framework is the foundation for credit derivatives pricing, allowing us to move from intuitive concepts to precise values.
Hazard Rate Model
The hazard rate represents the instantaneous conditional probability of default at time , given survival up to that point. This concept may seem abstract at first, but it captures something fundamentally important: the intensity of default risk at each moment, measured as a rate rather than a probability. For a small time interval :
where:
- : instantaneous hazard rate (default intensity) at time
- : small time interval
- : current time
Think of the hazard rate as measuring the "default pressure" at each instant. A hazard rate of 0.02 means that, conditional on having survived until now, the probability of defaulting in the next small time interval is approximately . If is one year, that corresponds to a 2% conditional probability; if is one day (roughly 1/365 of a year), the conditional probability is approximately 0.0055%.
The survival probability , which is the probability of no default by time , is related to the hazard rate by:
where:
- : probability of surviving (no default) until time
- : hazard rate at time
- : exponential function
- : time horizon
This formula emerges from the mathematical requirement that survival requires avoiding default at every instant from time 0 to time . The integral accumulates the hazard rate over time, and the exponential function converts this accumulated hazard into a probability. The negative sign in the exponent ensures that higher cumulative hazard translates to lower survival probability, as intuition demands.
For a constant hazard rate , this simplifies to:
where:
- : survival probability to time
- : constant hazard rate
- : time in years
This exponential survival function is analogous to radioactive decay: at any moment, there is a constant probability rate of default, regardless of how long the entity has survived. Just as a radioactive atom has no memory of how long it has existed, the constant hazard rate model assumes that a company's instantaneous default risk depends only on its current state, not on how long it has been solvent. While this "memoryless" property is a simplification of reality, it provides tractable mathematics and is a reasonable first approximation for many applications.
import numpy as np
# Model survival probability with constant hazard rate
hazard_rates = [
0.01,
0.02,
0.05,
0.10,
] # 1%, 2%, 5%, 10% annual default intensity
times = np.linspace(0, 10, 100)
def survival_probability(t, hazard_rate):
"""Calculate survival probability under constant hazard rate model."""
return np.exp(-hazard_rate * t)
Cumulative Default Probability
The cumulative default probability by time is simply one minus the survival probability. This relationship follows directly from the fact that an entity must either survive or default; there is no third possibility.
where:
- : survival probability
- : constant hazard rate
- : time horizon
For short time horizons and small hazard rates, this approximates to:
where:
- : hazard rate
- : time horizon
This linear approximation works because the exponential function when is small. For a hazard rate of 2% and a one-year horizon, the exact default probability is , while the approximation gives . The difference is negligible for practical purposes. However, as time extends or hazard rates increase, the approximation breaks down because it fails to account for the compounding effect: the longer the time horizon, the more the cumulative default probability accelerates relative to the linear approximation.
import numpy as np
# Calculate cumulative default probabilities for a hazard rate of 2%
hazard = 0.02
horizons = [1, 3, 5, 10]
# Calculate probabilities for different horizons
results = []
for t in horizons:
exact_prob = 1 - np.exp(-hazard * t)
approx_prob = hazard * t
results.append((t, exact_prob, approx_prob))Cumulative Default Probabilities (λ = 2%): --------------------------------------------- 1-year: 0.0198 (exact), 0.0200 (approximation) 3-year: 0.0582 (exact), 0.0600 (approximation) 5-year: 0.0952 (exact), 0.1000 (approximation) 10-year: 0.1813 (exact), 0.2000 (approximation)
The approximation works well for short horizons but diverges over longer periods where compounding effects matter.
