TY - GEN A1 - Gelß, Patrick A1 - Matera, Sebastian A1 - Schütte, Christof T1 - Solving the master equation without kinetic Monte Carlo: tensor train approximations for a CO oxidation model N2 - In multiscale models of heterogeneous catalysis, one crucial point is the solution of a Markovian master equation describing the stochastic reaction kinetics. This usually is too high-dimensional to be solved with standard numerical techniques and one has to rely on sampling approaches based on the kinetic Monte Carlo method. In this study we break the curse of dimensionality for the direct solution of the Markovian master equation by exploiting the Tensor Train Format for this purpose. The performance of the approach is demonstrated on a first principles based, reduced model for the CO oxidation on the RuO_2(110) surface. We investigate the complexity for increasing system size and for various reaction conditions. The advantage over the stochastic simulation approach is illustrated by a problem with increased stiffness. T3 - ZIB-Report - 15-41 KW - heterogeneous catalysis KW - master equation KW - kinetic Monte Carlo KW - tensor decompositions KW - tensor train format KW - alternating linear scheme Y1 - 2015 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-55743 SN - 1438-0064 ER - TY - JOUR A1 - Gelß, Patrick A1 - Matera, Sebastian A1 - Schütte, Christof T1 - Solving the master equation without kinetic Monte Carlo: Tensor train approximations for a CO oxidation model JF - Journal of Computational Physics N2 - In multiscale modeling of heterogeneous catalytic processes, one crucial point is the solution of a Markovian master equation describing the stochastic reaction kinetics. Usually, this is too high-dimensional to be solved with standard numerical techniques and one has to rely on sampling approaches based on the kinetic Monte Carlo method. In this study we break the curse of dimensionality for the direct solution of the Markovian master equation by exploiting the Tensor Train Format for this purpose. The performance of the approach is demonstrated on a first principles based, reduced model for the CO oxidation on the RuO2(110) surface. We investigate the complexity for increasing system size and for various reaction conditions. The advantage over the stochastic simulation approach is illustrated by a problem with increased Y1 - 2016 U6 - https://doi.org/10.1016/j.jcp.2016.03.025 VL - 314 SP - 489 EP - 502 ER - TY - JOUR A1 - Gelß, Patrick A1 - Klein, Rupert A1 - Matera, Sebastian A1 - Schmidt, Burkhard T1 - Quantum dynamics of coupled excitons and phonons in chain-like systems: tensor train approaches and higher-order propagators JF - The Journal of Chemical Physics Y1 - 2025 U6 - https://doi.org/10.1063/5.0258904 VL - 162 ER - TY - JOUR A1 - Gelß, Patrick A1 - Klein, Rupert A1 - Matera, Sebastian A1 - Schmidt, Burkhard T1 - Solving the time-independent Schrödinger equation for chains of coupled excitons and phonons using tensor trains JF - The Journal of Chemical Physics N2 - We demonstrate how to apply the tensor-train format to solve the time-independent Schrödinger equation for quasi-one-dimensional excitonic chain systems with and without periodic boundary conditions. The coupled excitons and phonons are modeled by Fröhlich–Holstein type Hamiltonians with on-site and nearest-neighbor interactions only. We reduce the memory consumption as well as the computational costs significantly by employing efficient decompositions to construct low-rank tensor-train representations, thus mitigating the curse of dimensionality. In order to compute also higher quantum states, we introduce an approach that directly incorporates the Wielandt deflation technique into the alternating linear scheme for the solution of eigenproblems. Besides systems with coupled excitons and phonons, we also investigate uncoupled problems for which (semi-)analytical results exist. There, we find that in the case of homogeneous systems, the tensor-train ranks of state vectors only marginally depend on the chain length, which results in a linear growth of the storage consumption. However, the central processing unit time increases slightly faster with the chain length than the storage consumption because the alternating linear scheme adopted in our work requires more iterations to achieve convergence for longer chains and a given rank. Finally, we demonstrate that the tensor-train approach to the quantum treatment of coupled excitons and phonons makes it possible to directly tackle the phenomenon of mutual self-trapping. We are able to confirm the main results of the Davydov theory, i.e., the dependence of the wave packet width and the corresponding stabilization energy on the exciton–phonon coupling strength, although only for a certain range of that parameter. In future work, our approach will allow calculations also beyond the validity regime of that theory and/or beyond the restrictions of the Fröhlich–Holstein type Hamiltonians. Y1 - 2022 U6 - https://doi.org/10.1063/5.0074948 VL - 156 IS - 2 SP - 024109 ER -