Credit Risk Modeling Approaches › Structural Models: The Merton Framework

Michael BrenndoerferDecember 22, 202525 min read
Quantitative FinanceAll chapters
Part 10 of 15Risk Management and Financial Risk Controls
Chapter 55Credit Risk Modeling Approaches › Structural Models: The Merton Framework
Reading section 2 of 6 · Structural Models: The Merton Framework · 25 minSection 2 of 6 · 25 min

Book contents

Section 2 of 6 · 25 min read

Credit Risk Modeling Approaches › Structural Models: The Merton FrameworkOne section · Sections in this chapter2 / 6 ↓
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

Structural models derive default probability from the fundamental economics of a firm's capital structure. The key insight is straightforward: shareholders hold a residual claim on the firm's assets after debt is paid. If assets exceed debt at maturity, shareholders receive the difference. If assets fall short, they walk away with nothing while creditors absorb the loss. This simple observation, when combined with the machinery of option pricing theory, yields a powerful framework for quantifying credit risk.

Equity as a Call Option on Assets

To understand the Merton framework, recognize the option-like nature of equity in a leveraged firm. Consider a firm with total asset value VV financed by equity EE and zero-coupon debt with face value DD maturing at time TT. At maturity, the firm's assets must be distributed between debt holders and equity holders according to the priority of claims established by corporate law. Debt holders have senior claims and must be paid first. Only after debt is fully satisfied do equity holders receive any value.

This priority structure means that at maturity, the payoff to equityholders is:

ET=max⁡(VT−D,0)E_T = \max(V_T - D, 0)

where:

  • ETE_T: equity value at maturity
  • VTV_T: total firm asset value at maturity
  • DD: face value of debt (the strike price of the option)
  • TT: time to maturity in years
  • max⁡(VT−D,0)\max(V_T - D, 0): equityholders receive the residual after paying debt, or nothing if assets are insufficient

This payoff has the same shape as a European call option on the firm's assets with strike price DD. The analogy explains how equity and debt behave. Equity holders enjoy limited downside. They can walk away with nothing if the firm fails without being forced to inject additional capital. At the same time, equity holders have unlimited upside potential if assets appreciate materially above the debt level. This asymmetry, participating in gains while limiting losses, is the defining characteristic of an option.

Debtholders' position is the complement of equity. Debtholders receive

BT=min⁡(VT,D)=D−max⁡(D−VT,0)B_T = \min(V_T, D) = D - \max(D - V_T, 0)

where:

  • BTB_T: debt value at maturity
  • VTV_T: total firm asset value at maturity
  • DD: face value of debt
  • min⁡(VT,D)\min(V_T, D): debtholders receive the lesser of asset value or promised payment
  • max⁡(D−VT,0)\max(D - V_T, 0): payoff of a put option on assets struck at DD (represents credit risk)
  • D−max⁡(D−VT,0)D - \max(D - V_T, 0): decomposition showing debt equals risk-free bond minus a put option

The second form decomposes the debt payoff into the face value minus a put option on assets struck at D. This decomposition reveals that bondholders are effectively short a put option: the put option represents the credit risk embedded in corporate debt. When a firm defaults (assets fall below debt), bondholders suffer losses exactly equal to the intrinsic value of this embedded put. This insight turns credit risk analysis into option pricing, allowing us to apply the tools developed for derivatives valuation.

Out[3]:
Visualization
Line chart of Merton-model payoffs against firm asset value from 0 to 150. Debt payoff rises to 80 and then stays flat, while equity remains zero until asset value exceeds the debt face value of 80 and then rises linearly. A shaded, hatched region below 80 marks default.
Equity and debt payoffs at maturity versus firm asset value under the Merton structural model with debt face value D=80. Equity exhibits a call option payoff, starting at zero below 80 and increasing linearly above it, while debt payoff rises with asset value below 80 and then flattens. The shaded, hatched default region marks asset values below the debt obligation, where equity expires worthless and creditors absorb the shortfall.

The shaded region in the figure shows where default occurs. When the terminal asset value falls below the debt face value, equity holders receive nothing, as their call option finishes out of the money. Meanwhile, debtholders recover only VTV_T instead of the promised DD and absorb losses proportional to the shortfall. The visual representation makes clear why leverage matters: as debt increases relative to assets, the default region expands, and the probability of equity finishing worthless increases correspondingly.

The Merton Model Derivation

Now we develop the mathematical framework for valuing equity and debt while extracting default probabilities. The key modeling assumption is that asset value follows geometric Brownian motion under the physical measure:

dVt=μVtdt+σVVtdWtdV_t = \mu V_t dt + \sigma_V V_t dW_t

where:

  • VtV_t: firm asset value at time tt
  • μ\mu: expected asset return (drift rate under the physical measure)
  • σV\sigma_V: asset volatility, the standard deviation of asset returns
  • WtW_t: standard Brownian motion
  • dWtdW_t: infinitesimal Brownian increment
  • dtdt: infinitesimal time increment
  • dVtdV_t: infinitesimal change in asset value over time dtdt

This stochastic differential equation describes asset dynamics with two components. The deterministic drift term μ V_t dt captures the expected appreciation of assets over time, which reflects the firm's earning power and investment returns. The stochastic term σVVtdWt\sigma_V V_t dW_t introduces randomness that scales with the current asset level, so percentage changes in asset value have constant volatility regardless of the firm's size.

As we derived in Part III when developing Black-Scholes, this stochastic differential equation implies that asset returns are normally distributed, which means the terminal asset value has a lognormal distribution. Applying Itô's lemma to ln⁡Vt\ln V_t and integrating over the interval from 0 to TT yields:

ln⁡VT=ln⁡V0+(μ−σV22)T+σVTZ\ln V_T = \ln V_0 + \left(\mu - \frac{\sigma_V^2}{2}\right)T + \sigma_V \sqrt{T} Z

where:

  • VTV_T: asset value at maturity time TT
  • V0V_0: current asset value
  • ln⁡VT\ln V_T and ln⁡V0\ln V_0: natural logarithms of asset values
  • μ\mu: expected asset return (physical drift)
  • σV\sigma_V: asset volatility
  • TT: time to maturity in years
  • ZZ: standard normal random variable, Z∼N(0,1)Z \sim N(0,1)
  • μ−σV22\mu - \frac{\sigma_V^2}{2}: drift adjustment accounting for Itô's lemma, the volatility drag that reduces geometric growth
  • σVT\sigma_V \sqrt{T}: total volatility (standard deviation) accumulated over period TT
  • The formula shows that log asset returns are normally distributed with mean (μ−σV2/2)T(\mu - \sigma_V^2/2)T and variance σV2T\sigma_V^2 T

The drift adjustment term μ−σV22\mu - \frac{\sigma_V^2}{2} is important to understand. This correction arises from Itô's lemma and reflects the difference between arithmetic and geometric average returns. When returns are volatile, the geometric mean (which determines wealth accumulation) falls below the arithmetic mean by approximately half the variance, a phenomenon sometimes called volatility drag. This explains why highly volatile assets can have high expected returns yet still disappoint long-term investors.

Physical vs. Risk-Neutral Default Probability

The Merton model produces two default probabilities depending on the measure used. The physical default probability uses the actual asset drift μ\mu and represents real-world default likelihood. The risk-neutral default probability uses the risk-free rate rr for pricing credit derivatives. The two differ by the market price of risk.

Now we can derive the default probability by working through the mathematical condition for default. Default occurs when VT<DV_T < D, meaning assets at maturity are insufficient to cover debt obligations. We derive the default condition by starting from this inequality and substituting the lognormal distribution:

VT<D(default condition)ln⁡VT<ln⁡D(take logarithms)ln⁡V0+(μ−σV22)T+σVTZ<ln⁡D(substitute lognormal distribution)σVTZ<ln⁡D−ln⁡V0−(μ−σV22)T(isolate terms with Z)Z<ln⁡(D/V0)−(μ−σV22)TσVT(divide both sides by σVT)\begin{aligned} V_T &< D && \text{(default condition)} \\ \ln V_T &< \ln D && \text{(take logarithms)} \\ \ln V_0 + \left(\mu - \frac{\sigma_V^2}{2}\right)T + \sigma_V \sqrt{T} Z &< \ln D && \text{(substitute lognormal distribution)} \\ \sigma_V \sqrt{T} Z &< \ln D - \ln V_0 - \left(\mu - \frac{\sigma_V^2}{2}\right)T && \text{(isolate terms with } Z \text{)} \\ Z &< \frac{\ln(D/V_0) - \left(\mu - \frac{\sigma_V^2}{2}\right)T}{\sigma_V \sqrt{T}} && \text{(divide both sides by } \sigma_V \sqrt{T} \text{)} \end{aligned}

The derivation proceeds by first taking the logarithm of both sides, which is valid since both VTV_T and DD are positive. Next, we substitute the expression for ln⁡VT\ln V_T from our lognormal distribution result. Rearranging to isolate the random component ZZ on one side yields a threshold that ZZ must fall below for default to occur. Since ZZ is standard normal, the probability that default occurs is the probability that ZZ falls below this threshold. Defining the distance to default d2d_2 as:

d2=ln⁡(V0/D)+(μ−σV22)TσVTd_2 = \frac{\ln(V_0/D) + \left(\mu - \frac{\sigma_V^2}{2}\right)T}{\sigma_V \sqrt{T}}

The physical default probability is:

P(default)=Φ(−d2)P(\text{default}) = \Phi(-d_2)

where:

  • d2d_2: distance to default measured in standard deviations, indicates how many standard deviations the asset value is above the default threshold
  • V0V_0: current asset value
  • DD: debt face value (default barrier)
  • ln⁡(V0/D)\ln(V_0/D): log moneyness, measuring the relative position of assets versus debt
  • μ−σV22\mu - \frac{\sigma_V^2}{2}: risk-adjusted drift
  • TT: time to maturity
  • σVT\sigma_V \sqrt{T}: total volatility over the horizon
  • Φ(⋅)\Phi(\cdot): cumulative distribution function of the standard normal distribution
  • The sign flip from ln⁡(D/V0)\ln(D/V_0) to ln⁡(V0/D)\ln(V_0/D) accounts for taking the probability of the complementary event

The quantity d_2 is called the distance to default when measured in standard deviations. This terminology captures its intuitive interpretation: the parameter d2d_2 tells us how many standard deviations the expected log asset value is above the log default threshold. A higher d_2 means the firm is further from default: a larger adverse shock is required to trigger insolvency. A firm with d2=3d_2 = 3 would need assets to fall by more than three standard deviations to default, an event with very low probability under the normal distribution.

Valuing Equity and Debt

Having established the default probability framework, we now turn to valuation. Using the Black-Scholes framework from Part III Chapter 6, we can value equity as a European call option on the firm's assets. Since equity holders receive max⁡(VT−D,0)\max(V_T - D, 0) at maturity, and the underlying asset VV follows geometric Brownian motion, the current equity value equals the call option value:

E0=V0Φ(d1)−De−rTΦ(d2Q)E_0 = V_0 \Phi(d_1) - D e^{-rT} \Phi(d_2^Q)

where:

  • E0E_0: current equity value
  • V0V_0: current asset value
  • DD: debt face value (strike price)
  • rr: risk-free rate
  • TT: time to maturity
  • Φ(⋅)\Phi(\cdot): cumulative standard normal distribution
  • d1d_1 and d2Qd_2^Q: risk-neutral distance to default parameters (defined below)
  • V0Φ(d1)V_0 \Phi(d_1): present value of the asset value in scenarios where the option is exercised, weighted by delta
  • De−rTΦ(d2Q)D e^{-rT} \Phi(d_2^Q): present value of the strike price payment, weighted by risk-neutral probability of exercise
  • This is the standard Black-Scholes call option formula applied to equity as a call on firm assets

This formula has a clear economic interpretation: the first term represents the expected asset value received by equity holders in favorable scenarios, discounted and weighted by the hedge ratio. The second term represents the present value of the debt payment, weighted by the probability that equity finishes in the money. The difference captures the value of equity's limited liability and upside potential.

The parameters d1d_1 and d2Qd_2^Q are given by:

d1=ln⁡(V0/D)+(r+σV22)TσVTd_1 = \frac{\ln(V_0/D) + \left(r + \frac{\sigma_V^2}{2}\right)T}{\sigma_V \sqrt{T}} d2Q=d1−σVTd_2^Q = d_1 - \sigma_V \sqrt{T}

where:

  • V0V_0: current asset value
  • DD: debt face value
  • rr: risk-free rate (replaces μ\mu in the risk-neutral measure)
  • σV\sigma_V: asset volatility
  • TT: time to maturity
  • d1d_1: measures moneyness adjusted for volatility and time
  • d2Qd_2^Q: risk-neutral distance to default
  • Superscript QQ indicates risk-neutral (pricing) measure
  • Φ(d2Q)\Phi(d_2^Q) represents the risk-neutral probability that equity finishes in-the-money

The superscript QQ indicates risk-neutral parameters, distinguishing them from the physical measure quantities: the key difference is that risk-neutral valuation replaces the physical drift μ\mu with the risk-free rate rr, which reflects the fundamental theorem of asset pricing. Debt value can be decomposed as risk-free debt minus the embedded put option:

B0=De−rT−P0=De−rTΦ(d2Q)+V0Φ(−d1)B_0 = D e^{-rT} - P_0 = D e^{-rT}\Phi(d_2^Q) + V_0 \Phi(-d_1)

where:

  • B0B_0: current value of risky debt
  • DD: debt face value
  • rr: risk-free rate
  • TT: time to maturity
  • V0V_0: current asset value
  • P0P_0: value of the put option on assets (credit risk component)
  • Φ(⋅)\Phi(\cdot): cumulative standard normal distribution
  • d1,d2Qd_1, d_2^Q: risk-neutral distance to default parameters
  • De−rTD e^{-rT}: present value of the debt face value (risk-free component)
  • De−rTΦ(d2Q)D e^{-rT}\Phi(d_2^Q): present value of debt payment in non-default scenarios
  • V0Φ(−d1)V_0 \Phi(-d_1): present value of asset recovery in default scenarios
  • Φ(−d1)\Phi(-d_1) and Φ(−d2Q)\Phi(-d_2^Q): probabilities associated with default outcomes

This decomposition reveals the fundamental nature of credit risk. Corporate debt can be valued as risk-free debt minus a put option, where the put option represents the creditors' exposure to default. When assets fall below the debt face value, creditors receive only the asset value rather than the full promised payment. The difference between risk-free and risky debt value represents the credit spread that compensates lenders for bearing default risk. This spread increases with leverage, volatility, and time to maturity, exactly as the option framework predicts.

Implementation: The Merton Model

Let's implement the complete Merton model with both physical and risk-neutral default probabilities:

In[4]:
Code
import numpy as np
from scipy.stats import norm


class MertonModel:
    """
    Merton structural model for credit risk.

    Parameters
    ----------
    V0 : float
        Current asset value
    D : float
        Face value of debt (default barrier)
    T : float
        Time to maturity in years
    sigma_V : float
        Asset volatility
    r : float
        Risk-free rate
    mu : float
        Expected asset return (physical measure)
    """

    def __init__(self, V0, D, T, sigma_V, r, mu=None):
        self.V0 = V0
        self.D = D
        self.T = T
        self.sigma_V = sigma_V
        self.r = r
        self.mu = mu if mu is not None else r  # Defaults to risk-neutral

    def _d1_d2(self, drift):
        """Calculate d1 and d2 for given drift rate."""
        d1 = (
            np.log(self.V0 / self.D) + (drift + 0.5 * self.sigma_V**2) * self.T
        ) / (self.sigma_V * np.sqrt(self.T))
        d2 = d1 - self.sigma_V * np.sqrt(self.T)
        return d1, d2

    def physical_default_prob(self):
        """Physical (real-world) probability of default."""
        _, d2 = self._d1_d2(self.mu)
        return norm.cdf(-d2)

    def risk_neutral_default_prob(self):
        """Risk-neutral probability of default."""
        _, d2 = self._d1_d2(self.r)
        return norm.cdf(-d2)

    def distance_to_default(self):
        """Distance to default in standard deviations (physical measure)."""
        _, d2 = self._d1_d2(self.mu)
        return d2

    def equity_value(self):
        """Black-Scholes equity value (call on assets)."""
        d1, d2 = self._d1_d2(self.r)
        return self.V0 * norm.cdf(d1) - self.D * np.exp(
            -self.r * self.T
        ) * norm.cdf(d2)

    def debt_value(self):
        """Value of risky debt."""
        d1, d2 = self._d1_d2(self.r)
        # Debt = Risk-free debt - Put option on assets
        return self.D * np.exp(-self.r * self.T) * norm.cdf(
            d2
        ) + self.V0 * norm.cdf(-d1)

    def credit_spread(self):
        """Credit spread over risk-free rate (annualized)."""
        B = self.debt_value()
        # Yield on risky debt: -ln(B/D)/T
        risky_yield = -np.log(B / self.D) / self.T
        return risky_yield - self.r


# Example: A firm with moderate leverage
model = MertonModel(
    V0=100,  # Asset value
    D=70,  # Debt face value (70% leverage)
    T=1,  # 1-year horizon
    sigma_V=0.25,  # 25% asset volatility
    r=0.05,  # 5% risk-free rate
    mu=0.10,  # 10% expected asset return
)
Out[5]:
Console
Merton Model Analysis
========================================
Asset Value:          $100.00
Debt Face Value:      $70.00
Leverage Ratio:       70.0%
Asset Volatility:     25.0%

Risk Metrics:
  Distance to Default:    1.702 std devs
  Physical Default Prob:  4.4406%
  Risk-Neutral PD:        6.6587%

Valuation:
  Equity Value:           $33.86
  Debt Value:             $66.14
  Total Firm Value:       $100.00
  Credit Spread:          66.7 bps

The distance to default of 1.70 standard deviations indicates the firm has moderate default risk. A higher distance indicates lower credit risk, as larger adverse shocks are required to trigger default. The risk-neutral default probability exceeds the physical probability because the risk-neutral drift (the risk-free rate) is lower than the physical drift (the expected asset return), making downward paths more likely under the pricing measure. This difference reflects the market price of risk. The credit spread compensates lenders for bearing this default risk.

Out[6]:
Visualization
Line chart showing default probability declining as distance to default rises from minus 1 to 4 standard deviations. Shaded, patterned regions identify high, moderate, and low risk. Dashed vertical and dotted horizontal references mark the current firm at distance 1.70 and default probability 4.44 percent.
Relationship between distance to default and default probability in the Merton framework. The curve falls from 50% at distance zero toward 0% at four standard deviations. Three shaded zones partition the credit spectrum: high risk (DD < 1.5), moderate risk (1.5 < DD < 2.5), and low risk (DD > 2.5). The example firm is marked at DD=1.70 and physical default probability 4.44%.

The figure illustrates how distance to default maps to default probability through the cumulative normal distribution. Firms with distance to default below 1.5 standard deviations face default risk exceeding 7%, placing them in the high-risk zone. The moderate risk zone spans 1.5 to 2.5 standard deviations, with default probabilities between 0.6% and 7%. Firms with distance to default above 2.5 standard deviations are considered low risk, with default probability below 0.6%. Our example firm has distance to default of 1.70 and physical default probability of 4.44%, placing it within the moderate-risk zone.

Key Parameters

These key parameters drive model behavior and credit risk dynamics.

  • V0: Current firm asset value. This represents the total value of the firm's assets, including both tangible and intangible assets. Higher asset values relative to debt reduce default probability by increasing the buffer between current value and the default barrier.
  • D: Face value of debt (the default barrier). This parameter determines the threshold that assets must exceed to avoid default. Higher debt increases leverage and default risk by raising the barrier that must be cleared.
  • T: Time to maturity in years. Longer horizons increase uncertainty about future asset values and option value. Time to maturity has competing effects on default probability. Longer time allows both more opportunity for assets to appreciate and more opportunity for adverse shocks.
  • sigma_V: Asset volatility. Higher volatility increases both equity value (through increased option value) and default probability (through increased probability of large adverse movements). This dual effect explains why equity holders may prefer higher volatility while debt holders prefer lower volatility.
  • r: Risk-free rate. Used in the risk-neutral pricing measure for valuation purposes.
  • mu: Expected asset return under the physical measure. Used for physical default probability calculations to assess real-world default risk.

Sensitivity Analysis

Credit risk depends critically on leverage and asset volatility. Let's examine how default probability varies with these parameters to build intuition for the model's behavior:

Out[8]:
Visualization
Line chart of physical default probability against leverage from 30 to 90 percent. Three curves represent asset volatilities of 15, 25, and 35 percent; each rises nonlinearly, with the 35 percent scenario highest and the 15 percent scenario lowest.
Physical default probability versus leverage ratio for three asset volatility levels (15%, 25%, 35%) under the one-year Merton model. Default probability increases nonlinearly with leverage, and the separation among volatility scenarios widens sharply above roughly 60% leverage. The distinct line styles show that both leverage and asset volatility are central to credit risk.

The nonlinear relationship between leverage and default probability reflects the option-like nature of credit risk. At low leverage (30-40%), even large changes in the leverage ratio have minimal impact on default probability, remaining below 5% even at 25% volatility. This corresponds to deeply in-the-money equity options where the probability of finishing out of the money is negligible. As leverage increases beyond 60%, credit risk accelerates because the equity option moves from deep in the money toward at the money, where option values are most sensitive to changes in moneyness. At 80% leverage, the 35% volatility scenario produces default probability above 20%, while the 15% scenario remains near 2%. The widening gap shows how leverage amplifies the effect of volatility on credit quality.

Estimating Unobservable Asset Value and Volatility

The Merton model's key practical challenge is that asset value VV and asset volatility σV\sigma_V are not directly observable. While equity prices trade on public exchanges, a firm's total asset value includes intangible assets, growth options, and synergies that cannot be directly measured. We observe equity value EE (market capitalization) and equity volatility σE\sigma_E from option prices or historical returns. We invert the model to recover asset parameters from these observables.

Applying Itô's lemma (as we did when deriving the Black-Scholes PDE in Part III), we can derive how equity volatility relates to asset volatility through the chain rule of stochastic calculus:

σE=∂E∂V⋅VE⋅σV=Φ(d1)VEσV\sigma_E = \frac{\partial E}{\partial V} \cdot \frac{V}{E} \cdot \sigma_V = \frac{\Phi(d_1) V}{E} \sigma_V

where:

  • σE\sigma_E: equity volatility (observable from market prices)
  • σV\sigma_V: asset volatility (unobservable, must be estimated)
  • EE: equity value
  • VV: asset value
  • ∂E∂V\frac{\partial E}{\partial V}: delta of the equity call option, equal to Φ(d1)\Phi(d_1)
  • VE\frac{V}{E}: leverage ratio in the option context
  • Φ(d1)\Phi(d_1): the option delta, measuring sensitivity of equity value to asset value changes
  • The product Φ(d1)VE\frac{\Phi(d_1) V}{E} represents the leverage multiplier: equity volatility exceeds asset volatility because equity is a levered claim

This relationship captures the leverage effect in credit markets. Since equity represents a call option on assets, small percentage changes in asset value are magnified into larger percentage changes in equity value. The leverage multiplier Φ(d1)VE\frac{\Phi(d_1) V}{E} exceeds one for leveraged firms, explaining why equity volatility typically exceeds asset volatility. A highly leveraged firm near default will have extremely high equity volatility even if asset volatility is moderate, because small asset fluctuations translate into large percentage changes in the thin equity cushion.

This gives us two equations to work with:

  1. Equity pricing: E=VΦ(d1)−De−rTΦ(d2)E = V\Phi(d_1) - De^{-rT}\Phi(d_2)
  2. Volatility relationship: σEE=Φ(d1)VσV\sigma_E E = \Phi(d_1) V \sigma_V

Both equations involve the same unknowns (VV and σV\sigma_V), and both can be expressed using market observables (EE, σE\sigma_E, DD, rr, TT). We solve this system numerically:

In[9]:
Code
import numpy as np
from scipy.optimize import fsolve
from scipy.stats import norm


def estimate_merton_params(E, sigma_E, D, T, r):
    """
    Estimate asset value and volatility from equity market data.

    Parameters
    ----------
    E : float
        Market value of equity
    sigma_E : float
        Equity volatility
    D : float
        Face value of debt
    T : float
        Time to maturity
    r : float
        Risk-free rate

    Returns
    -------
    V : float
        Estimated asset value
    sigma_V : float
        Estimated asset volatility
    """

    def equations(params):
        V, sigma_V = params
        d1 = (np.log(V / D) + (r + 0.5 * sigma_V**2) * T) / (
            sigma_V * np.sqrt(T)
        )
        d2 = d1 - sigma_V * np.sqrt(T)

        # Equation 1: Equity pricing equation
        eq1 = V * norm.cdf(d1) - D * np.exp(-r * T) * norm.cdf(d2) - E

        # Equation 2: Volatility relationship
        eq2 = norm.cdf(d1) * V * sigma_V - sigma_E * E

        return [eq1, eq2]

    # Initial guess: V = E + D, sigma_V = sigma_E
    initial_guess = [E + D * np.exp(-r * T), sigma_E]
    V, sigma_V = fsolve(equations, initial_guess)

    return V, sigma_V


# Example: Estimate from market data
E_market = 35  # Market cap
sigma_E_market = 0.50  # Equity volatility (50%)
D_face = 70  # Debt face value
T_horizon = 1  # 1 year
r_free = 0.05  # Risk-free rate

V_est, sigma_V_est = estimate_merton_params(
    E_market, sigma_E_market, D_face, T_horizon, r_free
)
Out[10]:
Console
Merton Model Parameter Estimation
========================================
Market Observables:
  Equity Value:       $35.00
  Equity Volatility:  50%
  Debt Face Value:    $70.00

Estimated Parameters:
  Asset Value:        $101.55
  Asset Volatility:   17.3%
  Leverage (D/V):     68.9%

Model Outputs:
  Distance to Default: 2.637
  Physical PD:         0.42%
  Credit Spread:       5 bps

The estimation successfully backs out the unobservable asset parameters from market data. The estimated asset value of approximately 101.6millionexceedstheequityvalueof101.6 million exceeds the equity value of 35 million by the present value of debt, confirming the accounting identity that assets equal equity plus debt. The estimated asset volatility of 17.3% is well below the observed equity volatility of 50%, which reflects the leverage effect: debt acts as a fixed obligation that amplifies equity movements. Think of a seesaw: if total assets move by 1%, equity (the smaller end) must move by proportionally more to maintain balance. The resulting 68.9% leverage ratio and distance to default of 2.64 standard deviations translate to a physical default probability of 0.42% and a model credit spread near 5 basis points.

Out[11]:
Visualization
Line chart showing equity volatility rising from about 35 percent at 30 percent leverage to nearly 100 percent at 85 percent leverage. A horizontal dashed line marks the assumed 25 percent asset volatility, and a vertical dotted line plus annotation marks the estimated firm at 68.9 percent leverage.
Equity volatility increases nonlinearly with leverage ratio (D/V) for a fixed illustrative asset volatility of 25%, accelerating above 60% leverage and approaching 90% near 80% leverage. The horizontal reference marks the assumed asset volatility, while the vertical reference and annotation identify the estimated firm's 68.9% leverage ratio.

The leverage effect visualization demonstrates why equity volatility systematically exceeds asset volatility for leveraged firms. At 30% leverage, equity volatility is only slightly above asset volatility because the equity cushion is substantial. As leverage increases, the equity cushion shrinks, and any percentage change in asset value translates into a proportionally larger percentage change in equity value. At 80% leverage, equity volatility reaches nearly three times asset volatility. This relationship has practical implications for credit risk assessment: when estimating asset volatility from equity volatility, failing to account for the leverage effect would substantially overestimate asset risk.