arrow
Volume 12, Issue 1
A Branching Random Walk Method for Many-Body Wigner Quantum Dynamics

Sihong Shao & Yunfeng Xiong

Numer. Math. Theor. Meth. Appl., 12 (2019), pp. 21-71.

Published online: 2018-09

Export citation
  • Abstract

A branching random walk algorithm for many-body Wigner equations and its numerical applications for quantum dynamics in phase space are proposed and analyzed in this paper. Using an auxiliary function, the truncated Wigner equation and its adjoint form are cast into integral formulations, which can be then reformulated into renewal-type equations with probabilistic interpretations. We prove that the first moment of a branching random walk is the solution for the adjoint equation. With the help of the additional degree of freedom offered by the auxiliary function, we are able to produce a weighted-particle implementation of the branching random walk. In contrast to existing signed-particle implementations, this weighted-particle one shows a key capacity of variance reduction by increasing the constant auxiliary function and has no time discretization errors. Several canonical numerical experiments on the 2D Gaussian barrier scattering and a 4D Helium-like system validate our theoretical findings, and demonstrate the accuracy, the efficiency, and thus the computability of the proposed weighted-particle Wigner branching random walk algorithm.

  • Keywords

  • AMS Subject Headings

  • Copyright

COPYRIGHT: © Global Science Press

  • Email address
  • BibTex
  • RIS
  • TXT
@Article{NMTMA-12-21, author = {Sihong Shao and Yunfeng Xiong}, title = {A Branching Random Walk Method for Many-Body Wigner Quantum Dynamics}, journal = {Numerical Mathematics: Theory, Methods and Applications}, year = {2018}, volume = {12}, number = {1}, pages = {21--71}, abstract = {

A branching random walk algorithm for many-body Wigner equations and its numerical applications for quantum dynamics in phase space are proposed and analyzed in this paper. Using an auxiliary function, the truncated Wigner equation and its adjoint form are cast into integral formulations, which can be then reformulated into renewal-type equations with probabilistic interpretations. We prove that the first moment of a branching random walk is the solution for the adjoint equation. With the help of the additional degree of freedom offered by the auxiliary function, we are able to produce a weighted-particle implementation of the branching random walk. In contrast to existing signed-particle implementations, this weighted-particle one shows a key capacity of variance reduction by increasing the constant auxiliary function and has no time discretization errors. Several canonical numerical experiments on the 2D Gaussian barrier scattering and a 4D Helium-like system validate our theoretical findings, and demonstrate the accuracy, the efficiency, and thus the computability of the proposed weighted-particle Wigner branching random walk algorithm.

}, issn = {2079-7338}, doi = {https://doi.org/10.4208/nmtma.OA-2018-0074}, url = {http://global-sci.org/intro/article_detail/nmtma/12690.html} }
TY - JOUR T1 - A Branching Random Walk Method for Many-Body Wigner Quantum Dynamics AU - Sihong Shao & Yunfeng Xiong JO - Numerical Mathematics: Theory, Methods and Applications VL - 1 SP - 21 EP - 71 PY - 2018 DA - 2018/09 SN - 12 DO - http://doi.org/10.4208/nmtma.OA-2018-0074 UR - https://global-sci.org/intro/article_detail/nmtma/12690.html KW - AB -

A branching random walk algorithm for many-body Wigner equations and its numerical applications for quantum dynamics in phase space are proposed and analyzed in this paper. Using an auxiliary function, the truncated Wigner equation and its adjoint form are cast into integral formulations, which can be then reformulated into renewal-type equations with probabilistic interpretations. We prove that the first moment of a branching random walk is the solution for the adjoint equation. With the help of the additional degree of freedom offered by the auxiliary function, we are able to produce a weighted-particle implementation of the branching random walk. In contrast to existing signed-particle implementations, this weighted-particle one shows a key capacity of variance reduction by increasing the constant auxiliary function and has no time discretization errors. Several canonical numerical experiments on the 2D Gaussian barrier scattering and a 4D Helium-like system validate our theoretical findings, and demonstrate the accuracy, the efficiency, and thus the computability of the proposed weighted-particle Wigner branching random walk algorithm.

Sihong Shao and Yunfeng Xiong. (2018). A Branching Random Walk Method for Many-Body Wigner Quantum Dynamics. Numerical Mathematics: Theory, Methods and Applications. 12 (1). 21-71. doi:10.4208/nmtma.OA-2018-0074
Copy to clipboard
The citation has been copied to your clipboard