Related papers: Fourth Order Algorithms for Solving the Multivaria…
Langevin simulation provides an effective way to study collisional effects in beams by reducing the six-dimensional Fokker-Planck equation to a group of stochastic ordinary differential equations. These resulting equations usually have…
We show that the method of factorizing the evolution operator to fourth order with purely positive coefficients, in conjunction with Suzuki's method of implementing time-ordering of operators, produces a new class of powerful algorithms for…
Formulated is a new systematic method for obtaining higher order corrections in numerical simulation of stochastic differential equations (SDEs), i.e., Langevin equations. Random walk step algorithms within a given order of finite $\Delta…
We show that the method of splitting the operator ${\rm e}^{\epsilon(T+V)}$ to fourth order with purely positive coefficients produces excellent algorithms for solving the time-dependent Schr\"odinger equation. These algorithms require…
A consequent approach is proposed to construct symplectic force-gradient algorithms of arbitrarily high orders in the time step for precise integration of motion in classical and quantum mechanics simulations. Within this approach the basic…
Starting from a forward--backward path integral of a point particle in a bath of harmonic oscillators, we derive the Fokker-Planck and Langevin equations with and without inertia. Special emphasis is placed upon the correct operator order…
We demonstrate a method which allows the stochastic modelling of quantum systems for which the generalised Fokker-Planck equation in the phase space contains derivatives of higher than second order. This generalises quantum stochastics far…
The Fokker-Planck equation models rare events across sciences, but its high-dimensional nature challenges classical computers. Quantum algorithms for such non-unitary dynamics often suffer from exponential {decay in} success probability. We…
In this paper, we introduce second order and fourth order space discretization via finite difference implementation of the finite element method for solving Fokker-Planck equations associated with irreversible processes. The proposed…
We propose a discrete lattice version of the Fokker-Planck kinetic equation along lines similar to the Lattice-Boltzmann scheme. Our work extends an earlier one-dimensional formulation to arbitrary spatial dimension $D$. A generalized…
We obtain exact results for fractional equations of Fokker-Planck type using evolution operator method. We employ exact forms of one-sided Levy stable distributions to generate a set of self-reproducing solutions. Explicit cases are…
In this work, we present a second-order numerical scheme to address the solution of optimal control problems constrained by the evolution of nonlinear Fokker-Planck equations arising from socio-economic dynamics. In order to design an…
Based on the continuous time random walk, we derive the Fokker-Planck equations with Caputo-Fabrizio fractional derivative, which can effectively model a variety of physical phenomena, especially, the material heterogeneities and structures…
Recently, the numerical schemes of the Fokker-Planck equations describing anomalous diffusion with two internal states have been proposed in [Nie, Sun and Deng, arXiv: 1811.04723], which use convolution quadrature to approximate the…
The solution of many physical evolution equations can be expressed as an exponential of two or more operators acting on initial data. Accurate solutions can be systematically derived by decomposing the exponential in a product form. For…
Our paper "Solving Third Order Linear Difference Equations in Terms of Second Order Equations" gave two algorithms for solving difference equations in terms of lower order equations: an algorithm for absolute factorization, and an algorithm…
We present a method for constructing numerical schemes with up to 3rd strong convergence order for solution of a class of stochastic differential equations, including equations of the Langevin type. The construction proceeds in two stages.…
By decomposing the important sampled imaginary time Schr\"odinger evolution operator to fourth order with positive coefficients, we derived a number of distinct fourth order Diffusion Monte Carlo algorithms. These sophisticated algorithms…
By implementing the exact density matrix for the rotating anisotropic harmonic trap, we derive a class of very fast and accurate fourth order algorithms for evolving the Gross-Pitaevskii equation in imaginary time. Such fourth order…
We present an algorithm which allows to solve analytically linear systems of differential equations which factorize to first order. The solution is given in terms of iterated integrals over an alphabet where its structure is implied by the…