The Wigner Branching Random Walk: Efficient Implementation and Performance Evaluation

The Wigner Branching Random Walk: Efficient Implementation and Performance Evaluation

Year:    2019

Communications in Computational Physics, Vol. 25 (2019), Iss. 3 : pp. 871–910

Abstract

To implement the Wigner branching random walk, the particle carrying a signed weight, either −1 or +1, is more friendly to data storage and arithmetic manipulations than that taking a real-valued weight continuously from −1 to +1. The former is called a signed particle and the latter a weighted particle. In this paper, we propose two efficient strategies to realize the signed-particle implementation. One is to interpret the multiplicative functional as the probability to generate pairs of particles instead of the incremental weight, and the other is to utilize a bootstrap filter to adjust the skewness of particle weights. Performance evaluations on the Gaussian barrier scattering (2D) and a Helium-like system (4D) demonstrate the feasibility of both strategies and the variance reduction property of the second approach. We provide an improvement of the first signed-particle implementation that partially alleviates the restriction on the time step and perform a thorough theoretical and numerical comparison among all the existing signed-particle implementations. Details on implementing the importance sampling according to the quasi-probability density and an efficient resampling or particle reduction are also provided.

You do not have full access to this article.

Already a Subscriber? Sign in as an individual or via your institution

Journal Article Details

Publisher Name:    Global Science Press

Language:    English

DOI:    https://doi.org/10.4208/cicp.OA-2018-0141

Communications in Computational Physics, Vol. 25 (2019), Iss. 3 : pp. 871–910

Published online:    2019-01

AMS Subject Headings:    Global Science Press

Copyright:    COPYRIGHT: © Global Science Press

Pages:    40

Keywords:    Wigner equation branching random walk signed particle bootstrapping weighted particle Monte Carlo method quantum dynamics importance sampling resampling particle reduction.

  1. Solving the Wigner equation with signed particle Monte Carlo for chemically relevant potentials

    Wang, Yu | Simine, Lena

    The Journal of Chemical Physics, Vol. 155 (2021), Iss. 3

    https://doi.org/10.1063/5.0055603 [Citations: 2]
  2. A higher-order accurate operator splitting spectral method for the Wigner–Poisson system

    Chen, Zhenzhu | Jiang, Haiyan | Shao, Sihong

    Journal of Computational Electronics, Vol. 21 (2022), Iss. 4 P.756

    https://doi.org/10.1007/s10825-022-01904-x [Citations: 3]
  3. Branching Random Walk Solutions to the Wigner Equation

    Shao, Sihong | Xiong, Yunfeng

    SIAM Journal on Numerical Analysis, Vol. 58 (2020), Iss. 5 P.2589

    https://doi.org/10.1137/19M1272408 [Citations: 2]