The Black-Scholes formula stands as one of the most significant achievements in quantitative finance, providing a theoretical framework for pricing European options. Aspiring quant researchers, especially those interviewing at top-tier firms like Citadel, are frequently tested on their understanding of this formula—not just its application, but its derivation and underlying assumptions. 

Citadel Quant Research Intern Interview Question: Deriving the Black-Scholes Formula

Table of Contents

• Introduction to the Black-Scholes Formula
• Assumptions of the Black-Scholes Model
• Mathematical Preliminaries
• Stochastic Processes and the SDE for Stock Prices
• Portfolio Construction and Hedging
• Derivation of the Black-Scholes Partial Differential Equation (PDE)
• Risk-Neutral Valuation
• Solving the Black-Scholes PDE
• The Final Black-Scholes Formula
• Implications and Limitations of the Model

Introduction to the Black-Scholes Formula

The Black-Scholes formula, developed by Fischer Black, Myron Scholes, and Robert Merton in the early 1970s, revolutionized the field of financial engineering. It provides an analytic solution for pricing European-style options—contracts that can only be exercised at expiration. The formula is a direct consequence of modeling the evolution of asset prices using stochastic calculus and constructing a riskless hedged portfolio.

Understanding the Black-Scholes formula is not only crucial for passing quant interviews at elite firms like Citadel, but also foundational for anyone pursuing a career in quantitative finance or derivatives trading.

Assumptions of the Black-Scholes Model

Before diving into the mathematical derivation, it is critical to understand the assumptions underlying the Black-Scholes model. These assumptions simplify the real financial markets and make the mathematics tractable.

• Frictionless Markets: No transaction costs or taxes; assets are infinitely divisible.
• Continuous Trading: Trading can occur continuously in time.
• Short Selling Allowed: Investors can short-sell securities with full use of proceeds.
• Constant Risk-Free Rate: The risk-free interest rate \( r \) is constant and known.
• Log-Normal Asset Prices: The price of the underlying asset follows a geometric Brownian motion with constant drift (\( \mu \)) and volatility (\( \sigma \)).
• No Dividends: The underlying asset does not pay dividends during the option's life.
• European Option: The option can only be exercised at maturity.

These idealized conditions are not fully met in the real world, but they are useful for building a tractable model.

Mathematical Preliminaries

The Black-Scholes model is built upon several mathematical concepts. Understanding these is essential for following the derivation.

Brownian Motion and Geometric Brownian Motion

• Brownian Motion (\( W_t \)): A stochastic process with continuous paths and stationary, independent increments that are normally distributed. Formally, \( W_0 = 0 \), and for all \( t \), \( W_t \sim N(0, t) \).
• Geometric Brownian Motion: If \( S_t \) denotes the asset price at time \( t \), then it follows the SDE:
  $$ dS_t = \mu S_t dt + \sigma S_t dW_t $$

Itô's Lemma

Itô's Lemma is a fundamental result in stochastic calculus. It provides the differential of a function of a stochastic process. For a function \( f(t, S_t) \):

$$ df = \left( \frac{\partial f}{\partial t} + \mu S_t \frac{\partial f}{\partial S} + \frac{1}{2} \sigma^2 S_t^2 \frac{\partial^2 f}{\partial S^2} \right) dt + \sigma S_t \frac{\partial f}{\partial S} dW_t $$

We will use Itô's Lemma to derive the dynamics of an option price on a stock following geometric Brownian motion.

Stochastic Processes and the SDE for Stock Prices

The Black-Scholes model assumes the following stochastic differential equation (SDE) for the asset price:

$$ dS_t = \mu S_t dt + \sigma S_t dW_t $$

• \( S_t \): Asset price at time \( t \)
• \( \mu \): Expected return (drift)
• \( \sigma \): Volatility of returns
• \( dW_t \): Increment of a Wiener process (Brownian motion)

This SDE indicates that the asset price evolves with a deterministic trend (\( \mu S_t dt \)), and a stochastic component (\( \sigma S_t dW_t \)), both proportional to the current price.

Solution to the SDE

The explicit solution to this SDE is:

$$ S_t = S_0 \exp \left( \left(\mu - \frac{1}{2}\sigma^2\right)t + \sigma W_t \right) $$

This demonstrates that the log-returns of the asset are normally distributed, hence the asset price is log-normally distributed.

Portfolio Construction and Hedging

A central idea in the Black-Scholes derivation is constructing a riskless hedged portfolio by holding a position in the option (whose price we want to determine) and a position in the underlying stock. By dynamically adjusting the portfolio to eliminate risk, we can argue that the portfolio should earn the risk-free rate.

The Replicating Portfolio

Let \( V(S, t) \) denote the price of an option as a function of the asset price \( S \) and time \( t \). Construct a portfolio \( \Pi \) as follows:

$$ \Pi = V(S, t) - \Delta S $$

• Long 1 option (value \( V \))
• Short \( \Delta \) shares of the underlying (value \( -\Delta S \))

Here, \( \Delta \) is chosen to eliminate risk.

Derivation of the Black-Scholes Partial Differential Equation (PDE)

By applying Itô's Lemma to \( V(S, t) \), we have:

$$ dV = \frac{\partial V}{\partial t} dt + \frac{\partial V}{\partial S} dS + \frac{1}{2} \frac{\partial^2 V}{\partial S^2} (dS)^2 $$

Recall that under the SDE:

$$ dS = \mu S dt + \sigma S dW $$

and

$$ (dS)^2 = (\sigma S dW)^2 = \sigma^2 S^2 (dW)^2 = \sigma^2 S^2 dt $$

because \( (dW)^2 = dt \) in stochastic calculus.

Therefore,

$$ dV = \frac{\partial V}{\partial t} dt + \frac{\partial V}{\partial S} (\mu S dt + \sigma S dW) + \frac{1}{2} \frac{\partial^2 V}{\partial S^2} \sigma^2 S^2 dt \\ = \left[ \frac{\partial V}{\partial t} + \mu S \frac{\partial V}{\partial S} + \frac{1}{2} \sigma^2 S^2 \frac{\partial^2 V}{\partial S^2} \right] dt + \sigma S \frac{\partial V}{\partial S} dW $$

Portfolio Dynamics

The change in the hedged portfolio is:

$$ d\Pi = dV - \Delta dS $$

Substitute the expressions for \( dV \) and \( dS \):

$$ d\Pi = \left[ \frac{\partial V}{\partial t} + \mu S \frac{\partial V}{\partial S} + \frac{1}{2} \sigma^2 S^2 \frac{\partial^2 V}{\partial S^2} \right] dt + \sigma S \frac{\partial V}{\partial S} dW - \Delta (\mu S dt + \sigma S dW) $$

Group the \( dt \) and \( dW \) terms:

$$ d\Pi = \left[ \frac{\partial V}{\partial t} + \mu S \frac{\partial V}{\partial S} - \Delta \mu S + \frac{1}{2} \sigma^2 S^2 \frac{\partial^2 V}{\partial S^2} \right] dt + (\sigma S \frac{\partial V}{\partial S} - \Delta \sigma S) dW $$

Eliminating Risk

To create a riskless portfolio, we set the coefficient of \( dW \) to zero:

$$ \sigma S \frac{\partial V}{\partial S} - \Delta \sigma S = 0 \implies \Delta = \frac{\partial V}{\partial S} $$

Thus, our portfolio will be riskless over an infinitesimal interval.

Portfolio Growth at Risk-Free Rate

Since the portfolio is riskless, it should grow at the risk-free rate \( r \):

$$ d\Pi = r \Pi dt $$

Recall:

$$ \Pi = V - \Delta S = V - \frac{\partial V}{\partial S} S $$

Plugging in the above expressions:

$$ \left[ \frac{\partial V}{\partial t} + \frac{1}{2} \sigma^2 S^2 \frac{\partial^2 V}{\partial S^2} \right] dt = r \left( V - S \frac{\partial V}{\partial S} \right) dt $$

Rearrange terms:

$$ \frac{\partial V}{\partial t} + \frac{1}{2}\sigma^2 S^2 \frac{\partial^2 V}{\partial S^2} + r S \frac{\partial V}{\partial S} - r V = 0 $$

This is the Black-Scholes Partial Differential Equation (PDE).

Risk-Neutral Valuation

A key insight is that in the Black-Scholes world, the expected return (\( \mu \)) of the underlying does not appear in the PDE. Instead, the only relevant rate is the risk-free rate (\( r \)). This is a consequence of risk-neutral valuation.

Risk-Neutral Measure (\( \mathbb{Q} \))

Under the risk-neutral measure, all assets are expected to grow at the risk-free rate. The SDE for the underlying under \( \mathbb{Q} \) is:

$$ dS_t = r S_t dt + \sigma S_t dW_t^{\mathbb{Q}} $$

where \( dW_t^{\mathbb{Q}} \) is a Brownian motion under the risk-neutral measure.

The value of a derivative (such as a European option) is the discounted expected payoff under the risk-neutral measure:

$$ V(S, t) = \mathbb{E}^{\mathbb{Q}}\left[ e^{-r(T-t)} \Phi(S_T) \mid S_t = S \right] $$

• \( \Phi(S_T) \): Option payoff at maturity
• \( T \): Expiry time

For a European call option,

$$ \Phi(S_T) = \max(S_T - K, 0) $$

where \( K \) is the strike price.

Solving the Black-Scholes PDE

Let’s derive an explicit formula for a European call option. We need to solve the Black-Scholes PDE with the terminal condition:

$$ V(S, T) = \max(S - K, 0) $$

Transformation to the Heat Equation

We can transform the Black-Scholes PDE into the classic heat/diffusion equation. Define:

• \( \tau = T - t \) (time to maturity)
• \( x = \ln\left( \frac{S}{K} \right) \) (log-moneyness)

After substitution and some algebra, the Black-Scholes PDE becomes a heat equation, which has a known solution. For brevity, we’ll state the final result:

Solution: The Black-Scholes Formula

The price of a European call option is:

$$ C(S, t) = S N(d_1) - K e^{-r(T-t)} N(d_

Where:

$$ d_1 = \frac{\ln\left(\frac{S}{K}\right) + \left(r + \frac{\sigma^2}{2}\right)(T-t)}{\sigma \sqrt{T-t}} $$ $$ d_2 = d_1 - \sigma \sqrt{T-t} $$

• \( S \): Current stock price
• \( K \): Strike price
• \( r \): Risk-free interest rate
• \( T-t \): Time to maturity
• \( \sigma \): Volatility of the asset
• \( N(\cdot) \): Cumulative distribution function (CDF) of the standard normal distribution

This formula gives the fair price of a European call option in a frictionless and idealized market.

Put Option (Put-Call Parity)

By using put-call parity, we can write the price of a European put option as:

$$ P(S, t) = K e^{-r (T-t)} N(-d_2) - S N(-d_1) $$

The Final Black-Scholes Formula

For clarity, here is the full Black-Scholes formula for a European call option:

$$ C(S, t) = S N(d_1) - K e^{-r (T-t)} N(d_2) $$

And for a European put option:

$$ P(S, t) = K e^{-r (T-t)} N(-d_2) - S N(-d_1) $$

Where:

$$ d_1 = \frac{\ln(S/K) + (r + \frac{\sigma^2}{2})(T-t)}{\sigma \sqrt{T-t}}, \quad d_2 = d_1 - \sigma \sqrt{T-t} $$

Implications and Limitations of the Model

The Black-Scholes formula was a breakthrough, enabling the widespread use of options and derivatives in modern finance. However, it is crucial to understand its limitations and where its assumptions may not hold in practice:

• No Dividends: The formula does not account for dividends. Extensions exist (e.g., Black-Scholes-Merton) for dividend-paying stocks.
• Constant Volatility: Real market volatility is stochastic and can change over time, leading to phenomena like volatility skew and smile.
• Continuous Trading: In reality, trading is not continuous, and market liquidity can dry up, especially during stress.
• No Transaction Costs: Real markets have bid-ask spreads, commissions, and other frictions.
• Lognormal Returns: Asset returns may display fat tails, jumps, and other non-Gaussian behaviors.

Despite these limitations, the Black-Scholes model remains the starting point for most option pricing models and is a core concept tested in quant research interviews at top firms like Citadel.

Python Implementation Example

Here is a simple Python function for calculating the Black-Scholes price for a European call option:

import math
from scipy.stats import norm

def black_scholes_call(S, K, T, r, sigma):
    """
    S: Spot price
    K: Strike price
    T: Time to maturity (in years)
    r: Risk-free rate
    sigma: Volatility
    """
    d1 = (math.log(S / K) + (r + 0.5 * sigma ** 2) * T) / (sigma * math.sqrt(T))
    d2 = d1 - sigma * math.sqrt(T)
    call = S * norm.cdf(d1) - K * math.exp(-r * T) * norm.cdf(d2)
    return call

Summary Table: Black-Scholes Formula Inputs and Interpretation

Frequently Asked Interview Questions

• What are the main assumptions behind the Black-Scholes model?
• Why does the expected return (\( \mu \)) of the stock not appear in the final pricing formula?
• How would you extend the Black-Scholes formula to account for dividends?
• What is the significance of the risk-neutral measure?
• How would stochastic volatility or jumps in asset prices affect the formula?

Conclusion

The Black-Scholes formula is a cornerstone of quantitative finance and a common topic in Citadel Quant Research Intern interviews. Its derivation elegantly combines stochastic calculus, risk-neutral valuation, and financial economics. Mastery of this derivation—and a deep understanding of the model's assumptions and limitations—will not only help you excel in interviews but also lay the groundwork for advanced study and real-world application in derivatives pricing and risk management.

For any aspiring quant, practicing the derivation of the Black-Scholes formula and being able to clearly articulate each step is an essential skill that showcases both mathematical rigor and financial intuition.

References

• Black, F., & Scholes, M. (1973). The Pricing of Options and Corporate Liabilities. Journal of Political Economy, 81(3), 637-654.
• Hull, J.C. Options, Futures, and Other Derivatives. Pearson.
• Shreve, S.E. Stochastic Calculus for Finance II: Continuous-Time Models.

Related Articles

• JP Morgan Quant Interview Questions: What to Expect and How to Prepare
• Interview Assessment Experience - ADIA
• Quant Analyst Interview Questions at Millennium
• Why Downside Correlation Rises in Market Crashes and Its Impact on Volatility
• Goldman Sachs Quant Interview: Probability and Statistics Questions

Related Articles

• Data Scientist Interview Questions Asked at Netflix, Apple & LinkedIn
• Top Data Scientist Interview Questions Asked at Google and Meta
• Quant Finance and Data Science Interview Prep: Key Strategies and Resources
• JP Morgan AI Interview Questions: What to Expect and Prepare For
• Best Books to Prepare for Quant Interviews: Essential Reading List