Modern Economy, 2017, 8, 878-887 
http://www.scirp.org/journal/me 
ISSN Online: 2152-7261 
ISSN Print: 2152-7245 
 
 
 
The Effects of Negative Nominal Rates on the 
Pricing of American Calls: Some Theoretical 
and Numerical Insights 
Alessia Cafferata1, Pier Giuseppe Giribone2, Marina Resta1 
1DIEC, University of Genova, Genova, Italy 
2Banca Carige, Genova, Italy 
 
 
 
How  to  cite  this  paper:  Cafferata,  A., 
Giribone,  P.G.  and  Resta,  M.  (2017)  The 
Effects  of  Negative  Nominal  Rates  on  the 
Pricing of American Calls: Some Theoreti-
cal  and  Numerical  Insights. Modern Eco- 
nomy, 8, 878-887. 
https://doi.org/10.4236/me.2017.87061   
   
Received: May 29, 2017 
Accepted: July 10, 2017 
Published: July 13, 2017 
 
Copyright © 2017 by authors and   
Scientific Research Publishing Inc. 
This work is licensed under the Creative 
Commons Attribution International   
License (CC BY 4.0). 
http://creativecommons.org/licenses/by/4.0/   
  
Open Access
 
Abstract 
The  article  investigates  the  effects  played  on  options  pricing  by  negative 
risk-free rates when the underlying is an equity with null dividends. In such 
anomalous conditions, in fact, the fair value at early exercise of the American 
Call would not match the value of the European Call with the same financial 
features. We originally motivate this assumption with theoretical arguments. 
We  then  move  to  an  empirical  investigation  where  we  put  at  work  some 
quasi-closed formulas for pricing an American option and the stochastic tri-
nomial trees algorithm. We then draw the conclusion that from a numerical 
viewpoint, the bias between the fair value of the American Call and the value 
of the corresponding. European Call is mainly due to approximation errors, 
which can be mitigated when Trinomial Stochastic Trees are used. 
 
Keywords 
American Options, Quasi-Closed Formulas, Negative Interest Rates, Stochastic 
Trinomial Trees 
 
1. Introduction 
As  outlined  in  a  recent  note  from  the  Actuarial  Association  of  Europe  [1], 
nowadays negative nominal interest rates for long term maturities are observable 
in both European and American financial markets. In addition to the economic 
effects, this led to several technical problems, as the existing pricing models do 
not give proper valuations, so that either the financial position cannot be cor-
rectly priced or the results can be questioned [2].   
From the financial standpoint, it is therefore necessary to check to what extent 
the  existing  pricing  models  can  be  adapted  to  incorporate  negative  nominal 
DOI: 10.4236/me.2017.87061    July 13, 2017 
A. Cafferata et al. 
 
rates. This aspect has been already investigated in some research papers: [3] and 
[4] discuss the issue for options written on interest rates, both from the practical 
and the theoretical viewpoint; [5], focusing on foreign exchange and index op-
tions investigate whether the use of models allowing for negative interest rates 
can  improve  option  pricing  and  implied  volatility  forecasting;  [6],  discusses  a 
new closed form  for option pricing that leads to sensitively lower the error in 
European options pricing. Besides, [7] adapts the Nelson-Siegel model [8] to in-
clude the negative interest. Finally, the Hull and White model [9], has been re-
cently adapted to calibrate in a more proper way when the underlying is a nega-
tive interest rate [10]; however, to the best of our knowledge, much less efforts 
have been devoted to model the effects of negative nominal interest rates in op-
tion pricing for other types of underlying. 
In particular, some issues might arise in the case of equity that does not pay 
dividends: finding the fair value at early exercise of an American Call might be 
tricky, as it could not match the value of a European Call option with the same 
parameters. The problem is relevant, because of its corporate implications, as the 
option evaluation could make the difference when valuing a firm.   
This paper aims to fill in this gap. Our research question is discussing whether 
the known approximation formulas can effectively bypass the above highlighted 
problems. We illustrate an empirical application, where we compare the estima-
tion of a number of quasi-closed formulas, with that provided by the stochastic 
trinomial trees algorithm, and we highlight how the bias between the fair value 
of the American Call at early exercise and the value of the corresponding Euro-
pean Call can be strongly mitigated using this latter methodology. 
The paper is structured as follows. Section 2 starts by providing a snapshot- 
wise demonstration of why the equivalence between the fair value of the Ameri-
can Call at early exercise and the European Call can be violated, for an underly-
ing  with  null  dividends.  The  section  then  contains  a  brief  overview  about  the 
approximation schemes employed in the work. Section 3 illustrates the numeri-
cal case study with discussion. Section 4 concludes. 
2. Theoretical Issues and Methodology   
2.1. On the Violation of the Equivalence between American   
and European Call Value 
,
C
=
f
(
)
S K T r q σ
,
,
,
,
A
E
E
A
  the value of  an  American Call with 
Let us denote by 
spot price  S , strike price  K , time to expiration  T , interest rate  r , dividend 
yield  q ,  and  volatility  σ.  We  also  denote  by 
  the 
corresponding  value  for  a  European  contingent  claim.  We  focus  on  the  well- 
known property [11] according to which in case of an underlying with null divi-
dends we have: 
)
S K T r q σ
,
C
=
(
f
,
,
,
,
 
f
A
(
S K T r
,
,
,
,0,
)
σ
=
f
E
(
S K T r
,
,
,
,0,
)
σ
.                            (1) 
In the case of negative interest rate, (1) might not be satisfied. The value of a 
European Call Option, in fact, cannot be lower than the difference between the 
879 
A. Cafferata et al. 
spot price and the actual value of the strike price: 
)
σ
S T K r
,
,
S T K r
,
,
)
σ
,0,
,0,
≥
(
(
f
f
,
,
A
E
 
≥
S K
t
−
exp
(
−
rT
)
.              (2) 
Besides: 
f
A
(
S T K r
,
,
,
,0,
)
σ ≥
[
max 0,
S K
t
−
]
,                              (3) 
Joining (2) and (3) we get: 
S K S K
t
−
<
−
t
exp
(
−
rT
)
.                                        (4) 
(4)  clearly  holds  if 
r < ,  (4)  is  no  more  consistent,  because  the 
term in the right-hand side might be either is negative or lower than the value in 
the left-hand side, as: 
r > .  If 
1
> .   
rT−
exp
0
0
(
)
2.2. Methodology 
Pricing American contingent claims has traditionally represented a stimulating 
field of analysis as, in contradistinction to European options, they can be exer-
cised at any time before or at maturity. In this case, the Black-Scholes method-
ology cannot be applied, and it is necessary to use approximations schemes. Re-
viewing the related literature requires paramount efforts, besides it is out of the 
scope of this work: the interested reader can refer to [11].   
Nevertheless,  we  are  mainly  concerned  with  two  sub-groups  of  the  above 
methods. In the first group, we consider three quasi-closed formulas that con- 
veys in different ways the original idea discussed in [12]. In particular, the Bar-
one-Adesi and Whaley-BAW-model [13] is a quadratic approximation method 
for pricing exchange-traded American call and put options on commodities and 
commodity  futures.  Using  the  same  notational  conventions  as  in  Sec.2.1,  we 
consider an American Call option whose underlying has a cost of carry equal to 
r≥ , ceteris paribus the value of the American Call is equal to 
b
that of the European Call so that the Generalized Black-Scholes-GBS-formula for 
European contingent claims applies:   
= − . When  b
r q
=
)
C
(
GBS
E
+
(
A S S
2
*
y
) 2
,
S
<
*
S
; and 
AC
= −
S K S
,
≥
*
S
, 
  is the value of the European Call according to the GBS formula, 
A
2
=
*
S
{
1
−
(
N d S
1
*
)
exp
(
−
)
b r T
}
y
2
; 
C
A
)GBS
(
EC
where 
and: 
with: 
y
2
Finally, 
*
S
(
(
(
+
−
b
2
b
2
2
σ
2
σ
)
1
)
1
− +
= −
*S   is the price level such that: 
{
1 exp
−
K C S K T r q
,
)
,
σ
−
=
+
S
(
8
r
,
*
,
,
E
2
*
2
σ
)
1 exp
−
(
rT
)
2
. 
(
−
(
qT N d S
)
1
*
)
}
y
2
.        (5) 
The  Newton-Raphson  algorithm  can  be  then  used  to  solve  (5)  with  initial 
value: 
880 
A. Cafferata et al. 
 
S
*
START
=
K
+
∞
*
S
−
K
1 exp
−
(
h
2
)
, 
h
2
= −
(
K bT
+
2
σ
T
)
∞
*
S
−
K
 , 
S
*
∞
=
K
1 2
−
−
(
b
2
2
σ
)
1
− +
(
b
2
2
σ
−
2
)
1
+
8
r
2
σ
1
−
, 
1iS +   estimator is given by: 
)
−
(
K RHS S
=
+
S
i
1
+
i
b S
i
i
(
1
−
b
i
)
, 
where: 
and: 
so that the best 
where 
RHS S   is the right-hand side of (5) at the i-th step.   
(
)i
The second method is due to Bjerksund and Stensland—BS1993—and it is more 
general than the BAW, as the underlying can be a stock, a future or an exchange 
rate, and it is based on a feasible but non-optimal exercise strategy correspond-
ing to a trigger price I [14]. If S > I, it is optimal to exercise the option immedi-
ately, and the value must be equal to the intrinsic value S-K. On the other hand, 
if S ≤ I, it will never be optimal to exercise the American call option before expi-
ration, and the value can be found using the Black-Scholes formula. Finally, the 
third  approximation  method  is  due  to  Bjerksund  and  Stensland,  again  [15]— 
BS2002—, and it is based on the extension of the flat boundary concept by di-
viding  the  time  to  maturity  into  two  parts,  and  allowing  two  separate  flat 
boundaries in each of them. 
An alternative to the above-mentioned approximation methods is represented 
by stochastic binomial and trinomial trees. Assuming the stock price to follow a 
discrete time process, in the binomial tree scheme [16] the life of the option until 
the maturity T is decomposed into N time steps of equal length. At each time 
step,  the  underlying  will  move  either  up  or  down  by  a  specific  factor 
∆ ,  with  probability  p   and  1 p− ,  respec-
u
tively. The value of the America Call exercised at expiration is: 
(
N i N i
−
(
tσ
∆   or 
tσ
−
Su d
max
(
BIN
E
exp
exp
exp
i N i
−
∑
(
1
rT
i
p
u
p
u
,0
C
K
N i
−
(
)
)
{
)
}
!
=
=
(
−
−
!
−
=
d
)
)
)
)
 
N
i
=
0
)
(
!
} (
where: 
up
{
exp
(
To properly assess 
=
t
d
r q
− ∆ −
EC   in case of early exercise, at each node of the three the 
u d
−
. 
)
following pay-off must be applied: 
(
,exp
max
K
−
S
i
,
j
{
r t
− ∆
)
p f
u i
,
j
1
+
(
1
+ −
p
u
)
f
i
,
j
}
 , 
jf
,i
  is the value of the Call for the node of position (i,j) in the tree. The 
where 
initial value of the option can be then derived by way of the standard backward 
induction technique. A straightforward extension of this procedure is given by 
the trinomial scheme algorithm [17], with the underlying that can now assume 
three  different  states:  up,  down  or  unchanged.  The  increase  in  the  number  of 
possible states allows to lower the number of necessary steps for the convergence 
of the  procedure, without  any loss in the estimation  accuracy. The size of the 
881 
A. Cafferata et al. 
 
(
tσ
2
)
u
=
exp
jumps is usually set to: 
2
probability of reaching upward/downward branches is given by: 
(
∆ , and 
(
σ
exp
exp
exp
exp
exp
b t
∆
σ
−
t
∆
t
∆
)
(
)
=
−
−
=
2
2
2
d
(
{
σ
−
(
tσ
−
∆ , so that the 
)
up
{
)
)
(
σ
=
exp
dp
the probability 
(
2
−
b t
∆
exp
t
exp
∆
mp   of reaching the intermediate node is: 
exp
t
∆
−
2
2
)
)
(
−
σ
t
∆
p
m
1
= −
(
σ
)
2
t
∆
}2
−
)
p
u
2
p
d
 
}2
  while 
.   
3. Examples and Discussion 
We consider three scenarios (A, B, and C), and for each of them we compute the 
value of the American Call with the approximation schemes illustrated in Sec. 
2.2. In detail: 
•  A  represents  a  typical  market  situation,  with  a  positive  risk-free  rate,  and 
with a dividend-paying stock as underlying; 
•  B simulates a market situation with a positive risk-free rate, and with a null 
dividend stock as underlying: in this case, as 
r > , Equation (1) holds;   
0
•  C considers an atypical situation, with a negative risk-free rate. The underly-
ing stock, likewise in the B case, does not pay any dividend. 
The parameters employed in the simulation are reported in Table 1: we have 
used the annualized value of the volatility, while  T   is expressed as a fraction of 
the  year;  finally,  the  value  of  r   in  the  third  scenario  corresponds  to  the 
3-months value of the Euribor at 9 September 2016, as provided by Bloomberg. 
The simulation results are shown in Table 2, where we employed the follow-
ing  abbreviations:  BAW  to  indicate  the  Barone-Adesi  and  Whaley  model, 
BS1993  and  BS2002 referring  to  1993  and  2002  Bjerksund  and  Stensland  ap-
proximation formulas, respectively, and T-TREE for the trinomial tree. In this 
latter case, the discretization steps were set to N = 9000. 
Looking  at  Table  2,  several  remarks  come  out.  First,  in  the  scenario  A,  by 
construction, the early exercise of the American call is sometimes optimal, and 
this is duly taken into consideration by every approximation scheme. In the sce-
nario B, as it replicates a situation where the early exercise is never optimal and 
(1)  holds,  all  the  examined  schemes  have  properly  applied  the  Black-Scholes   
Table 1. Parameters employed in the three scenarios simulation. 
Parameters 
S 
K 
r 
q 
b 
σ 
T 
A 
100 
100 
10% 
10% 
0 
25% 
0.25 
Scenarios 
B 
100 
100 
10% 
0 
10% 
25% 
0.25 
C 
100 
100 
−0.301% 
0 
−0.301% 
25% 
0.25 
 
882 
 
A. Cafferata et al. 
 
Table 2. Simulation results for the three scenarios under different estimation models. 
Scheme 
BAW 
BS1993 
BS2002 
T-TREE 
A 
4.8908 
4.8765 
4.8802 
4.8801 
Scenarios 
B 
6.2545 
6.2545 
6.2545 
6.2544 
C 
4.9479 
4.9479 
4.9479 
5.0461 
formula  for  the  European  Call.  In  the  third  case,  the  methods  relying  on 
quasi-closed approximation formulas (BAW, BS1993 and BS2002) have still ex-
ploited (1) which is no more verified, so that they all incorrectly estimated the 
American Call value. On the other hand, the T-TREE scheme generated a more 
robust estimation, because the convenience for the early exercise was checked on 
each node of the tree. As preliminary conclusion, we can therefore state that us-
ing the trinomial trees rather than other approximation schemes might be pref-
erable, as this methodology seems being more robust to anomalous parameters 
values. 
We  then  moved  one  step  further,  giving  additional  instruments  to  evaluate 
such robustness. To such aim, we focused on the scenario C (i.e. the one where 
critical issues arose) and we studied the behaviour of the estimation errors using 
the T-TREE (ErrorT-TREE) and the BS2002 (Error BS2002) schemes, varying once per 
time S, K, T, σ and r. The choice of BS2002 is motivated as it is generally ac-
knowledged to be the more accurate among the examined quasi-closed formulas. 
Figure 1 shows the behaviour of the variable Error = ErrorT-TREE − Error BS2002. 
Looking  at  the  results,  from  Figure  1(a)  we  observe  that,  varying  the  spot 
value S, Error lies within the interval [0.02,0.35], and tends to increase, originally 
in a more than proportional fashion. This suggests the existence of model risk, 
raising as the option’s moneyness increases. Similar considerations apply also to 
the behaviour of Error with respect to the strike price K, shown in Figure 1(b). 
In this second case, in fact, the lower K (high moneyness) the higher Error is, i.e. 
the higher the gap between the T-TREE and the best quasi-closed approximation 
method. In the case of the time to maturity T, observable in Figure 1(c), the di-
vergence between T-TREE and BS2002 is very evident, with Error varying in the 
range [0.05, 0.5]: the longer the hedging period, the worst the performance of 
conventional methods is. For what is concerning the behaviour of Error varying 
r, Figure 1(d) examines only the case of negative risk-free rates. In this case, we 
can observe an elbow-like curve, with higher Error values (more than 0.1) con-
centrated around lowest (and quite unrealistic) negative nominal rates. In every 
case, as r < 0, Error never falls under the 0.09 threshold. Finally, from Figure 
1(e) we can state that there is a positive correlation between the behaviour of 
Error and σ, with the former monotonically growing as the annualized volatility 
increases.   
We  then  examined  the  impact  of  different  approximation  schemes  on  the   
883 
A. Cafferata et al. 
 
 
Figure 1. From top to bottom and from left to right: behaviour of error varying S (a), K (b), T (c), r (d), and σ (e). 
 
value of the most used Greeks [11], because of the paramount role that they play 
in the hedging activity. We therefore evaluated Delta (Δ), Vega (ν), and Theta 
( Θ ), being: 
884 
A. Cafferata et al. 
 
C
∆ = ∂
A
S
∂
;
ν
C
= ∂
A
;
σ
∂
C
Θ = ∂
A
∂ , 
T
AC   is the option value, and  S ,  σ  and  T   are as usual. Table 3 con-
where 
tains the estimated values. 
From the results in Table 3, we look at replicating the situations already dis-
cussed in the first sensitivity analysis, with all the methods generating the same 
values for each Greek in the Scenarios A and B, and with the T-TREE scheme 
providing different results in the case C.   
4. Conclusion 
In  this  paper,  we  examined  how  the  existing  numerical  schemes  react  in  the 
pricing of an American Call option, in presence of anomalous conditions. We 
focused on the case of negative risk-free rate and zero dividends stock as under-
lying  and  we  put  at  work  three  quasi-closed  approximation  formulas  and  the 
Trinomial  Trees  technique.  We  then  analysed  three  toy  scenarios,  replicating 
different market conditions, to conclude that in the case of negative risk-free rate 
it should be preferable pricing the American Calls by way of the Trinomial tree 
(T-TREE) scheme. This is because unlike the other techniques, T-TREE does not 
price the American Call using the equivalence between its fair value at early ex-
ercise and the corresponding value of the European Call with the same financial 
features,  but  rather  the  convenience  for  the  early  exercise  is  checked  on  each 
node  of  the  tree.  In  this  way,  the  T-TREE  is  protected  from  the  risk  that  this 
property  is  no  longer  valid,  as  it  happens  in  case  of  negative  nominal  rates. 
Moreover, in such anomalous conditions, the accuracy of the T-TREE with re-
spect  to  the  other  methods  is  very  robust  to  both  hard  negative  values  of  the 
risk-free rate, and to increases with respect to the moneyness of the underlying,   
Table 3. The impact of different approximation schemes on the value of Delta, Vega and 
Theta Greeks. 
 
Greeks 
ΔBAW 
ΔBS1993 
ΔBS2002 
ΔT-TREE 
ν
BAW 
BS1993 
BS2002 
T-TREE 
Θ
ν
ν
ν
BAW 
BS1993 
BS2002 
T-TREE 
Θ
Θ
Θ
A 
0.5156 
0.5151 
0.5154 
0.5152 
19.5379 
19.4847 
19.4951 
19.5035 
−9.4356 
−9.3404 
−9.3612 
−9.3662 
Scenarios 
B 
0.6035 
0.6035 
0.6035 
0.6033 
19.2716 
19.2716 
19.2716 
19.2714 
−15.0457 
−15.0457 
−15.0457 
−15.0424 
C 
0.5225 
0.5225 
0.5225 
0.5294 
19.9153 
19.9153 
19.9153 
20.3289 
−9.8153 
−9.8153 
−9.8153 
−10.0203 
885