Efficient sampling for the conditional time integrated variance process in the Heston stochastic volatility model is key to the simulation of the stock price based on its exact distribution. We construct a new series expansion for this integral in te
rms of double infinite weighted sums of particular independent random variables through a change of measure and the decomposition of squared Bessel bridges. When approximated by series truncations, this representation has exponentially decaying truncation errors. We propose feasible strategies to largely reduce the implementation of the new series to simulations of simple random variables that are independent of any model parameters. We further develop direct inversion algorithms to generate samples for such random variables based on Chebyshev polynomial approximations for their inverse distribution functions. These approximations can be used under any market conditions. Thus, we establish a strong, efficient and almost exact sampling scheme for the Heston model.
The time evolution problem for non-self adjoint second order differential operators is studied by means of the path integral formulation. Explicit computation of the path integral via the use of certain underlying stochastic differential equations, w
hich naturally emerge when computing the path integral, leads to a universal expression for the associated measure regardless of the form of the differential operators. The discrete non-linear hierarchy (DNLS) is then considered and the corresponding hierarchy of solvable, in principle, SDEs is extracted. The first couple members of the hierarchy correspond to the discrete stochastic transport and heat equations. The discrete stochastic Burgers equation is also obtained through the analogue of the Cole-Hopf transformation. The continuum limit is also discussed.
We explore the connections between the theories of stochastic analysis and discrete quantum mechanical systems. Naturally these connections include the Feynman-Kac formula, and the Cameron-Martin-Girsanov theorem. More precisely, the notion of the qu
antum canonical transformation is employed for computing the time propagator, in the case of generic dynamical diffusion coefficients. Explicit computation of the path integral leads to a universal expression for the associated measure regardless of the form of the diffusion coefficient and the drift. This computation also reveals that the drift plays the role of a super potential in the usual super-symmetric quantum mechanics sense. Some simple illustrative examples such as the Ornstein-Uhlenbeck process and the multidimensional Black-Scholes model are also discussed. Basic examples of quantum integrable systems such as the quantum discrete non-linear hierarchy (DNLS) and the XXZ spin chain are presented providing specific connections between quantum (integrable) systems and stochastic differential equations (SDEs). The continuum limits of the SDEs for the first two members of the NLS hierarchy turn out to be the stochastic transport and the stochastic heat equations respectively. The quantum Darboux matrix for the discrete NLS is also computed as a defect matrix and the relevant SDEs are derived.
The transition probability of a Cox-Ingersoll-Ross process can be represented by a non-central chi-square density. First we prove a new representation for the central chi-square density based on sums of powers of generalized Gaussian random variables
. Second we prove Marsaglias polar method extends to this distribution, providing a simple, exact, robust and efficient acceptance-rejection method for generalized Gaussian sampling and thus central chi-square sampling. Third we derive a simple, high-accuracy, robust and efficient direct inversion method for generalized Gaussian sampling based on the Beasley-Springer-Moro method. Indeed the accuracy of the approximation to the inverse cumulative distribution function is to the tenth decimal place. We then apply our methods to non-central chi-square variance sampling in the Heston model. We focus on the case when the number of degrees of freedom is small and the zero boundary is attracting and attainable, typical in foreign exchange markets. Using the additivity property of the chi-square distribution, our methods apply in all parameter regimes.
We present a new method for sampling the Levy area for a two-dimensional Wiener process conditioned on its endpoints. An efficient sampler for the Levy area is required to implement a strong Milstein numerical scheme to approximate the solution of a
stochastic differential equation driven by a two-dimensional Wiener process whose diffusion vector fields do not commute. Our method is simple and complementary to those of Gaines-Lyons and Wiktorsson, and amenable to quasi-Monte--Carlo implementation. It is based on representing the Levy area by an infinite weighted sum of independent Logistic random variables. We use Chebychev polynomials to approximate the inverse distribution function of sums of independent Logistic random variables in three characteristic regimes. The error is controlled by the degree of the polynomials, we set the error to be uniformly 10^(-12). We thus establish a strong almost-exact Levy area sampling method. The complexity of our method is square logarithmic. We indicate how our method can contribute to efficient sampling in higher dimensions.
We study solutions to nonlinear stochastic differential systems driven by a multi-dimensional Wiener process. A useful algorithm for strongly simulating such stochastic systems is the Castell--Gaines method, which is based on the exponential Lie seri
es. When the diffusion vector fields commute, it has been proved that at low orders this method is more accurate in the mean-square error than corresponding stochastic Taylor methods. However it has also been shown that when the diffusion vector fields do not commute, this is not true for strong order one methods. Here we prove that when there is no drift, and the diffusion vector fields do not commute, the exponential Lie series is usurped by the sinh-log series. In other words, the mean-square error associated with a numerical method based on the sinh-log series, is always smaller than the corresponding stochastic Taylor error, in fact to all orders. Our proof utilizes the underlying Hopf algebra structure of these series, and a two-alphabet associative algebra of shuffle and concatenation operations. We illustrate the benefits of the proposed series in numerical studies.
We present Lie group integrators for nonlinear stochastic differential equations with non-commutative vector fields whose solution evolves on a smooth finite dimensional manifold. Given a Lie group action that generates transport along the manifold,
we pull back the stochastic flow on the manifold to the Lie group via the action, and subsequently pull back the flow to the corresponding Lie algebra via the exponential map. We construct an approximation to the stochastic flow in the Lie algebra via closed operations and then push back to the Lie group and then to the manifold, thus ensuring our approximation lies in the manifold. We call such schemes stochastic Munthe-Kaas methods after their deterministic counterparts. We also present stochastic Lie group integration schemes based on Castell--Gaines methods. These involve using an underlying ordinary differential integrator to approximate the flow generated by a truncated stochastic exponential Lie series. They become stochastic Lie group integrator schemes if we use Munthe-Kaas methods as the underlying ordinary differential integrator. Further, we show that some Castell--Gaines methods are uniformly more accurate than the corresponding stochastic Taylor schemes. Lastly we demonstrate our methods by simulating the dynamics of a free rigid body such as a satellite and an autonomous underwater vehicle both perturbed by two independent multiplicative stochastic noise processes.
