Quantum effects in an expanded Black–Scholes model

Abstract The limitations of the classical Black–Scholes model are examined by comparing calculated and actual historical prices of European call options on stocks from several sectors of the S &P 500. Persistent differences between the two prices point to an expanded model proposed by Segal and Segal (PNAS 95:4072–4075, 1988) in which information not simultaneously observable or actionable with public information can be represented by an additional pseudo-Wiener process. A real linear combination of the original and added processes leads to a commutation relation analogous to that between a boson field and its canonical momentum in quantum field theory. The resulting pricing formula for a European call option replaces the classical volatility with the norm of a complex quantity, whose imaginary part is shown to compensate for the disparity between prices obtained from the classical Black–Scholes model and actual prices of the test call options. This provides market evidence for the influence of a non-classical process on the price of a security based on non-commuting operators. Graphic Abstract Supplementary Information The online version contains supplementary material available at 10.1140/epjb/s10051-022-00402-0.


Introduction
The Black-Scholes-Merton model [1,2] (usually referred to as the Black-Scholes model) was proposed as a rational basis for determining the fair price of options. The key idea behind this model is to hedge an asset, such as a stock, by buying and selling the asset and a call option on that asset in a way that eliminates risk. For a European call option, the Black-Scholes model yields a pricing formula based on five variables: the strike price, the current stock price, the time to expiration, the risk-free interest rate, and the volatility, with the risk-free rate and volatility assumed to be constant [3][4][5]. Because stock prices fluctuate, the value of an option fluctuates accordingly. The Black-Scholes model assumes that stock prices follow a stationary lognormal distribution or, equivalently, that the fractional stock prices follow a normal distribution. Despite the known inadequacy of the normal distribution for market movements, the resulting closed-form formulae for the price of an option at any time prior to expiry provide post hoc justification for this assumption. With many refinements and adjustments, the Black-Scholes model and its variations have become the de-facto standard for The mathematical form of the Black-Scholes equation has led to parallels being drawn with quantum mechanics. As a linear second-order partial differential equation, the Black-Scholes equation can be interpreted as an imaginary-time Schrödinger equation, whose solution can then be formulated as a path integral [6][7][8][9][10][11]. This provides a natural setting for relaxing some of the assumptions of the original Black-Scholes model, such as allowing for stochastic volatility [12,13] and stochastic interest rates [13], and extending its applicability to other types of financial instruments, such as barrier and Asian options [14,15]. Evaluating such path integrals benefits from the vast array of established methods for carrying out exact, approximate, and numerical evaluations [6][7][8][9][10][11]16]. Nevertheless, these studies represent an interpretation of the Black-Scholes equation based on a particular formulation of solutions to Schrödinger's equation in imaginary time, rather than providing a link to the fundamental principles of quantum mechanics.
Here, our focus is on the work of Segal and Segal [17], who introduced quantum effects into the Black-Scholes model to incorporate market features such as the impossibility of simultaneously measuring prices and their instantaneous forward time derivatives. They argued that such effects provide a natural explanation for the extreme irregularities in the short-term movements of market prices.
The mathematical framework for the Segal-Segal model is built on the formalization of non-commuting operators, which, through an accompanying uncertainty principle, is a pillar of quantum mechanics. Segal and Segal modified the Black-Scholes model by adding a new stochastic process to the Black-Scholes stochastic differential equation to account for the influence of factors not simultaneously measurable with those described by the original Wiener process. They determined a calculus for such processes from which a modification of the Black-Scholes pricing formula was obtained.
In this paper, we examine the efficacy of the nonclassical, or 'quantum' Black-Scholes model proposed by Segal and Segal [17] for the valuation of options by considering options prices based on several types of underlying stock in various sectors of the S&P 500. We will first compute the price based on the original Black-Scholes model with estimates of the interest rate and volatility from market data. The quantum Black-Scholes pricing formula [17] is then used to find the magnitude of the additional volatility to associate discrepancies in the original calculations with the presence of the non-classical effects predicted by Segal and Segal. To our knowledge, this is the first time that the Segal-Segal model has been tested against actual market data.
The outline of our paper is as follows. In Sect. 3, we briefly review the salient points of the original Black-Scholes model [1,2] for European options, including the modified model with time-dependent interest rates and volatility. The time-dependent Black-Scholes pricing formulae are used to calculate the prices of call options based on several stocks by using publicly available market data for the strike price, the current stock price, and the time to expiration, and to estimate the riskfree interest rate and the volatility. The basic principles behind the quantum Black-Scholes model are reviewed in Sect. 4. The quantum pricing formulae obtained from this model are used to calculate the imaginary contributions to the generalized (implied) volatility to identify the origins of discrepancies between actual and calculated option prices in Sect. 3. A summary and our conclusions are provided in Sect. 5. Additional comparisons between the Black-Scholes and expanded Black-Scholes models with market data for European call options are provided in the Supplementary Material.

Types of option
Options come in two main forms: calls, which give the owner of the option the right (but not the obligation) to buy the underlying security at a fixed strike price before or on a specified date (the expiry date), depending on the type of option, and puts, which give the owner the right (but not the obligation) to sell the underlying security at a fixed strike price before or on a specified date. These are called long positions: the investor owns the right to buy or sell to the writer of the option at the strike price. Conversely, writing a call or put option, whereby the writer must sell to or buy from the long position holder or buyer of the option, is known as holding a short position.
The most common option styles are European, which can be exercised only at expiry, and American, which can be exercised at any time prior to expiry (and so are worth more). For example, if an investor buys shares in Apple, Inc. (AAPL) at $200, but wants to protect the investment from the downside, that investor can buy puts or sell calls at a strike price of, say, $200. If a short call position is taken, the investor pockets the premium if the share price is below $200 at expiry of the option. Likewise, taking a long position on a put, the investor can buy AAPL shares at the market price at option expiry and sell using the right obtained from the put ownership, thereby making a profit. Although the maximum amount at risk for a long position is the premium paid, the maximum risk is the strike price for a short put and unlimited for a short call.

Background
The Black-Scholes model provided a rational basis for the pricing of options, originally for European options, but later extended to other types of option. This helped to spearhead the expansion of options markets, such as the Chicago Board Options Exchange (CBOE), to the point where these are the largest and most active securities markets. The notional value of outstanding derivatives contracts stood at $600 trillion at the end of June 2021 [27]. The Black-Scholes model is widely used by traders, often with modifications, in the knowledge that, despite well-known deficiencies, the analytic pricing formulas provide considerable insight into the actual price of an option.
The Black-Scholes model was the culmination of work by many authors over many years [28][29][30]. The widespread influence of the Black-Scholes model has spawned formulations [31] that relax several assumptions of the original model and leads to more general expressions for the minimum risk, which is the basis of the Black-Scholes differential equation. More modern approaches [29,32,33] focus on the collective outcome of individual traders with different behaviors to ascertain the origins of price movements and risk. The calculations reported here examine an extension of the Black-Scholes model in an altogether different direction using an analogy with quantum field theory to incorporate quantities that cannot be simultaneously known.

Black-Scholes pricing formulae
This section provides a brief summary of the main steps for obtaining the Black-Scholes formulae for a European call option. More detailed derivations may be found in Refs. [3][4][5]. There are two types of assets in the Black-Scholes model: a risky (i.e., fluctuating) asset S(t), such as a company stock, which is assumed to follow geometric Brownian motion with drift, and a risk-free asset B(t), such as a bank account or a Treasury Bill (T-Bill) paying an interest rate r. The stock price at time t ∈ [0, T m ], where T m is the fixed duration between the issuance of the option and its maturity, is where μ is the percentage drift, σ the percentage variance, both assumed constant, and W is a standard Wiener process, also known as standard Brownian motion. The differential equation for B(t) is whose solution is B(t) = B 0 e rt , where B 0 is the initial balance of a bank account or the initial value of a Treasury Bill. We consider an option on a stock with strike price K and time to maturity T = T m − t. As the stock fluctuates, so does the value V (S, t) of the option. The stochastic differential equation for V (S, t) is, from Itô's lemma, Consider a portfolio that contains the option, which has been sold (short position), and Δ shares of the underlying asset (long position). The value Π of this portfolio is According to Itô's lemma, the stochastic differential equation for Π is Black and Scholes [1] observed that the coefficient of dS, which corresponds to the random part of the Brownian motion and leads to volatility and, therefore, to risk, can be eliminated by setting Δ = ∂V /∂S. The resulting differential equation (5) is deterministic. Once we obtain a solution for V , this relation provides a hedging prescription for maintaining an instantaneously riskless portfolio. We now invoke the assumption that there is no arbitrage in the market, so any riskless portfolio must earn the risk-free interest rate r of the market: dΠ = r Π dt. Using this equation with (4) and (5) leads to the celebrated Black-Scholes equation Explicit formulae for the values of options are obtained by transforming this linear parabolic equation into the backward heat equation (also linear and parabolic). The solution of the heat equation is then obtained in terms of its fundamental solution by applying the appropriate boundary conditions. The details of the passage from (6) to the solution of the heat equation may be found in Refs. [4,5]. For European call options C(S, T ), after transforming back to the original variables, we have that for T m ≥ T ≥ 0. Therefore, the solution to the Black-Scholes equation is in which N (x) is the cumulative normal distribution function, and

Stock option pricing from market data
The pricing formulae (10) and (11) for a European call option are specified by five variables: the strike price K, the current stock price S, the time T to expiration, the risk-free interest rate r, and the volatility σ. The strike price and the time to expiration are set by the writer of the option. The current stock price is available from several sources to whatever frequency desired. We have used daily stock price returns. The specification of interest rates and the volatility merit further discussion. For the interest rate, we take the continuously compounded yield on a 3-month T-bill whose maturity date is closest to the expiry date of the option. T-bills are guaranteed by the government of the United States and are, therefore, considered to be free of default risk. Because the interest rates were so low, and the changes were essentially negligible over the 6-week period of our options window, we took the constant interest rate of 0.08% [34] for all our calculations. An altogether different approach is used for the instantaneous volatility of a stock. We consider companies listed on the S&P 500, which is an index determined by the 500 largest companies in the United States listed on the New York Stock Exchange or NASDAQ ranked and weighted by total market capitalization (the product of the share price and the number of shares of a corporation held by stock holders). We determine the percentage variance for each stock from Here, VIX is the ticker symbol for the volatility index of the Chicago Board Options Exchange, which is a real-time measure of volatility based on S&P 500 index options with near-term expiration dates [35,36]. VIX is a forward-looking index obtained from the implied volatilities of S&P index options and, therefore, represents the market's expectation of the 30-day future volatility of the S&P index, while incorporating the information content of historical (i.e., backwardlooking) volatility. The factor β is a backward-looking measure of the volatility of a stock compared to the volatility of all other stocks in a particular index, in our case, the S&P 500 index [37]. The regimes of interest are: 0 < β < 1 for stocks with a lower volatility than the S&P 500, β = 1 for stocks with the same volatility, and 1 < β for stocks with a higher volatility than this index. The value β = 0 means that the stock is uncorrelated with the S&P 500, and β < 0 that the stock is negatively correlated. For the calculations reported here, we used the following ( In the Black-Scholes pricing formulae (10) and (11), the interest rate and volatility are constant. Allowing the interest rates and the volatility to be timedependent (but not stochastic) leads to modified pricing formulae [4] where the new functions d 1 and d 2 are given by As we are taking the interest rate as constant, the factors in these pricing formulas that include the interest rate reduce to those in (11). Prices of call options were calculated from time series for the volatility. For each time T , this necessitated replacing the integrals in (15) with the corresponding Riemann sums. Historical data [38] were used for the closing price of the call option on each trading day over the 6-week period from October 8, 2020 to November 20, 2020. This was the time period for which the most complete call option data were available. VIX was obtained from Yahoo Finance [39] and β from CNBC [34]. Figure 1 shows the comparison between the actual price of a European call option and the prices calculated from the pricing formulae (14) and (15) with the interest rates and percentage volatilities as described above. There are several noteworthy trends in these comparisons. Most apparent is that the Black-Scholes prices reproduce the gross qualitative trends of the actual prices for options based on each of the stocks in Fig. 1. Looking more closely shows that the Black-Scholes prices are typically lower than the actual prices, with the two prices sometimes showing substantial differences. Nevertheless, there are fleeting moments where the Black-Scholes price exceeds the actual price.
The degree to which the Black-Scholes model reproduces the actual prices depends on the stock used for the option. For example, the price of an option based on stock of the American Airlines Group (Fig. 1a) is consistently underestimated by the Black-Scholes model, while options based on J. P. Morgan Chase & Co (Fig. 1c) and Royal Caribbean Cruises (Fig. 1e) are reproduced quite well by the Black-Scholes price over the entire 6-week period. However, the Black-Scholes model shows large discrepancies for options based on stock from Berkshire Hathaway Class B, Nike, and Tesla (Fig. 1b, d, f). Comparisons between actual and calculated option prices based on other stocks may be found in the Supplementary Material.
While it is possible in principle to extend our analysis to intra-day resolutions, this is unlikely to yield significant improvements, due primarily to the behavior of participants in the options market. Options tend to be traded much less often compared to their underlying stocks, so intra-day price histories tend to be incomplete or unequally spaced in time, with the extreme case being tick data. This means that standardized comparisons between various options lose their noteworthiness. The main conclusion is that, while the agreement in the broad trends of the price trajectories (and, in some cases, the accuracy of the prices) provides validation for our method, the discrepancies also point towards the possibility of factors that are not taken into account by the usual Wiener process in (1). Figure 2 shows the variations of stock prices and volatilities for the same options as in Fig. 1. Most apparent from this figure is the sharp drop in stock prices with a concomitant increase in the volatility midway between weeks 2 and 3, followed by a reversal between weeks 4 and 5. The first period corresponds to the time when the S&P500 fell by some 200 points due  (10) and (11), is then T = Tm − t, where Tm is the fixed duration between the issuance and maturity dates to the rise in coronavirus cases, the reintroduction of lockdowns and the lack of fiscal stimulus, compounded by the uncertainty of the presidential election in the United States. The VIX generally trades in opposite directions to stocks, so the opposite movement is in line with the expected behavior (more fear amongst investors manifesting as increased volatility). Options tend to be more forward-looking than stocks, so the market sentiment reflected a more positive outlook over the longer term, given the earnings beats that were reported in the third quarter of 2020. A comparison of Figs. 1 and 2 shows that the options react in substantially different ways to the variations in the underlying stock. Options based on stocks in the financials sector (BRK-B and JPM) showed marked increases in prices that coincide with the corresponding 13% and 17% increases, respectively, in their stock prices. In the consumer discretionary sector (NKE, RCL, and TSLA), only RCL showed evidence of the steep increase; the stock increase was near 20%. The TSLA and NKE options showed no unusual behavior, which is likely due to increases of less than 10%.

Quantum Black-Scholes model 4.1 Volatility in the classical Black-Scholes model
The comparisons in Fig. 1 used the 'classical' Black-Scholes model [1,2] of a European call option based on publicly available information. These comparisons reveal several types of trend: (i) The two prices track each other, with the simulated price below the actual price and a time-varying difference that may approach 50%, but with no crossing. Figure 1a, d shows this type of behavior. (ii) Same as (i), but with crossing. Figure 1b, e, f shows this type of behavior. (iii) The two prices track each other, with a small timevarying difference that does not exceed a few percent, and with crossing. Figure 1c shows this type of behavior. (iv) The simulated price exceeds the actual price for most or all times during the lifetime of the call. Examples are provided in the Supplementary Material. Figure 2 shows that there are large swings in the volatility with proportionately smaller swings in the prices of the underlying stocks. Of the five quantities that determine the formulae (14) and (15), the volatility of the stock is the most difficult to estimate, because, unlike the other four parameters, volatility cannot be observed directly and, therefore, has been open to interpretation. According to the pricing formulae (14) and (15), the price of a European call option increases with the volatility. Thus, an accurate estimate of the price of an option depends on an accurate volatility. The most common estimates of the volatility are based on the his-torical volatility of the stock (backward-looking) and the implied volatility of the option (forward-looking) [3][4][5].
Historical volatility depends on past performance of the stock and, in this regard, is not necessarily reflective of future volatility. This gives it limited use in the calculation of option prices, which by their very nature are inherently forward-looking. Implied volatility involves using the prices of near-the-money options, which are largely set by market forces of supply and demand, to calculate the resulting volatility using the Black-Scholes formulae (10) and (11). This implies that one assumes the Black-Scholes price to fully reflect the market price, even taking into account the assumptions that are known to make it an idealized model, rather than a true reflection of market dynamics.
We propose VIX as a better estimate for the following reasons. Firstly, its calculation does not depend on the Black-Scholes formula, but instead is modelindependent, and is found using [40,41] F = Strike + e RT (Call Price − Put Price) , (16) from which we obtain and Here, T is the time to expiration, F is the forward index level derived from index option prices (i.e., determined by identifying the strike price at which the absolute difference between the call and put prices is the smallest), K 0 is the first strike below the forward index level F , K i is the strike price of the i th out of the money option, ΔK i = 1 2 (K i+1 − K i−1 ) is the interval between the strike prices, R is the risk-free interest rate, and Q(K i ) is the midpoint of the bid-ask spread for each option with strike K i .
Thus, the calculation of VIX is unambiguously formulated and (more significantly) not restricted by the usual Black-Scholes assumptions, making it more reflective of market conditions than an implied volatility based on the latter. Second, the formula (18) uses the prices of 30-day forward index options, and so is forward-looking. This makes it a more reliable indicator than historical volatility, since it is a measure of the future expected volatility of the index.
The central conclusion of Sect. 3 is that option prices determined by (14) and (15) are typically below the actual prices, with exceptions resulting from particular events that cause pessimism amongst investors. The particular events and the effect on the stock and option prices are best explained on a case-by-case basis. The reason that such events lead to calculated option prices exceeding actual prices is that our volatility, which is obtained from and reflects public information, may not (e) (f) Fig. 2 The variations of the stock prices (blue lines) and volatilities (red lines) over the same 6-week period for the same call options with the same strike prices as in Fig. 1. The left and right vertical scales are for the stock prices and volatilities, respectively change quickly enough to account for the reaction of investors to such news (over small time scales). Indeed, there are several instances in Fig. 1 where the trends in calculated and actual prices differ significantly. In general, a falling stock price implies higher volatility (due to added fear) (Fig. 2), but a greater volatility in (14) and (15), with all other quantities held constant, yields a higher call price. Hence, we must look to the interplay between the changes in stock price and volatility, and which has a larger effect on the option price. In the resolution of this, the approach of Segal and Segal [17] provides a viable alternative to the classical Black-Scholes model. One may ask whether it is valid to simply attribute a discrepancy between market volatility and that obtained from any model to an inaccurate estimation (i.e., the inevitable experimental error), and take the market price as reflecting the true volatility σ of a given stock. However, this does not take into account the fact that the information content of volatility is derived from multiple sources, all of which may not be simultaneously actionable independent of the observer. In particular, we note that non-public information here holds a different meaning to the usual connotation of additional information that might be known to, and thus simultaneously usable by, a company insider (which merely increases the existing information content of the classical volatility σ). Hence, while using VIX does correspond to an accurate estimate of the usual parameter σ, simply taking a different value does not adequately reflect a stock's true volatility. In essence, we bring into focus the limit to which the meaning of the usual volatility is applicable. As we show in the next section, considering public and non-public information (for example, the prices of stocks and their instantaneous forward time derivatives, respectively) as nonsimultaneous observables in the physical sense implied by quantum field theory (i.e., a nontrivial equal time commutation relation) capture this aspect well, yielding a generalized volatility parameter.
The time period chosen certainly happened to be unusual, given the impact of external factors such as those mentioned in Sect. 3.3; this is reflected in the larger volatility σ, as depicted in Fig. 2. However, one would expect this to only reduce the discrepancy between the classically simulated and actual option prices. Hence, any remnant discrepancy can truly be attributed to a (simultaneously) non-actionable source. Therefore, while we expect the Segal-Segal model to be applicable in all situations, it is rendered best in two cases-during such unusual periods, and during periods when the original volatility is within 'normal' ranges but when the option price discrepancy is large.

Quantum Black-Scholes model and pricing formulae
Segal and Segal [17] introduced non-classical effects into the Black-Scholes model as a way of incorporating market features, such as the impossibility of simultaneous measurement of prices and their instantaneous derivatives. The basic idea is to add to the Wiener process W in (1) for the evolution of public information, a process X that represents processes that are not simultaneously observable with those in W . By invoking the formal structure of bosonic quantum field theory [42], Segal and Segal showed that the linear combination aW + bX can be represented as Φ((a + ib)c t ) = Φ(f (t)), where Φ is a mapping from vectors in a complex Hilbert space H to Hermitian operators in the quantum field Hilbert space K and c t is the characteristic function of the interval [0, t]. The extension from real to complex functions means that the generalized process is referred to as a pseudo-Wiener process. For any element f of a Hilbert space H, e iΦ(f ) is the corresponding Weyl operator, whose definition is restricted to the interval [0, t] because of the characteristic function. Hence, the basic equation (1) becomes To understand the significance of the mapping Φ, we invoke the duality transform between the particle and wave representations of a quantum particle (Theorem 3 in Ref. [43]). The fundamental connection to the principles of quantum mechanics is readily established, where the operators representing public and private information are generalizations [44] of the mappings in Refs. [45,46]. Applying the duality transform to e iΦ(f ) (Theorem 4 in Ref. [43]), we obtain, for arbitrary f, g in a real Hilbert space, the usual canonical commutation (Weyl) relations (Corollary 4.1 in Ref. [43]), which take the form e iΦ(f) e iΦ(g) = e iΦ(f+g) e 1 2 iIm( f,g ) (20) in a complex Hilbert space H [47].
The operators Φ(f ) mutually commute if f ∈ R in H, and likewise for Φ(if ). However for f, g ∈ R, where the right-hand side is an inner product defined on H. This is shown to be formally equivalent to (20) by taking the closure of the conventional boson field creation and annihilation operators and using the appropriate infinitesimal form of (20), as in Ref. [47]. The non-vanishing commutator (21) makes apparent the underlying reason why the combination of the public and private information processes cannot be modelled as two-dimensional Brownian motion.
Returning to (19), the real part of f (t) can be regarded as the process representing public information, whereas the imaginary part corresponds to a pro-cess that cannot be observed simultaneously with public information (for example, the price of a stock and its instantaneous forward time derivative, respectively). The Feynman-Kac formula can then be extended in the non-commutative form (see Refs. [17,48,49] for a detailed derivation) and the normal distribution of Φ(x) used, with so that the stock price follows a lognormal distribution over [0, T ] with The price of a call option is then expressed as where Here, s = |f (T )|, γ ≥ 1, and k = s 2 or k = 0 for Wiener or serial correlated processes, respectively. Comparing the denominators for this pricing formula with the corresponding denominators in (10) and (11) for the classical Black-Scholes model suggests that γs generalizes the classical volatility σ. Indeed, by setting γ = 1, s = σ, and k = σ 2 in (25) and (26), we recover (10) and (11). Finally, (25) and (26) can be generalized to accommodate time-dependent volatilities and interest rates in analogy to (14) and (15), to yield withd 1 andd 2 given bỹ The results (25) and (26) obtained by Segal and Segal [17] are based on mathematical development of noncommuting simultaneous processes, each governed by a Wiener process in the stochastic differential equation for the stock price (19). We present here a test of these ideas against market data by determining the values of Im(f (T )), representing the additional part of the volatility for which the calculated and market prices in Fig. 1 coincide. We have set γ = 1, since an increasing γ was found to increase the volatility of the option price towards either extremum in the time interval, i.e., towards the initial time or the time of expiry.
Results are shown in Fig. 3 for the stock options in Fig. 1. As is immediately evident from a comparison of the option price trajectories and the corresponding bars representing the magnitudes of the additional volatility Im(f (T )) at times well before expiry, the value of Im(f (T )) increases with the difference between the actual and calculated option prices, and its magnitude is of order Re(f (T )) = σ. It is also interesting to note that Im(f (T )) is particularly sensitive to discrepancies at times closer to expiry, i.e., at times when the expectation of the classical simulation converging to the actual price is greater (for example, in Fig. 3d, the price discrepancy just before expiry is very small, while the value of Im(f (T )) is in fact at its highest across the 6 weeks). In other words, the divergence of option prices in the securities market obtained from the Black-Scholes model is compensated by the inclusion of a pseudo-Wiener process incorporating information that is not simultaneously observable with public information. This additional information may take the form of speculation or the reaction of large market participants to one or more events. The evidence we have presented points directly to the presence of inefficiencies in market pricing, which lends further credibility to the proposition that the efficient market hypothesis does not offer a complete description of the stock market at the level of individual securities or over short time scales, and that such inefficiencies have the potential to be exploited to make greater returns.
There is a noteworthy issue that results from an analysis of the under-performance of the option. Such instances yield entirely imaginary solutions for the value of Im(f (T )) and, hence, cannot be plotted on the same axis with instances that yield real solutions. One can interpret them as reflecting the situations where the pressure of the falling stock price outweighs that of the increasing classical volatility σ. It was found that the moduli of these imaginary solutions were, in general, much smaller than those of the real solutions. From a physical point of view, however, an imaginary value of b in the coefficient of the generalized quantum Wiener process, Φ((a+ib)c t ) = Φ(f (t)), would be indistinguishable from a purely real contribution a of the appropriate smaller magnitude, with b = 0. Since the deviation from the classical pricing model is greater when the quantum process has a non-vanishing contribution, this leads us to the conclusion that, while the probabilities of over-performance and under-performance may be equal at any given instant, the amount by which an option over-performs is greater than that by which it under-performs. If an option over-performs at a given instant in time, the discrepancy with respect to the clas- (e) (f) Fig. 3 Comparison between actual prices of European call options (shown in red) and prices calculated from the pricing formulae shown in Fig. 1 (blue). The vertical bars represent the values of Im(f (T )). Gaps in these bars indicate that Im(f (T )) is itself imaginary, i.e., that the optimized value of f (T ) is real sical model is greater than if it had under-performed at that instant, since s = |f (T )| is naturally greater in the former case than in the latter.

Summary and conclusions
We began by calculating option prices based on stocks from several sectors of the S&P 500 from the classical Black-Scholes formulae with time-dependent coefficients. By comparing these calculations with actual prices, we observed several trends that suggested the possibility of additional factors that determine the price of the option, but which cannot be observed (or actioned) concurrently with the publicly available information used in the classical theory. Following the suggestion of Segal and Segal, the Black-Scholes model was expanded using a free boson field with infinite degrees of freedom to represent this additional information. This statement was shown to take its usual meaning in the quantum sense, in the form of a non-zero commutation relation between the mappings denoting the two Wiener processes.
The motivation behind using the underlying canonical commutation relations (in the Weyl form) was outlined by showing that the essential general mathematical structure for both systems is the same. The price of a call in the generalized development was determined, which yielded a generalized complex-valued volatility. These formulae were used to price the aforementioned options, with the imaginary part of the volatility shown to compensate for the difference between the market and classical Black-Scholes prices. This provides direct proof of a non-commuting correction to the classical model, providing market evidence that, over short enough time intervals, a securities market reflects essential quantum features, just like any other physical system.

Author contributions
AB designed the study, carried out the computations, analyzed the results, and wrote the first draft. DDV supervised the project, contributed to the mathematical formulation and analysis of the results, and edited the first draft.
Data availability statement All of the market data were obtained from publicly available sources, which are cited at appropriate places in the paper. All of the computed data are presented in the figures in the paper and Supplementary Information. This manuscript has data included as electronic supplementary material. The online version of this article contains supplementary material, which is available to authorized users.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/.