arXiv:cond-mat/9507128

Quantum Mechanics Simulated as Branching Process

Fricke, Thomas

Published Jul 28, 1995• 0 citations

Abstract

Diffusion processes with branching play an important role in statistical dynamics. They are a common approach to the computing of quantum mechanical groundstates, and serve as models for population dynamics and as physical pictures for biological evolution. On a computer the efficiency of this simulation method is limited by the approach to the infinitesimal time step, which is necessary to perform alternating diffusion and branching steps. In this paper, a method is described, which eliminates the infinitesimal time step for a certain class of branching processes, if the process of interest can be ``embedded'' into another process, which is solvable by other analytic and/or numerical methods. The simplest choice for the embbeding process is given by a process with a constant branching rate, which dominates the rate of the embedded process.

DOI: 10.48550/arxiv.cond-mat/9507128

Related papers

Simulating Stochastic Diffusions by Quantum Walks

Yan Wang

Stochastic differential equation (SDE) and Fokker-Planck equation (FPE) are two general approaches to describe the stochastic drift-diffusion processes. Solving SDEs relies on the Monte Carlo samplings of individual system trajectory, whereas FPEs describe the time evolution of overall distributions via path integral alike methods. The large state space and required small step size are the major challenges of computational efficiency in solving FPE numerically. In this paper, a generic continuous-time quantum walk formulation is developed to simulate stochastic diffusion processes. Stochastic diffusion in one-dimensional state space is modeled as the dynamics of an imaginary-time quantum system. The proposed quantum computational approach also drastically accelerates the path integrals with large step sizes. The new approach is compared with the traditional path integral method and the Monte Carlo trajectory sampling.

ASME International Design Engineering Technical Conferences and Computers and Information in Engineering Conference (IDETC/CIE)Aug 4, 20136 citations

Accelerating Stochastic Dynamics Simulation With Continuous-Time Quantum Walks

Yan Wang

Stochastic diffusion is a general phenomenon observed in various national and engineering systems. It is typically modeled by either stochastic differential equation (SDE) or Fokker-Planck equation (FPE), which are equivalent approaches. Path integral is an accurate and effective method to solve FPEs. Yet, computational efficiency is the common challenge for path integral and other numerical methods, include time and space complexities. Previously, one-dimensional continuous-time quantum walk was used to simulate diffusion. By combining quantum diffusion and random diffusion, the new approach can accelerate the simulation with longer time steps than those in path integral. It was demonstrated that simulation can be dozens or even hundreds of times faster. In this paper, a new generic quantum operator is proposed to simulate drift-diffusion processes in high-dimensional space, which combines quantum walks on graphs with traditional path integral approaches. Probability amplitudes are computed efficiently by spectral analysis. The efficiency of the new method is demonstrated with stochastic resonance problems.

ASME International Design Engineering Technical Conferences and Computers and Information in Engineering Conference (IDETC/CIE)Aug 21, 20165 citations

Quantum Computing Through Multi-Scale Modeling and Simulation

Jincheng Zhang

This paper presents a novel approach to quantum computing leveraging multi-scale modeling and simulation. The core claim is to develop a method capable of simultaneously simulating quantum systems at both microscopic and macroscopic levels, thereby enhancing the precision and efficiency of quantum computations. This is achieved through the construction of a multi-scale model integrating quantum mechanics and classical mechanics, coupled with advanced numerical simulation techniques. The proposed framework addresses the limitations of current quantum simulation methods by offering a more holistic representation of quantum systems. This allows for a deeper understanding of quantum phenomena and ultimately contributes to the advancement of quantum computing technologies. The key mathematical components utilized within this framework include the Schrödinger equation, perturbation theory, and various numerical integration schemes. Specifically, the time evolution of a quantum state is governed by the equation: (iħ ∂/∂t) |ψ(t)> = H |ψ(t)> where |ψ(t)> is the quantum state at time *t*, *ħ* is the reduced Planck constant, and *H* is the Hamiltonian operator. Perturbation theory, used to approximate the solution when the Hamiltonian is not solvable exactly, relies on the expansion: |ψperturbed>(t) = |ψunperturbed>(t) + ∑n=1∞ (-1)n+1 *An* |ψn>(t) where *An* are the perturbation coefficients. Numerical integration techniques, such as the Runge-Kutta method, are employed to approximate the time evolution operator. The accuracy of the simulation is critically dependent on the chosen time step, Δt, and the order of the integration method. A smaller Δt generally leads to higher accuracy but requires more computational resources. The overall approach seeks to bridge the gap between theoretical quantum mechanics and practical quantum computation, offering a viable path towards building more sophisticated and reliable quantum computers.

Sep 5, 20260 citations

Optimization of the Split-Operator Method for Modeling of Quantum System Time Evolution

D Vlasov

A continuing and important problem in modeling physical processes deals with creating optimum calculation algorithms. These algorithms should provide the necessary precision yet not demand a lot of storage and time for evaluations. In the present paper the numerical integration of the non-stationary Schrödinger equation (NSSE) is considered. There are many numerical methods of integration of this equation to the present moment. In this paper, we explore the split-operator method (SOM) and consider a way of lessening time computing expenses in view of specificity of the given method. We also consider the choice of optimum discretization steps in coordinate and momentum space. The computer application that we outline is applicable, in method and preparation, to the solution more difficult problems, and it has been applied to simple chemical reactions. We use standard atomic units to simplify equations in the analysis that follows.

American Journal of Undergraduate ResearchSep 17, 20050 citations

Stationary quantum Markov process for the Wigner function

Abstract. Many stochastic models have been investigated for quantum mechanics because of its stochastic nature. In 1988, Cohendet et al. introduced a dichotomic variable to quantum phase space and proposed a background Markov process for the time evolution of the Wigner function. However, in their method the whole distribution function is required to determine the next step of a constituent particle. In this paper, we discuss a stationary quantum Markov process which enables us to treat each particle independently introducing U(1) extension for the phase space. The process has branching and vanishing points but we can keep a finite time interval between the branchings. The procedure to make the simulation of the process is also discussed. Stationary quantum Markov process for the Wigner function 2 1.

Jan 1, 20070 citations

A computable branching process for the Wigner quantum dynamics

Sihong Shao, Yunfeng Xiong

A branching process treatment for the nonlocal Wigner pseudo-differential operator and its numerical applications in quantum dynamics is proposed and analyzed. We start from the discussion on two typical truncations of the nonlocal term, i.e., the k-truncated and y-truncated models. After introducing an auxiliary function γ(x), the (truncated) Wigner equation is reformulated into the integral formulation as well as its adjoint correspondence, both of which can be regarded as the renewal-type equations and have transparent stochastic interpretation. We prove that the moment of a branching process happens to be the solution for the adjoint equation, which connects rigorously the Wigner quantum dynamics to the stochastic branching process, and thus a sound mathematical framework for the Wigner Monte Carlo methods is established. Within the framework, the branching process for the y-truncated model recovers the popular signed particle Monte Carlo method which needs a discretization of the momentum space beforehand, whereas the k-truncated branching method allows a continuous momentum sampling. More interestingly, with the help of the rigorous connection, we find that a constant auxiliary function γ(x) ≡ γ0 performs better from the point of view of both theoretical and numerical aspects. In theory, we obtain a simple but analytical calculation formula for the growth rate of particle number, i.e., e2γ0t with t being the time. In numerics, the constant auxiliary function may facilitate the random sampling from the exponential distribution and improve the accuracy. Typical numerical experiments on the Gaussian barrier scattering validate our theoretical findings, as well as demonstrate the accuracy, the efficiency and thus the computability of the Wigner branching process.

Mar 1, 20160 citations