- Journal Home
- Volume 42 - 2024
- Volume 41 - 2023
- Volume 40 - 2022
- Volume 39 - 2021
- Volume 38 - 2020
- Volume 37 - 2019
- Volume 36 - 2018
- Volume 35 - 2017
- Volume 34 - 2016
- Volume 33 - 2015
- Volume 32 - 2014
- Volume 31 - 2013
- Volume 30 - 2012
- Volume 29 - 2011
- Volume 28 - 2010
- Volume 27 - 2009
- Volume 26 - 2008
- Volume 25 - 2007
- Volume 24 - 2006
- Volume 23 - 2005
- Volume 22 - 2004
- Volume 21 - 2003
- Volume 20 - 2002
- Volume 19 - 2001
- Volume 18 - 2000
- Volume 17 - 1999
- Volume 16 - 1998
- Volume 15 - 1997
- Volume 14 - 1996
- Volume 13 - 1995
- Volume 12 - 1994
- Volume 11 - 1993
- Volume 10 - 1992
- Volume 9 - 1991
- Volume 8 - 1990
- Volume 7 - 1989
- Volume 6 - 1988
- Volume 5 - 1987
- Volume 4 - 1986
- Volume 3 - 1985
- Volume 2 - 1984
- Volume 1 - 1983
Cited by
- BibTex
- RIS
- TXT
This paper discusses a numerical method for computing the evolution of large interacting system of quantum particles. The idea of the random batch method is to replace the total interaction of each particle with the $N-1$ other particles by the interaction with $p\ll N$ particles chosen at random at each time step, multiplied by $(N-1)/p$. This reduces the computational cost of computing the interaction potential per time step from $O(N^2)$ to $O(N)$. For simplicity, we consider only in this work the case $p=1$ — in other words, we assume that $N$ is even, and that at each time step, the $N$ particles are organized in $N/2$ pairs, with a random reshuffling of the pairs at the beginning of each time step. We obtain a convergence estimate for the Wigner transform of the single-particle reduced density matrix of the particle system at time $t$ that is both uniform in $N>1$ and independent of the Planck constant $\hbar$. The key idea is to use a new type of distance on the set of quantum states that is reminiscent of the Wasserstein distance of exponent $1$ (or Monge-Kantorovich-Rubinstein distance) on the set of Borel probability measures on $\mathbf{R}^d$ used in the context of optimal transport.
}, issn = {1991-7139}, doi = {https://doi.org/10.4208/jcm.2107-m2020-0306}, url = {http://global-sci.org/intro/article_detail/jcm/19946.html} }This paper discusses a numerical method for computing the evolution of large interacting system of quantum particles. The idea of the random batch method is to replace the total interaction of each particle with the $N-1$ other particles by the interaction with $p\ll N$ particles chosen at random at each time step, multiplied by $(N-1)/p$. This reduces the computational cost of computing the interaction potential per time step from $O(N^2)$ to $O(N)$. For simplicity, we consider only in this work the case $p=1$ — in other words, we assume that $N$ is even, and that at each time step, the $N$ particles are organized in $N/2$ pairs, with a random reshuffling of the pairs at the beginning of each time step. We obtain a convergence estimate for the Wigner transform of the single-particle reduced density matrix of the particle system at time $t$ that is both uniform in $N>1$ and independent of the Planck constant $\hbar$. The key idea is to use a new type of distance on the set of quantum states that is reminiscent of the Wasserstein distance of exponent $1$ (or Monge-Kantorovich-Rubinstein distance) on the set of Borel probability measures on $\mathbf{R}^d$ used in the context of optimal transport.