Repository: Freie Universität Berlin, Math Department

WavePacket: A Matlab package for numerical quantum dynamics. II: Open quantum systems, optimal control, and model reduction

Schmidt, B. and Hartmann, C. (2018) WavePacket: A Matlab package for numerical quantum dynamics. II: Open quantum systems, optimal control, and model reduction. Comp. Phys. Comm., 228 . pp. 229-244.

Full text not available from this repository.

Official URL:


WavePacket is an open-source program package for numeric simulations in quantum dynamics. It can solve time-independent or time-dependent linear Schrödinger and Liouville-von Neumann-equations in one or more dimensions. Also coupled equations can be treated, which allows, e.g., to simulate molecular quantum dynamics beyond the Born-Oppenheimer approximation. Optionally accounting for the interaction with external electric fields within the semi-classical dipole approximation, WavePacket can be used to simulate experiments involving tailored light pulses in photo-induced physics or chemistry. Being highly versatile and offering visualization of quantum dynamics 'on the fly', WavePacket is well suited for teaching or research projects in atomic, molecular and optical physics as well as in physical or theoretical chemistry. Building on the previous Part I which dealt with closed quantum systems and discrete variable representations, the present Part II focuses on the dynamics of open quantum systems, with Lindblad operators modeling dissipation and dephasing. This part also describes the WavePacket function for optimal control of quantum dynamics, building on rapid monotonically convergent iteration methods. Furthermore, two different approaches to dimension reduction implemented in WavePacket are documented here. In the first one, a balancing transformation based on the concepts of controllability and observability Gramians is used to identify states that are neither well controllable nor well observable. Those states are either truncated or averaged out. In the other approach, the H2-error for a given reduced dimensionality is minimized by H2 optimal model reduction techniques, utilizing a bilinear iterative rational Krylov algorithm. The present work describes the MATLAB version of WavePacket 5.3.0 which is hosted and further developed at the Sourceforge platform, where also extensive Wiki-documentation as well as numerous worked-out demonstration examples with animated graphics can be found.

Item Type:Article
Subjects:Physical Sciences > Chemistry > Physical Chemistry
Physical Sciences > Physics > Chemical Physics
Physical Sciences > Physics > Mathematical & Theoretical Physics > Quantum Mechanics
Mathematical and Computer Sciences > Mathematics > Mathematical Modelling
Divisions:Department of Mathematics and Computer Science > Institute of Mathematics > BioComputing Group
ID Code:2090
Deposited By: BioComp Admin
Deposited On:19 Jul 2017 09:34
Last Modified:02 May 2018 13:12

Repository Staff Only: item control page