From quantum mechanics to finance: Microfoundations for jumps, spikes and high volatility phases in diffusion price processes
Abstract
We present an agent behavior based microscopic model that induces jumps, spikes and high volatility phases in the price process of a traded asset. We transfer dynamics of thermally activated jumps of an unexcited/ excited two state system discussed in the context of quantum mechanics to agent socioeconomic behavior and provide microfoundations. After we link the endogenous agent behavior to price dynamics we establish the circumstances under which the dynamics converge to an Itôdiffusion price processes in the large market limit.
Key Words: behavioral finance, diffusion process, microscopic foundations, agent based model, econophysics, jumps, spikes
MSC2010: 60F05, 60J60, 60K35, 91D30, 91B69, 91B80
10mm
1 Introduction and Methods
1.1 Introduction
In 1900 the french mathematician Louis Bachelier suggested to use brownian motion to model price fluctuations at the paris stock market (see Bachelier [4]) and so laid the foundations of modern financial mathematics. Through the last century, especially since the introduction of the BlackScholes model in 1970 (see Black and Scholes [8]), Itô diffusion processes became the standard tool for modeling and pricing of financial assets. On the other hand people are aware that the single asset price process is a macroscopic result of many microscopic factors like for example individual traders behavior. So it is not surprising that many and various probabilistic models have been invented to study not only the individual agents trading behavior but also their interaction.
In Föllmer and Schweizer [12] as well as in Horst [17] stock prices are modeled in discrete time as sequence of temporary equilibria. Those emerge as a consequence of simultaneous matching of supply and demand of several agents. Bayraktar, Horst and Sircar [6],[7], to account for this asynchronous order arrival, as well as Horst and Rothe [18] use the mathematical framework of queuing theory earlier examined by Mandelbaum, Pats et al. [26] and Mandelbaum, Massey and Reiman [25] for their models. Also the model explained in Lux [23] and Lux [24] takes asynchronous order arrivals into consideration by using a so called market maker who matches supply and demand and alters the price accordingly. Thereby the individual agents behavior is dependend on his opinion. Further Opinionbased models range from binary (e.g. Föllmer [11], Arthur [3], Orléan [27], Weisbuch and Boudjema [31] and SznajdWeron and Sznajd [29]) to opinions from a continuous spectrum, which are used, for example, to describe large social networks or ratings (see Duffant et al. [9], GómezSerrano, Graham and LeBoudec [14] or Weisbuch, Deffuant and Amblard [32]). However the characteristics and the interaction of the agents are described in the respective model, it mostly can be classified, in a wider sense, as interacting objects with assigned states forming a link to other scientific fields. Especially the application of methods and dynamics derived from physics, captured under the term of econophysics, is becoming increasingly popular and also this paper builds a bridge from physics to socio economics and finance.
We use an extented version of the microscopic market model for diffusion price processes of Pakannen [28], which is presented in Henkel [16] and apply dynamics of a quantum system from Bauer, Bernard and Tilloy [5], respectively Tilloy, Bauer and Bernard [30], to agents social behavior.
We transfer statistical properties observed in thermally activated jumps in a quantum system to the average dynamics of market participants and induce the same in diffusion price processes in the large market limit. As such we provide a microscopic explanation for jumps in asset price processes as oberserved in Aït et al. [1] without using a jump process, as for example employed by Deng [10]. Additionally our dynamics induce spikes and high volatility phases that have also present in various price processes (e.g. Ham [15]).
We structure the content as following. In the first section we set the general market framework as a pool of interacting agents. To describe their interactive behavior we assign each agent heterogeneously an excitement state similar to the excitement of a two level quantum system as in Bauer, Bernard and Tilloy [5] and express the overall market excitement as the distribution of those excitement states. Furthermore we state conditions under which, as a meanfield like result, the overall market excitement can be expressed as a single diffusion process in the large market limit. Then, with defining the agents propensity to trade and by specifying the impact on the asset price we link the endogenous market dynamics to the asset price movement. We specify conditions under which, in a large market, the asset price development can be approximated by a diffusion price process and conclude with a proposition summarizing the diffusion approximation.
1.2 Methods
The model presented in this paper was implemented in the statistical programming language R (Ihaka and Gentleman [19]) in order to simulate the related distributions and stochastic processes and to illustrate trajectories. Especially functions, which are presented as solutions of stochastic differential equations without provision of a closed analytic form, are illustrated using an implementation of the EulerMaruyamascheme (see e.g. Glasserman [13]). All Rscripts are available on request.
2 Endogenous dynamics
Before we link the agents behavior to an asset price, we specify all model components, that directly affect the interaction between market participants. Let be a finite set of agents.
Following Bauer, Bernard and Tilloy [5] we consider a state space , where represents an ”unexcited” state and an ”excited” state. The vector of all individual states takes values in the configuration space .
During the time each agent can consider a state transition. The time of the kth consideration is nominated by , . The action times are described later in detail, however we use the terminology to describe the development of the states within discrete time in the following definition.
The state of agent a at time is defined as . We capture the development of agent a’s state by the process and the development of all agents states by the ndimensional process . We assume that the vector of initial states is distributed following some ndimensional distribution function, in particular , and assign to each agent an initial state .
Definition 2.1 (Market excitement).
We measure the proportion of excited agents in the market at time by the market excitement
(1) 
Additionally, we denote the initial distribution of the market excitement resulting from as . Note that, by construction, is the average excitement of all agents and a probability measure on .
Definition 2.2 (Endogenous market history).
We capture all endogenous information up to in the endogenous market history, which is given by the sigma algebra . Here the tupel represents agent who acts at and the resulting market excitement . We assume that only one agent changes his state at a specific point in time. Although this assumption seems rather strong, it is reasonable as transitions are performed in continuous time and are unlikely to happen at the same time.
Definition 2.3 (Transition intensity).
In order to heterogeneously specify the agents tendency to consider a state transition we assign to each agent a transition intensity , that is an agent dependent constant times the number of agents participating in the market. Next, we define the probability that it is agent a who wants to reevaluate his state. Heuristically we weight the individual transition intensity by the sum of all trading intensities, which we call aggregated transition intensity and denote by . In summary,
(2) 
Next, we characterize the state transition laws, i.e. the probability that agent a changes from excited to unexcited and vice versa, given that he is the one that considers a state change, before we consequentially derive the dynamics of the market excitement.
Definition 2.4 (Transition probabilities).
We use the following notation for the two individual state transition probabilities.
(3) 
(4) 
where are agent dependent constants and
(5) 
We capture all state transition probabilities per agent in a transition matrix, i.e. we define
(6) 
Remark 2.5.
We choose this explicit form of transition probabilities presented in Equations (3) and (4) for the following reasons. The first part of the sum models individual intrinsic disposition for excitement. Thereby , respectively , captures the distance from agents ’s actual excitement state and his individual preference.
Since the average opinion can from time to time either change by or stay unchanged, it has values on the valued lattice from 0 to 1, viz.
(7) 
In summary, is a Markov chain on with value dependent transition probabilities, which are stated in the next Lemma.
Lemma 2.6 (Discrete market excitement dynamics).
The probability that any excited agent becomes unexcited and therefore that decreases by is given by
(8) 
and similarly the probability that the market excitement increases by is given by
(9) 
In order to embed the Markov chain homogeneously in continuous time and thus describing the market excitement by a time homogeneous Markov process, we further characterize the points in times at which the agents decide to make a transition.
Definition 2.7 (Transition times).
The transition times are defined as
Since we want the transition times to be memoryless for the sake of simplicity, i.e.
(10) 
the transition times are assumed to be exponentially distributed. Heuristically we assume that the rate of the exponential distribution is given by the aggregated transition intensity, i.e.
(11) 
Definition 2.8 (Market excitement index).
After fixing we can define the market excitement index via
(12) 
Note that, by construction, is cádlág and a well defined time homogeneous pure jump type Markov process. Its existence is stated in the following lemma. The basis of the lemma builds the synthesis theorem (e.g. Theorem 12.18 of Kallenberg [21]), which embeds a discrete Markov chain into continuous time using exponentially distributed waiting times.
Lemma 2.9 (Existence).
There exists a probability space in which is a time homogeneous pure jump Markov process with rate kernel
(13) 
The transition kernel is a regular version of the conditional distribution , which is given by Lemma 2.6.
Proof.
By the construction of the Markov chain and the assumption made in Definition 2.7 the synthesis theorem (e.g. Theorem 12.18 of Kallenberg [21]) states that is a pure jumptype Markov process and also gives the rate kernel. Time homogeneity is given by the recursive definition of (see e.g. Proposition 8.6 of Kallenberg [21]) as the transition matrix is independent of time. ∎
Although the heterogenous agents are allowed to have individual parameters, in order to ensure a convergence to a meanfield like single equation in the large market limit, the scaled parameters should tend to their mean when the number of market participants goes to infinity. This we summarize in the next Assumption.
Assumption 2.10.
We assume



and
for some constants and being a probability distribution.
Now, we are ready to state the large market limit for the market excitement index.
Proposition 2.11 (Large market approximation).
If Assumption 2.10 holds, then
(14) 
with being the unique strong solution of the stochastic differential equation (SDE)
(15) 
where is a one dimensional standard Brownian motion, which is independent of .
Proof.
See Appendix 5.1. ∎
To illustrate properties of the large market limit we show two trajectories of the solution of Equation (15) for , below in Figure 1 and Figure 2. The appearance of the process strongly depends on the value of . For a small , as shown in Figure 1 with , the market excitement index moves towards an equilibrium at the constant , which is given as the mean of the individual preference level (see Assumption 2.10 3.). Setting a high value of (see Figure 2, where ) the market excitement index is pulled towards the two states and with jumps and spikes in between. Although the separation into two cases (i.e. one or two equilibria) by the value of is not obvious from the underlying SDE (15), it is expected from the microscopic modeling. In the individual transition probabilities (see Equation (3) and (4)) scales down the agents individual autonomy and hence represents the exposure to herding. Respectively, reflects the average herding intensity. It was already shown in Lux [23], with a similar setup, that minor herding behavior results in a single equilibrium, while strong herd behavior results in two temporary equilibria with phase transitions. In our model the two states and serve as the two temporary equilibria, where jumps represent phase transitions and the spikes imply unsucessful jump attempts. Note that the probability to be in the equilibrium is equal to (see Tilloy, Bauer and Bernard [30]) and that the behavior is similar to Kramer’s double well potential (See Kramer [22]).
Remark 2.12.
Although the SDE presented in Equation (15) is exactly the same as in Bauer, Bernard and Tilloy [5], it arises differently. In Bauer, Bernard and Tilloy [5] the SDE is the result of a transition from discrete to continuous time in the measurement of a single quantum system, while in our model the SDE is induced by the number of interacting objects tending to infinity. Also the origin of the spikes and jumps in the large market limit is quite different. While in Bauer, Bernard and Tilloy [5] the jumps and spikes are a result of tight monitoring, in our model both is a result of agents social behavior. Here the root cause of the jumps as well as of the spikes is the agent’s individual exposure to mutual interference , that is his herding behavior. When the average herding behavior is strong (i.e. is large), the fast contagion of agents excitement leads to hypes. We show in Section 3 that these hypes result in phases of high price volatility, where transition to a hype is indicated by a jump in the price process and spikes are unsuccessful jump attempts.
3 Price dynamics
In this section we link the endogenous dynamics of the previous section with dynamics of an asset price. We assign to each trader an individual trading intensity quantifying his propensity to trade and a excess demand function which characterizes the quantity of shares bought or sold. We then define a pricing rule according to which the number of bought or sold shares impacts the price. To further show the flexibility of our model, we introduce an additional group of traders called fundamentalists. Last are characterized by basing their behavior on the difference between actual price and a fundamental value . In particular, when the price is below (above) , they consider the asset cheap (expensive) and want to buy (sell). We assume the fundamentalists are homogeneous, viz. is common, and the fundamental value is timeinvariant.
Let be the set of agents of which a fixed subset with are fundamentalists and the rest are noise traders. We denote the portion of fundamentalists with . We assume that fundamentalists are unexcited and have no desire to change their state, i.e. . Note that, alternatively we could have introduced fundamentalists as an additional state. However, the current setup illustrates the flexibility arising from heterogeneous transition probabilities.
We assume that any change of the market is a direct consequence of agents behavior. The behavior is given by actions that can either be the change of state or a trade of the asset. We index each of these actions by and extent the endogenous market history defined in Definition 2.2 with the trading behavior.
Definition 3.1 (kth action, market history).
The kth action is characterized by the tupel , where is the time when the action occurs, is the acting agent at time and is an action indicator whether the agent trades () or changes his state (). is the price per share and the above mentioned market excitement. All information is captured in the market history, which is given by .
Similar to the transition intensity (Defintion 2.3), we set a function, that captures each agents propensity to trade the asset and call the sum of both action rate.
Definition 3.2 (Trading intensity, action rate).
We assume that the agents propensity to trade is given by the trading intensity
(16) 
with being an agent dependent basic trading intensity and a positive constant, which reflects the positive impact of excitement on the propensity to trade.
Moreover we introduce the action rate for each agent with
(17) 
The aggregated action rate is then given as
(18) 
Analogue to Equation (2) we identify the acting agent and the related action by weighting the respective intensity functions.
Definition 3.3.
(Acting probabilities)
The probability, that agent a trades at is defined as
(19) 
Similarly, we define the probability, that agent a changes his state by
(20) 
Moreover the probability that the kth action is a state transition is set as
(21) 
and analogously the probability that the kth action is a trade is given by
(22) 
Next we define the traded quantity per agent once he decided to trade. Thereby we differ between noise traders and fundamentalists. While fundamentalists base their excess demand on the difference between the last known price and the fundamental value, noise traders trade according to random signals which are assumed to be i.i.d. with and . Thereby the variance of the traded quantity is determined by the variance of the market excitement .
Definition 3.4 (Excess demand function).
In summary we set the following excess demand function
(23) 
After an agent decides to trade, the new price at time will be set by a market maker, who is assumed to handle all trades, and is defined by a pricing rule depending on the excess demand of the acting agent and the old price. We assume that the difference of old and new price is linear in the quantity traded rootscaled by the number of market participants. That is
(24) 
where
(25) 
By construction and are now two interacting Markov chains. In order to embed them homogeneously in continuous time and thus describing the price as well as the character by a time homogeneous Markov process, we further characterize the points in times at which the agents decide to act.
Definition 3.5 (Intraaction times).
The intraaction times are defined as
Analogously to Definition 2.7 we assume that the rate of the exponential distribution is given by the aggregated action rate, i.e.
(26) 
More precisely, to ensure a sufficient level of independence between the source of randomization and the price as well as market character we need to assume that the Intraaction times are given by
(27) 
where are i.i.d. random variables independent of with .
Definition 3.6 (Price process).
After setting an initial price and fixing we can define the price process as
(28) 
We extent Lemma 2.9 with the Price process in the following.
Lemma 3.7 (Existence).
If the preceding Assumptions hold true, then a probability space exists, which carries the model in a way that is a time homogeneous pure jump Markov process with rate kernel
(29) 
where the transition kernel is a regular version of the conditional distribution
,
Proof.
Analogous to Lemma 2.9. ∎
Before we can state the large market limit for the market excitement index and the price process, we assume some stability of the proportion of fundamentalists. Additionally we require a meanconvergence of the trading intensities of the fundamentalists as well as of the noise traders.
Assumption 3.8.
We assume
for some constants and .
Next, we state the SDE whose solution approximates the endogenous dynamics and the price process in a large market, that is a market with many participants.
Proposition 3.9 (Diffusion approximation).
(30) 
where is the unique strong solution of the SDEs
(31) 
where and are independent one dimensional standard Brownian motion, independent of , and independent of .
Proof.
See Appendix 5.2. ∎
Equation (31) summarizes our model in the large market limit. The endogenous behavior is described by given by the first SDE, which is not depending on the price process and is the same as in the previous section. On the contrary depends on . Not only the volatility coefficient of reappears in the SDE defining , but also scales the volatility of with the factor . Last leads to high volatility phases when the majority of agents is excited. To illustrate the properties of we show two trajectories. Thereby we repeat the figures of from the previous section for readers convenience.
In Figure 3 and 4 we show the first case with a trajectory of with , , and endogenous dynamics, that is , with parameters , and . Driven by 20% of the agents being fundamentalists,
drifts to the fundamental value. Thereby the volatility is rather stable, since has a single equilibrium at and the rest of the volatility coefficient of consists of constants. We illustrate the second case in Figure 5 and Figure 6 with the same parameters but setting . There the spikes and phase transitions from are transferred to the price process and result in spikes and jumps.
Moreover, s explained above, the phases of temporary equilibrium of at comply with high volatility phases of , since the factor increases the volatility coefficient of when is large. The intensity of the effect of last is specifically steered by the constant .
Remark 3.10 (Approximation with Poisson Jump Process).
Since the SDE of in Equation (31) is exactly the same as in Tilloy, Bauer and Bernard [30], we can leverage from a result on statistical properties presented there in Proposition 2. The position of spikes and jumps (that is positions of local maxima and minima) can be approximated, when is large, by two Poisson point processes and on with intensities
(32) 
where is the delta function.
Remark 3.11.
While the statistics of are discussed in Bauer, Bernard and Tilloy [5], respectively Tilloy, Bauer and Bernard [30], the statistics of (especially the structure of the spikes) were not studied yet although they are a direct consequence of . One might use the Fokker Planck approximation to approximate the stationary distribution of , respectively , to get more insight. Additionally, in order to study , it might also be worth to study the underlying Markov chain . However both is out of this articles scope.
4 Conclusion and Outlook
We have proposed a microscopic agentbased model to explain jumps, spikes and high volatility phases in diffusion price processes. Using the mathematical framework of Henkel [16] we set a finite network of heterogenous agents interacting in continuous time. The agents behavior is thereby inspired by the dynamics of excited particles in a quantum system (see Bauer, Bernard and Tilloy [5]). In a second step we link the endogenous dynamics to an asset price process by specifying agents individual trading propensity and excess demand functions together with an overall pricing rule. Furthermore, we showed the conditions under which the average agent excitement as well as the price process converge to a diffusion process when the number of market participants tends to infinity. Since our model induces large market dynamics that are likewise present in the discussion of quantum systems coupled to a thermal bath with continuous monitoring (see Bauer, Bernard and Tilloy [5]) we build a bridge between quantum mechanics and financial mathematics. So we could leverage from the statistical properties of quantum trajectories and apply a result of Tilloy, Bauer and Bernard [30] to our asset price model by which the occurring jumps and spikes can be approximated by two Poisson processes.
For the sake of simplicity several assumptions have been made, that also show limitations of the model. For instance, the missing feedback of the price process on the endogenous dynamics as well as the strong Markov property of the model seem unrealistic. Although the first can be addressed rather easy with more complicated state transition probabilities which consider also the asset price, the Markov property is critical for the convergence to a diffusion process. More complicated microscopic models leading to nonmarkovian limits, e.g. solutions of stochastic differential delay equations (see Arriojas et al. [2]), could be investigated in the future, although far less literature in form of limit theorems is available for this case.
5 Appendix
5.1 Proof of Proposition 2.11
We apply Theorem 3.6 of Henkel [16]. Note that our model fits in the framework of Henkel [16] with the market character and market charecter index, which are defined in Definition 2.3 and 2.16 in Henkel [16], given by , and . The expected aggregated transition and transition volume (Henkel [16] Definition 3.3 and 3.4) necessary to determine the limit of , respectively , when are given by
(33) 
(34) 
Now, by Assumption 2.10
(35) 
and
(36) 
So, by Theorem 3.6 of Henkel [16] , where is the solution of
(37) 
and hence , where is the solution of
(38) 
since .
5.2 Proof of Proposition 3.9
Also here we apply Theorem 3.6 of Henkel [16]. Since the dynamics of do not depend on , the convergence of and its limit is given by Proposition 2.11. To also show the convergence of and to determine the limit when , we calculate the expected aggregated excess demand and the trading volume (see Henkel [16] Definition 3.1 and 3.3).
(39) 
(40) 
(41) 
and
(42) 
After realizing that Assumption 3.5 of Henkel [16] is fulfilled, we apply Theorem 3.6 of Henkel [16] and get Proposition 3.9 as a result.
Footnotes
References
 Y. AïtSahalia, J. CachoDiaz, and R. J. Laeven. Modeling financial contagion using mutually exciting jump processes. Journal of Financial Economics, 117(3):585–606, 2015.
 M. Arriojas, Y. Hu, S.E. Mohammed, and G. Pap. A delayed black and scholes formula. Stochastic Analysis and Applications, 25(2):471–492, 2007.
 W. B. Arthur. Increasing returns and path dependence in the economy. University of Michigan Press, 1994.
 L. Bachelier. Théorie de la spéculation. GauthierVillars, 1900.
 M. Bauer, D. Bernard, and A. Tilloy. Computing the rates of measurementinduced quantum jumps. Journal of Physics A: Mathematical and Theoretical, 48(25):25FT02, 2015.
 E. Bayraktar, U. Horst, and R. Sircar. A limit theorem for financial markets with inert investors. Mathematics of Operations Research, 31(4):789–810, 2006.
 E. Bayraktar, U. Horst, and R. Sircar. Queuing theoretic approaches to financial price fluctuations. Handbooks in Operations Research and Management Science, 15:637–677, 2007.
 F. Black and M. Scholes. The pricing of options and corporate liabilities. The journal of political economy, pages 637–654, 1973.
 G. Deffuant, D. Neau, F. Amblard, and G. Weisbuch. Mixing beliefs among interacting agents. Advances in Complex Systems, 3(01n04):87–98, 2000.
 S. Deng. Stochastic models of energy commodity prices and their applications: Meanreversion with jumps and spikes. University of California Energy Institute Berkeley, 2000.
 H. Föllmer. Random economies with many interacting agents. Journal of mathematical economics, 1(1):51–62, 1974.
 H. Föllmer and M. Schweizer. A microeconomic approach to diffusion models for stock prices. Mathematical Finance, 3(1):1–23, 1993.
 P. Glasserman. Monte Carlo methods in financial engineering, volume 53. Springer Science & Business Media, 2003.
 J. GómezSerrano, C. Graham, and J.Y. Le Boudec. The bounded confidence model of opinion dynamics. Mathematical Models and Methods in Applied Sciences, 22(02), 2012.
 J. D. Hamilton and G. Lin. Stock market volatility and the business cycle. Journal of Applied Econometrics, 11(5):573–593, 1996.
 C. Henkel. An agent behavior based model for diffusion price processes with application to phase transition and oscillations. arXiv preprint arXiv:1606.08269, 2016.
 U. Horst. Financial price fluctuations in a stock market model with many interacting agents. Economic Theory, 25(4):917–932, 2005.
 U. Horst and C. Rothe. Queuing, social interactions, and the microstructure of financial markets. Macroeconomic Dynamics, 12(02):211–233, 2008.
 R. Ihaka and R. Gentleman. R: a language for data analysis and graphics. Journal of computational and graphical statistics, 5(3):299–314, 1996.
 J. Jacod and A. N. Shiryaev. Limit theorems for stochastic processes. Springer, 2003.
 O. Kallenberg. Foundations of modern probability. Springer, 2002.
 H. A. Kramers. Brownian motion in a field of force and the diffusion model of chemical reactions. Physica, 7(4):284–304, 1940.
 T. Lux. Herd behaviour, bubbles and crashes. The economic journal, pages 881–896, 1995.
 T. Lux. Time variation of second moments from a noise trader/infection model. Journal of Economic Dynamics and Control, 22(1):1–38, 1997.
 A. Mandelbaum, W. A. Massey, and M. I. Reiman. Strong approximations for markovian service networks. Queueing Systems, 30(12):149–201, 1998.
 A. Mandelbaum, G. Pats, et al. Statedependent stochastic networks. part i. approximations and applications with continuous diffusion limits. The Annals of Applied Probability, 8(2):569–646, 1998.
 A. Orléan. Bayesian interactions and collective dynamics of opinion: Herd behavior and mimetic contagion. Journal of Economic Behavior & Organization, 28(2):257–274, 1995.
 M. S. Pakkanen. Microfoundations for diffusion price processes. Mathematics and Financial Economics, 3(2):89–114, 2010.
 K. SznajdWeron and J. Sznajd. Opinion evolution in closed community. International Journal of Modern Physics C, 11(06):1157–1165, 2000.
 A. Tilloy, M. Bauer, and D. Bernard. Spikes in quantum trajectories. Physical Review A, 92(5):052111, 2015.
 G. Weisbuch and G. Boudjema. Dynamical aspects in the adoption of agrienvironmental measures. Advances in Complex Systems, 2(01):11–36, 1999.
 G. Weisbuch, G. Deffuant, and F. Amblard. Persuasion dynamics. Physica A: Statistical Mechanics and its Applications, 353:555–575, 2005.