Nonlinear geometric optics based multiscale stochastic Galerkin methods for highly oscillatory transport equations with random inputs

We develop generalized polynomial chaos (gPC) based stochastic Galerkin (SG) methods for a class of highly oscillatory transport equations that arise in semiclassical modeling of non-adiabatic quantum dynamics. These models contain uncertainties, particularly in coefficients that correspond to the potentials of the molecular system. We first focus on a highly oscillatory scalar model with random uncertainty. Our method is built upon the nonlinear geometrical optics (NGO) based method, developed in \cite{NGO} for numerical approximations of deterministic equations, which can obtain accurate pointwise solution even without numerically resolving spatially and temporally the oscillations. With the random uncertainty, we show that such a method has oscillatory higher order derivatives in the random space, thus requires a frequency dependent discretization in the random space. We modify this method by introducing a new "time" variable based on the phase, which is shown to be non-oscillatory in the random space, based on which we develop a gPC-SG method that can capture oscillations with the frequency-independent time step, mesh size as well as the degree of polynomial chaos. A similar approach is then extended to a semiclassical surface hopping model system with a similar numerical conclusion. Various numerical examples attest that these methods indeed capture accurately the solution statistics {\em pointwisely} even though none of the numerical parameters resolve the high frequencies of the solution.

[1]  Branko Ruscic,et al.  Uncertainty quantification in thermochemistry, benchmarking electronic structure computations, and Active Thermochemical Tables , 2014 .

[2]  Nicolas Crouseilles,et al.  Uniformly accurate numerical schemes for highly oscillatory Klein–Gordon and nonlinear Schrödinger equations , 2013, Numerische Mathematik.

[3]  Jian-Guo Liu,et al.  The Reconstruction of Upwind Fluxes for Conservation Laws , 1998 .

[4]  O. Morandi,et al.  Effective classical Liouville-like evolution equation for the quantum phase-space dynamics , 2010 .

[5]  Shi Jin,et al.  Semi-classical models for the Schrödinger equation with periodic potentials and band crossings , 2013 .

[6]  Jie Shen,et al.  Spectral Methods: Algorithms, Analysis and Applications , 2011 .

[7]  Dongbin Xiu,et al.  Galerkin method for wave equations with uncertain coefficients , 2008 .

[8]  Dongbin Xiu,et al.  The Wiener-Askey Polynomial Chaos for Stochastic Differential Equations , 2002, SIAM J. Sci. Comput..

[9]  Markus Reiher,et al.  Uncertainty quantification for quantum chemical models of complex reaction networks. , 2016, Faraday discussions.

[10]  R. Ghanem,et al.  Stochastic Finite Elements: A Spectral Approach , 1990 .

[11]  Omar M. Knio,et al.  Spectral Methods for Uncertainty Quantification , 2010 .

[12]  Nicolas Crouseilles,et al.  Asymptotic Preserving schemes for highly oscillatory Vlasov-Poisson equations , 2012, J. Comput. Phys..

[13]  J. Tully,et al.  Trajectory Surface Hopping Approach to Nonadiabatic Molecular Collisions: The Reaction of H+ with D2 , 1971 .

[14]  J. Tully.  Molecular dynamics with electronic transitions , 1990 .

[15]  R. Parr.  Density-functional theory of atoms and molecules , 1989 .

[16]  Michael I. Weinstein,et al.  Honeycomb Lattice Potentials and Dirac Points , 2012, 1202.3839.

[17]  Tanja van Mourik,et al.  A critical note on density functional theory studies on rare-gas dimers , 2002 .

[18]  Nicolas Crouseilles,et al.  Nonlinear Geometric Optics method based multi-scale numerical schemes for highly-oscillatory transport equations , 2016, 1605.09676.

[19]  Qin Li,et al.  A Multiband Semiclassical Model for Surface Hopping Quantum Dynamics , 2014, Multiscale Model. Simul..

[20]  Guannan Zhang,et al.  Stochastic finite element methods for partial differential equations with random input data* , 2014, Acta Numerica.

[21]  O. Morandi,et al.  Wigner model for Klein tunneling in graphene , 2011 .

[22]  F. Guinea,et al.  The electronic properties of graphene , 2007, Reviews of Modern Physics.

[23]  O. Morandi,et al.  Wigner model for quantum transport in graphene , 2011, 1102.2416.

[24]  W. Heisenberg,et al.  Zur Quantentheorie der Molekeln , 1924 .

[25]  Dongbin Xiu,et al.  High-Order Collocation Methods for Differential Equations with Random Inputs , 2005, SIAM J. Sci. Comput..

[26]  Biao Wu,et al.  Nonlinear Landau-Zener tunneling , 2000 .

[27]  Pavel Hobza,et al.  Unexpectedly strong energy stabilization inside the hydrophobic core of small protein rubredoxin mediated by aromatic residues: correlated ab initio quantum chemical calculations. , 2005 .