Modeling Asset Prices in Financial Markets: Arithmetic and Geometric Brownian Motions

Saral BINDAL

In this article, Saral BINDAL (Indian Institute of Technology Kharagpur, Metallurgical and Materials Engineering, 2024-2028 & Research assistant at ESSEC Business School) presents two statistical models used in finance to describe the time behavior of asset prices: the arithmetic Brownian motion (ABM) and the geometric Brownian motion (GBM).

Introduction

In financial markets, performance over time is governed by three fundamental variables: the drift (μ), volatility (σ), and maybe most importantly time (T). The drift represents the expected growth rate of the price and corresponds to the expected return of assets or portfolios. Volatility measures the uncertainty or risk associated with price fluctuations around this expected growth and corresponds to the standard deviation of returns. The relationship between these variables reflects the trade-off between risk and return. Time, which is related to the investment horizon set by the investor, determines how both performance and risk accumulate. Together, these variables form the foundation of asset pricing to model the behavior of market price over time, and in fine the performance of the investor at their investment horizon.

Modeling asset prices

Asset price modeling is used to understand the expected return and risk in asset management, risk management, and the pricing of complex financial products such as options and structured products. Although asset prices are influenced by countless unpredictable risk factors, quants in finance always try to find a parsimonious way to model asset prices (using a few parameters only).

The first study of asset price modelling dates from Louis Bachelier in 1900, in his doctoral thesis Théorie de la Spéculation (The Theory of Speculation), where he modelled stock prices as a random walk and applied this framework to option valuation. Later, in 1923, the mathematician Norbert Wiener formalized these ideas as the Wiener process, providing the rigorous stochastic foundation that underpins modern finance.

In the 1960s, Paul Samuelson refined Bachelier’s model by introducing the geometric Brownian motion, which ensures positive stock prices following a lognormal statistical distribution. His 1965 paper “Rational Theory of Warrant Pricing” laid the groundwork for modern asset price modelling, showing that discounted stock prices follow a martingale.

We detail below the two models usually used in finance to model the evolution of asset prices over time: the arithmetic Brownian motion (ABM) and the geometric Brownian motion (GBM). We will then use these models to simulate the evolution of asset prices over time with the Monte Carlo simulation method.

Arithmetic Brownian motion (ABM)

Theory

One of the most widely used stochastic processes in financial modeling is the arithmetic Brownian motion, also known as the Wiener process. It is a continuous stochastic process with normally distributed increments. Using the Wiener process notation, an asset price model in continuous time based on an ABM can be expressed as the following stochastic differential equation (SDE):


SDE for the arithmetic Brownian motion

where:

  • dSt = infinitesimal change in asset price at time t t
  • μ = drift (growth rate of the asset price)
  • σ = volatility (standard deviation)
  • dWt = infinitesimal increment of wiener process (N(0,dt))

Note that the standard Brownian motion is a special case of the arithmetic Brownian motion with a mean equal to zero and a variance equal to one.

In this model, both μ and σ are assumed to be constant over time. It can be shown that the probability distribution function of the future price is a normal distribution implying a strictly positive (although negligible in most cases) probability for the price to be negative.

Integrating the SDE for dSt over a finite interval (from time 0 to time t), we get:


Integrated SDE for the arithmetic Brownian motion

Here, Wt is defined as Wt = √t · Zt, where Zt is a normal random variable drawn from the standard distribution N(0, 1) with mean equal to 0 and variance equal to 1.

At any date t, we can also compute the expected value and a confidence interval such that the asset price St lies between the lower and upper bound of the interval with probability equal to 1-α.


Theoritical formulas for mean, upper and lower limits of ABM model

Where S0 is the initial asset price and zα.

The z-score for a confidence level of (1 – α) can be calculated as:


z-score formula

where Φ-1 denotes the inverse cumulative distribution function (CDF) of the standard normal distribution.

For example the statistical z-score (zα) values for 66%, 95%, and 99% confidence intervals are as the following:


z-score examples

Monte Carlo simulations with ABM

Since Monte Carlo simulations are performed in discrete time, the underlying continuous-time asset price process (ABM) is approximated using the Euler–Maruyama discretization of SDEs (see Maruyama, 1955), as shown below.


Discretization formula for the arithmetic Brownian motion (ABM)

where Δt denotes the time step, expressed in the same time units as the drift parameter μ and the volatility parameter σ (usually the annual unit). For example, Δt may be equal to one day (=1/252) or one month (=1/12).

Figure 1 below illustrates a single simulated asset price path under an arithmetic Brownian motion (ABM), sampled at monthly intervals (Δt = 1/12) over a 10-year horizon (T = 10). Alongside the simulated path, the figure shows the expected (mean) price trajectory and the corresponding upper and lower bounds of a 66% confidence interval. In this example, the model assumes an annual drift (μ) of $8, representing the expected growth rate, and an annual volatility (σ) of $15, capturing random price fluctuations. The initial asset price (S0) is equal to $100.

Figure 1. Single Monte Carlo–simulated asset price path under an Arithmetic Brownian Motion model.
A Monte Carlo–simulated price path under an arithmetic Brownian motion model
Source: computation by the author (with Excel).

Figure 2 below illustrates 1,000 simulated asset price paths generated under an arithmetic Brownian motion (ABM). In addition to the simulated paths, the figure displays the expected (mean) price trajectory along with the corresponding upper and lower bounds of a 66% confidence interval, using the same parameter settings as in Figure 1.

Figure 2. Monte Carlo–simulated asset price paths under an Arithmetic Brownian Motion model.
Monte Carlo–simulated price paths under an arithmetic Brownian motion model.
Source: computation by the author (with R).

Geometric Brownian motion (GBM)

Theory

Since an arithmetic Brownian motion (ABM) can take negative values, it is unsuitable for directly modeling stock prices if we assume limited liability for investors. Under limited liability, an investor’s maximum possible loss is indeed confined to their initial investment, implying that asset prices cannot fall below zero. To address this limitation, financial models instead use geometric Brownian motion (GBM), a non-negative stochastic process that is widely employed to describe the evolution of asset prices. Using the Wiener process notation, an asset price model in continuous time based on a GBM can be expressed as the following stochastic differential equation (SDE):


SDE for the geometric Brownian motion (GBM)

where:

  • St = asset price at time t t
  • μ = drift (growth rate of the asset price)
  • σ = volatility (standard deviation)
  • dWt = infinitesimal increment of wiener process (N(0,dt))

Integrating the SDE for dSt/St over a finite interval, we get:


Integrated SDE for the geometric Brownian motion (GBM)

The theoretical expected value and confidence intervals are given analytically by the following expressions:


Theoritical formulas for mean, upper and lower limits of GBM model

Monte Carlo simulations with GBM

To implement Monte Carlo simulations, we approximate the underlying continuous-time process in discrete time, yielding:


Asset price under discrete GBM

where Zt is a standard normal random variable drawn from the distribution N(0, 1) and Δt denotes the time step, chosen so that it is expressed in the same time units as the drift parameter μ and the volatility parameter σ.

Figure 3 below illustrates a single simulated asset price path under a geometric Brownian motion (GBM), sampled at monthly intervals (Δt = 1/12) over a 10-year horizon (T = 10). Alongside the simulated path, the figure shows the expected (mean) price trajectory and the corresponding upper and lower bounds of a 66% confidence interval. In this example, the model assumes an annual drift (μ) of 8%, representing the expected growth rate, and an annual volatility (σ) of 15%, capturing random price fluctuations. The initial asset price is S0 €100.

Figure 3. Monte Carlo–simulated asset price path under a Geometric Brownian Motion model.
Monte Carlo–simulated asset price path under a GBM model.
Source: computation by the author (with Excel).

Figure 4 below illustrates 1,000 simulated asset price paths generated under a geometric Brownian motion (GBM). In addition to the simulated paths, the figure displays the expected (mean) price trajectory along with the corresponding upper and lower bounds of a 66% confidence interval, using the same parameter settings as in Figure 3.

Figure 4. Monte Carlo–simulated asset price paths under a Standard Brownian Motion model.
 Monte Carlo–simulated asset price paths under a Geometric Brownian Motion model.
Source: computation by the author (with R).

Discussion

The drift μ represents the expected rate of growth of asset prices, so its cumulative contribution increases linearly with time as μT. In contrast, volatility σ captures investment risk, and its cumulative impact scales with the square root of time as σ√T. As a result, over short horizons stochastic shocks tend to dominate the deterministic drift, whereas over longer horizons the expected growth component becomes increasingly prominent.

When many paths for the asset price are simulated and plotted over time, the resulting trajectories form a cone-shaped region, commonly referred to as a fan chart. The center of this fan traces the smooth expected path governed by the drift μ, while the widening envelope reflects the growing dispersion of outcomes induced by volatility σ.

This representation underscores a key implication for long-term investing and risk management: uncertainty expands with the investment horizon even when model parameters remain constant. While the expected value evolves predictably and linearly through time, the range of plausible outcomes broadens at a slower, square-root rate, shaping the risk–return trade-off across different time scales.

You can download the Excel file provided below for generating Monte Carlo Simulations for asset prices modeled on arithmetic and geometric Brownian motion.

Download the Excel file.

You can download the Python code provided below, for generating Monte Carlo Simulations for asset prices modeled on arithmetic and geometric Brownian motion.

Download the Python code.

Alternatively, you can download the R code below with the same functionality as in the Python file.

 Download the R code.

Link between the ABM and the GBM

The ABM and GBM models are fundamentally different: the drift for the ABM is additive while the drift for the GBM is multiplicative. Moreover, the statistical distribution for the price for the ABM is a normal distribution while the statistical distribution for the GBM is a log-normal distribution. However, we can study the relationship between the two models as they are both used to model the same phenomenon, the evolution of asset prices over time in our case.

We can especially study the relationship between the two parameters of the two models, μ and σ. In the presentation above, we used the same notations for μ and σ for the two models, but the values of these parameters for the two models will be different when we apply these models to the same phenomenon. There is no mapping of the ABM and GBM in the price space such that we get the same results as the two models are fundamentally different.

Let us rewrite the two models (in terms of SDE) by differentiating the parameters for each model:


SDE for the ABM and GBM

To model the same phenomenon, we can use the following relationship between the parameters of the ABM and GBM models:


Link between the ABM and GBM parameters.

To make the two models comparable in terms of price behavior, an ABM can locally approximate GBM by matching instantaneous drift and volatility such that:


Local link between the ABM and GBM parameters.

This local correspondence is state-dependent and time-varying, and therefore not a true parameter equivalence.

Figure 5 below compares the asset price path for an ABM, monthly adjusted ABM and a GBM.


Simulated asset price paths for ABM, adjusted ABM and GBM.

Why should I be interested in this post?

Understanding how asset prices are modeled, and in particular the difference between additive and multiplicative price dynamics, is essential for building strong intuition about how prices evolve over time under uncertainty. This understanding forms the foundation of modern risk management, as it directly informs concepts such as capital protection, downside risk, and the long-term behavior of investment portfolios.

Related posts on the SimTrade blog

   ▶ Saral BINDAL Historical Volatility

   ▶ Saral BINDAL Implied Volatility and Option Prices

   ▶ Jayati WALIA Brownian Motion in Finance

   ▶ Jayati WALIA Monte Carlo simulation method

Useful resources

Academic research

Bachelier L. (1900) Théorie de la spéculation. Annales scientifiques de l’École Normale Supérieure, 3e série, 17, 21–86.

Kataoka S. (1963) A stochastic programming model. Econometrica, 31, 181–196.

Lawler G.F. (2006) Introduction to Stochastic Processes, 2nd Edition, Chapman & Hall/CRC, Chapter “Brownian Motion”, 201–224.

Maruyama G. (1955) Continuous Markov processes and stochastic equations. Rendiconti del Circolo Matematico di Palermo, 4, 48–90.

Samuelson P.A. (1965) Rational theory of warrant pricing. Industrial Management Review, 6(2), 13–39.

Telser L. G. (1955) Safety-first and hedging. Review of Economic Studies, 23, 1–16.

Wiener N. (1923) Differential-space. Journal of Mathematics and Physics, 2, 131–174.

Other

H. Hamedani, Brownian Motion as the Limit of a Symmetric Random Walk, ProbabilityCourse.com Online chapter section.

About the author

The article was written in January 2026 by Saral BINDAL (Indian Institute of Technology Kharagpur, Metallurgical and Materials Engineering, 2024-2028 & Research assistant at ESSEC Business School).

   ▶ Read all posts written by Saral BINDAL.

Implied Volatility

Jayati WALIA

In this article, Jayati WALIA (ESSEC Business School, Grande Ecole Program – Master in Management, 2019-2022) explains how implied volatility is computed from option market prices and a option pricing model.

Introduction

Volatility is a measure of fluctuations observed in an asset’s returns over a period of time. The standard deviation of historical asset returns is one of the measures of volatility. In option pricing models like the Black-Scholes-Merton model, volatility corresponds to the volatility of the underlying asset’s return. It is a key component of the model because it is not directly observed in the market and cannot be directly computed. Moreover, volatility has a strong impact on the option value.

Mathematically, in a reverse way, implied volatility is the volatility of the underlying asset which gives the theoretical value of an option (as computed by Black-Scholes-Merton model) equal to the market price of that option.

Implied volatility is a forward-looking measure because it is a representation of expected price movements in an underlying asset in the future.

Computation methods for implied volatility

The Black-Scholes-Merton (BSM) model provides an analytical formula for the price of both a call option and a put option.

The value for a call option at time t is given by:

 Call option value

The value for a put option at time t is given by:

Put option value

where the parameters d1 and d2 are given by:,

call option d1 d2

with the following notations:

St : Price of the underlying asset at time t
t: Current date
T: Expiry date of the option
K: Strike price of the option
r: Risk-free interest rate
σ: Volatility of the underlying asset
N(.): Cumulative distribution function for a normal (Gaussian) distribution. It is the probability that a random variable is less or equal to its input (i.e. d₁ and d₂) for a normal distribution. Thus, 0 ≤ N(.) ≤ 1

From the BSM model, both for a call option and a put option, the option price is an increasing function of the volatility of the underlying asset: an increase in volatility will cause an increase in the option price.

Figures 1 and 2 below illustrate the relationship between the value of a call option and a put option and the level of volatility of the underlying asset according to the BSM model.

Figure 1. Call option value as a function of volatility.
Call option value as a function of volatility
Source: computation by the author (BSM model)

Figure 2. Put option value as a function of volatility.
Put option value as a function of volatility
Source: computation by the author (BSM model)

You can download below the Excel file for the computation of the value of a call option and a put option for different levels of volatility of the underlying asset according to the BSM model.

Excel file to compute the option value as a function of volatility

We can observe that the call and put option values are a monotonically increasing function of the volatility of the underlying asset. Then, for a given level of volatility, there is a unique value for the call option and a unique value for the put option. This implies that this function can be reversed; for a given value for the call option, there is a unique level of volatility, and similarly, for a given value for the put option, there is a unique level of volatility.

The BSM formula can be reverse-engineered to compute the implied volatility i.e., if we have the market price of the option, the market price of the underlying asset, the market risk-free rate, and the characteristics of the option (the expiration date and strike price), we can obtain the implied volatility of the underlying asset by inverting the BSM formula.

Example

Consider a call option with a strike price of 50 € and a time to maturity of 0.25 years. The market risk-free interest rate is 2% and the current price of the underlying asset is 50 €. Thus, the call option is ‘at-the-money’. If the market price of the call option is equal to 2 €, then the associated level of volatility (implied volatility) is equal to 18.83%.

You can download below the Excel file below to compute the implied volatility given the market price of a call option. The computation uses the Excel solver.

Excel file to compute implied volatility of an option

Volatility smile

Volatility smile is the name given to the plot of implied volatility against different strikes for options with the same time to maturity. According to the BSM model, it is a horizontal straight line as the model assumes that the volatility is constant (it does not depend on the option strike). However, in practice, we do not observe a horizontal straight line. The curve may be in the shape of the alphabet ‘U’ or a ‘smile’ which is the usual term used to refer to the observed function of implied volatility.

Figure 3 below depicts the volatility smile for call options on the Apple stock on May 13, 2022.

Figure 3. Volatility smile for call options on Apple stock.
Apple volatility smile
Source: Computation by author.

Excel file for implied volatility from Apple stock option

We can also observe that the for a specific time to maturity, the implied volatility is minimum when the option is at-the-money.

Volatility surface

An essential assumption of the BSM model is that the returns of the underlying asset follow geometric Brownian motion (corresponding to log-normal distribution for the price at a given point in time) and the volatility of the underlying asset price remains constant over time until the expiration date. Thus theoretically, for a constant time to maturity, the plot of implied volatility and strike price would be a horizontal straight line corresponding to a constant value for volatility.

Volatility surface is obtained when values for implied volatilities are calculated for options with different strike prices and times to maturity.

CBOE Volatility Index

The Chicago Board Options Exchange publishes the renowned Volatility Index (also known as VIX) which is an index based on the implied volatility of 30-day option contracts on the S&P 500 index. It is also called the ‘fear gauge’ and it is a representation of the market outlook for volatility for the next 30 days.

Related posts on the SimTrade blog

   ▶ All posts about Options

   ▶ Akshit GUPTA Options

   ▶ Jayati WALIA Brownian Motion in Finance

   ▶ Jayati WALIA Brownian Motion in Finance

   ▶ Youssef LOURAOUI Minimum Volatility Factor

   ▶ Youssef LOURAOUI VIX index

Useful resources

Academic articles

Black F. and M. Scholes (1973) “The Pricing of Options and Corporate Liabilities” The Journal of Political Economy, 81, 637-654.

Dupire B. (1994). “Pricing with a Smile” Risk Magazine 7, 18-20.

Merton R.C. (1973) “Theory of Rational Option Pricing” Bell Journal of Economics, 4, 141–183.

Business

CBOE Volatility Index (VIX)

CBOE VIX tradable products

About the author

The article was written in May 2022 by Jayati WALIA (ESSEC Business School, Grande Ecole Program – Master in Management, 2019-2022).

Monte Carlo simulation method

Jayati WALIA

In this article, Jayati WALIA (ESSEC Business School, Grande Ecole – Master in Management, 2019-2022) explains the Monte Carlo simulation method and its applications in finance.

Introduction

Monte Carlo simulations are a broad class of computational algorithms that rely majorly on repeated random sampling to obtain numerical results. The underlying concept is to model the multiple possible outcomes of an uncertain event. It is a technique used to understand the impact of risk and uncertainty in prediction and forecasting models.

The Monte Carlo method was invented by John von Neumann (Hungarian-American mathematician and computer scientist) and Stanislaw Ulam (Polish mathematician) during World War II to improve decision making under uncertain conditions. It is named after the popular gambling destination Monte Carlo, located in Monaco and home to many famous casinos. This is because the random outcomes in the Monte Carlo modeling technique can be compared to games like roulette, dice and slot machines. In his autobiography, ‘Adventures of a Mathematician’, Ulam mentions that the method was named in honor of his uncle, who was a gambler.

How Monte Carlo simulation works

The main idea is to repeatedly run a large number of simulations of a random process for a variable of interest (such as an asset price in finance) covering a wide range of possible situations. The outcomes of this variables are drawn from a pre-specified probability distribution that is assumed to be known, including the analytical function and its parameters. Thus, Monte Carlo simulations inherently try to recreate the entire distribution of asset prices.

Example: Apple stock

Consider the Apple stock as our asset of interest for which we will generate stock prices according to the Monte Carlo simulation method.

The first step in the simulation is choosing a stochastic model for the behavior of our random variable (the Apple stock price in our case). A commonly used model is the geometric Brownian motion (GBM) model. The model assumes that future asset price changes are uncorrelated over time and the probability distribution function of the future price is a log-normal distribution. The movements in price in GBM process can be expressed as:

img_SimTrade_GBM_process

with dS being the change in asset price in continuous time dt. dW is the Wiener process (Wt+1 – Wt is a random variable from the normal distribution N(0, 1)). σ represents the price volatility considering the unexpected changes that can result from external effects (σ is assumed to be constant over time). μdt together represents the deterministic return within the time interval with μ representing the growth rate of the asset price or the ‘drift’.

Integrating dS/S over a finite interval, we get :

img_SimTrade_simulated_asset_price

Where ε is a random number generated from a normal distribution N(0,1).

This equation thus gives us the evolution of the asset price from a simulated model from day t-1 to day t.

We can now generate a simulation path for 100 days using the above formula.

The figure below shows five simulations for the price of the Apple stock over 100 days with Δt = 1 day. The initial price for Apple stock (i.e, price at t=0) is $146.52.

Figure 1. Simulated Apple stock prices according to the Monte Carlo simulation method.
img_SimTrade_Apple_MonteCarloSim
Source: computation by author.

Thus, we can observe that the prices obtained by just these five simulations range from $100 to over $220.

You can download below the Excel file for generating Monte Carlo Simulations for Apple stock.

 Download the Excel file for generating Monte Carlo Simulations for Apple stock

Applications in finance

The Monte Carlo simulation method is widely used in finance for valuation and risk analysis purposes.

One popular application is option pricing. For option contracts with complicated features (such as Asian options) or those with a combination of assets as their underlying, Monte Carlo simulations help generate multiple potential payoff scenarios for the option which are averaged out to determine the option price at the issuance date.

The Monte Carlo method is also used to assess potential risks by generating simulations of market variables affecting portfolios such as asset returns, interest rates, macroeconomic factors, etc. over different time periods. These simulations are then assessed as required for risk modelling and to compute risk metrics such as the value at Risk (VaR) of a position.

Other applications include personal finance planning and corporate project finance where simulations are generated to construct stochastic financial models for sensitivity analysis and net present value (NPV) projections.

Related posts on the SimTrade blog

   ▶ Jayati WALIA Quantitative Risk Management

   ▶ Jayati WALIA Brownian Motion in Finance

   ▶ Jayati WALIA The Monte Carlo simulation method for VaR calculation

   ▶ Shengyu ZHENG Pricing barrier options with simulations and sensitivity analysis with Greeks

Useful resources

Hull, J.(2008) Risk Management and Financial Institutions, Fifth Edition, Chapter 7 – Valuation and Scenario Analysis.

About the author

The article was written in March 2022 by Jayati WALIA (ESSEC Business School, Grande Ecole Program – Master in Management, 2019-2022).

Black-Scholes-Merton option pricing model

Jayati WALIA

In this article, Jayati WALIA (ESSEC Business School, Grande Ecole Program – Master in Management, 2019-2022) explains the Black-Scholes-Merton model to price options.

The Black-Scholes-Merton model (or the BSM model) is the world’s most popular option pricing model. Developed in the beginning of the 1970s, this model introduced to the world, a mathematical way of pricing options. Its success was essentially a starting point for new forms of financial derivatives in the knowledge that they could be priced accurately using the ideas and analyses pioneered by Black, Scholes and Merton and it set the foundation for the flourishing of modern quantitative finance. Myron Scholes and Robert Merton were awarded the Nobel Prize for their work on option pricing in 1997. Unfortunately, Fischer Black had died several years earlier but would certainly have been included in the prize had he been alive, and he was also listed as a contributor by Scholes and Merton.

Today, the Black-Scholes-Merton formula is widely used by traders in investment banks to price and hedge option contracts. Options are used by investors to hedge their portfolios to manage their risks.

Assumptions of the BSM Model

As any model, the BSM model relies on a set of assumptions:

  • The model considers European options, which we can only be exercised at their expiration date.
  • The price of the underlying asset follows a geometric Brownian motion (corresponding to log-normal distribution for the price at a given point in time).
  • The risk-free rate remains constant over time until the expiration date.
  • The volatility of the underlying asset price remains constant over time until the expiration date.
  • There are no dividend payments on the underlying asset.
  • There are no transaction costs on the underlying asset.
  • There are no arbitrage opportunities.

The BSM equation

The value of an option is a function of the price of the underlying stock and its statistical behavior over the life of the option.

A commonly used model is Geometric Brownian Motion (GBM). GBM assumes that future asset price differences are uncorrelated over time and the probability distribution function of the future prices is a log-normal distribution (or equivalently the probability distribution function of the future returns is a normal distribution). The price movements in a GBM process can be expressed as:

GBM equation

with dS being the change in the underlying asset price in continuous time dt and dX the random variable from the normal distribution (N(0, 1) or Wiener process). σ is the volatility of the underlying asset price (it is assumed to be constant). μdt represents the deterministic return within the time interval with μ representing the growth rate of asset price or the ‘drift’.

Therefore, option price is determined by these parameters that describe the process followed by the asset price over a period of time. The Black-Scholes-Merton equation governs the price evolution of European stock options in financial markets. It is a linear parabolic partial differential equation (PDE) and is expressed as:

BSM model equation

Where V is the value of the option (as a function of two variables: the price of the underlying asset S and time t), r is the risk-free interest rate (think of it as the interest rate which you would receive from a government debt or similar debt securities) and σ is the volatility of the log returns of the underlying security (say stocks).

The key idea behind the equation is to hedge the option and limit exposure to market risk posed by the asset. This is achieved by a strategy known as ‘delta hedging’ and it involves replicating the option through an equivalent portfolio with positions in the underlying asset and a risk-free asset in the right way so as to eliminate risk.

Thus, from the BSM equation we can derive the BSM formulae that describe the price of call and put options over their life time.

The BSM formulae

Note that the type of option we are valuing (call or put), the strike price and the maturity date do not appear in the above BSM equation. These elements only appear in the ‘final condition’ i.e., the option value at maturity, called the payoff function.

For a call option, the payoff C is given by:

CT = max⁡(ST – K; 0)

For a put option, the payoff is given by:

PT = max⁡(K – ST; 0)

The BSM formula is a solution to the BSM equation, given the boundary conditions (given by the payoff equations above). It calculates the price at time t for both a call and a put option.

The value for a call option at time t is given by:

Call option value equation

The value for a put option at time t is given by:

Put option value equation

where

With the notations:
St: Price of the underlying asset at time t
t: Current date
T: Expiry date of the option
K: Strike price of the option
r: Risk-free interest rate
σ: Volatility (the standard deviation of the return on the underlying asset)
N(.): Cumulative distribution function for a normal (Gaussian) distribution. It is the probability that a random variable is less or equal to its input (i.e. d₁ and d₂) for a normal distribution. Thus, 0 ≤ N(.) ≤ 1

Figure 1 gives the graphical representation of the value of a call option at time t as a function of the price of the underlying asset at time t as given by the BSM formula. The strike price for the call option is 50€ with a maturity of 0.25 years and volatility of 50% in the underlying.

Figure 1. Call option value
Call option value
Source: computation by author.

Figure 2 gives the graphical representation of the value of a put option at time t as a function of the price of the underlying asset at time t as given by the BSM formula. The strike price for the put option is 50€ with a maturity of 0.25 years and volatility of 50% in the underlying.

Figure 2. Put option valuePut option value
Source: computation by author.

You can download below the Excel file for option pricing with the BSM Model.

Download the Excel file for option pricing with the BSM Model

Some Criticisms and Limitations

American options

The Black-Scholes-Merton model was initially developed for European options. This is a limitation of the equation for American options which can be exercised at any time before the expiry date. The BSM model would then not accurately determine the option value (an important case when the underlying asset pays a discrete dividend).

Stocks paying dividends

Also, in reality, most stocks pay dividends, and no dividends was an assumption in the initial BSM model, which analysts now eliminated by accommodating the dividend yield in the formula if required.

Constant volatility

Another limitation is the use of constant volatility. Volatility is the measure of risk based on the standard deviation of the return on the underlying asset. In reality the value of an asset will change randomly, not with a specific constant pattern regarding the way it can change.

Finally, the assumption of no transaction cost neglects the liquidity risk in the market since transaction costs are clearly incurred in the real world and there exists a bid-offer spread on most underlying assets. For the most heavily traded stocks, this cost may be low but for others it may lead to an inaccuracy.

Related posts on the SimTrade blog

All posts about Options

▶ Jayati WALIA Brownian Motion in Finance

▶ Akshit GUPTA Options

▶ Akshit GUPTA The Black-Scholes-Merton model

▶ Akshit GUPTA History of options market

Useful resources

Black F. and M. Scholes (1973) The Pricing of Options and Corporate Liabilities The Journal of Political Economy 81, 637-654.

Merton R.C. (1973) Theory of Rational Option Pricing Bell Journal of Economics 4, 141–183.

About the author

The article was written in March 2022 by Jayati WALIA (ESSEC Business School, Grande Ecole Program – Master in Management, 2019-2022).

Brownian Motion in Finance

Jayati WALIA

In this article, Jayati WALIA (ESSEC Business School, Grande Ecole Program – Master in Management, 2019-2022) explains the Brownian motion and its applications in finance to model asset prices like stocks traded in financial markets.

Introduction

Stock price movements form a random pattern. The prices fluctuate everyday resulting from market forces like supply and demand, company valuation and earnings, and economic factors like inflation, liquidity, demographics of country and investors, political developments, etc. Market participants try to anticipate stock prices using all these factors and contribute to make price movements random by their trading activities as the financial and economics worlds are constantly changing.

What is a Brownian Motion?

The Brownian motion was first introduced by botanist Robert Brown who observed the random movement of pollen particles due to water molecules under a microscope. It was in the 1900s that the French mathematician Louis Bachelier applied the concept of Brownian motion to asset price behavior for the first time, and this led to Brownian motion becoming one of the most important fundamental of modern quantitative finance. In Bachelier’s theory, price fluctuations observed over a small time period are independent of the current price along with historical behavior of price movements. Combining his assumptions with the Central Limit Theorem, he also deduces that the random behavior of prices can be said to be represented by a normal distribution (Gaussian distribution).

This led to the development of the Random Walk Hypothesis or Random Walk Theory, as it is known today in modern finance. A random walk is a statistical phenomenon wherein stock prices move randomly.

When the time step of a random walk is made infinitesimally small, the random walk becomes a Brownian motion.

Standard Brownian Motion

In context of financial stochastic processes, the Brownian motion is also described as the Wiener Process that is a continuous stochastic process with normally distributed increments. Using the Wiener process notation, an asset price model in continuous time can be expressed as:

brownian motion equation

with dS being the change in asset price in continuous time dt. dX is the random variable from the normal distribution (N(0, 1) or Wiener process). σ is assumed to be constant and represents the price volatility considering the unexpected changes that can result from external effects. μdt together represents the deterministic return within the time interval with μ representing the growth rate of asset price or the ‘drift’.

When the market is modeled with a standard Brownian Motion, the probability distribution function of the future price is a normal distribution.

Geometric Brownian Motion

weiner notation

with dS being the change in asset price in continuous time dt. dX is the random variable from the normal distribution (N(0, 1) or Wiener process). σ is assumed to be constant and represents the price volatility considering the unexpected changes that can result from external effects. μdt together represents the deterministic return within the time interval with μ representing the growth rate of asset price or the ‘drift’.

When the market is modeled with a geometric Brownian Motion, the probability distribution function of the future price is a log-normal distribution.

Properties of a Brownian Motion

  • Continuity: Brownian motion is the continuous time-limit of the discrete time random walk. It thus, has no discontinuities and is non-differential everywhere.
  • Finite: The time increments are scaled with the square root of the times steps such that the Brownian motion is finite and non-zero always.
  • Normality: Brownian motion is normally distributed with zero mean and non-zero standard deviation.
  • Martingale and Markov Property: Martingale property states that the conditional expectation of the future value of a stochastic process depends on the current value, given information about previous events. The Markov property instead focusses on the ‘no memory’ theory that the expected future value of a stochastic process does not depend on any past values except the current value. Brownian motion follows both these properties.

Simulating Random Walks for Stock Prices

In quantitative finance, a random walk can be simulated programmatically through coding languages. This is essential because these simulations can be used to represent potential future prices of assets and securities and work out problems like derivatives pricing and portfolio risk evaluation.

A very popular mathematical technique of doing this is through the Monte Carlo simulations. In option pricing, the Monte Carlo simulation method is used to generate multiple random walks depicting the price movements of the underlying, each with an associated simulated payoff for the option. These payoffs are discounted back to the present value and the average of these discounted values is set as the option price. Similarly, it can be used for pricing other derivatives, but the Monte Carlo simulation method is more commonly used in portfolio and risk management.

For instance, consider Microsoft stock that has a current price of $258.65 with a growth trend of 55.2% and a volatility of 35.92%.

A plot of daily returns represented as a random normal distribution is:

Normal Distribution

The above figure represents the simulated price path according to the Geometric Brownian motion for the Microsoft stock price. Similarly, a plot of 10 such simulations would be like this:

Microsoft GBM Simulations

Thus, we can see that with just 10 simulations, the prices range from $100 to over $600. We can increase the number of simulations to expand the data set for analysis and use the results for derivatives pricing and many other financial applications.

Brownian motion and the efficient market hypothesis

If the market is efficient in the weak sense (as introduced by Fama (1970)), the current price incorporates all information contained in past prices and the best forecast of the future price is the current price. This is the case when the market price is modelled by a Brownian motion.

Related Posts

▶ Jayati WALIA Black-Scholes-Merton option pricing model

▶ Jayati WALIA Plain Vanilla Options

▶ Jayati WALIA Derivatives Market

▶ Saral BINDAL Modeling Asset Prices in Financial Markets: Arithmetic and Geometric Brownian Motions

Useful Resources

Academic articles

Fama E. (1970) Efficient Capital Markets: A Review of Theory and Empirical Work, Journal of Finance, 25, 383-417.

Fama E. (1991) Efficient Capital Markets: II Journal of Finance, 46, 1575-617.

Books

Malkiel B.G. (2020) A Random Walk Down Wall Street: The Time-tested Strategy for Successful Investing, WW Norton & Co.

Code

Python code for graphs and simulations

Brownian Motion

What is the random walk theory?

About the author

The article was written in August 2021 by Jayati WALIA (ESSEC Business School, Grande Ecole Program – Master in Management, 2019-2022).