The Role of the Volatility in the Option Market

: We review some general aspects about the Black–Scholes equation, which is used for predicting the fair price of an option inside the stock market. Our analysis includes the symmetry properties of the equation and its solutions. We use the Hamiltonian formulation for this purpose. Taking into account that the volatility inside the Black–Scholes equation is a parameter, we then introduce the Merton–Garman equation, where the volatility is stochastic, and then it can be perceived as a ﬁeld. We then show how the Black–Scholes equation and the Merton–Garman one are locally equivalent by imposing a gauge symmetry under changes in the prices over the Black–Scholes equation. This demonstrates that the stochastic volatility emerges naturally from symmetry arguments. Finally, we analyze the role of the volatility on the decisions taken by the holders of the options when they use the solution of the Black–Scholes equation as a tool for making investment decisions.


Introduction
A financial market is defined as the place where two or more individuals buy and sell financial assets and make general transactions.The stock market is a well known financial market, and it consists of trading shares in the stocks of some company or organization [1].Inside the stock market, some derivative instruments are also traded.One important derivative instrument traded in the market is the option.Options were traded for the first time in 1973 [2].Interestingly, the same year, Black and Scholes introduced their famous Black-Scholes (BS) formula [3].This formula gives the fair price of an option as a function of the stock price, time to maturity, strike price, interest rate and volatility [2].Since the formula is dynamical, then the time as a variable also appears inside the solution.The BS formula is by itself a differential equation containing two free-parameters, namely, the volatility and the interest rate [4].The brilliance of the BS formulation is in the fact that it proposed a way to cancel random fluctuations by mixing up inside a single portfolio the prices of a security (stock price) with the prices of a derivative instrument.This cancellation of the random fluctuations is what makes the BS equation a predictive one [3,4].Different methods for analyzing the option prices have been developed previously.Among them, we have the fractional-order option pricing [5] and the Laplace homotopy analysis method (LHAM) using the Caputo-Fabrizio (CF) fractional derivative operator [6].Additionally, some ideas trying to modify the BS equation and compare it with variance models were proposed in [7].In general, the BS equation, as a differential equation, can be expressed as a Hamiltonian (eigenvalue) equation after some variable transformations [4].This formulation gives us significant advantages with respect to the standard differential equation formulation.Among the numerous advantages of the Hamiltonian formulation, we have the following: (1).The problem becomes an eigenvalue problem such as in quantum mechanics [8].(2).We can identify the free-parameters that control the dynamic of the system, which, in the case of the BS equation, are the volatility and interest rate.(3).We can analyze the symmetries of the system and then identify symmetry breaking patterns [9][10][11].(4).We can identify the operators generating transformations in certain quantities, such as the prices of the stocks, for example, or even the volatility if we consider the Merton-Garman equation [12,13].(5).We can simply use analogies with quantum mechanics (QM) in order to develop several interesting new tools inside the financial market [4,[8][9][10][11][14][15][16][17].This does not mean that the financial market is quantum; this rather means that we can use certain analogies, with the advantage of knowing in advance the form of the solutions for certain equations.However, we must be careful at the moment of working over these analogies because there are certain differences between the Hamiltonians appearing in ordinary quantum mechanics [8] and the financial Hamiltonians, such as in the BS case.One strong difference is that in ordinary QM, the Hamiltonians are normally Hermitian [8], the Hamiltonians in finance, including the BS Hamiltonian, are generically non-Hermitian [4,[9][10][11].Then we have loss of unitarity in the stock market, or equivalently, the information is not preserved inside the stock market [10,17].Although, in QM, non-unitarity is considered to be a serious problem [8], inside the stock market, the loss of information is, in some sense, expected.However, even in ordinary physics, there are cases where non-Hermiticity is valid, as far as the PT symmetry is satisfied [18][19][20][21][22][23][24].However, this is not the case inside the stock market.In this paper, we focus our analysis on the role of the volatility in the option market.For this purpose, we show the solutions for the BS equation and we explain how the investors use them to decide whether they should buy an option inside the market or not.From this perspective, we explain how the investors estimate the volatility value and how they compare it with the value obtained from the implied volatility [2].Finally, given the importance of the volatility inside the stock market, in this paper, we revise the seminal result demonstrating that the BS equation is locally equivalent to the Merton-Garman (MG) equation [11].This result is a natural consequence of imposing the local symmetry condition under changes in the prices of the stocks over the BS Hamiltonian.In this way, the stochastic volatility emerges as a field necessary for satisfying a local symmetry over the BS equation.The local equivalence of the BS and the MG equations is a very important result because it means that the volatility as a variable emerges from a series expansion around the BS solution.The paper is organized as follows: In Section 2, we revise the concept of the option market.We then explain the difference between call option and put option.In Section 3, we explain some generalities about the BS equation.We then derive its Hamiltonian form and then use the corresponding Hamiltonian for calculating certain quantities.In Section 4, we analyze the Merton-Garman equation (MG), which is an extension of the BS case in order to include the volatility field as a stochastic variable.We then explain how to obtain the corresponding Hamiltonian, and then we use it for calculating certain quantities.
In Section 5, we analyze the interesting mechanism of spontaneous symmetry breaking inside the scenario of quantum finance by using the BS equation as an example where the mechanism appears.In Section 6, we demonstrate the local equivalence between the BS equation and the MG one.In Section 7, we explain the Higgs mechanism describing the dynamical origin of the volatility.In Section 8, we analyze the solutions of the BS equation, which is used by the investors in general for making decisions in the market.Finally, in Section 9, we conclude.

The Option Market
The options were created and traded for the first time in 1973.Interestingly, during the same year, the famous Black-Scholes (BS) formula was introduced [2].This equation has been used since then by investors and academics for determining the fair price of an option.Before exploring the details of the BS equation, we have to understand how the option market operates.The option market usually has a holder of the option and the writer of the same option [4].There are two types of options: (1).The call option, where the holder has the right to buy the shares of a stock at some predetermined price.(2).The put option, where the holder has the right to sell the shares of some stock at some predetermined price.
The stock market charts, such as those appearing in Yahoo Finance [25], give us the price S(t) of different stocks corresponding to different companies.An option is nothing else than a contract signed between the holder (buyer of the option) and a writer (seller of the option).For the call options, the holder expects the prices of the stock to increase (bullish), and the writer expects the same prices to decrease (bearish).This is the case because for a call option, the holder has the right to buy a stock in future at the same price as the stock is priced today (at the moment of signing the contract).Then the holder is interested in a tendency to increase the value of the option that he holds.
The holder and the writer of an option, at the most fundamental level, are playing a zero-sum game, as it can be perceived when we compare the Figures 1 and 2. Indeed, Figure 1 is the mirror image of Figure 2, showing in this way that the profits of the holder are the losses of the writer and vice versa.In this particular example, the holder buys a call option at 5 USD per share.This option gives him the right to buy some specific stock of some company at the strike price of 42 USD.For the holder to exercise the option, he/she needs it to have at least a stock price larger than 47 USD per share.This is the case because the holder bought the option at 5 USD, starting then with a negative balance of 5 USD per share, as is illustrated in Figure 1.He/she has this negative balance when he/she does not exercise the option, which is the case when the stock price is inferior to 47 USD per share.If by the expiration date of the Option, the price is larger than 47 USD per share, then the holder earns some positive profits after exercising the option.It is at this point where the writer starts to lose profits, as can be seen from Figure 2, which analyzes the profits of the writer for the same option and stock price.Note that the writer earns positive profits if the option is not exercised by the holder.In summary, the holder earns money when he exercises the option after verifying that the price of the same option is larger than 47 USD (for closing prices of the stock between 47 USD and 42 USD, the European option could be also exercised in order to minimize the losses suffered by the holder).What he earns (the holder) is the difference between the price of the stock minus what he/she paid for the option, multiplied by the number of shares involved in the process.On the other hand, the writer makes money when the holder does not exercise the option after verifying that the stock price closes at a value inferior to 47 USD (more strictly for stock values inferior to 42 USD because for stock prices between 42 USD and 47 USD, the holder could still exercise the European call option for minimizing his/her losses).Here, we want to remark that while the maximum possible profit for the writer is 5 USD per share (see Figure 2), the maximum profit for the holder is just limited by the price of the stock involved in the process.We can now analyze the situation with the put option.A put option is one where the holder, once he/she buys it, he/she has the right to sell a stock at some specific prearranged price.This means that the holder of the option pays to the writer of the option some amount in order to obtain the right to sell a stock at a pre-specified price, which is usually the price of the stock on the day of buying the option.In this case, the holder of the option wants the price of the stock to fall (bearish), while the writer of the option wishes for the prices to increase (bullish).Figure 3 illustrates an example for the possible profits earned by the holder of a European call option.We can notice that the lower the price of the stock is, the larger the earnings of the holder.In this particular example, the cost of the option is 3 USD.Then naturally, the holder loses 3 USD per share and then the writer of the same option earns 3 USD per share, as can be seen from Figure 4. Figures 3 and 4 are mirror images of each other.This is a typical relation for the pay-offs between two players in a zero-sum game.Inside the European call option, for this example, the holder receives positive pay-offs when the price of the stock is inferior to 47 USD.This is the case because the holder has the right to sell the stock at 50 USD, no matter what happens with the prices of the stock.Since he/she has to recover the investment of 3 USD, then he/she only exercises the option if the price of the involved stock falls below 47 USD (The holder could also decide to exercise the European put option for prices of the stock between 47 USD and 50 USD, if by the time of closing, he/she wishes to minimize his/her losses).For prices of the stock larger than this value (47 USD), the holder will still possibly exercise the option if the price of the stock is inferior to 50 USD, if only to minimize losses.However, if the prices of the stock increase over 50 USD, exercising the option brings losses larger than 3 USD per share to the holder.Then the option in this case is never exercised and the writer earns 3 USD per share, while the holder loses the same amount.For this particular example, the highest pay-off for the holder is 47 USD, while the writer also loses the same amount in this situation.Here again, the highest possible pay-off received by the writer is just the price of the option per share, and in the same situation, the holder would lose the same amount.

The Black-Scholes Equation
In this paper, we consider the stock price as S(t), which is normally taken as a random stochastic variable evolving in agreement to a stochastic differential equation, given by Here, φ is the expected return of the security, R(t) is the Gaussian white noise with zero mean and σ is the volatility.If we exclude the volatility, then we can guarantee the evolution of the price of the stock with certainty [4].In this way, by imposing σ = 0, we obtain a simple solution for Equation (1) as Additionally, the possibility of arbitrage is excluded if we can make a perfect hedged portfolio.The key ingredient inside the analysis of Black and Scholes was to eliminate any possibility of uncertainty, such that we can analyze the evolution of the price free of any white noise [4].We can then consider the following portfolio: This is a portfolio where an investor holds the option and then short sells the amount ∂ψ ∂S for security S. The price of the option C and its derivative ∂C ∂S are naturally correlated, and then, any random fluctuation can be eliminated through this correlation.By using the Ito calculus (stochastic calculus) [4], it is possible to demonstrate that Here, the change in the value of Π does not have uncertainty associated with it [4].Any random fluctuation has vanished because the random term has disappeared due to the choice of portfolio.Since, in this case, we have a risk-free rate of return (no arbitrage) [26,27], then the following equation is satisfied If we use the results (3) and ( 4) together with the previous equation, then we obtain This is the Black-Scholes equation and its details and assumptions can be found in [4,12,13].

Black-Scholes Hamiltonian Formulation
Here, we explain how Equation ( 6) can be expressed as an eigenvalue problem after a change in variable.The resulting equation is the Schrödinger equation with a non-Hermitian Hamiltonian [8].For Equation (6), consider the change of variable S = e x , where −∞ < x < ∞.In this way, the BS equation becomes where we have defined the operator as the BS Hamiltonian [4,[9][10][11].Note that the resulting Hamiltonian is non-Hermitian since Ĥ = Ĥ+ [14,15].If we compare with the standard Schrödinger equation, defined as [8] Ĥψ then it is evident that the BS Hamiltonian is defined in Equation ( 8), and that this Hamiltonian can help us to understand the evolution of the prices of the options.The result (8) has important consequences when we analyze the probability flow inside the scenario of the BS equation [10].In addition, note that since the spot interest rate r is constant, then the potential term is simply a constant term.This means that the vacuum condition is trivial for this case.Under the BS Hamiltonian, the evolution in time of the option is non-unitary in general (in addition, the Hamiltonian non-necessarily obeys the PT symmetry).This means that the probability is not necessarily preserved in time inside the BS model, although it is certainly well-defined, and its total value is equal to one, if the normalization factor is updated from time to time.In general, as was mentioned in the introduction, there are some cases in ordinary quantum mechanics, as well as in quantum field theory, where it is interesting to explore non-Hermitian Hamiltonians (Lagrangians) [18,20,28].Based on the previous explanations, we cannot expect the financial market to obey unitarity.This is precisely the point that we want to analyze in deep detail in this paper.

The Evolution of Probability in the Black-Scholes Equation
The solutions for the BS equation are well-known [2,4], and we explore them subsequently.Before analyzing the solutions, we study the evolution of the prices of the options, taking advantage of the Hamiltonian formulation.We can then find the kernel for the evolution of the prices in the stock market.The kernel naturally depends on the Hamiltonian under analysis.To understand this issue, we can solve the Schrödinger equation defined in Equation ( 9), obtaining then Here E is the eigenvalue for the Hamiltonian Ĥ.The same result can be expressed in the bra-ket notation as Considering the final condition at t = T, we have where g can be interpreted as the pay-off function if we assume |ψ >= |C >.Under this equality, the initial option price is then defined as and then with the time defined as τ = T − t and then running backwards.This means that in this formulation, the initial condition t = 0, is equivalent to τ = T, which, inside the stock market scenario, corresponds to the maturity time for the option.The solution (14) represents an exponential decaying behavior.If we project the result ( 14) over < x|, then we obtain the functional behavior, as follows: Note that in the last step, we have used the completeness relation for the basis of the stock prices, represented by the variable x.Then, we define the pricing kernel as

Black-Scholes Pricing Kernel
By introducing the Hamiltonian (8) inside Equation ( 16), we should be able to find the explicit form of the pricing kernel for the BS case.However, we have to develop a few definitions initially, as follows [4]: The completeness relation, is satisfied and this means that we are taking the basis of the prices as a complete basis.
With the previous results, we can define the product as The BS kernel can be obtained by using Equation ( 16) with the BS Hamiltonian Ĥ = ĤBS .We can re-express Equation ( 16) as If we want to find the explicit result, we need to evaluate the following matrix element explicitly: Here, E BS is the eigenvalue of the Hamiltonian operator corresponding to the state |p .We can find the explicit value of E BS if we evaluate the following matrix element: Here, we are using p|p = p|p because in this case, the operator ĤBS is acting over the momentum space.In position space, on the other hand, the momentum operator is a derivative operator.If we use the result (22), then the kernel defined by Equation ( 20) becomes Here again, τ = T − t.The kernel (23) allows us to analyze the probability flow in the market, and then we can analyze the non-conservation of probability in this case.First, we can evaluate the time-derivative of the kernel as [10] ∂p BS ∂t Since ∂τ/∂t = −1, then we obtain in general.The previous derivative is equal to zero for some relations between the variables and between parameters.In fact, we obtain We also have the additional condition It is interesting to notice that in Equation ( 26), the risk-free interest rate takes negative values under the imposed limit.In fact, r can be either positive or negative [29].However, in order to satisfy the demands of the BS equation, we impose the condition r ≥ 0, avoiding negative risk-free interest rates.From Equation ( 26), we can see that an infinite maturity time is necessary, under the condition r = −σ 2 /2, in order to preserve the probability in the market, something which can only be approximated in some circumstances.In Equation ( 27), the condition for preserving the probability is the logical fact that for small scales of maturity time, we can know with certainty the market behavior.This, in addition, implies that x → x , which means that the stock prices do not change so much under the conditions proposed in Equation (27).On the other hand, since Equation ( 26) allows the possibility of having negative values for r, there must be a limit value for the time τ where r = 0.This important limit can be found if we impose r → 0 in the time-derivative of the kernel (23).Assuming again the conservation of probability ( ∂p BS ∂t = 0), we obtain the polynomial If we solve for τ, then we obtain Here, the lowest sign is the correct one in order to keep τ positive.The limit value τ r=0 min = 0 (t = T, final value considering that the evolution of the system is analyzed from the final condition toward the initial one) again corresponds to the condition x → x , which is consistent with the fact that x corresponds to the final price.Note that when x >> x , namely, when (x − x ) → ∞ such that the price at maturity is much larger than the initial price of the stock, then we have the result Another important limit to be considered is r = σ 2 2 because it represents the limit when the BS system preserves the information, because in this limit, the BS Hamiltonian becomes Hermitian.For this special limit, the condition Since we know again that the condition τ ≥ 0 must be satisfied, we select the lowest sign in Equation (31).Note that for the minimal value of τ r=σ 2 /2 min = 0, we obtain This limit corresponds to the case where the stock is close to maturity.If, additionally, we assume r = 0 and σ = ∞ (finite volatility), then we obtain which again shows the consistency of the present formulation.Here again, we can predict the largest value for τ (the largest maturity time) by considering the limit |x − x | → ∞.
For this case, we obtain where we have used r = σ 2 /2 for this special case.The limit r = σ 2 /2 is extremely important because it corresponds to the case where the BS Hamiltonian becomes Hermitian, and for this case, the information (unitarity) should be preserved in principle.The currents of probability associated with the BS equation were analyzed in [10].

The Influence of the Non-Conservation of the Probability in the Prices of the Stocks under the Black-Scholes Formulation
The prices in the stock are represented by the variable x.This variable then affects the behavior of the option prices, here represented by C(x, t).The way this function evolves is contained in the kernel defined in Equation (23).It is evident that in the stock market, the probability is not conserved.For this case, the conservation of probability only appears with respect to the analytically extended time µ for the case r = σ 2 2 [10].This special condition appears when the non-Hermitian portion of the Hamiltonian ĤBS vanishes.Then, under the Hermiticity condition, the kernel (23) becomes This result corresponds to a standard Gaussian distribution with center x and with timedependent height.If we compare with Equation ( 23), we notice that the effect of the non-Hermitian portion of the BS Hamiltonian is to create a shift in the location of the center of the Gaussian.In such a case, the center of the Gaussian is also time-dependent.Under analytical extension, we can perceive from Equation ( 35) that the Kernel corresponding to the Hermitian portion of the BS Hamiltonian is just a plane-wave with a time-dependent amplitude (decreasing amplitude actually).Under the same analytical extension, the full kernel (23) has some additional time-dependence on the amplitude of the function.
Then, we can conclude that in fact, the loss of information due to the presence of the non-Hermitian contributions in the BS Hamiltonian affects the location of the most probable price of the stock to be observed, as well as the total effect of the volatility

The Merton-Garman Equation
In this section, we consider the MG equation.Inside the MG scenario, the volatility and the price of the stock are both stochastic.In this situation, the market is incomplete because the volatility cannot be traded in the market [4].Although several stochastic processes have been considered for modeling the case with stochastic volatility [30][31][32][33][34][35][36], here, we consider the generic case, defined by the set of equations [4] Here, the volatility is defined by the stochastic variable V = σ 2 , and S is the stochastic variable representing the price of the stock.Additionally, φ, λ, µ and ζ are free-parameters of the system [26].R 1 is the Gaussian noise for the price of the stock, and R 2 is the same but for the case of volatility.R 1 and R 2 are correlated in the following way: Here, −1 ≤ ρ ≤ 1, and the bra-kets < AB > correspond to the correlation between A and B. If we consider a function f , depending on S(t), explicitly on the time, as well as on the white noises R 1 and R 2 , and with the help of the Ito calculus, it is possible to derive the total derivative in time of this function as The compact expression for this equation is with the corresponding definitions This is the same notation used in [3].

Derivation of the Merton-Garman Equation
If we consider two different options defined as C 1 and C 2 on the same underlying security (stock price) with strike prices and maturities given by K 1 , K 2 , T 1 and T 2 , respectively, then it is always possible to create a portfolio, that combines both Options.If we consider the result (39), then the total derivative for the folio with respect to time is This result is obtained after defining f (t) = C 1 or f (t) = C 2 in Equation (39).It has been demonstrated that even in this particular case, where we have stochastic volatility, it is still possible to generate a hedged folio, and then at the end, we arrive again to the condition dΠ dt = rΠ, which is the typical one satisfying the requirements of complete markets.This important result is obtained after removing the white noises R 1 and R 2 .This is implemented after imposing the following conditions: and then solving for Γ 1 and Γ 2 .The solution for Π is non-trivial for this case, and then it requires the definition of the parameter This parameter does not appear for the case of the Black-Scholes equation for example.This is the case because β in the MG equation is the market price volatility risk.This can be understood after observing the behavior of the parameter β.In this way, the higher the value of β is, the lower the intention of the investors to risk.This is a natural consequence of the fact that in the MG equation, the volatility is a stochastic variable.Since the volatility is not traded in the market, then it is not possible to make a direct hedging process over this quantity [4].Then, when we have stochastic volatility, the expectations of the investors are important and they are reflected on the parameter β.Previously, it was demonstrated in [37] that the value of β in agreement with Equation ( 44) is always a non-vanishing result.
In general, it is always assumed that the risk of the market (in price) has been included inside the MG equation.The MG equation is then obtained by rewriting the Equation (44) in the form In this equation, the effects related to the intentions of the investors, which were originally stored in β, now appear contained inside the parameter λ.In other words, we have shifted the parameter λ → λ − β in Equation (45).In what follows, we transform the MG equation to its Hamiltonian form, which is very convenient for analyzing some interesting properties of the MG equation.

Hamiltonian Form of the Merton-Garman Equation
The MG equation can be also formulated as a Hamiltonian (eigenvalue) equation after performing the appropriate redefinition of variables.The changes in variables in this case are and then, the MG Equation ( 45) becomes [4,16,38] We can interpret this equation here as an eigenvalue equation.In such a case, we have with the MG Hamiltonian defined as Some exact solutions for the MG equation have been found for the special cases where α = 1, as was implemented in [16] by using path-integral techniques.The same equation has also been solved in [30][31][32][33][34][35][36] for the case α = 1/2 by using standard techniques of differential equations.The equation has two degrees of freedom, and they are transmitted to its Hamiltonian.The MG Hamiltonian, as is the case of the BS Hamiltonian, is non-Hermitian, and the non-conservation of probability is expected.

The Merton-Garman Kernel
The flow of probability for the MG case can be analyzed when we again calculate its kernel in the same way it was calculated for the BS case [4].The kernel describes the evolution of the system between two well-defined prices, which are represented by a variable x, as follows with T representing the maturity price.By following the same techniques explained in [4], we find Here, p • x = p x x + p y y, with x connected with the prices and y related to the volatility in the market.In Equation ( 49), we can see that a portion of the MG Hamiltonian looks like the BS one with σ 2 = e y .Then, one portion of the kernel looks similar to the BS kernel.We can now define the MG kernel in agreement with Equation (50) as Introducing Equation ( 51) inside the previous result, we obtain e −τρζe y(α−1/2) p x p y .( This kernel contains all the time-dependence of the prices of the stock plus the fluctuations in volatility (variations in time of volatility).

The Flow of Probability in the MG Case
The evolution of the prices and all the dynamics in general is concentrated in the kernel.Then, it is useful to evaluate the possible variations in the kernel with respect to time.Here again, since τ = T − t, then we have ∂p MG /∂t = −∂p MG /∂τ.In [17], it was demonstrated that for the case of the stock market, the following relation is valid: This equation expresses the evolution of the probability for some price of the option to occur.This probability is defined as ρ = |C(x, t)| 2 .Here, C(x, t) is a functional distribution of the prices of an option.If we use the Schrödinger equation defined in (48), we obtain If we integrate over the whole range of stock prices, as well as over the whole range of possible volatility values (expressed through the variable y), then we obtain the following expression: This equation suggests that in general, the probability is not conserved in the stock market.
Even under asymptotically well-behaved solutions for the prices of the option, we would obtain approximately the solution which is equal to the Black-Scholes case under the same conditions [17].This part is the generic non-conservation of probability because it appears in any Financial equation, no matter how the volatility behaves.Equation ( 56) has the solution up to some constant, normalizing the result.The information is, however, not preserved here because the normalization factor changes in time.This means that the total probability at one instant is normalized to be equal to one; however, at a subsequent time, the same probability is not normalized anymore under the same normalization factor.Then, a new normalization factor would be required for keeping consistency.This problem is called loss of unitarity in quantum mechanics [8], and in general, it represents a big problem in fundamental science.However, at the stock market level, it is expected that the information is not preserved, and we should not worry so much about the loss of unitarity.The unitarity of any system is restored when the Hermiticity of the Hamiltonian is restored [8].Then, if we analyze the Hamiltonian of the MG equation defined in Equation ( 49), we understand the the non-Hermitian contributions disappear when the conditions are satisfied.These conditions then guarantee the conservation of probability but not with respect to the ordinary time but with respect to the imaginary time-coordinate, here defined as t → −iω, as it can be easily proven [10,17].The result (58) can be summarized in the following expression: For the special case where α = 1, we can then solve the result for r, as follows: Additionally, for the special case with α = 1/2, we obtain The results (60) and (61) correspond to the special values that the risk-free interest rate has to take for some special values of α such that the market can still preserve unitarity or, in other words, preserve the probability.For more general values of α, we only have to solve the Equation (59) with the corresponding value taken by α.

Spontaneous Symmetry Breaking from the Black-Scholes Hamiltonian
Spontaneous symmetry breaking is a phenomenon taken from physics [39,40].It consists of the fact that under certain circumstances, the ground state of a system can violate certain symmetries that are still satisfied by the Hamiltonian.The typical example is the Mexican hat case as it appears in the Figure 5.For illustrating this example, we can consider the following Hamiltonian corresponding to the linear σ-model [39,41]: From this Hamiltonian, the potential term corresponds to the expression ignoring the kinetic terms, as follows: The Hamiltonian (62), and then the potential (63), are invariant under the symmetry transformation φ i → R ij φ j .The ground state for the potential (63) is The field φ i in this example has N-components.Then, the vacuum state condition (64) does not specify the direction taken by the ground state.It only specifies its magnitude.For example, we could select a ground state consistent with (64), as follows: Then, the ground state selects a preferred direction, and it does not respect the symmetry of the original Hamiltonian (62).From the ground state (65), we can see that the system could perfectly select N-different components for the ground state.All these ground states have the same energy, and we call them degenerate ground states [39,40,[42][43][44][45][46][47]].
An additional remark is that the ground state (64) would be trivial if µ 2 < 0. However, when µ 2 > 0, the non-trivial result (64) remains.This means that whether the symmetry is spontaneously broken or not depends on the values taken by the free-parameters of the system.Further details about the mechanism of spontaneous symmetry breaking can be found in [39].In this paper, we review how to apply this principle when we deal with the financial Hamiltonians already mentioned.We start with the BS Hamiltonian because all the financial Hamiltonians are locally equivalent to the BS Hamiltonian [11].Subsequently, we demonstrate that the BS Hamiltonian and the MG one are equivalent locally after imposing local symmetry conditions over the BS equation.

The Martingale Condition as a Vacuum State
The fundamental theorem of finance establishes the existence of a martingale measure as far as the market is complete and the condition of no-arbitrage holds [4].The martingale condition is a risk-neutral measure, and it suggests that the conditional probability of the discounted value of an equity at time t = 0 is equal to its present value at t = 0 [1,[48][49][50].This can be expressed mathematically as In this equation, S(t) is the stock price, and r is the short-term risk-free interest rate.In the same equation, t is a dummy variable representing the time, integrated from zero up to a certain value as the limits of the integral suggest.In ordinary situations, the martingale condition is unique.Here, we can formulate this condition in its Hamiltonian form.For this purpose, we can consider the standard change in variable S = e x , which allows us to express the martingale condition as If we use the complete base condition ∞ −∞ dx|x >< x| = Î, with Î being the identity, and in addition S(x) =< x|S >, then we obtain This is equivalent to This is the ground state of the financial system, and it coincides with the martingale condition.The martingale condition, expressed in the form (69), can be also verified by applying the explicit form of the financial Hamiltonian over the state |S > [4].

Broken Symmetries in the Financial Equations
The symmetries of a system are defined by the operators (generators of symmetries) that commute with the Hamiltonian of the system in the following way: [ Ĥ, Â] = 0.If the same symmetry is not respected by the ground state of the system, then we have the additional condition Â|S > = 0.Under these circumstances, we say that the symmetry of the system is spontaneously broken [9,11].If we translate this language to quantum finance, this means that the martingale state does not obey the same symmetries of the financial Hamiltonian.In order to illustrate this phenomena in finance, we can start with the BS Hamiltonian [4].The Schrödinger equation is with the BS Hamiltonian defined as in Equation ( 8).Here, < x| p|C >= ∂C(x, t)/∂x.It can be proven that if we replace the option C(x, t) in Equation ( 70) by the stock price S(x, t) = e x , then we obtain the vacuum condition (69).Note, however, that although the operator defined as pC(x, t) = ∂C(x, t)/∂x is a conserved quantity (symmetry of the Hamiltonian), the same operator does not annihilate the ground state (martingale state), as can be seen from the following expression: This means that the ground state does not respect the symmetry under changes in prices, even if the BS Hamiltonian does.Then, this symmetry is spontaneously broken and the martingale state is represented by a degenerate vacuum.This is in some sense expected because this means that any value taken by the stock price S(x, t) = e x represents a potential martingale state, as far as the market is complete and there is no chance of arbitrage.Then, the market can reach the equilibrium for any value of S(x, t).Equivalently, for any ket |S >, we have potentially a martingale condition if the market is in equilibrium.If p is extended to the complex plane such that it becomes p → i p, then spontaneously broken symmetry would correspond to rotations on the complex plane.These rotations are of the form U = e −i pα .

Spontaneous Symmetry Breaking: Symmetries under Changes of Prices
We can use the result (71) to derive some expressions by defining the martingale state as < x|S >= e x = ∑ n φ n , such that < x|S >= S(x, t).The effect of the broken generator p over the field φ defined in this way is to map it toward another field φ with different components.In the same way, applying p over φ appropriately maps it back to the original field φ.Without loss of generality, we can select φ to satisfy < S| φ|S >= 0. With these definitions, we have which is obtained from and taking into account the completeness relation ∞ −∞ dx|x >< x| = Î, together with the standard expression < x| p|S >= ∂S/∂x.In addition, we can define < x| φ|S >= φ(x, t).These results are general and they do not depend on how we define the field φ.With these arguments, we can demonstrate that < S|φ|S > = 0 because < S| ∑ n φ n |S >= e x , which never vanishes unless x → −∞.Note, however, that here we have fixed φ as a quantum field satisfying < S| φ|S >= 0, which precisely corresponds to x → This is the Nambu-Goldstone field as it is defined equivalently in [39].In addition, it is well-known that the effect of p is to move the system from one vacuum to another one.Then for example, the system can be translated from the vacuum defined by φ toward a new vacuum defined by φ.Following the standard terminology, we can define the order parameter as taking into account that in general.The consistency of this expression with the results (72) and (74) suggests the validity of the condition (75).Then, what the broken symmetry generator p does is to map one vacuum toward another one.This is typical from spontaneous symmetry breaking.We can have a better understanding of this phenomena if we define the potential of the Hamiltonian (8), as follows: If we analyze the minimum for this potential, then the state C(x, t) approaches to S(x, t) = e x near the ground state (martingale state).Considering the expansion S(x, t , we can define the potential (77) as a function of S(x, t) taken as a field as Here, we have used the definition of φ n (x, t) and the fact that At the moment of comparing terms in the potential, ideally we perform it order by order.However, at any order in the expansion, the relative exponent for the terms in the expansion do not change.Then, without loss of generality, we can select any order in the expansion of n to perform the corresponding calculations.Then, we can focus on the second order terms in the expansion (n = 2), as follows: The relative exponent between the derivative and the non-derivative term is the most important point to compare.The relation between exponents does not change, no matter which order in the expansion we decide to compute.In this way, for simplicity, we can consider Equation (79) as the appropriate potential for studying the vacuum conditions of the system.The minimum of the potential can be expressed as a function of the freeparameters of the system, namely, σ and r.The potential defined in Equation ( 79) contains a term that comes from the derivative-term in Equation (77) ( pS(x, t) = ∂S(x, t)/∂x).This term is the source of non-Hermiticity, and it is responsible for the non-conservation of information.In other words, this term is related to the flow of information through the boundaries of the system [17].The minimum for the potential (79) occurs when ∂V(S) ∂φ = 0. Solving for φ(x), we find This vacuum condition clearly depends on the free-parameters of the system.The freeparameters are, naturally, the interest rate and volatility.If r = σ 2 /2, then the vacuum is trivial since φ vac = 0.This naturally corresponds to a trivial value of the security S(x, t).
The same relation guarantees the no-flow of information through the boundaries of the system (Hermiticity).Here, we take σ 2 ≤ 2r conventionally to avoid the theory being unstable.Then, the largest possible value for the volatility is constrained by the value of the interest rate r.In addition, when r >> σ 2 , the vacuum converges to a constant value.The relation (104) is non-unique because its value depends on the values taken by the free-parameters of the system.

Local Equivalence between the Black-Scholes and the Merton-Garman Equation
In several cases, it comes out that two apparently different theories are equivalent locally.The same occurs with certain groups of equations belonging to the corresponding theories.One important example happens when we go from special (SR) to general relativity (GR), for example [51].In such a case, both theories are locally equivalent, and their differences appear when we compare global aspects of both theories.Such global aspects are mathematically equivalent to the introduction of gauge fields [52].In GR, the Christoffel connections are equivalent to these gauge-fields, as can be seen from several different theories of gravity [53].The presence of gauge fields modifies the definition of derivative.Then for example, when we are working in SR, the covariant derivative for a vector V ν is defined as However, when we go to GR, the definition of derivative is modified to Globally, Equations (81) and (82) are different.However, locally they are equivalent because in such a case, the factors Γ ν µγ vanish.In GR this is called the equivalence principle and this is the reason why, in a free-falling frame of reference, an observer does not feel gravity, for example.Interestingly, we can translate similar arguments to other theories.We call these sets of theories gauge theories, and the symmetries involved are called gauge symmetries.In general terms, a gauge theory introduces physical (gauge) fields that restore the local symmetry of the Hamiltonian (Lagrangian) of the theory under analysis.The new physical (gauge) field normally appears to modify the definition of derivative inside the Hamiltonian.In this section, we show that the MG equation emerges naturally from the BS equation after imposing gauge symmetries over the BS equation [11].

The Merton-Garman Equation Emerging from the Black-Scholes Equation
We study how the BS Hamiltonian changes when a local transformation under some shift on the prices emerges.We can consider a local transformation under the changes in prices as U = e ωθ(x) .Here, θ(x) is a variable depending on x, which also depends on the price of the stock S(t), as Equation ( 46) suggests.If the operator U were a symmetry of the system, then it would satisfy the condition [ ĤBS , U] = 0.However, it is possible to demonstrate that this is not the case here after using the Hamiltonian given in Equation ( 8), together with the definition of the local changes in price.In this way, we obtain If we want to obtain an exact symmetry under local changes in the prices (U = e ωθ(x) ), then we need to add additional terms inside the Hamiltonian ( 8).This can be perceived after applying the transformation U = e ωθ(x) over the BS Hamiltonian, as follows: If we want to keep the invariance under the action of U, then we have to add some additional terms to the Hamiltonian such that the symmetry is restored.These additional terms are not added arbitrarily.The best way to add them is by modifying the definition of derivative as it was exemplified in Equation (82) for the case of general relativity.Here, for the BS equation, we redefine the derivative to become a covariant derivative: defined as Looking into the analogy with the gauge theories, we interpret py as the gauge field that restores the symmetry of the system.It comes out that this gauge field is related to the stochastic volatility of the system.More exactly, py can be interpreted as the momentum associated with the volatility.After replacing the ordinary derivative with the covariant derivative in Equation ( 8), we obtain The minus sign difference in the products appearing in the first term is a natural consequence of the fact that the momentum associated with the changes in prices, as well as the momentum associated with the changes in the stochastic volatility are both non-Hermitian quantities, satisfying then the conditions Here, the index + means Hermitian conjugate operation.Equation (86) can be expressed as The gauge invariance under a general transformation of the form U = e ωθ(x,y) for the new financial Hamiltonian defined in Equation ( 88 These conditions are obtained after checking the invariance of Equation (88).Then, the changes due to local transformations of the new terms appearing in Equation ( 88) exactly balances the additional terms appearing in Equation (84).Interestingly, when σ 2 = 2r, then ∂ 2 θ ∂x∂y = 0.This is also the same Hermiticity condition for the BS Hamiltonian, which would guarantee conservation of information.The symmetry is then completely restored after including the volatility as a gauge field inside the BS Hamiltonian, and this is independent from the free-parameters of the theory.The only consideration we need to make in order to restore the local symmetry of the Hamiltonian is to select the appropriate values for the function θ(x, y).In summary, the gauge invariance of the Hamiltonian (88), under local changes of the stock prices, is fully satisfied.Interestingly, it comes out that the Hamiltonian defined in Equation ( 88) is just the Merton-Garman Hamiltonian if we redefine the parameters appropriately, as follows: These expressions guarantee the equivalence of the Hamiltonian in Equation (88) and the MG Hamiltonian [4].Interestingly, the relations in Equation (90) give us the conditions ρ = ±1, which correspond to the extreme conditions for the parameter ρ.We have to remember that in the standard MG equation, the parameter ρ satisfies the following condition [4]: Then, as far as we define the covariant derivative in the form (85), the gauge invariance of the BS Hamiltonian suggests that the MG parameter ρ can only take the extreme values ±1 from the condition (91).If this is true, then the white noises related to the time evolution of the stock price and volatility satisfy the following conditions: under the mentioned circumstances where the gauge invariance connects the BS and the MG equations [4,10].There is an interesting connection between the interest rate r and the two volatility coefficients λ and µ inside Equation (90).This means that the local equivalence between the BS equation and the MG equation gives a natural connection between the interest rate and volatility of the market.

The Higgs Mechanism in Quantum Finance: The Dynamical Origin of the Volatility
It should not be a surprise that the stochastic volatility in the market, inside the scenario of gauge theories, behaves as a massive term inside the financial Hamiltonian.Here, we remark once again that the Hamiltonian obtained in Equation ( 88) is the Merton-Garman Hamiltonian as long as the conditions (90) are satisfied.Before, in [9,10], the ground (martingale) state for the MG Hamiltonian was analyzed.It was found to be the most general martingale condition for the MG equation, which is also consistent with Equation (88) and is given by ĤMG e x+y = ĤMG S(x, y, t) = 0. ( We can consider this as the martingale condition for the Hamiltonian (88).In [9], the martingale condition (93) is a real vacuum condition, only when the constraint, is satisfied.With the constraints defined in Equation (90), Equation (94) becomes equivalent to the condition If we solve this equation, we obtain This equation suggests an interesting relation between the stochastic volatility and the parameters λ and µ, a consequence of the martingale condition (93).This can be seen explicitly if we introduce the last relation in Equation (90) inside Equation ( 95).In such a case, we obtain e y = −r, which connects the interest rate r with the volatility (defined through y).Here, however, a negative interest rate emerges in order to have the general market equilibrium defined in Equation (93).However, we must remark that the martingale state suggested in Equation (93) is not the only possible definition of a martingale state inside the MG scenario.Still, Equation ( 93) is the only definition for martingale, considering simultaneously the prices of the stock plus the volatility.Although the martingale condition has some arbitrariness, it is still universal to define it as the ground state of the MG Hamiltonian defined in Equation (88).The potential to minimize is [9,10].
Here we consider the terms linear in px and py as potential terms.We can compare Equation (97) with the standard MG Hamiltonian defined in [9], and repeated here as It is a trivial task to show that the Equations ( 97) and ( 98) are the same under the conditions (90).
In the neighborhood of the minimal defined by the condition (93), we can consider that < x, y|S >= S(x, y, t) = e x+y = ∑ ∞ n=0 (x + y) n /n! = ∑ ∞ n=0 φ n x φ n y .Additionally, we know that ∂S(x, y, t)/∂x = ∂S(x, y, t)/∂y = . If, here again, we limit the series expansion to second order, then the potential term to analyze is The conditions defined in Equation (90) make the Equations ( 98) and (99) equivalent.Indeed, for this case, we have This only suggests that the coefficients involving φ x φ 2 y and φ 2 x φ y are the same under the conditions (90).This occurs because of the symmetric character of the covariant derivative defined in Equation (85).In principle we could have defined the covariant derivative in different ways.If, instead of the definition (85), we use introducing then another parameter γ.In this general case, the symmetric character of the terms linear in px and py in Equation ( 98) is not valid anymore, and then, the coefficients involving φ x φ 2 y and φ 2 x φ y would not be necessarily the same.Then, the introduction of γ in general modifies the possible values taken by the parameters of the MG Hamiltonian.No matter which value γ takes in Equation (101); at the most fundamental level, the MG equation emerges from the BS equation after imposing local symmetry under changes in the prices over the BS equation.Still, in the previous calculations, we took γ = 1 for the sake of simplicity.

The Dynamical Origin of the Volatility
The volatility, entering as a massive term in the Hamiltonian, has a dynamical origin through a mechanism similar to the Higgs one [11,54,55].In [9], from the general vacuum condition, emerged the relation which, under the condition γ = 1 in Equation (101), or equivalently, under the condition (100), gives the result φ xvac = φ yvac .The result (102) is a natural consequence of the general martingale state, defined in Equation (93) [9].Visualizing the dynamical origin of volatility through the general definition of martingale given in Equation ( 102) is difficult.Then, for illustrating the generic mechanism, we use a simpler definition of martingale, which is the standard one [4].In such a case, the martingale state would only depend on the prices of the stock.Considering this, then Equation (99) becomes (ignoring the derivatives with respect to volatility for the martingale state) This result ignores the term < x| py |S >= ∂S(x, t)/∂y = 0, since in this special case, we are taking S(x, t) (martingale state) as a state independent of y.From Equation (103), we can find that the minimum of the potential gives us a martingale state that only depends on the prices of the stock.This is the ordinary martingale condition that is the same for the BS and MG cases.The ground state in Equation (103) (∂V/∂φ x = 0) gives us the result [9] The expansion of the field φ x around the ground state (104) gives This is a shift on the field definition such that the ground state for φ(x) trivially vanishes.Introducing (105) inside Equation (103) gives us the potential term as a function of the shifted field φ(x), as follows: After simplification, the result (106) becomes In Equation (107), we have only focused on the massive terms for the volatility field.We have not written explicitly the other terms because we want to emphasize the existence of the massive term for the volatility field.From Equation (107), we observe that the dynamical mass for the volatility field is If φ vac vanishes (trivial result), then the massive term corresponding to the volatility field vanishes.Then the dynamical origin of the volatility mass emerges from the relation between the parameters σ and r in Equation ( 104).The dynamical mechanism for generating the mass for the volatility field, is a non-linear mechanism in this case because of the appearance of σ inside the definition of φ vac (Equation ( 104)) and because of the existence of the term e y in the mass definition in Equation (108) (remember that e y = σ 2 for the MG case).Then, the volatility field interacts with itself, generating its own mass.For analyzing the Higgs mechanism in quantum finance, we have used only second-order terms on the polynomial expansion.If we consider higher-order terms in the expansion, still the same arguments used for obtaining the result (107) work.The behavior of the kinetic terms for the MG equation, at all orders, was analyzed in [10].

Solutions for the Black-Scholes Equation
After analyzing all the details behind the BS equation, it is time to look into its solutions.We follow the solutions as they are illustrated in [2] for this purpose.We denote the prices for the European call option as C and those for the European put option as P. The solutions for the BS equation are for the call option, and for the put option.In Equations ( 109) and (110), N(x) corresponds to the standard normal distribution on the x-variable [56].In Equations ( 109) and (110), K is the strike price, T is the maturity period for the European option, S is the price of the stock, r is the interest rate and t is the time-variable.The standard normal distributions appear in Equations ( 109) and (110), depending on the variables d 1 and d 2 , here defined as [2] Then, d 1 and d 2 are variables depending explicitly on the parameters σ and r.It is clear that the volatility cannot be solved explicitly in the BS solution.The best we can do is to make an iterative process in order to find the value of volatility corresponding to the price of an option, given the values of the interest rate and the variables of the problem [2].

Implied Volatility and Making Decisions Based on the Volatility Estimation
Implied volatility is just volatility, such that when it is substituted into the Black-Scholes formula, the output is equal to the market price of the option [2,26].This is the volatility taken from the data available in Yahoo Finance and other websites [25].The investors then compare the implied value of volatility, with the value of volatility, which they estimate based on the period of time they consider.For example, they could make a daily, monthly or annual estimation of the volatility value.No matter which period of time they select for estimating the volatility, they have to be aware of the fact that the implied volatility and the estimated volatility have to be compared in the same scale of time.From these arguments, it is clear that at the moment of buying an option, the investors normally do not know the value of the volatility.Then, what the investors do is to analyze the historical behavior of the volatility; in other words, they try to see how the volatility changed in time during the previous days, months or years.With such information, the investors then can estimate the volatility, subsequently convert it to some specific scale of time, such that it can be compared with the value taken by the implied volatility.If the estimated volatility is higher than the implied volatility, then the investors should definitely buy the option.Otherwise, if the estimated value is lower than the implied value, then the investors decide not to buy the option.For closing this paper, we illustrate an example taken from [2], and here repeated for clarifying how the investors decide how whether or not they buy an option.Suppose that the value of a European call option on a non-dividend paying stock is 3.67 USD when S = 33, K = 30, r = 0.05 and T − t = 0.25 years (three months).The implied volatility would be the value of σ, such that when we introduce it inside the BS solution in Equation ( 109), together with the given values for S, K, r and T − t, it gives us as a result the value for the call option C = 3.67 USD per share.It is not difficult to realize that the process is iterative.The result for this example gives a volatility of 22% per annum.This result is taken from data in the charts, after solving σ in Equation (109).On the other hand, let us assume that the investor thinks that the volatility will be 1.5% per day.Then we have to convert this daily estimation to an annual value.The formula to apply is where D year is the number of trading days in a year.If we use this formula, then we obtain σ annual = 1.5 where we have assumed that there are 256 trading days in a year (Taking into account that the trading cannot occur everyday).From the standard theory we know that the higher the volatility of an option is, the higher is its value.Since the investor estimates a volatility of σ estimated = 24%, which is a larger value than the estimated volatility, then he will certainly invest in the option value.

Conclusions
In this paper, we have revised some fundamental properties of the BS equation.We have also expressed its Hamiltonian formulation, which facilitates the analysis of the dynamics of the option market.We then explained some fundamental mechanisms, suggesting the spontaneous symmetry breaking and Higgs mechanism inside the scenario of quantum finance, by considering symmetries under changes of the prices of stocks.If we impose the local symmetry with respect to the changes in the prices over the BS equation, then the MG equation emerges naturally as a gauge-theory.This seminal result is a natural consequence of extending the concept of gauge theories to the scenario of quantum finance.In this paper, we also analyzed the flow of information in the stock market, and finally, we explained how the investors use the solutions of the BS equation for making investment decisions by comparing the estimated volatility with the implied value of volatility.Implied volatility is just the value of volatility obtained from the BS equation after replacing all the other variables and parameters observed from the stock charts.On the other hand, estimated volatility is the value of volatility, which the investor believes should be the right one.Since the value of an option increases with the volatility; when an investor estimates a value of volatility larger than the implied value of volatility, then he/she should invest in the option.However, if the estimated value is smaller than the implied value of volatility, then the investor should decide not to invest at all in the option.In general, the option market is a zero-sum game between the holder of an option and the writer of the same option.This interesting game is modeled through the BS equation, which helps the investors to make smart decisions.We must remark that although the Black-Scholes has certain limitations, such as, for example, the fact that the volatility appears as a parameter instead of appearing as a field or variable; still, several investors use it for making investment decisions.The key power of the equation is its capability of canceling random fluctuations on the prices of the Options and its simplicity.Still, it can be complemented with the Merton-Garman equation in order to develop models for analyzing the volatility.In this paper we simply limited ourselves to mention the Merton-Garman equation, and we then explained its local equivalence with the Black-Scholes equation.This part is important because it touches on the most difficult parameter to predict from the BS equation, namely, the volatility.Finally, we must remark that the BS equation predicts the fair price of options for both call options and put options.For the call options, the holders want the prices of the same options to increase in time, because the price they agreed on for buying the option is predetermined.This is the case because the holder of a call option has the right to buy the option at the initially agreed price.For the put option, the situation is different, namely, the holder has the right to sell it at an initially determined price.Then the holder of the put option wishes the prices of the option would decrease in time, such that he/she can sell the option at a higher price than the price at maturity.It is for this reason that the Holders of a call option are called bullish and the holders of the put option are called bearish.

Figure 1 .
Figure 1.Example of the possible profits received by the holder of a European call option.The option price is 5 USD, and the strike price is 42 USD.This is an example of a zero-sum game. Figure taken from [2].

Figure 2 .
Figure 2. Example of the possible profits received by the writer of a European call option.The option price is 5 USD and the strike price is 42 USD.This is an example of a zero-sum game. Figure taken from [2].

Figure 3 .
Figure 3. Example of the possible profits received by the holder of a European put option.The option price is 3 USD and the strike price is 50 USD.This is an example of a zero-sum game. Figure taken from [2].

Figure 4 .
Figure 4. Example of the possible profits received by the writer of a European put option.The option price is 3 USD and the strike price is 50 USD.This is an example of a zero-sum game. Figure taken from [2].

Figure 5 .
Figure 5. Mexican hat potential, typical for potentials of the form defined in Equation (63). Figure taken from [41].