The concept of risk-neutral pricing of European call options is investigated from a mathematical approach. It is found that risk-neutral pricing used in the pricing of European call options is a specious concept    that is only approximately correct and that ignores terms of , where is the risk premium and σ2 is the variance of the one-day return of the asset that underlies the call option. The risk premium equals with α the true drift rate of the underlying option and r the risk-free rate.
First some notation and background information. Let the probability density function (pdf) for a random variable be . The probability that a measurement of the random variable takes a value between and , , is given by . The cumulative density function (CDF) .
Let be the value of an asset at time t. Let the 1-day return be and the n-day return be . Note that
and, for simplicity in notation, that .
When returns over non-lapping time periods are independent, the pdf for the n-day return is the n-fold convolution of the pdf of the 1-day return since the n-day return is the sum of n independent 1-day returns. Variances add under convolution, and hence the variance of n-day returns is n times the variance that describes the distribution of the 1-day returns. Since the normal distribution is stable under self-convolution, if the 1-day return distribution is a normal distribution with a mean μ and a variance σ2, then the pdf for an n-day return is a normal distribution with a mean and with a variance  .
The returns of assets are better described by Student’s t-distributions than by normal distributions. However, Student’s t-distributions are not stable under self-convolution and have fat tails. The fat tails lead to integrals that diverge and these integrals are needed to price options. Both of these characteristics make pricing with Student’s t-distributions difficult      . For simplicity and to be specific, the pricing of European call options for returns with normal distributions is discussed in this note. In particular, attention is paid to the mathematical basis for risk-neutral pricing of European call options. The results apply to fat tailed distributions as well  .
Following   , the price of a European call option at time T is . In the expression for , is the expectation of the maximum value of 0 or the difference of the value of the asset and the strike price at time T. At time , when the option is purchased, the value of the option is as the sale of the option is a cash transaction and therefore the value of the option at time T, , is discounted by the time value of money, with r the risk-free rate.
The expression for follows from the arbitrage theorem      .
2.1. Black-Scholes Option Pricing Formula
The Black-Scholes option pricing formula gives prices for European call options and is obtained under the constraints of: 1) no arbitrage; 2) the price of an asset is described by a geometric Brownian motion process with support ; 3) risk-neutral pricing; and, 4) the future price is a martingale. Constraint 1) requires, from the arbitrage theorem    , that . Constraints 3) and 4) require where is the price of the asset at time t, is the price at time , and r is the risk-free rate. Given that is a geometric Brownian motion process derived from a stochastic process with mean and variance , i.e., the 1-day returns are normally distributed with drift rate μ and variance σ2, then the expectation of is given by
for the time origin chosen such that .
Since by assumption is geometric Brownian motion process and since the distribution of is an n-fold convolution of the 1-day return distribution, then the constraint specifies the drift rate for and the n-fold convolution dictates the shape and the variance of the distribution of , which then dictates the shape and variance of . From Equation (2) and constraints 3) and 4), must have a drift rate and from the n-fold convolution, must be normally distributed (by assumption the 1-day return is normally distributed) with a variance of , and when the time origin is chosen such that .
The value of a stock is expected to drift at the rate of where is called the risk premium and r is the risk-free rate.
2.2. Risk-Neutral Pricing
Assume that where is the T-day return and that the T-day return is normally distributed with mean and variance . Thus where .
The price of a European call option with strike price is then   , assuming normal statistics for the returns with support ,
The difference between the Black-Scholes option pricing formula and Equation (3) rests in the α in Equation (3). If one sets in Equation (3), then one obtains the Black-Scholes formula. The Black-Scholes formula is obtained by using the concept of risk-neutral pricing, wherein one essentially sets the risk premium, , equal to zero. Mathematically this approach can not be correct. Setting violates the no arbitrage condition that . The distribution for is centred about (i.e., the mean of the pdf for the n-day return is ). Arbitrarily setting the mean of the distribution will change the value of the expectation and hence mis-price European call options. In general, for small T, and the error introduced by arbitrarily setting the mean of the distribution will be small. Application of Girsanov’s theorem in risk-neutral pricing is discussed in Sec. 3.
A series expansion in α of about the risk-neutral value shows that risk-neutral pricing underestimates the price of a call option when . is the cost at time of a European call option whereas is the cost of a European call option using risk-neutral pricing (i.e., using the same assumptions that yield the Black-Scholes option pricing formula).
and and the pdf and CDF for a normal distribution with a mean of zero and a variance of σ2. See Equation (8) for definitions of and .
when since the expansion is only valid for , which follows from the definition . Thus the risk-neutral pricing underestimates the value of the call option and gives, on average, an advantage to the option buyer when . The converse holds. The seller has, on average, an advantage under risk-neutral pricing when . Neither party has, on average, an advantage when the true value of α is used to price an option. The best estimate of a true value is typically the sample mean.
Table 1 gives pricing of European call options for years, per annum, as measured over one-year, and for geometric Brownian motion with Gaussian increments. Note that risk-neutral pricing (i.e., the Black-Scholes formula) underestimates the price of a European call option when and overestimates the price of a European call option when .
The “success” of risk-neutral pricing owes to the fact that the magnitude of the random fluctuations are typically significantly greater than the magnitude of the risk premium, i.e., , and thus the random fluctuations
Table 1. The costs of European options, , for various strike prices and values of α, given per annum, , years, and . The column labelled corresponds to prices obtained via risk-neutral pricing, i.e., by the Black-Scholes formula.
obscure the drift. Presumably fluctuations about the mean value are described well by the shape and scale parameters of the distribution and one should use the best available estimate of the location parameter of the distribution (i.e., the mean drift rate, α) to price an option.
Consider a normal distribution with mean μ and variance σ2. If , then to an error of , μ can be ignored. This can be verified from a series expansion of the CDF
Note that has a factor of , that is the cumulative density function (CDF), and that is the zero mean pdf with variance σ2:
The standard deviation σ is a measure of the width of the distribution whereas the mean μ shifts the curve left or right. For a broad curve, small shifts left or right make a small difference.
In an Ito calculus formalism, risk-neutral pricing is explained as
where is an increment of Brownian motion (or Weiner process (  , p. 79) with , , and . is for a transformed process with non-zero mean such that . Provided that is negligible, then the process that underlies is to a good approximation equal to the process that underlies and the values of as calculated under the processes that yield and are essentially the same. If is not negligible, then the risk premium needs to be known to calculate the value of an option and risk-neutral pricing is not possible. It should be noted that strictly speaking is not a Brownian motion: and . Girsanov’s theorem is discussed in Sec. 3.
The solution to Equation (9) is not since . The solution to Equation (9) is, with ,
which follows from the solution, in a Langevin formalism   , of the equivalent stochastic differential equation to Equation (9):
where and is a zero mean stochastic process that in a limit is delta function correlated. In the limit, is a white noise and the Wiener process , (  , p. 79). Equation (11) is equivalent to the first line of Equation (9). In the Langevin formalism, the equation for the development in time of the average value of , , is
In the event that α is a constant, . The mean value of (or , ) drifts at the rate α when the time development of (or ) is described by Equation (9) in an Ito formulation or equivalently by Equation (11) in a Langevin formulation.
Note that the contribution to the drift arises from averaging over an ensemble of realizations of the stochastic process (c.f. Equations (2), (10), and (13)). Care must be employed in obtaining and in interpreting results within the Ito formalism. One could attempt, in the Langevin picture, to hide the risk premium in , as was attempted in Equation (9) to justify risk-neutral pricing. However, one would experience a similar difficulty. The transformed would not be zero mean and would not be delta function correlated in the limit, and knowledge of would still be required to price the option unless is negligible.
Langevin equations are first order differential equations with noise driving terms. The Langevin equations should be interpreted as integral equations (  , p. 172;  , Ch. 10.2;  , p. 79). Average values found by the Langevin approach are identical to solutions found by Ito’s calculus (  , p. 189). The Langevin approach has the advantage that transformations obey the usual rules of calculus. Ito’s lemma is not applied for each transformation in the Langevin approach (  , pp. 189, 282).
3. Girsanov’s Theorem
Girsanov’s theorem  provides for a multiplicative transformation that alters the shape of the pdf. This allows the mean of a random variable to be set to an arbitrary value, without using an additive correction to the original problem or pdf (  , pp. 36-39). In option pricing, the goal behind the multiplicative transformation is an economic one: to find a risk-neutral measure to ensure no risk-free arbitrages in pricing (  , p. 39). The multiplicative transformation must decrease probabilities for values greater than the desired mean and increase probabilities for values that are less than the desired mean.
3.1. Probability Measures
A pdf defines a probability measure P. If the pdf for a random variable is , then
If is a small neighbourhood of a specific outcome x, then . For a function of the random variable , then the expectation of over is
Two probability measures P and Q are said to be equivalent if, for any set A in the probability space, and (  , p. 245). This means that all events that are possible under the measure P are possible under the measure Q, and vice versa. If for some set B, and , then the event B is impossible in measure P whereas event B is possible in measure Q, and the two measures P and Q are not equivalent. The possibilities that not-equivalent measures describe are different.
3.2. Transformed Stochastic Process 
Consider a stochastic process (for a stochastic process, is a random variable for each point in time t) that has a drift rate μ and is driven by a Wiener process , such that, e.g., . Define (  , pp. 36-39 and pp. 210-214)
and P are equivalent probability measures and are related by (  , pp. 36-39 and pp. 210-214)
where is the equivalent risk-neutral measure. is called the Radon-Nikodým derivative.
If dP is the unit-normal distribution (i.e., the underlying pdf is normally distributed with zero mean and standard deviation of unity, ) then
and the equivalent risk-neutral measure for this example, for a small neighbourhood near , is
Clearly Z is selected to give the desired mean to the equivalent measure. In the example given here, is a unit normal distribution with a mean of . A change of variable starts the transformation of the process to a zero mean process .
The transformed process for the example given is where could be the risk-free rate and is a Gaussian increment at time t of a Brownian motion. In the equivalent measure , is a geometric Brownian motion that increases on average at the risk-free rate : .
Gardiner (  , p. 245) writes “can show that two stochastic differential equations can be considered equivalent if their noise terms are the same even if their drift terms are different”. See Equations (20) below for examples of stochastic differential equations with different drift and noise terms.
Gardiner (  , p. 246) writes “the possible sample paths from the two equations are identical, but depending on the choice of measure for the underlying driving process the relative frequency of the paths is different”. Emphasis added. All paths are available in the two measures P and , but different paths and hence end points are emphasized in the two measures.
Gardiner (  , p. 246) writes that “the stochastic differential equations
are equivalent if . This result is Girsanov’s theorem.”
Gardiner (  , p. 247) shows simulations of a stochastic differential equation. The simulations for the same noise but different drifts look qualitatively indistinguishable whereas the simulation with different noise looks different than the other simulations with same noise but different drifts: “it is quite credible that either could be a simulation of the other equation.” Qualitatively indistinguishable appears to mean the rms noises about the local trend lines appear to be similar.
This result is not surprising. The solution to Equation (9) or to Equation (11) in a Langevin approach is given by Equation (10) and can be recast as
is a scaled version of and is determined solely by the noise , . Provided that the drift is not too large compared to σ, then any with the same σ but different drift should look similar. In this case, one is using the higher frequency noise as a fiducial to compare observations of for for different and/or .
Figure 1 presents simulations of for 200 time steps (200 days) for and for , with per annum. The smooth lines are the expected values. All simulations used . Note the similarity between the simulations with different drifts. Note also the difference between the simulations and the expected means for the different drift rates: a simulation with (red curves) gives values that are less than the values obtained with (blue curves).
Figure 2 presents simulations of for 1000 time steps (1000 days) for and for , with per annum. The smooth lines are the
Figure 1. Simulations of for per annum (red curves) and for , (blue curves) with given by Equation (10). The smooth curves are the expected values , Equation (13). All calculations in this figure used .
Figure 2. Simulations of for per annum (red curves) and for , (blue curves) with given by Equation (10). The smooth curves are the expected values , Equation (13). All calculations in this figure used .
expected values. All simulations used . Figure 1 shows the first 201 points of Figure 2. The sharp rise in owes to the sequence of random draws. A different seed for the pseudo random number generator gives drastically different results.
Figure 2 is included to show the drift possible (observe the trend from to ) owing to random draws. The noise had a mean of 0.06 over the interval to ; 0.026 over to ; 0.30 from to ; and, −0.007 from to .
The simulations for the same noise but different drifts look qualitatively indistinguishable, as pointed out by Gardiner (  , p. 247).
Gardiner (  , p. 247) writes:
“Girsanov’s theorem is now the justification for use of the drift rate r instead of μ in the valuation of options using the risk-neutral procedure. The noise term is identical for both cases, and in the case we can say that the two processes can be seen as arising from the choice of a different probability measure to the same set of sample paths. In some sense it can be shown that this is a rigorously justifiable procedure [10.11], although not everyone would accept that. However, the use of change of measure is now an accepted part of the procedure for valuing options and other derivatives when one goes beyond the simple geometric Brownian motion picture.”
Gardiner does not seem to be a true believer, and seems resigned to a deeply engrained status quo.
From  :
“The relation (5) has been called by Cox-Ingersoll-Ross (1981) the “Local Expectation Hypothesis”, a terminology which has led to some confusion. Note that the equilibrium process has not been changed, it is the same under both measures P and . Girsanov’s Theorem allows us to replace the relation (4) through the equivalent simpler relation (5). In particular, no assumption has been made about the existence of risk-neutral investors. In a real economy neither a “representative” nor a “risk-neutral” investor will exist, since both assumptions would prevent the existence of a (stable) equilibrium. The great advantage of the representation (5) under the (martingale) measure is that we do not have to know anything about the individual expectations P and the investors’ attitude towards risk. In summary: has not been changed. It is the same equilibrium price process as under P, but in simpler representation under . is called the “equivalent risk-neutral measure” or the “P-equivalent martingale measure”. (  , p. 72)
It would appear that has been changed. The possible paths are the same, but more weight (probability) has been placed on lower yielding paths (assuming the risk premium is greater than zero). From  :
“After the change of measure, we are still considering the same set of stock price paths, but we have shifted the probability on them. If , as it normally is, then the change of measure puts more probability on the paths with lower return so that the overall mean rate of return is reduced from α to r.” (  , p. 217) “In finance, the change from the actual to the risk-neutral probability measure changes the distribution of asset prices without changing the asset prices themselves, ...” (  , p. 37)
Altering the distribution changes the problem, unless the alteration is undone by an inverse transformation to return to the original frame of reference. Essentially, the risk-neutral approach appears to be to multiply one term on one side of an equation by Z. This is not a valid mathematical approach. Consider with solution . Now solve . The solutions for are not the same unless .
In the development of the risk-neutral measure in Equation (9), it is considered that is a Weiner process (Brownian motion). Strictly speaking, Brownian motion is a zero mean process. is not a zero mean process. One might wish to apply a coordinate transformation such that is a zero mean Brownian motion in the transformed frame of reference. However, one must remember that one is working in a transformed coordinate system, and provide a reverse transformation at the end to obtain the answer in the original frame of reference. This is similar to the problem of relative motion.
Consider an airplane that can cruise at v km/hour with respect to the air mass in which the airplane is embedded (i.e., the local air) and assume that the local air is moving relative to the ground. If one is interested in the location of the plane after a given time, one can solve the problem by using a coordinate system that is embedded in the moving air. Relative to this coordinate system the plane has travelled a distance in time t. The choice of coordinate system makes the problem look simpler. However, to know the location of the plane relative to a reference coordinate system on the ground such as an airport, one must know the relationship between the reference coordinate system and the moving coordinate system. In a similar manner, to find the value of an asset or an option, one must know the risk premium.
3.4. A Thought Experiment
Consider pricing a very long lived option, one so long lived that the mean value is . With 99.99% certainty (normal statistics are assumed in this work), the value of the underlying will be . Does it make sense to use risk-neutral pricing and force ? In this case, the approximation that does not hold, and it appears that risk-neutral pricing would not be a reasonable approach.
3.5. Apply Girsanov’s Theorem
Let us examine the justification for risk-neutral pricing, which is presented as Equation (9). As before, define a risk premium, move it to the noise term , and absorb the risk premium in .
Now is a non-zero mean noise term, except in the special case that . Multiply the noise term by , c.f. Equation (16), to obtain
is chosen such that the mean value of is zero, i.e., .
One could redefine on the right hand side of Equation (23) to obtain
which is a simple-looking equation but it must be remembered that is no longer a constant. The dependence of on undoes the work to create a zero mean noise term . Alternatively, one could multiply both sides of Equation (23) by Z, follow the rules for transformation of stochastic differential equations and variables (  , Sec. 10.2-10.4, p. 275), only to find the same result that the work undid the desired result of setting the mean of to zero.
If in Equation (9) or in Equation (22) can be ignored, then under this approximation risk-neutral pricing would be accurate.
3.6. A Scaled Approach
Equation (21) suggests an approach to understand risk-neutral pricing of a European call option. Start with the expression for the value of the call option and manipulate to remove the drift owing to the risk premium in :
where both and have been scaled by the same factor to remove drift in owing to the risk premium. The expectation in the last line might be anticipated at first look to be risk neutral. However, is a function of the
risk premium: . In addition, the inverse scaling
factor remains in the expression for as an overall multiplicative factor. The risk premium is thus required to price the option. If the risk premium is negligible, then the scaling factor is approximately unity and risk-neutral pricing of the European call option would be sufficiently accurate. The approach presented in this section is silent on the magnitude of the risk premium that is negligible―one would need to examine the expectation, as was done in Sec. 2.2, to determine the magnitude. The rms fluctuations of the one day returns, σ, is the relevant metric. It is the magnitude of the ratio that determines what is negligible or not.
Risk-neutral pricing of European call options is mathematically an approximation. Provided that is small, then the drift rate is obscured over short time intervals by random fluctuations and one is justified in ignoring the risk premium in the drift rate. It would seem honest to state this as the rationale behind risk-neutral pricing, rather than appealing to theorems and pretending that risk-neutral pricing is exact.
Risk-neutral pricing underestimates the price of a European call option when and overestimates the price of a European call option when .
It is interesting to note that risk-neutral pricing of European call options ignores but takes great care to include in the average drift rate. If the risk premium can be ignored, then likely can also be ignored. For the S & P 500, on average whereas the risk premium .
1Specious “refers to something that appears at first encounter to be genuine or to be soundly argued or reasoned”  . A specious argument is “an argument that seems correct only if you do not think about it carefully”  . Something specious is “seemingly well-reasoned, plausible or true, but actually fallacious”  .
 Bouchaud, J.-P. and Sornette, D. (1994) The Black-Scholes Option Pricing Problem in Mathematical Finance: Generalization and Extensions for a Large Class of Stochastic Processes. Journal de Physique I France, 4, 863-881.
 Cassidy, D.T. (2012) Effective Truncation of a Student’s t-Distribution by Truncation of the Chi Distribution in a Mixing Integral. Open Journal of Statistics, 2, 519-525.
 Cassidy, D.T., Hamp, M.J. and Ouyed, R. (2010) Pricing European Options with a Log Student’s t-Distribution: A Gosset Formula. Physica A, 389, 5736-5748.
 Coffey, W.T., Kalmykov, Yu.P. and Waldron, J.T. (2004) The Langevin Equation: With Applications to Stochastic Problems in Physics, Chemistry, and Electrical Engineering. 2nd Edition, World Scientific Publishing, Singapore.
 Girsanov, I.V. (1960) On Transforming a Certain Class of Stochastic Processes by Absolutely Continuous Substitution of Measures. Theory of Probability and Its Applications, 5, 285-301.