Multistate electron transfer dynamics in the condensed phase: Exact calculations from the reduced hierarchy equations of motion approach

Entanglement dynamics of two qubits coupled individually to Ohmic baths

Developed originally for the Holstein polaron, the Davydov D1 ansatz is an efficient, yet extremely accurate trial state for time-dependent variation of the spin-boson model [N. Wu, L. Duan, X. Li, and Y. Zhao, J. Chem. Phys. 138, 084111 (2013)]. In this work, the Dirac-Frenkel time-dependent variational procedure utilizing the Davydov D1 ansatz is implemented to study entanglement dynamics of two qubits under the influence of two independent baths. The Ohmic spectral density is used without the Born-Markov approximation or the rotating-wave approximation. In the strong coupling regime finite-time disentanglement is always found to exist, while at the intermediate coupling regime, the entanglement dynamics calculated by Davydov D1 ansatz displays oscillatory behavior in addition to entanglement disappearance and revival.

Optimizing hierarchical equations of motion for quantum dissipation and quantifying quantum bath effects on quantum transfer mechanisms

We present an optimized hierarchical equations of motion theory for quantum dissipation in multiple Brownian oscillators bath environment, followed by a mechanistic study on a model donor-bridge-acceptor system. We show that the optimal hierarchy construction, via the memory-frequency decomposition for any specified Brownian oscillators bath, is generally achievable through a universal pre-screening search. The algorithm goes by identifying the candidates for the best be just some selected Padé spectrum decomposition based schemes, together with a priori accuracy control criterions on the sole approximation, the white-noise residue ansatz, involved in the hierarchical construction. Beside the universal screening search, we also analytically identify the best for the case of Drude dissipation and that for the Brownian oscillators environment without strongly underdamped bath vibrations. For the mechanistic study, we quantify the quantum nature of bath influence and further address the issue of localization versus delocalization. Proposed are a reduced system entropy measure and a state-resolved constructive versus destructive interference measure. Their performances on quantifying the correlated system-environment coherence are exemplified in conjunction with the optimized hierarchical equations of motion evaluation of the model system dynamics, at some representing bath parameters and temperatures. Analysis also reveals the localization to delocalization transition as temperature decreases.

Open Quantum Dynamics Calculations with the Hierarchy Equations of Motion on Parallel Computers

Calculating the evolution of an open quantum system, i.e., a system in contact with a thermal environment, has presented a theoretical and computational challenge for many years. With the advent of supercomputers containing large amounts of memory and many processors, the computational challenge posed by the previously intractable theoretical models can now be addressed. The hierarchy equations of motion present one such model and offer a powerful method that remained under-utilized so far due to its considerable computational expense. By exploiting concurrent processing on parallel computers the hierarchy equations of motion can be applied to biological-scale systems. Herein we introduce the quantum dynamics software PHI, that solves the hierarchical equations of motion. We describe the integrator employed by PHI and demonstrate PHI's scaling and efficiency running on large parallel computers by applying the software to the calculation of inter-complex excitation transfer between the light harvesting complexes 1 and 2 of purple photosynthetic bacteria, a 50 pigment system.