Probability Distributions in Finance: Normal, Lognormal

Probability distributions are a core modeling tool in quantitative finance. Many portfolio-risk estimates, option-pricing models, and default-likelihood models require assumptions about how uncertain outcomes are distributed. A misspecified distribution can materially change estimated tail probabilities, option values, and risk limits.

The choice becomes especially visible in the tails. Later in this chapter, the same simulated heavy-tailed sample produces a lower 95% VaR but higher expected shortfall and higher 99% risk measures than a fitted normal model. That mixed result is the practical reason to inspect a distribution's assumptions rather than relying on a blanket rule about normal models.

This chapter covers several common models in quantitative finance: the normal distribution as a baseline for returns and mean-variance analysis, the lognormal distribution for positive quantities such as modeled asset prices, discrete distributions like the binomial and Poisson for counts and event models, and heavier-tailed alternatives. Whether any of them fits a particular dataset or decision well must be checked against the relevant horizon and evidence. For each distribution, we examine both the mathematical properties and the financial intuition behind its use.

The normal distribution, also called the Gaussian distribution, is a common baseline in financial modeling. One mathematical reason is the classical Central Limit Theorem: standardized sums of independent, identically distributed random variables with finite variance converge in distribution to a normal law. Thinking of a return as the aggregate effect of many small shocks can motivate a normal approximation, but market shocks need not be independent or identically distributed, so this is a modeling heuristic rather than a theorem about returns.

To understand what it means for a random variable to follow a normal distribution, we must first grasp the concept of a probability density function. Unlike discrete random variables that take on specific values with certain probabilities, a continuous random variable like one following a normal distribution can take any value along the real line. Probability over a region is the integral of the density over that region; density at one point is not itself a probability. A higher density at a point means that sufficiently small, equal-width neighborhoods around that point have greater probability than neighborhoods around a lower-density point.

A random variable $X$ follows a normal distribution with mean $\mu$ and variance $\sigma^2$, written $X \sim N(\mu, \sigma^2)$, if its probability density function is:

Each component serves a specific purpose. Let's examine each component:

• $x$: the value at which we evaluate the density
• $\mu$: the mean (center of the distribution)
• $\sigma$: the standard deviation (controls the spread)
• $\frac{1}{\sigma\sqrt{2\pi}}$: a normalization constant ensuring the density integrates to 1
• $(x-\mu)^2/(2\sigma^2)$: one half of the squared standardized distance from $x$ to the mean

The heart of the normal distribution lies in the exponential term. The quantity $(x-\mu)^2$ measures the squared distance between the point $x$ and the center of the distribution $\mu$. Dividing by $2\sigma^2$ scales this distance relative to the spread of the distribution. A value two units away from the mean is much more unusual when $\sigma = 0.5$ than when $\sigma = 10$, because the standardized distance is larger. The negative sign in the exponent ensures that the function decreases as we move away from the mean, and the exponential function turns this squared distance into a probability density that decays smoothly and symmetrically.

The exponential of a negative quadratic creates the characteristic bell shape: values near the mean have high probability density, while values far from the mean have exponentially decreasing density. For a standard normal variable, the two-sided probability beyond four standard deviations is about $6.3 \times 10^{-5}$, and beyond five it is about $5.7 \times 10^{-7}$. Whether those benchmark probabilities adequately describe a financial return series is an empirical question tied to the asset, horizon, and sample. The quadratic term in the exponent is the key: because squaring grows faster than linear growth, the tails of the normal distribution thin out quickly.

The normalization constant $\frac{1}{\sigma\sqrt{2\pi}}$ deserves special attention. This factor ensures that when we integrate the density function over all possible values from negative infinity to positive infinity, we obtain exactly 1. This is the total probability. The $\sqrt{2\pi}$ arises from the Gaussian integral, one of the most famous results in calculus, while the $\sigma$ in the denominator accounts for how wider distributions (larger $\sigma$) must have lower peak heights to maintain the same total area under the curve.

Key parameters:

• Mean ($\mu$): The center of the distribution and the expected value of $X$. In finance, this represents the expected return.
• Variance ($\sigma^2$): Measures the spread of the distribution. The square root, $\sigma$, is the standard deviation, which in finance we call volatility.
• Skewness: Zero for the normal distribution, meaning it is perfectly symmetric around the mean.
• Kurtosis: Equal to 3 for the normal distribution (or excess kurtosis of 0). Kurtosis is the fourth standardized moment; by itself, it does not uniquely determine how probability is divided among the center, shoulders, and tails.

The symmetry of the normal distribution is both its elegance and its limitation. Symmetry is a modeling assumption: it assigns equal probabilities to equal-sized deviations above and below the mean. Some financial return series exhibit asymmetry, with the direction and magnitude depending on the asset, sampling horizon, and market regime. This possibility is one reason to treat the normal distribution as a baseline rather than a universal model of returns.

The standard normal distribution is the special case where $\mu = 0$ and $\sigma^2 = 1$. Any normal variable can be standardized:

This standardization process turns any normal random variable into the canonical standard normal form. The operation $\frac{X - \mu}{\sigma}$ first centers the distribution by subtracting the mean (so the new center is at zero), then scales it by dividing by the standard deviation (so the new spread equals one). This step preserves the essential shape of the distribution while placing it in a standard coordinate system.

This standardization is central because it allows us to use a single table of probabilities for any normal distribution. The cumulative distribution function of the standard normal, denoted $\Phi(z)$, gives the probability that a standard normal variable is less than $z$. Before computers, statisticians relied on printed tables of $\Phi(z)$ values, and the standardization formula allowed them to answer questions about any normal distribution using this single reference. Even today, the standard normal is the universal reference point: when we speak of a "3-sigma event," we mean an outcome more than three standard deviations from the mean, regardless of what the actual mean and standard deviation are.

One of the most useful properties of the normal distribution is how probability mass is concentrated around the mean:

• Approximately 68% of observations fall within one standard deviation of the mean
• Approximately 95% fall within two standard deviations
• Approximately 99.7% fall within three standard deviations

These percentages emerge directly from integrating the normal density function, but their approximate values (68, 95, and 99.7) are worth committing to memory because they provide instant intuition for assessing how unusual any particular observation is. If someone tells you that an event was a "two-sigma move," you immediately know that such events occur roughly 5% of the time, or about one trading day in twenty.

This rule provides quick mental estimates for risk. If daily stock returns have a standard deviation of 1%, then under normality assumptions, the probability of an absolute move larger than 3% is about 0.27%. That is about once every 370 trading days, or roughly once every 1.5 trading years. This calculation reveals both the power and the danger of the normal assumption: it gives us precise quantitative predictions, but those predictions can be spectacularly wrong if the true distribution differs from normal.

The normal distribution appears throughout quantitative finance:

Portfolio theoryLink Copied

Portfolio theory often summarizes portfolios by expected return and variance. Normality is sufficient for those two moments to determine the full distribution of portfolio returns, but the mean-variance framework itself does not require normally distributed returns. When candidate portfolio returns are normal, expected utility for a fixed utility function can be evaluated from the mean and variance because those parameters determine the distribution.

Value at Risk (VaR)Link Copied

VaR calculations often assume normal returns. Under that model, the signed VaR loss quantile at confidence level $\alpha$ is:

This formula translates the abstract concept of a tail quantile into a concrete loss estimate. Under the assumption that returns are normally distributed with mean $\mu$ and standard deviation $\sigma$, we seek the threshold such that returns fall below this level only $(1-\alpha)$ percent of the time. The standard normal quantile $z_\alpha$ tells us how many standard deviations below the mean this threshold lies.

Each component of this formula plays a specific role:

• $\text{VaR}_\alpha$: the signed loss quantile at confidence level $\alpha$; under a continuous model, losses exceed this threshold with probability $1-\alpha$
• $\mu$: the expected return over the time horizon
• $\sigma$: the volatility (standard deviation of returns) over the time horizon
• $z_\alpha$: the $\alpha$-quantile of the standard normal distribution (e.g., $z_{0.95} \approx 1.645$ for 95% confidence)

The formula works because under normality, the $(1-\alpha)$-percentile of returns is exactly $\mu - z_\alpha \sigma$. For example, the 95% VaR uses $z_{0.95} \approx 1.645$, meaning we expect losses to exceed this threshold only 5% of the time. Taking the negative of this return quantile yields the signed loss quantile $z_\alpha \sigma - \mu$. It is positive when $z_\alpha\sigma>\mu$, zero at equality, and negative when the modeled lower-tail return quantile is still a gain. A reporting policy that disallows negative VaR can use $\max(0, z_\alpha\sigma-\mu)$, but that floor is a reporting convention rather than the raw quantile.

These VaR figures are loss thresholds, not maximum or expected losses. The 1-day 95% VaR of approximately 2.4% means that, under the fitted normal model, losses exceed 2.4% with 5% probability. The 10-day VaR is higher here because the volatility term grows with $\sqrt{n}$ and outweighs the mean term, which grows with $n$. Under the iid normal model, $\text{VaR}_n=z_\alpha\sigma\sqrt{n}-n\mu$, so VaR need not increase with the horizon for every parameter choice. The square-root scaling of the volatility term follows from variance addition: if daily variances are $\sigma^2$, then the variance over $n$ days is $n\sigma^2$, and the standard deviation is $\sqrt{n}\sigma$.

Returns may be approximately normal, but a nonnegative stock-price level cannot follow an exact nondegenerate normal law because that law assigns some probability below zero. A local normal approximation may still be useful. A stock price can rise 50% but cannot fall more than 100%. This asymmetry is captured by the lognormal distribution, which constrains values to be positive. The lognormal distribution arises when multiplicative returns produce an accumulated log return that is normally distributed: exponentiating that normal log return gives a positive, lognormally distributed price. Multiplicative compounding by itself does not determine the distribution.

The connection between normal returns and lognormal prices is one of the most important relationships in quantitative finance. If log-returns are normally distributed, then prices are lognormally distributed. Consider a stock with initial price $S_0$ and continuously compounded return $r$ over some period:

This equation captures a key relationship in how prices evolve. The exponential function $e^r$ converts the additive world of log-returns into the multiplicative world of price changes. When $r = 0.10$, the stock has grown by a factor of $e^{0.10} \approx 1.105$, or about 10.5%. The slight difference between the log-return and the percentage return becomes more pronounced for larger moves.

The components of this price equation are:

• $S_T$: stock price at time $T$
• $S_0$: initial stock price
• $r$: continuously compounded return over the period
• $e^r$: the growth factor, converting log-return to price ratio

If $r \sim N(\mu, \sigma^2)$, then $S_T$ follows a lognormal distribution. This relationship is fundamental: the exponential function turns the symmetric, unbounded normal distribution into a right-skewed distribution bounded below by zero, which matches the behavior of asset prices that can rise without limit but cannot fall below zero. The step preserves the elegant properties of the normal distribution while mapping them onto the economically sensible domain of positive prices.

The probability density function of the lognormal distribution takes a form that reveals its connection to the normal:

Compare this to the normal density to see the main differences. The argument of the exponential now involves $\ln x$ rather than $x$, which reflects that we are measuring distance in the log-space. The presence of $1/x$ in front of the exponential arises from the chain rule when transforming from the normal to the lognormal through the change of variables.

Each element of this formula serves a distinct purpose:

• $x$: the value at which we evaluate the density (must be positive)
• $\mu$: the mean of $\ln(X)$, not the mean of $X$ itself
• $\sigma$: the standard deviation of $\ln(X)$, not of $X$ itself
• $1/x$: this factor (compared to the normal PDF) accounts for the change of variables from $\ln(x)$ to $x$

The key insight is that $\mu$ and $\sigma$ are parameters of the underlying normal distribution of log-values, not the lognormal variable itself. This distinction often causes confusion: when we say a stock has lognormal volatility of 20%, we mean the standard deviation of its log-returns, not of the returns themselves. The actual mean and variance of a lognormal variable are:

These formulas reveal how the parameters of the underlying normal distribution translate into the moments of the lognormal:

• $E[X]$: the expected value (mean) of the lognormal variable $X$
• $\text{Var}(X)$: the variance of $X$
• $\mu$: the mean of $\ln(X)$, the underlying normal distribution
• $\sigma^2$: the variance of $\ln(X)$
• $e^{\mu + \sigma^2/2}$: the mean is greater than $e^\mu$ when $\sigma > 0$, and equal to it in the degenerate case $\sigma = 0$

The $\sigma^2/2$ term in the mean formula is a convexity adjustment: because the exponential function is convex, Jensen's inequality tells us that $E[e^Y] > e^{E[Y]}$ for any non-constant random variable $Y$. Specifically, the curvature of the exponential means that high values of $Y$ contribute disproportionately to the mean of $e^Y$. This is why the expected value of a lognormal variable includes a variance correction rather than equaling $e^\mu$, since higher variance increases this convexity effect. Holding the mean log-return $\mu$ fixed, greater $\sigma$ therefore raises $E[X]$. In the GBM notation below, however, $\mu$ denotes arithmetic drift: holding that drift fixed, the log-mean falls by $\sigma^2T/2$ as variance rises, and $E[S_T]=S_0e^{\mu T}$ does not depend on $\sigma$.

A standard baseline model for stock price dynamics is geometric Brownian motion (GBM), which produces lognormally distributed prices. This model has been central to quantitative finance since Samuelson's work in the 1960s and forms the foundation for the Black-Scholes option pricing formula. Understanding GBM requires familiarity with stochastic calculus, but the core idea is intuitive. Prices evolve continuously with both a deterministic trend and random fluctuations proportional to the current price.

Under GBM, the stock price evolves according to:

This stochastic differential equation describes infinitesimal changes in the stock price. Each term has a clear interpretation:

• $dS$: infinitesimal change in stock price
• $S$: current stock price
• $\mu$: drift rate (expected return per unit time)
• $\sigma$: diffusion volatility (annualized when time is measured in years); over an interval $\Delta t$, the log-return standard deviation is $\sigma\sqrt{\Delta t}$
• $dt$: infinitesimal time increment
• $dW$: increment of a Wiener process (Brownian motion), with $dW \sim N(0, dt)$

The first term $\mu S \, dt$ represents the deterministic trend, while the second term $\sigma S \, dW$ captures random fluctuations proportional to the current price. The key feature is that both terms are proportional to $S$: a $100 stock and a $10 stock with the same parameters will have the same percentage return distribution, but the absolute dollar changes will be ten times larger for the more expensive stock. This proportionality is what makes GBM appropriate for modeling asset prices, which tend to grow (or shrink) in percentage terms.

Solving this stochastic differential equation yields:

This solution tells us exactly where the stock price will be at any future time, given the realized path of the Brownian motion. The formula has a clear structure:

• $S_T$: stock price at time $T$
• $S_0$: initial stock price at time 0
• $\mu$: drift rate (annualized expected return)
• $\sigma$: volatility (annualized standard deviation)
• $T$: time horizon in years
• $W_T$: value of the Wiener process at time $T$, with $W_T \sim N(0, T)$
• $\sigma^2/2$: Itô correction term arising from the quadratic variation of Brownian motion

The Itô correction $\sigma^2/2$ appears because we apply Itô's lemma to $\ln(S)$: the second derivative term $\frac{1}{2}\frac{\partial^2 \ln S}{\partial S^2}(\sigma S)^2 = -\frac{\sigma^2}{2}$ contributes to the drift. This correction is one of the key differences between ordinary calculus and stochastic calculus. In ordinary calculus, if we know $dS/S$, we would conclude that $d(\ln S) = dS/S$. However, when $S$ follows a diffusion process, the additional $-\sigma^2/2$ term arises from the non-zero quadratic variation of Brownian motion. For $\sigma>0$, the expected GBM log-growth rate $\mu-\sigma^2/2$ is below its arithmetic drift $\mu$; at $\sigma=0$, the two are equal.

This model forms the foundation of the Black-Scholes option pricing formula.

Monte Carlo simulation of stock prices under GBM is simple. The approach below draws exact GBM transition increments at each point on the chosen time grid and cumulatively sums them. Unlike an Euler approximation, this produces the exact GBM distribution at the grid points for constant parameters; it does not simulate between-grid extrema.

Not all financial phenomena are continuous. Events like defaults, dividend announcements, or large jumps in prices are discrete occurrences. The binomial and Poisson distributions model these situations and give a framework for counting how many times something happens rather than measuring continuous quantities. These distributions are particularly important in credit risk, where we count defaults in a portfolio, and in derivatives pricing, where discrete lattice models approximate continuous price processes.

The binomial distribution models the number of successes in a fixed number of independent trials, each with the same probability of success. The classic example is flipping a coin multiple times and counting heads, but the same mathematics applies to questions such as: How many stocks in a portfolio will outperform? How many corporate bonds will default? How many days this month will have positive returns?

If $X \sim \text{Binomial}(n, p)$, the probability mass function is:

This formula answers the question: if we have $n$ independent trials, each with success probability $p$, what is the probability of getting exactly $k$ successes? The answer decomposes into two multiplicative components: the number of ways the successes can be arranged, and the probability of any particular arrangement.

Each component of this formula serves a specific purpose:

• $k$: the number of successes we're computing the probability for
• $n$: total number of independent trials
• $p$: probability of success on each trial
• $\binom{n}{k} = \frac{n!}{k!(n-k)!}$: the binomial coefficient, representing the number of ways to choose $k$ successes among $n$ trials
• $p^k$: success-probability factor for $k$ specified success positions
• $(1-p)^{n-k}$: failure-probability factor for $n-k$ specified failure positions

The formula multiplies the number of arrangements by the probability of each specific arrangement. Intuitively, if you flip a coin $n$ times with probability $p$ of heads, the PMF tells you the probability of getting exactly $k$ heads: you need $k$ successes (probability $p^k$), $n-k$ failures (probability $(1-p)^{n-k}$), and the binomial coefficient counts how many orderings produce that outcome. The key insight is that because the trials are independent, the probability of any specific sequence of outcomes is simply the product of the individual probabilities, and the binomial coefficient tells us how many such sequences lead to exactly $k$ successes.

The mean and variance have elegant formulas that reveal the distribution's structure:

These formulas encode intuitive relationships:

• $E[X]$: the expected number of successes
• $\text{Var}(X)$: the variance in the number of successes
• $n$: total number of trials
• $p$: probability of success on each trial
• $(1-p)$: probability of failure, which appears in variance because maximum variance occurs when $p = 0.5$ (maximum uncertainty)

The expected number of successes, $np$, is simply the number of trials times the success probability, exactly what intuition suggests. The variance formula $np(1-p)$ is more subtle: it shows that variance is maximized when $p = 0.5$ (complete uncertainty about each trial's outcome) and decreases toward zero as $p$ approaches 0 or 1 (near certainty). This makes sense because if $p$ is very close to 1, almost all trials succeed, leaving little room for variability in the total count.

The Cox-Ross-Rubinstein (CRR) binomial model prices options by assuming that at each time step, the stock price either moves up by a factor $u$ or down by a factor $d$. This discretization turns the continuous option pricing problem into a tree of possible outcomes that can be analyzed one step at a time. While this may seem like a crude approximation to the continuous world of actual markets, the model provides deep insights into option pricing and converges to the continuous Black-Scholes formula as the number of time steps increases.

After $n$ steps, the stock price has experienced $k$ up moves and $n-k$ down moves, giving:

This formula shows how the final price depends entirely on the number of up moves, not their specific timing. The components are:

• $S_T$: stock price at maturity after $n$ steps
• $S_0$: initial stock price
• $u$: up factor (price multiplier for an up move)
• $d$: down factor (price multiplier for a down move)
• $k$: number of up moves
• $n-k$: number of down moves

Whether the stock rises first and then falls, or falls first and then rises, the final price is the same as long as the total number of up and down moves matches. This recombination supports efficient pricing for European options, whose payoff depends only on the final price. Path-dependent claims such as Asian options require the tree state to be augmented with path information; the simple terminal-price state is not sufficient.

The number of up moves follows a binomial distribution under the risk-neutral measure, with probability $q$ of an up move:

This risk-neutral probability is central to derivatives pricing. Notice that $q$ is not the actual probability that the stock will go up. It is the artificial probability that makes the expected return on the stock equal to the risk-free rate. The components are:

• $q$: risk-neutral probability of an up move
• $r$: risk-free interest rate (annualized)
• $\Delta t$: length of each time step in years
• $u$: up factor (stock price multiplier for an up move)
• $d$: down factor (stock price multiplier for a down move)
• $e^{r\Delta t}$: the growth factor for a risk-free investment over one time step

This formula is derived from the no-arbitrage condition: the expected stock return under the risk-neutral measure must equal the risk-free rate. Setting $q \cdot u + (1-q) \cdot d = e^{r\Delta t}$ and solving for $q$ yields this result. In this complete one-period binomial model, a replicable option is priced from the risk-neutral probability rather than the physical up-move probability. This conclusion depends on no arbitrage, market completeness, and the ability to form the replicating portfolio; it is not a universal statement about every derivative market.

Assuming $d<u$, this $q$ is a strict probability and the one-period stock-and-bond model is arbitrage-free only when

which is equivalent to $0<q<1$. Inputs outside this range do not define valid CRR risk-neutral probabilities.

The Poisson distribution models the number of events occurring in a fixed interval when events happen at a constant average rate and independently of each other. Unlike the binomial distribution, which assumes a fixed number of trials, the Poisson distribution allows for any number of events, making it useful for modeling rare occurrences over a continuous interval.

If $X \sim \text{Poisson}(\lambda)$, the probability mass function is:

This formula has a clear structure when we trace its derivation from first principles. Imagine dividing the interval into many tiny subintervals, each with a small probability of containing an event. As the number of subintervals goes to infinity while keeping the expected total number of events fixed at $\lambda$, we obtain the Poisson distribution.

The components of this formula are:

• $k$: the number of events we're computing the probability for
• $\lambda$: the average rate of events (expected number of events in the interval)
• $e^{-\lambda}$: the probability of zero events, which is a normalization factor
• $\lambda^k$: the intensity factor for $k$ events; together with $e^{-\lambda}$, it makes $P(X=k)$ largest at $\lambda=k$ for $k>0$
• $k!$: accounts for the indistinguishability of event orderings

The Poisson distribution arises as the limit of a binomial distribution when $n \to \infty$ and $p \to 0$ while $np = \lambda$ remains constant. More generally, a Poisson approximation requires independent rare-event opportunities whose largest individual probability vanishes while the sum of their probabilities approaches $\lambda$. A finite expected count alone is not enough: dependence or a non-negligible individual probability can produce a different count distribution. Credit portfolios with small independent default probabilities and jump models with independent small-interval arrivals can satisfy the approximation; clustered defaults and jumps require a different model.

Both the mean and variance equal $\lambda$:

This property deserves emphasis:

• $E[X]$: the expected number of events
• $\text{Var}(X)$: the variance in the number of events
• $\lambda$: the average rate parameter

The equality of mean and variance is a distinctive property of the Poisson distribution. When analyzing count data, if the observed variance significantly exceeds the mean (overdispersion), this suggests the Poisson model may be inadequate. Dependence is one possible cause, but rate heterogeneity, omitted covariates, and excess zeros can also produce overdispersion. For instance, defaults may cluster because borrowers share economic conditions even when individual default mechanisms differ.

The Poisson distribution is particularly useful for modeling:

• Jump processes in asset prices: The Merton jump-diffusion model adds jumps whose arrivals follow a Poisson process and whose sizes follow a separate distribution to geometric Brownian motion. The stock price follows GBM between jumps and moves discontinuously at jump times. This combination represents continuous diffusion-like fluctuations together with occasional discontinuous price moves.

• Credit defaults in a portfolio: If defaults occur independently with a small probability, the number of defaults in a large portfolio is approximately Poisson distributed. This approximation underlies many credit risk models, although the assumption of independence often breaks down during systemic crises, when many borrowers default simultaneously.

• Order arrivals: In market microstructure, the arrival of orders is often modeled as a Poisson process. This assumption implies that the time between consecutive orders follows an exponential distribution, without clustering or periodicity in order flow. While real markets show more complex arrival patterns, the Poisson model provides a useful baseline.

A normal distribution is only one possible model for returns. When an observed or simulated return distribution places more probability in sufficiently remote tails than a fitted normal model, normal tail calculations can misstate extreme-event frequencies and downstream risk measures. The examples below use a seeded Student's t sample to show how to diagnose that mismatch without presenting synthetic observations as market data.

The term "fat tails" refers to distributions where extreme events occur more frequently than a normal distribution would suggest. A distribution with fat tails places more probability mass in the regions far from the mean, meaning that large positive or negative values are more common than the normal distribution's rapid exponential decay would predict. Several ways to measure this phenomenon:

Kurtosis: The normal distribution has a kurtosis of 3 (or excess kurtosis of 0). Kurtosis is the fourth central moment divided by the variance squared. A higher value signals more weight on large standardized deviations overall, but it does not by itself order every tail probability or identify where probability moved between the center and shoulders.

Tail exponents: For power-law-tailed distributions, the probability of extreme events decays polynomially rather than exponentially. The tail exponent $\alpha$ characterizes how fast this decay occurs. A smaller $\alpha$ means slower polynomial decay and fatter tails, while a larger $\alpha$ means faster polynomial decay and thinner tails within the power-law family. For every finite $\alpha$, however, the asymptotic tail remains heavier than a normal tail.

The seeded sample was generated from a heavy-tailed distribution, so its departure from the normal benchmark is expected rather than empirical evidence about markets. For this seed, absolute moves beyond four sample standard deviations occur 30 times, compared with about 0.38 expected in 6,000 independent normal draws, a ratio of about 79. Nine observations exceed five sample standard deviations, compared with about 0.0034 under that same normal benchmark. The exercise demonstrates how a known heavy-tailed data-generating process changes extreme-event counts; it does not estimate their frequency in an observed market.

The Student's t distribution provides a simple way to model fat tails while retaining many convenient properties of the normal distribution. William Sealy Gosset originally developed the t distribution (writing under the pseudonym "Student") for quality control at Guinness brewery. It is now widely used in finance for modeling the heavy tails observed in return data.

The probability density function of the standard t distribution with $\nu$ degrees of freedom is:

This formula may look complex, but its structure shows the key difference between the t distribution and the normal distribution. While the normal distribution uses the exponential function $e^{-x^2/2}$ which decays very rapidly, the t distribution uses a power function $(1 + x^2/\nu)^{-(\nu+1)/2}$ which decays much more slowly. This slower decay is precisely what produces the heavier tails.

Each component:

• $x$: the value at which we evaluate the density
• $\nu$: degrees of freedom parameter (controls tail thickness; lower $\nu$ means fatter tails)
• $\Gamma(\cdot)$: the gamma function, generalizing the factorial to non-integer values ($\Gamma(n) = (n-1)!$ for positive integers)
• $\frac{\Gamma\left(\frac{\nu+1}{2}\right)}{\sqrt{\nu\pi}\,\Gamma\left(\frac{\nu}{2}\right)}$: a normalization constant ensuring the density integrates to 1
• $(1 + x^2/\nu)^{-(\nu+1)/2}$: the core term that produces polynomial rather than exponential tail decay

The polynomial decay $(1 + x^2/\nu)^{-(\nu+1)/2}$ versus the normal's exponential decay $e^{-x^2/2}$ is why the t distribution has heavier tails. To understand this intuitively, consider what happens as $x$ becomes large. For the normal distribution, doubling $x$ quadruples the magnitude of the exponent because $(2x)^2 = 4x^2$, so the density decays extremely rapidly. For the t distribution, the density is asymptotically proportional to $|x|^{-(\nu+1)}$, so doubling a large threshold reduces the density by about $2^{-(\nu+1)}$. Its one-sided survival probability is proportional to $x^{-\nu}$, so doubling a large threshold reduces that probability by about $2^{-\nu}$. As $\nu \to \infty$, the polynomial term converges to the exponential, and the t distribution approaches the normal.

Key properties of the t distribution:

• Degrees of freedom ($\nu$): Controls tail thickness. Lower values mean fatter tails.
• Convergence: As $\nu \to \infty$, the t distribution approaches the standard normal.
• Variance: For $\nu > 2$, the variance is $\nu/(\nu-2)$, which is always greater than 1. For $1<\nu\leq2$, the mean exists but the variance is infinite. For $0<\nu\leq1$, the mean and therefore the variance are undefined. The second raw moment diverges throughout $0<\nu\leq2$.
• Kurtosis: For $\nu > 4$, excess kurtosis is $6/(\nu-4)$, always positive and decreasing in $\nu$. For $\nu \leq 4$, the kurtosis is undefined because the fourth moment does not exist.

Degrees of freedom must be estimated for the dataset and model at hand; there is no universal range for financial returns. Values near 4 are nevertheless useful for understanding the boundary behavior of the family: for $\nu \leq 4$, the fourth moment is undefined, while for $\nu > 4$ it is finite. The seeded example below uses $\nu=4$ deliberately to demonstrate a particularly heavy-tailed case.

The log-scale plot on the right reveals how dramatically the distributions differ in their tails. At 5 standard deviations, the t distribution with 3 degrees of freedom assigns probability thousands of times higher than the normal distribution. This difference, while invisible on a linear scale where both probabilities appear to be essentially zero, has enormous practical significance for risk management. The one-sided probability of a 5-sigma event under a normal distribution is about 1 in 3.5 million observations; under a t distribution with 3 degrees of freedom scaled to have unit variance, it is about 1 in 600 observations. If each observation is treated as an independent daily return and a year contains 252 trading days, those probabilities imply expected waiting times of roughly 13,900 trading years and 2.4 trading years, respectively. These are model-based average waiting times, not schedules for when events will occur.

We can fit a scaled t distribution to the seeded heavy-tailed sample and compare its maximized log-likelihood with that of a fitted normal distribution. Because the sample was generated from a t distribution, this is a reproducibility check and model-comparison demonstration, not evidence about observed market returns.

Power-law notation provides a general way to describe polynomial tail decay. Student's t distributions already have power-law tails, but other regularly varying models can use a different exponent or slowly varying factor:

This asymptotic relationship describes the right-tail survival function for large thresholds. It includes the polynomial right tail of a Student's t distribution as a special case. For a nonnegative variable, the exponent controls positive moments; for a signed variable, existence of a full moment also depends on the left tail, which this one-sided relationship does not specify.

The components of this relationship are:

• $P(X > x)$: the probability that the random variable exceeds value $x$ (the survival function)
• $x$: a large threshold value
• $\alpha$: the tail exponent, which controls how fast probabilities decay. Some authors also call $\alpha$ a tail index; in another common extreme-value-theory convention, the tail index is $\xi=1/\alpha$
• $\sim$: denotes asymptotic equivalence. The ratio of the left side to $Cx^{-\alpha}$ approaches 1 as $x \to \infty$

A larger $\alpha$ means faster polynomial decay and a smaller $\alpha$ means slower decay. As an illustration, if $\alpha=3$ and the asymptotic approximation is appropriate at both thresholds, doubling a sufficiently large threshold reduces its exceedance probability by approximately $2^{-3}=1/8$. This is much slower than the exponential falloff of normal tails.

For comparison, the one-sided standard-normal survival probability falls by a factor of about $5.1\times 10^{10}$ from 4σ to 8σ, roughly eleven orders of magnitude. Under the illustrative cubic law, the corresponding asymptotic reduction is only a factor of 8. Which approximation is appropriate must be established from a specified model or dataset rather than inferred from this contrast alone.

The implications of power-law tails are serious for risk management:

• Moment thresholds require both tails: For $p>0$ and a nonnegative variable with $P(X>x)\sim Cx^{-\alpha}$, the right-tail contribution to $E[X^p]$ is finite for $p<\alpha$ and diverges for $p\geq\alpha$. For a signed return, the displayed right-tail relation alone cannot guarantee that a full mean or variance exists because the left tail is unspecified. If both tails have comparable exponents near 3, absolute moments of orders below 3, including the mean and variance, are finite, while the fourth moment is not.
• Slow convergence: Under i.i.d. sampling with finite variance ($\alpha > 2$ in this tail model), the Central Limit Theorem applies, but its normal approximation can remain poor at practically available sample sizes.
• Tail concentration and dependence: With power-law marginal tails, a small number of observations can dominate a sample risk estimate. Diversification may be weaker when those tails combine with strong common or tail dependence, but marginal tail thickness alone does not determine diversification benefits.

Fat tails change how we should think about risk:

VaR model error: A fitted normal model can underestimate sufficiently extreme heavy left-tail quantiles, but it need not underestimate VaR at every confidence level. The direction and magnitude depend on the fitted distributions, skewness, and confidence level. The seeded example below slightly overestimates 95% VaR and underestimates 99% VaR.

Expected shortfall: Under the continuous-loss convention used here, expected shortfall is the mean loss conditional on being at or beyond the VaR threshold. Terminology such as "Conditional VaR" varies across sources, and distributions with atoms require a quantile-integral definition to handle the boundary consistently. For continuous heavy-tailed models, ES can be substantially larger than VaR.

Distribution choice begins with the variable being modeled. A normal or Student's t model can describe signed continuous returns, a lognormal model can describe a positive level generated by exponentiated normal log changes, a binomial model counts successes in a fixed number of trials, and a Poisson model counts events over an interval. A power-law tail specification is a statement about remote-tail decay rather than a complete model for the center of a distribution.

The horizon and decision also matter. A model that is adequate for central daily-return behavior may be poor for multi-day dependence, extreme quantiles, or expected shortfall. Fit diagnostics should therefore match the intended use: inspect residual dependence and stability across regimes, compare empirical and fitted quantiles, and check the particular risk measure or payoff rather than relying on one summary statistic.

The normal distribution provides a tractable baseline, while the lognormal distribution connects normal log changes to positive modeled prices. Binomial and Poisson distributions address different discrete questions: a fixed number of trials versus event counts over an interval. Student's t and power-law models make slower tail decay explicit, but their parameters and moment conditions must be checked rather than assumed.

No distribution is selected by its name alone. The relevant variable, horizon, dependence structure, evidence, and decision determine whether its assumptions are useful.

Ready to test your understanding? Take this quick quiz to reinforce what you've learned about probability distributions in finance.