TY - GEN A1 - Lasser, Caroline A1 - Röblitz, Susanna T1 - Computing expectation values for molecular quantum dynamics N2 - We compute expectation values for the solution of the nuclear Schrödinger equation. The proposed particle method consists of three steps: sampling of the initial Wigner function, classical transport of the sampling points, weighted phase space summation for the final computation of the expectation values. The Egorov theorem guarantees that the algorithm is second order accurate with respect to the semiclassical parameter. We present numerical experiments for a two-dimensional torsional potential with three different sets of initial data and for a six-dimensional Henon-Heiles potential. By construction, the computing times scale linearly with the number of initial sampling points and range between three seconds and one hour. T3 - ZIB-Report - 09-30 KW - time dependent Schrödinger equation KW - molecular quantum dynamics KW - expectation values KW - Wigner functions Y1 - 2009 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-11461 SN - 1438-0064 ER - TY - GEN A1 - Kube, Susanna A1 - Lasser, Caroline A1 - Weber, Marcus T1 - Monte Carlo sampling of Wigner functions and surface hopping quantum dynamics N2 - Wigner transformation provides a one-to-one correspondence between functions on position space (wave functions) and functions on phase space (Wigner functions). Weighted integrals of Wigner functions yield quadratic quantities of wave functions like position and momentum densities or expectation values. For molecular quantum systems, suitably modified classical transport of Wigner functions provides an asymptotic approximation of the dynamics in the high energy regime. The article addresses the computation of Wigner functions by Monte Carlo quadrature. An ad aption of the Metropolis algorithm for the approximation of signed measures with disconnected support is systematically tested in combination with a surface hopping algorithm for non-adiabatic quantum dynamics. The numerical experiments give expectation values and level populations with an error of two to three percent, which agrees with the theoretically expected accuracy. T3 - ZIB-Report - 07-17 KW - Metropolis Monte Carlo KW - approximation KW - quadrature KW - oscillating functions Y1 - 2007 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-9604 ER - TY - JOUR A1 - Lasser, Caroline A1 - Röblitz, Susanna T1 - Computing expectation values for molecular quantum dynamics JF - SIAM Journal on Scientific Computing Y1 - 2010 U6 - https://doi.org//10.1137/090770461 VL - 32 IS - 3 SP - 1465 EP - 1483 ER - TY - JOUR A1 - Kube, Susanna A1 - Lasser, Caroline A1 - Weber, Marcus T1 - Monte Carlo sampling of Wigner functions and surface hopping quantum dynamics JF - Journal of Computational Physics Y1 - 2008 U6 - https://doi.org/10.1016/j.jcp.2008.11.016 VL - 228 IS - 6 SP - 1947 EP - 1962 ER - TY - JOUR A1 - Kenfack, Anatole A1 - Marquardt, Falko A1 - K. Paramonov, Gennadij A1 - Barth, Ingo A1 - Lasser, Caroline A1 - Paulus, Beate T1 - Initial-state dependence of coupled electronic and nuclear fluxes in molecules JF - Phys. Rev. A Y1 - 2010 U6 - https://doi.org/10.1103/PhysRevA.81.052502 VL - 81 IS - 5 SP - 052502 PB - American Physical Society ER -