Quantum Dynamics in Phase Space using Projected von Neumann Bases

Control of concerted back-to-back double ionization dynamics in helium

Double ionization (DI) is a fundamental process that despite its apparent simplicity provides rich opportunities for probing and controlling the electronic motion. Even for the simplest multielectron atom, helium, new DI mechanisms are still being found. To first order in the field strength, a strong external field doubly ionizes the electrons in helium such that they are ejected into the same direction (front-to-back motion). The ejection into opposite directions (back-to-back motion) cannot be described to first order, making it a challenging target for control. Here, we address this challenge and optimize the field with the objective of back-to-back double ionization using a (1 + 1)-dimensional model. The optimization is performed using four different control procedures: (1) short-time control, (2) derivative-free optimization of basis expansions of the field, (3) the Krotov method, and (4) control of the classical equations of motion. All four procedures lead to fields with dominant back-to-back motion. All the fields obtained exploit essentially the same two-step mechanism leading to back-to-back motion: first, the electrons are displaced by the field into the same direction. Second, after the field turns off, the nuclear attraction and the electron-electron repulsion combine to generate the final motion into opposite directions for each electron. By performing quasi-classical calculations, we confirm that this mechanism is essentially classical.

A-VCI: A flexible method to efficiently compute vibrational spectra

The adaptive vibrational configuration interaction algorithm has been introduced as a new method to efficiently reduce the dimension of the set of basis functions used in a vibrational configuration interaction process. It is based on the construction of nested bases for the discretization of the Hamiltonian operator according to a theoretical criterion that ensures the convergence of the method. In the present work, the Hamiltonian is written as a sum of products of operators. The purpose of this paper is to study the properties and outline the performance details of the main steps of the algorithm. New parameters have been incorporated to increase flexibility, and their influence has been thoroughly investigated. The robustness and reliability of the method are demonstrated for the computation of the vibrational spectrum up to 3000 cm^-1 of a widely studied 6-atom molecule (acetonitrile). Our results are compared to the most accurate up to date computation; we also give a new reference calculation for future work on this system. The algorithm has also been applied to a more challenging 7-atom molecule (ethylene oxide). The computed spectrum up to 3200 cm^-1 is the most accurate computation that exists today on such systems.

Efficient computations of quantum canonical Gibbs state in phase space

The Gibbs canonical state, as a maximum entropy density matrix, represents a quantum system in equilibrium with a thermostat. This state plays an essential role in thermodynamics and serves as the initial condition for nonequilibrium dynamical simulations. We solve a long standing problem for computing the Gibbs state Wigner function with nearly machine accuracy by solving the Bloch equation directly in the phase space. Furthermore, the algorithms are provided yielding high quality Wigner distributions for pure stationary states as well as for Thomas-Fermi and Bose-Einstein distributions. The developed numerical methods furnish a long-sought efficient computation framework for nonequilibrium quantum simulations directly in the Wigner representation.