Expected 1DTE Option Returns
Expected 1DTE Option Returns
Abstract
∗
Johannes and Shah are at the Graduate School of Business, Columbia University ([email protected]
and [email protected]). Kaeck is with the University of Sussex Business School ([email protected]),
and Seeger is with VU Amsterdam ([email protected]).
where BtP is a Brownian motion under P, and µ and σ are the drift and volatility coefficients.
Consider a call option on St at time t, struck at K, and maturing at time T with the
price denoted as C = Ct,T (St , K). By combining the option price evolution under P using
Ito’s lemma and the pricing PDE, which captures Q-measure dynamics, the instantaneous
change in the option’s price is
dB P .
∂C ∂C
dC = rC + (µ − r) St dt + σSt
∂St ∂St t
EPt
dC 1 ∂C
− rdt = (µ − r) dt,
C C ∂St
which implies that expected call/put returns are always positive/negative, with magnitudes
determined by the equity premium and the option’s elasticity. While these representations
are useful for intuition, the conditional linearity in “dBtP ” and normality of instantaneous
option returns doesn’t hold generally, even in continuous-time models with crash-risk and
certainly doesn’t hold over the non-infinitesimal holding periods, even for 1DTE options.
Because of this, the literature on option returns often focuses on hold-to-maturity
(HTM) returns. For short-dated expiries, this is particularly the case; see, e.g., Beckmeyer
et al. (2023), Vilkov (2024) and Almeida et al. (2024). In the next section, we develop new
theoretical tools to analyze conditional moments of HTM option returns, not only for very
short-dated call or put options but also for HTM option returns for any maturity and also
for common option strategies.
where Rt,T = ST /St is the gross returns of the underlying St and [x]+ = max [x, 0].13 The
second equality holds given the homogeneity of option prices and highlights the nonlinear re-
lationship between realized gross asset returns and option returns. HTM put option returns
p
are defined as rt,T = [K − ST ]+ /Pt,T (St , K) − 1 where Pt,T (St , K) is the put option price.
HTM option returns are generally highly non-normal, with the degree of non-normality de-
termined by the underlying asset return characteristics (e.g., volatility, jump risk, etc.) and
contract characteristics like maturity and moneyness.
As mentioned earlier, characterizing option return distributions is difficult for many
reasons, including the idiosyncrasies of option samples and the impact of extreme outliers.
Compared to historical analyses of stock returns, option samples are always significantly
smaller as stock options were only introduced in the 1970s (e.g., no observations during the
Great Depression), and common databases like OptionMetrics started in 1996. This implies
that index option return analyses omit the largest historical outlier, the Crash of 1987.14
Even accounting for events like the Crash of 1987, Broadie et al. (2009) show that average
1-month HTM option returns are extremely noisy, especially for OTM strikes, due to a small
number of outlier returns combined with the convexity of option payouts.
It is important to note that these sampling issues have persisted and even been ex-
acerbated by events in the past few years. For example, the largest observed SPX index
OTM put returns occurred during the Covid outbreak, at least for longer-dated options,
which have dramatically changed average put returns.15 These issues are particularly true
for short-dated options, which not only have shorter samples but also had a staggered intro-
duction of different daily expiries, with Fridays introduced in 2005, Mondays/Wednesdays
in 2016, and Tuesdays/Thursdays in 2022. As mentioned above, the staggered expiry in-
troduction of daily expiries implies that important dates are missing from 1DTE or 0DTE
option samples, creating important biases. The next section describes new model-based
13
See Nawalkha and Zhuo (2022) for finite holding-period option returns.
14
See, e.g., Dew-Becker and Giglio (2023) for an analysis of historical synthetic option prices.
15
During March 2020, the 5% OTM SPX put had a 1-month HTM return of greater than 4500% and 10%
OTM put had a return of more than 10,000%, both multiples higher than any previously observed.
M eu log ST ,
t,T (u, St , Lt ) = Et
ΨM
where M ∈ {Q, P}, Lt are latent state variables (e.g., stochastic volatility), u ∈ C (the set
of complex numbers), and EM
t [·] denotes the time-t conditional expectation under measure
M. The characteristic function also depends on structural parameters (which may differ
under the two probability measures) and whose dependence is suppressed for simplicity.
The key to deriving expected option returns is the fact that the numerator and denom-
inator in the expected option return expression are similar. The denominator, the option
16
See, e.g., Broadie et al. (2009), Branger et al. (2010), Hu and Jacobs (2020), Israelov and Kelly (2017),
or Almeida and Freire (2022).
17
For the special case of Black-Scholes, Boyer and Vorkink (2014) calculate second and third-order mo-
ments using the properties of conditional log-normal random variables.
M M [ST − K]+
t,T (K, St , Lt ) := Et
Λc,
Z∞ " −iu log K M
#
= ΨM
e Ψ (iu, S , L )
t,T (1, St , Lt ) − K
1 + 1 ℜ t,T t t
du ,
2 π iu(1 − iu)
0
√
where ℜ [z] denoted the real part of a complex number z and i = −1. Given these
moments, expected call option returns are given by
P
EPt rt,T
Λc,
t,T (K, St , Lt )
Q
c
= − 1,
e−r(T −t) Λc,
t,T (K, St , Lt )
where r is the risk-free rate (assumed constant). Similar expressions for expected put option
returns are available.
To analyze higher moments, we follow a slightly different approach, which easily gener-
alizes also to portfolios of options. The key step is the fact that higher-order return moments
can be decomposed using the Binomial formula. The following proposition shows how to
calculate the n-th raw moment of put and call option returns using the mild assumption
that the corresponding moment of the stock price process exists. The result is general and
applies to any pricing model with a known characteristic function under P and Q.
Proposition 1. Suppose the underlying S has finite moments up to n-th order, and the
characteristic function of log S is known in closed form under the risk-neutral measure Q
and under the real-world measure P. Then, the n-th raw moment of the hold-to-maturity
call return is given by
EPt
n n
n X X k n
c n
Γct,T (ℓ) Ct,T (St , K)−k (−1)n−k (−K)k−ℓ
rt,T = (−1) + ×
ℓ=0
ℓ k
k=max[ℓ,1]
where
Z∞ −iu log K
Γct,T (ℓ) = ΨP P
1 1 e
t,T (ℓ, St , Lt ) + ℜ Ψt,T (iu + ℓ, St , Lt ) du.
2 π iu
0
EPt
n n
p
n X X k n
n
Γpt,T (ℓ) Pt,T (St , K)−k (−1)n−k K k−ℓ
rt,T = (−1) + ×
ℓ=0
ℓ k
k=max[ℓ,1]
where
There are a few notables at this stage. First, calculating the n-th moment requires
n + 2 numerical integrations. For the first moment, the calculation requires three numerical
integrations (one for the option price and one for ℓ = 0 and ℓ = 1 each). This compares
to only two numerical integrations using the results provided above. This approach mirrors
the early derivations of option prices as in Heston (1993) and others, which also lead to
an additional integral. The advantage of this formula is that when calculating higher-order
moments, as Γct,T (ℓ) is not a function of n, these calculations can be used for all the moments
and hence the computational increase for every additional moment is only one additional
integration for every additional moment. Second, the proof of the proposition can be easily
adjusted to other quantities, such as the moments of the P&L distribution, or when option
payoffs are normalized by other quantities, such as the underlying stock price rather than
the option price. Third, as shown in the following paragraph, our approach can be extended
to option portfolios.
To extend these results to returns on option portfolios, consider a portfolio consisting
of a risk-free zero bond investment (with fixed payoff B ∈ R at maturity T ) and call options
with different strikes, allowing for a zero strike price. This assumption is without loss of
generality, as with this extension, we can form portfolios of the bond, stock index (using a
strike zero call) and call options, which allows us to include put options by put-call parity.
We use U to denote the normalization constant in the option portfolio return; this is usually
the price of the portfolio but could be set to any other positive constant, especially in the
case of zero-cost portfolios in which the return is not defined. More specifically, we derive
where Ki ≥ 0 are the strikes of call options i = 1, . . . , m, ci denotes the holding of the i-th
option and U > 0 is known at time t.
Proposition 2. Suppose the underlying S has finite moments up to n-th order. The n-th
conditional return moment of a maturity T portfolio consisting of a bond investment (payoff
B ∈ R at time T ) and call options with strikes K = [K1 , . . . , Km ] ≥ 0 and option holdings
c = [c1 , . . . , cm ], with normalization constant U > 0, is given by
where
EPt [PCn ] =
m
!
X X X
α(n, k, j) × Γ K(k),
e ji
k∈K(n) j∈J (k) i=1
Z∞ −iu log K
1 P P
1 e
Ψt,T (ℓ, Vt , Lt ) +
ℜ × Ψt,T (iu + ℓ, Vt , Lt ) du if K > 0
Γ(K, ℓ) = 2 π iu
0
ΨP (ℓ, V , L )
t,T t t if K = 0
K(k)
e = max Ki
1≤i≤m, ki >0
m
Y m
Y
α(n, k, j) = γ(n, k)β(k, j) cki i (−Ki )ki −ji .
i=1 i=1
The sets K(n) and J (k) and the functions γ(n, k) and β(k, j) are defined in Appendix A.2.
Proof. The proof is provided in Appendix A.2 and follows from two Lemmas, which are also
proofed in Appendix A.2.
We remark that the main computational complexity for our higher-order moment for-
mula for portfolios is similar to option moments of vanilla options. The numerical integra-
tions required are of the same form, although there are now more of these depending on the
number of different strikes in the option portfolio.
10
To understand these option moments, we start with HTM returns in the Black-Scholes
model and then consider option return moments in a model augmented with event risk.
Expected option returns in the Black-Scholes model can be calculated using variants of the
Black-Scholes formula. For a call option, we have
e−r(T −t) EQ
c
= +
− 1.
t [ST − K]
EPt [ST − K]+ = eµ(T −t) EPt e−µ(T −t) [ST − K]+ ,
with the expectation given by a modification of the Black-Scholes formula, replacing r with
µ:
EPt e−µ(T −t) [ST − K]+ = St N (d1,µ,σ ) − e−µ(T −t) KN (d2,µ,σ ),
where
EPt rt,T
µ(T −t) −µ(T −t)
c
e St N (d 1,µ,,σ ) − e KN (d2,µ,σ )
= − 1,
St N (d1,r,σ ) − e−r(T −t) KN (d2,r,σ ),
which is a convenient formula. Similar formulas are available for puts as well as static
strategies like straddles, strangles, and spreads, all of which are easy to calculate.
To understand the quantitative implications for 1DTEs, Figure 1 shows expected Black-
Scholes option returns for a range of parameters. Expected call (put) option returns are
increasing in µ and generally decreasing in σ, and both can be quantitatively important.
Call (put) options are leveraged long (short) positions in the underlying, which explains
their dependence on expected returns.
Quantitatively, expected put returns for ATM 1-day options are small, approximately
-2% for the case with µ = 10% and σ = 25%, falling to about −5% when σ = 10%. These
11
expected returns are magnified for OTM calls and puts, as well as in low volatility and high
expected parameterizations. While some of the parameters in Figure 1 may seem extreme,
very high values of implied volatility and a range of expected returns are reasonable as both,
in reality, vary significantly with time.
The Black-Scholes model is also useful for thinking about volatility risk premiums,
albeit in an ad-hoc manner. To get a sense of the quantitative implications of a variance risk
premium, consider a Black-Scholes model assuming that the diffusion coefficient is different
under the P and Q measures. Historical comparisons of the VIX index to realized variance
(on a volatility scale) imply that the difference between average Q and P-measure volatility
is approximately 2.5%, which would be generated by, for example, an average historical VIX
around 20% and realized volatility about 17.5%. We provide evidence below for substantial
12
13
These population moments highlight the high sensitivity of option return volatility to
both moneyness and underlying asset volatility, as option return volatility is driven by the
combination of tail probabilities on index returns and the convexity of OTM option payouts.
As mentioned earlier, measuring option return volatilities in observed samples is challenging,
due to a combination of historical samples with modest lengths, the peculiarities of when
daily expiries were introduced, and the infrequency of outlier events. The impact on outliers
is well known and discussed in Broadie et al. (2009). Suppose the largest one-day index re-
turn is 2.5% in a given sample. This implies that all calls further OTM have -100% returns
(no zero observed volatility), clearly a Peso problem. These model-based return volatilities
provide benchmarks for difficult-to-estimate expected option return volatilities, which, in
turn, can be used for investment purposes (understanding risk/return tradeoffs or for con-
structing optimal portfolios) but also for calculating standard errors under the assumption
of a given model such as the Black-Scholes model or Merton’s jump-diffusion. Given the
14
obvious issues with small samples, these analytical moments are also a computationally effi-
cient alternative to Monte Carlo simulations. The advantage of using Black-Scholes is that
the inputs are observed – the realized average index returns and volatilities.
Short-dated options load heavily on event risk, as the relative importance of diffusive risk,
which scales linearly with time to maturity, is small. We now provide the main intuition
for analyzing option returns with event risk using a simple extension of the Black-Scholes
model. The model builds on Dubinsky et al. (2019) and Johannes et al. (2023), but the
predictions that we derive for option returns require modeling assumptions both under P
and Q since option returns depend crucially on P-measure dynamics. The option pricing
models in Dubinsky et al. (2019) and Johannes et al. (2023) are restricted to Q-measure
15
P P
1 2
ST = St × exp µ − σ (T − t) + σ BT − Bt + Zτz ,
2
where BtP is a standard Brownian motion under P and µ and σ are the diffusive drift and
standard deviation, respectively. This evolution implies that
1 Q 2 Q
2
Zτz ∼ N − σz , σz .
2
Q Q
1 2
ST = St × exp r − σ (T − t) + σ BT − Bt + Zτz .
2
16
−r(T −t) Q
! h
var Xt,T i+
BS St , K, , r, T − t = e Et St er(T −t)+Xt,T − 21 var(Xt,T )
−K
T −t
where Xt,T is a mean-zero normal variable with variance var(Xt,T ) under Q. Given this,
Q Q
we define Xt,T = σ BT − Bt + σzQ εz and the shocks with εz ∼ N (0, 1) and var (Xt,T ) =
σ 2 (T − t) + σzQ . This implies that the Black-Scholes formula applies with a modified
2
(σzQ )
2
Figure 4 shows the significant impact of event risk on expected put and call option
returns. Call and put option returns are extremely sensitive to µz , with call (put) returns
increasing (decreasing) in the macroeconomic return premium. The economic magnitudes
are quite striking, and point towards a first-order effect of this premium for explaining option
returns. The impact of macroeconomic announcement volatility and the announcement
volatility premium is also significant, especially for OTM options.
Merton’s model extends Black-Scholes to incorporate crash risk, and in addition, we incor-
porate macroeconomic event risk. We assume that under P stock prices evolve via
P
Nt Ntm
where λP is the (constant) intensity of the Poisson process NtP and τj denote the random
jump times of NtP . The jump sizes are given by ZjP ∼ N µP P
2
y , σy with jump com-
pensator µP = exp µP P 2 − 1. Event risk is modeled by Zj , which represents the
h i
1
y + 2 σy
additional return on the j-th announcement at time τz,j . Note that the counting process
Ntm is deterministic, which is a reasonable assumption for short option maturities we con-
17
Sτ Q (eZj − 1) − λQ µQ St dt + d
!
X Q X
dSt = rSt dt + St σdBt + d Sτz,j− (eZj − 1) ,
j−
j=1 j=1
where λQ is the (constant) intensity of the Poisson process NtQ and τjQ denote the ran-
dom jump times of the risk-neutral Poisson process. The jumps are driven by ZjQ ∼
18
N µQ Q Q µQ σyQ
2 h 2 i
1
y , σy with jump compensator µ = exp y + 2
− 1. The macroeco-
nomic announcement returns under the risk-neutral measure are again assumed normally
distributed: Zj ∼ N − 21 σzQ , σzQ , which imposes the martingale restriction as in the
2 2
µP u+ 1 P2 2
P(u, St) log(St )+(µ− 21 σ 2 )τ u+ 12 u2 σ 2 τ +λP τ e y 2 (σy ) u −1−µP u +(N m −N m )(µz u+ 1 (σ P )2 (u2 −u))
ΨM,
T t 2 z
t,T = e
Q 1 Q2 2
Q(u, St) log(St )+(r− 21 σ 2 )τ u+ 21 u2 σ 2 τ +λQ τ eµy u+ 2 (σy ) u −1−µQ u +(NTm −Ntm )( 12 (σzQ )2 (u2 −u))
ΨM,
t,T = e
where τ = T − t.
To provide some intuition on the model assumptions, Figure 5 plots expected option
returns for put and call options of different moneyness levels and their model-based standard
19
3.1 Data
We obtain intradaily option quote data on the SPX index from the CBOE.19 The data
include the national best bid/offer (NBBO) at the end of each 5-minute interval during
the regular trading session from 9:30 am ET until 4:15 pm ET, as well as corresponding
S&P 500 index levels. We discard most option maturities in our data set and focus on
contracts with 1-day to maturity. To avoid prices at the close (see Goyenko and Zhang,
2019), we record option mid-quotes at 3:55 pm ET one trading day before maturity. This
initial selection includes option contracts recorded on Fridays (with expiration the following
Monday), as well as options expiring after a public holiday. As weekends have distinct
seasonality patterns (see Jones and Shemesh, 2018), we mainly work with a second sample
for which we restrict the option maturities to exactly one calendar day. Furthermore, we
omit options if trading stops early at 1:15 pm ET on the maturity date (which usually
occurs the day before Independence Day and Thanksgiving and on Christmas Eve) and
options that are not weeklies as they expire in the morning rather than at the close.
In order to maximize the sample size, we focus on the period from January 2, 2012
until December 5, 2023. This is an unbalanced sample, as 1DTE options in the early part of
the sample are only available on Thursdays (for Friday expirations). From February 2016,
the addition of Monday and Wednesday options provides a denser maturity grid, and from
2022, 1DTE options exist for every weekday.
Using the raw data, we calculate a number of option market variables that are required
19
The CBOE does not provide data for three trading days in our sample: January 8, 2013, January 18,
2013, and February 13, 2013. None of these missing data affect our analysis. In addition, the data for March
6, 2020 mainly records unrealistic zero prices and is excluded.
20
for our analysis below. First, we use one-month treasury yields from FRED (St. Louis
Fed) as a risk-free rate proxy, and we closely follow the CBOE VIX methodology and
calculate on each day t in the sample the implied forward price Ft,T for 1DTE options
(maturity T ). To do so, we imply the futures price from put-call parity using the put/call
pair whose absolute price difference is the smallest. We then calculate standard moneyness
as M = K/Ft,T and also calculate Black-Scholes implied volatilities. We then linearly
interpolate the implied volatilities for options close to at-the-money to obtain the at-the-
atm
money implied volatility IVt,T for each day. We use this ATM volatility to define log
atm
√
moneyness as m = log(M )/ IVt,T T − t . Log moneyness is more informative in our
option dataset as it adjusts for the volatility of the underlying asset and allows us to more
cleanly compare contracts with different moneyness levels over time.
Our second main dataset consists of macroeconomic announcements and contains the
four main announcement categories: Consumer Price Index (CPI), Gross Domestic Product
(GDP), FOMC announcements (FOMC), and Employment Situation (Monthly Payrolls).
We collect dates and timestamps for these four categories for our sample period from January
2012 to December 2023. Announcement times are important as some announcements are
during the trading day (FOMC) while others usually occur before the market opens at
8:30 am ET. Our data is hand-collected from original sources (e.g., the Bureau of Labor
Statistics). Table 1 summarizes the announcements, including the number of observations
and their most frequent announcement time. There are no announcements on Monday
before the close of the market, which is helpful as weekend option returns behave differently
from the returns during the week.
21
Our first set of empirical results are for call and put options that are held until maturity. For
each option in the sample, we calculate the return of an option using the mid price at 3:55
pm ET one calendar day before maturity (hence we initially omit Monday expirations and
expirations after a public holiday). The payoff at maturity is calculated using the closing
price of the S&P 500 index on the maturity date.
Table 2 provides average returns and higher-order return statistics for puts and calls
for a number of log moneyness targets: -2, -1, 0, 1, and 2. To construct these returns, on
each day in the sample, we assign the option closest to the target that remains within a
tolerance of 0.5. SPX options close to the money usually trade in strike increments of 5,
with larger gaps further out of the money. For robustness, we also provide the same analysis
for the alternative moneyness definition (K/Ft,T ) in the Appendix (see Table A.1). Panel A
summarizes the return distributions for the entire sample period. Average call option returns
are small for in-the-money calls (2.84% for a log moneyness of -2) and monotonically increase
until a moneyness level of 1. The returns turn negative (-52.17%) for a log moneyness level
of 2. We provide standard HAC-corrected t-statistics, which may be severely biased given
the high kurtosis of out-of-the-money options (call returns for log moneyness of -2 have a
kurtosis of 181.78). According to this benchmark, only out-of-the-money calls have returns
that are significantly different from zero. Put option returns also, on average, show a
monotonic increase from -45.43% for out-of-the-money puts (log moneyness -2) to -3.78%
for in-the-money puts (log moneyness 2). Standard t-statistics show that most of the put
returns are significantly negative. Skewness and kurtosis are reasonably similar to those
of the call options, with high excess kurtosis for OTM puts. Panels B and C show that
our results are highly robust across different sample periods, which is important given the
unbalanced nature of our option data set.
In addition to standard t-statistics, Table 2 also contains Black-Scholes-based t-statistics
(labeled t-stat BS). To calculate these, for each return time series, we subtract from the real-
ized average option return the expected option return in the Black-Scholes model and divide
this quantity by the Black-Scholes standard error. This accounts for the sample size and
22
23
and Kurt are the skewness and the kurtosis of the returns. Obs are the number of observations.
Moneyness Call Option Contract Return Statistics Put Option Contract Return Statistics
t-stat Loss t-stat Loss
mean t-stat Skew Kurt mean t-stat Skew Kurt
BS (%) BS (%)
1.000 17.71 (1.79) (1.14) 82 3.60 14.47 -5.38 (-2.03) (-1.18) 18 1.19 2.84
2.000 -52.17 (-3.71) (-2.19) 98 12.89 181.78 -3.78 (-2.37) (-1.49) 2 0.71 2.00
the strike closest to the Moneyness level in column 1 (with a maximum tolerance of 1%) and record the return of the call and put option, respectively.
Call AvgRet (Put AvgRet) is the average call (put) option return in percentage. t-stat provides the HAC-corrected t-statistic for the mean and and
we also provide t-statistics based on Black-Scholes moments (t-stat BS ). Loss % gives the likelihood of complete capital loss (-100%, i.e. the option is
out-of-the-money at maturity). Skew and Kurt are the skewness and the kurtosis of the returns. Obs are the number of observations.
Panel A: Calls
25
-2.000 2.84 (1.90) (0.80) 1109 8.51 (3.25) (2.45) 314 0.59 (0.36) (-0.60) 795
-1.000 3.98 (1.68) (0.55) 1109 12.79 (2.99) (2.21) 314 0.50 (0.20) (-0.74) 795
0.000 8.15 (1.88) (0.86) 1109 22.15 (2.75) (2.10) 314 2.62 (0.57) (-0.31) 795
1.000 17.71 (1.79) (1.14) 1109 59.72 (2.76) (2.88) 314 1.11 (0.11) (-0.47) 795
2.000 -52.17 (-3.71) (-2.19) 1109 -10.05 (-0.24) (-0.37) 314 -68.81 (-5.68) (-2.36) 795
Panel B: Puts
-2.000 -45.43 (-3.39) (-1.45) 1109 -74.01 (-5.01) (-1.37) 314 -34.15 (-1.93) (-0.85) 795
-1.000 -22.97 (-2.70) (-1.84) 1109 -52.06 (-4.52) (-2.68) 314 -11.48 (-1.04) (-0.49) 795
0.000 -9.14 (-1.88) (-1.16) 1109 -28.22 (-3.60) (-2.97) 314 -1.60 (-0.27) (0.50) 795
1.000 -5.38 (-2.03) (-1.18) 1109 -15.08 (-3.38) (-2.79) 314 -1.55 (-0.51) (0.36) 795
2.000 -3.78 (-2.37) (-1.49) 1109 -10.00 (-3.77) (-3.06) 314 -1.32 (-0.74) (0.17) 795
Figure 6: Average SPX hold-to-maturity option returns
This figure provides average SPX option 1-day to maturity holding period returns for our sample period
January 2012 to December 2023. We separately show holding period returns for calls and puts, with and
without macroeconomic announcements occurring during the holding period. Additionally, the figure
presents the holding period returns for call options combined and the holding period returns for all
put combined. Shaded areas represent HAC-adjusted standard error bounds. Moneyness is defined as
log-moneyness and ranges from -2 to 2.
Average HTM Returns (Call Options) Average HTM Returns (Put Options)
75 0
50 20
25
40
0
25 60
50 Without Macro Announcement Without Macro Announcement
All Observations 80 All Observations
75 With Macro Announcement With Macro Announcement
2.0 1.5 1.0 0.5 0.0 0.5 1.0 1.5 2.0 2.0 1.5 1.0 0.5 0.0 0.5 1.0 1.5 2.0
Log Moneyness Log Moneyness
Similarly, put option returns also show very different patterns for macro and non-macro
days. Returns on macroeconomic announcement days are far more negative and again only
deviate from the BS benchmark on announcement days where all but one category shows
clear deviations from the model-based predictions.
26
mP P P
σ2
1 = µ− − λµ̄ + λµy ∆,
2
P
2
h
m 2 = σ + λ σ y + µy P 2 P 2 i
∆,
mP P2 + µPy 3 ∆,
P
3 = λ 3µy σy
mP P P
4 2 P 2 P 4
4 = λ 3 σ y + 6 σ y (µ y ) + µ y ∆
P P P P
2
h 2 2 i2 2
h 2 2 i
+ 3λ σy + µy + 6λσ σ y + µy 4
+ 3σ ∆2 .
where ∆ = 1/252 is the time between observations and µ̄P is as defined above.
Since the system of moment conditions is overidentified, we obtain Θ∗ by minimizing
a quadratic form of moment conditions. The GMM estimator is given by:
b Θ)′ W ψ(x;
ΘGM M = argmin ψ(x; b Θ).
Θ
with W38×38 required to be a positive definite weighing matrix. Appendix A.3 includes
further details on the estimation methodology.
Table 4 provides our parameter estimates. Our estimate for λP = λQ is 0.947, which
translates into slightly less than one jump per year, on average. The estimates in the
literature for this parameter vary, ranging from values lower than ours to substantially higher
values. Our estimate imposes infrequent jump events, and on occurrence, the jump impact
is high. On average, the options price a negative jump of -1.4% with a standard deviation of
8.5%. Again, these values are economically reasonable, especially considering that our jump
27
µ σ λ µQ
y σyQ σzQ µP
y σyP µz σzP
0.089 0.145 0.947 -0.014 0.085 0.005 -0.028 0.026 0.002 0.005
(0.005) (0.002) (0.525) (0.015) (0.067) (0.000) (0.024) (0.019) (0.000) (0.001)
frequency is relatively low. Our estimate for the macro-economic announcement volatility is
roughly 50 bps. This is comparable to findings in Johannes et al. (2023), although marginally
lower than their estimates. One major difference to Johannes et al. (2023) is that their
estimates are time-varying, and this may explain some of the differences we find.20
The P-measure estimates are also economically reasonable. The estimated expected
return is 8.9%, and σ is estimated to be 14.5%. More importantly, as expected and in line
with Broadie et al. (2007), we estimate a substantial jump risk premium for the random
jump component. Average jump sizes under P are -2.8% with a standard deviation of 2.6%,
and therefore, the impact of jumps under P is, on average, more negative but has a much
lower standard deviation than under Q. Of particular interest are the jump parameters of
the macroeconomic jumps. Our estimate for σzP is 50 bps, the same value as under Q, and
implies, contrary to the findings in Londono and Samadi (2023), that there is no significant
wedge between the two measures. But we find a significant drift premium of 20 bps, which
aligns with the findings of Savor and Wilson (2013) that macro-economic announcement
days have higher returns than ordinary trading days.
Figure 7 provides expected option return plots for the estimated model parameters,
using the Merton model with and without macroeconomic announcement events. The shapes
for the option returns are very similar to the empirical results presented earlier in Figure (6).
In particular, our estimated parameters imply a large difference between expected option
returns on announcement vs. non-announcement days. The shape of the expected return
distribution for calls on announcement days has the distinctive hump shape reported earlier,
20
It is also highly likely that jump distributions change over time. We do not explicitly allow for this, as
the exact variation of the jump distribution is not of interest to our purposes. Our estimated parameters
have to be interpreted as average values of a potentially changing jump distribution.
28
and there is a substantial difference to the shape observed on non-announcement days. Our
model achieves this mainly by imposing a return premium for announcement days. Average
put option returns are also similar to the empirical values, with lower option returns on
macroeconomic announcement days.
In this section, we provide further empirical results for a number of other popular option
trading strategies. We start with variance risk premium strategies, as these have been consid-
29
ered in the literature. Londono and Samadi (2023) find a significant variance risk premium
for macroeconomic announcements, whereas our model estimation in the last section does
not imply an economically meaningful difference between macroeconomic volatility under P
and Q.
We define the variance risk premium as the difference between the observed realized
variance and its expectation under the risk-neutral measure Q. Realized variance between
time-t and time-T is defined as
n
X Fti ,T Fti ,T
RVt,T =2× − 1 − log
i=1
Fti−1 ,T Fti−1 ,T
where Ft,T is the time-t forward contract with maturity T and t = t0 < t1 < · · · < tn = T (see
Neuberger, 2012, and Bondarenko, 2014). With this definition of variance, the expectation
of realized variance under the risk-neutral measure is (due to the martingale property of
Ft,T and the telescoping sum)
EQ Q
t [RVt,T ] = −2 × Et [log (ST )] − log (Ft,T ) ,
30
Overall, our results show that in the entire sample the variance risk premium was -26.22%
(t-statistic -17.66, Panel A), a value that is similar in magnitude to the monthly variance
risk premium reported in the literature (see Carr and Wu, 2006, Neuberger, 2012, Bon-
darenko, 2014). The VRP is slightly larger in the sample with macroeconomic announce-
ments (-28.48%, t-statistic -9.28, Panel B) compared to the sample without macroeconomic
announcements (-25.33%, t-statistic -15.09, Panel C). But overall, unlike our results for
the raw put and call returns, these risk premiums are economically similar across the two
samples with and without macroeconomic announcements.22
21
To do so, we first fill in option prices for strikes that are not observed by using values from an interpolated
implied volatility curve, and extend option prices beyond the traded range by fixing the implied volatility
to the last observed value. The option price integral required for the calculation is then approximated using
the trapezoidal rule.
22
As a robustness check, we calculate the Variance Risk Premium using an alternative methodology.
Specifically, the risk-neutral variance is calculated from option prices following Bakshi et al. (2003), and
realized variance is calculated using intraday and overnight returns. We then augment Corsi (2009)’s HAR
model with an announcement indicator to project the variance forward. The results using this approach are
qualitatively and quantitatively very similar to what we report above.
31
n
X
a,b
RVt,T = h(Fti ,T ) − h(Fti−1 ,T ) − h′ (Fti−1 ,T ) × (Fti ,T − Fti−1 ,T ),
i=1
where
2 × − log b − xb + 1 if x > b
h(x) = −2 log x if x ∈ [a, b]
2 × − log a − x + 1 if x < a
a
and h′ (x) denotes the partial derivative of h(x). As before, the expectation of the realized
volatility under the risk-neutral measure is an option portfolio related to the VIX index,
with options included only if their strikes are between a and b. This corridor realized
variance is a natural extension of the realized variance defined above and has the same
theoretical advantages (the estimation of the risk premium is model-free and independent of
the monitoring frequency).23 To be able to compare realized risk premia over time, rather
than fixing a and b, we calculate the bounds in log moneyness levels α and β. For instance,
if α = −2 and β = 2, on each trading day, we calculate the strikes a and b such that they
coincide with a two and minus two log moneyness level. This allows us to cleanly compare
our estimates over time. The corridor variance premium is a generalization of our earlier
Q h i
definition: V RPt,T = RVt,T /Et RVt,T − 1.
a,b a,b a,b
23
Setting a = 0 and b = ∞ yields the definition of realized variance used above.
32
mean
Risk Premium α β t-stat Skew Kurtosis Obs
HPR
33
4 Conclusions
This paper provides theoretical and empirical contributions to understanding short-dated
SPX options, in particular those with one day to maturity. Our theoretical contributions
34
35
Almeida, C., Freire, G., and Hizmeri, R. (2024). 0DTE asset pricing. SSRN working paper.
Bakshi, G., Cao, C., and Chen, Z. (1997). Empirical performance of alternative option
pricing models. Journal of Finance, 52(5):2003–2049.
Bakshi, G., Kapadia, N., and Madan, D. (2003). Stock return characteristics, skew laws,
and the differential pricing of individual equity options. Review of Financial Studies,
16(1):101–143.
Bandi, F. M., Fusari, N., and Renò, R. (2024). 0DTE option pricing. SSRN working paper.
Bates, D. S. (2000). Post-’87 crash fears in the S&P 500 futures option market. Journal of
Econometrics, 94:181–238.
Beckmeyer, H., Branger, N., and Gayda, L. (2023). Retail traders love 0DTE options... but
should they? SSRN working paper.
Bondarenko, O. (2014). Variance trading and market price of variance risk. Journal of
Econometrics, 180(1):81–97.
Branger, N., Hansis, A., and Schlag, C. (2010). Expected option returns and the structure
of jump risk premia. SSRN working paper.
Broadie, M., Chernov, M., and Johannes, M. (2007). Model specification and risk premia:
Evidence from futures options. Journal of Finance, 62(3):1453–1490.
Broadie, M., Chernov, M., and Johannes, M. (2009). Understanding index option returns.
Review of Financial Studies, 22(11):4493–4529.
36
Constantinides, G. M., Jackwerth, J. C., and Savov, A. (2013). The puzzle of index option
returns. Review of Asset Pricing Studies, 3(2):229–257.
Dew-Becker, I. and Giglio, S. (2023). Risk preferences implied by synthetic options. Tech-
nical report, NBER.
Dim, C., Eraker, B., and Vilkov, G. (2024). 0DTEs: Trading, gamma risk and volatility
propagation. SSRN working paper.
Dubinsky, A., Johannes, M., Kaeck, A., and Seeger, N. J. (2019). Option pricing of earnings
announcement risks. Review of Financial Studies, 32(2):646–687.
Goyenko, R. and Zhang, C. (2019). Option returns: Closing prices are not what you pay.
Working paper.
Heston, S. L. (1993). A closed-form solution for options with stochastic volatility with
applications to bond and currency options. Review of Financial Studies, 6(2):327–343.
Hu, G. and Jacobs, K. (2020). Volatility and expected option returns. Journal of Financial
and Quantitative Analysis, 55(3):1025–1060.
Israelov, R. and Kelly, B. T. (2017). Forecasting the distribution of option returns. SSRN
working paper.
Johannes, M. S., Kaeck, A., and Seeger, N. (2023). FOMC announcement event risk. SSRN
working paper.
37
Kelly, B., Pastor, L., and Veronesi, P. (2016). The price of political uncertainty: Theory
and evidence from the option market. Journal of Finance, 71(5):2417–2480.
Muravyev, D. and Ni, X. C. (2020). Why do option returns change sign from day to night?
Journal of Financial Economics, 136(1):219–238.
Pan, J. (2002). The jump-risk premia implicit in options: Evidence from an integrated
time-series study. Journal of Financial Economics, 63(1):3–50.
Rhoads, R. (2014). Trading weekly Options: Pricing characteristics and short-term trading
Strategies. John Wiley & Sons.
Savor, P. and Wilson, M. (2013). How much do investors care about macroeconomic risk?
Evidence from scheduled economic announcements. Journal of Financial and Quantitative
Analysis, 48(02):343–375.
Storn, R. and Price, K. (1997). Differential evolution - A simple and efficient heuristic for
global optimization over continuous spaces. Journal of Global Optimization, 11(4):341–
359.
38
Proof. Assume the processes are defined on a probability space (Ω, F, P) where ω ∈ Ω are
the elementary events. F = σ (∪t≤T Ft ) where Ft = σ (Ss , Ls : s ≤ t) (conditional expecta-
tions with respect to Ft are indicated by a subscript t). For the call, we can write
n
[ST − K]+
c
n
rt,T = −1
Ct,T (St , K)
where Ct,T (St , K) is the call price. Applying the binomial formula twice, this can be written
as
n X k
k n
Ct,T (St , K)−k (−1)n−k (−K)k−ℓ (1ST ≥K )k STℓ .
c
n X
rt,T =
k=0 ℓ=0
ℓ k
Nikodym density
ST (ω)ℓ
ΨP
L(ω) = .
0,T (ℓ, V0 , L0 )
L is equal to one in P-expectation, strictly positive (since stock prices are strictly positive)
and well-defined by assumption, hence a valid density.24 The new measure Sℓ is defined by
Z
S (A) =
ℓ
L(ω)P(dω) ∀A ∈ F.
A
For an integrable random variable Y the conditional expectations under P and Sℓ are related
by Bayes formula for conditional expectations (EP S P
t [Y L] = Et [Y ] × Et [L]). Here choose
ℓ
24
Note that for ℓ = 0 the density is equal to one.
A.1
P = ΨP S P
" #
STℓ
Et 1ST ≥K ST 0,T (ℓ, S0 , L0 )Et [1ST ≥K ]Et
ΨP
ℓ
ℓ
0,T (ℓ, S0 , L0 )
= ΨP S
t,T (ℓ, St , Lt )Et [1ST ≥K ] .
ℓ
To derive the expectation under the new measure, we use the characteristic function
of log(ST ) under Sℓ :
Therefore, putting these together and using the standard inversion formula for characteristic
functions yields
This concludes the proof for the call (after changing the order of summation). The proof
for the put returns follows the same logic.
We define the payoff of a call option portfolio consisting of m call options, where ci denotes
the holding of the i-th call option, as PC = m +
P
i=1 ci [ST − Ki ] . We allow Ki ≥ 0, which
implies that one of the call options may be the underlying asset, in which case Ki = 0.
We define the following set
( m
)
(k1 , . . . , km ) ∈ N0 × · · · × N0 :
X
K(n) := ki = n
i=1
which is the set of all non-negative integer-valued m-tuples for which the elements sum to
n. We suppress the dependence on m in our notation. Also define, for each k ∈ K(n),
( m
)
(j1 , . . . , jm ) ∈ N0 × · · · × N0 :
X
J (k) := ji ≤ n, ∀ 1 ≤ i ≤ m : ji ≤ ki
i=1
A.2
n!
γ(n, k) :=
k ! × · · · × km !
1
k1 km
β(k, j) := × ··· × .
j1 jm
Lemma 1. The n-th power of the payoff of a call option portfolio can be written as
Pm
i=1 ji
1ST ≥K(k)
X X
PCn = α(n, k, j) ST e
k∈K(n) j∈J (k)
where
K(k)
e := max Ki
1≤i≤m, ki >0
m
Y m
Y
α(n, k, j) := γ(n, k)β(k, j) cki i (−Ki )ki −ji .
i=1 i=1
Proof. Writing the i-th call payoff as ci ST 1ST ≥Ki + ci (−Ki )1ST ≥Ki and applying the multi-
nomial theorem to the n-th power of the call option payoff yields:
m
cji i STji 1jSiT ≥Ki cki i −ji (−Ki )ki −ji 1SkiT−j
X X Y
PCn = γ(n, k) β(k, j) i
≥Ki .
k∈K(n) j∈J (k) i=1
Lemma 2. Suppose the underlying S has finite moments up to n-th order. Then the n-th
conditional moment of a call option portfolio payoff is given by
EPt [PCn ] =
m
!
X X X
α(n, k, j) × Γ K(k),
e ji
k∈K(n) j∈J (k) i=1
where
Z∞ −iu log K
1 P P
1 e
Γ(K, ℓ) = Ψt,T (ℓ, St , Lt ) + ℜ × Ψt,T (iu + ℓ, St , Lt ) du
2 π iu
0
A.3
Proof. The results follow from the previous lemma, the linearity of the conditional expecta-
tion and the same line of argument as in Proposition 1, since Γ(K, ℓ) = EP 1
ℓ
t ST ≥K S T .
To extend these results to general option portfolios, we add a constant (bond invest-
ment) to the call portfolio. With this addition, we can form portfolios of the bond, stock
and call options, which provides full flexibility to include put options by put-call parity. We
define this extended payoff by PC,B = B + m +
P
i=1 ci [ST − Ki ] . This case may be of interest
in its own right, especially when the option portfolio is a zero-cost strategy. Finally, we can
derive the n-th conditional moment of the option portfolio return, provided its price is pos-
itive. Allowing for normalizations with quantities other than the option portfolio value, we
o PC,B
write rt,T = U
− 1 where U > 0 is either the option portfolio value or another normalizing
quantity such as the current stock price (Ft measurable). Putting these two steps together
yields Proposition 2 stated in the main text:
Proposition. Suppose the underlying S has finite moments up to n-th order. The n-th
conditional return moment of a maturity T portfolio consisting of a bond investment (payoff
B at time T ) and call options with strikes K = [K1 , . . . , Km ] ≥ 0 and option holdings
c = [c1 , . . . , cm ], with normalization constant U > 0, is given by
where
EPt [PCn ]
m
!
X X X
= α(n, k, j) × Γ K(k),
e ji
k∈K(n) j∈J (k) i=1
Z∞ −iu log K
P P
1 1 e
Ψt,T (ℓ, Vt , Lt ) + ℜ × Ψt,T (iu + ℓ, Vt , Lt ) du if K > 0
Γ(K, ℓ) = 2 π iu
0
ΨP (ℓ, V , L )
t,T t t if K = 0
K(k)
e = max Ki
1≤i≤m, ki >0
m
Y m
Y
α(n, k, j) = γ(n, k)β(k, j) cki i (−Ki )ki −ji .
i=1 i=1
A.4
Proof. This again follows directly from the Binomial theorem, and an application of the
previous Lemma.
To jointly estimate model parameters of the Merton model with macroeconomic announce-
ment event risk under P and Q we use a GMM approach as laid out in Section (3.2.2). To
do so, we minimize the quadratic form using differential evolution as proposed in Storn and
Price (1997). Following Hansen (1982), the optimal weighing matrix for the second step
is obtained by inverting the variance-covariance matrix of ψ(x;
b Θ1 ) where Θ1 refers to the
consistent estimates obtained from the first step of GMM. Since our theoretical moment con-
ditions are essentially point estimates (for announcement and non-announcement days), the
variance-covariance matrix essentially reflects the variation in empirical moments. Further,
the moment conditions on announcement days are independent of the moment conditions
on non-announcement days. The second stage optimal moment vector Θ2 = Θ∗ is used to
calculate standard errors using
r
b′ Ω
diag (G b −1 G)
b −1
h√ √ i ∗)
where Ω
b = E
b T ψ (x; Θ∗ ) T ψ (x; Θ∗ )′ and the gradient G = E ∂ψ(x;Θ
∂Θ ∗ is approxi-
mated numerically using simultaneous perturbation. Finally, since the parameters are over-
identified we test for the validity of the moment specification as proposed by Newey (1985).
b Θ∗ ) = 0.003 and hence fail to reject the null hypothesis that all moment
We obtain T ψ(x;
conditions have zero expectation at Θ∗ .
A.5
Moneyness Call Option Contract Return Statistics Put Option Contract Return Statistics
t-stat Loss t-stat Loss
mean t-stat Skew Kurt mean t-stat Skew Kurt
BS (%) BS (%)
1.000 8.12 (1.87) (0.85) 46 1.27 0.75 -9.14 (-1.88) (-1.16) 55 3.06 16.02
1.010 -3.57 (-0.30) (-0.99) 85 10.82 188.95 -3.97 (-1.67) (-0.65) 15 1.00 1.36
1.020 -68.94 (-7.58) (-2.32) 97 15.83 300.07 -2.89 (-1.96) (-1.01) 3 0.55 1.19
1.030 -87.62 (-11.46) (-0.59) 99 28.89 893.43 -1.93 (-1.87) (-0.98) 1 0.39 1.61
skewness and the kurtosis of the returns. Obs are the number of observations.
Panel A: Calls
0.970 1.34 (1.28) (0.33) 1109 4.43 (2.19) (1.95) 314 0.12 (0.11) (-0.84) 795
0.980 1.94 (1.36) (0.30) 1109 6.14 (2.23) (1.79) 314 0.27 (0.18) (-0.77) 795
0.990 3.27 (1.51) (0.32) 1109 10.34 (2.56) (1.79) 314 0.48 (0.21) (-0.74) 795
A.7
1.000 8.12 (1.87) (0.85) 1109 22.15 (2.75) (2.10) 314 2.57 (0.56) (-0.32) 795
1.010 -3.57 (-0.30) (-0.99) 1109 18.48 (0.92) (0.58) 314 -12.28 (-0.87) (-1.53) 795
1.020 -68.94 (-7.58) (-2.32) 1109 -31.07 (-1.03) (-0.64) 314 -83.90 (-14.39) (-2.33) 795
1.030 -87.62 (-11.46) (-0.59) 1108 -64.50 (-2.12) (-0.24) 314 -96.76 (-45.54) (-0.54) 794
Panel B: Puts
0.970 -79.98 (-8.19) (-0.38) 1109 -88.42 (-9.97) (-0.23) 314 -76.64 (-5.80) (-0.31) 795
0.980 -46.89 (-2.37) (-1.16) 1109 -71.56 (-6.81) (-1.02) 314 -37.15 (-1.36) (-0.73) 795
0.990 -24.55 (-2.26) (-1.84) 1109 -49.10 (-3.88) (-2.30) 314 -14.85 (-1.04) (-0.73) 795
1.000 -9.14 (-1.88) (-1.16) 1109 -28.22 (-3.60) (-2.97) 314 -1.60 (-0.27) (0.50) 795
1.010 -3.97 (-1.67) (-0.65) 1109 -10.60 (-2.39) (-1.89) 314 -1.36 (-0.51) (0.41) 795
1.020 -2.89 (-1.96) (-1.01) 1109 -6.97 (-2.49) (-2.12) 314 -1.28 (-0.80) (0.14) 795
1.030 -1.93 (-1.87) (-0.98) 1108 -4.84 (-2.40) (-2.19) 314 -0.79 (-0.71) (0.22) 794
Figure A.1: Average hold-to-maturity returns: Standard moneyness
This figure provides average SPX option 1-day to maturity holding period returns for our sample period
January 2012 to December 2023. We separately show holding period returns for calls and puts, with and
without macroeconomic announcements occurring during the holding period. We show returns for up to 3%
out of the money, sampled at intervals of 0.5%. Shaded areas represent HAC-adjusted standard error bounds.
Average HTM Returns (Call Options) Average HTM Returns (Put Options)
40 0
20
20
0
20 40
40 60
60
Without Macro Announcement 80 Without Macro Announcement
80 All Observations All Observations
100 With Macro Announcement With Macro Announcement
100
0.97 0.98 0.99 1.00 1.01 1.02 1.03 0.97 0.98 0.99 1.00 1.01 1.02 1.03
Moneyness Moneyness
A.8