When do perturbative approaches accurately capture the dynamics of complex quantum systems?

Driving of an open quantum system at finite temperature across first- and second-order quantum phase transitions

An open fully connected system of qubits at nonzero temperature is driven within a finite time interval along various paths in the space of its control parameters. The driving leads across finite-size precursors of first- and second-order quantum phase transition from factorized to entangled ground-state phases, aiming at the preparation of the complex ground state of the system at the final parameter point with maximal fidelity. During the drive, the system is coupled to a heat bath at a constant temperature, the dynamics being determined in a nonperturbative way by the method of hierarchical equations of motion. It is shown that the presence of the heat bath in combination with specific patterns of avoided crossings affecting the ground and excited states in the parameter region around the quantum phase transition may considerably improve the fidelity of preparation of the target ground state.

Self-consistent approach to the dynamics of excitation energy transfer in multichromophoric systems

Computationally tractable and reliable, albeit approximate, methods for studying exciton transport in molecular aggregates immersed in structured bosonic environments have been actively developed. Going beyond the lowest-order (Born) approximation for the memory kernel of the quantum master equation typically results in complicated and possibly divergent expressions. Starting from the memory kernel in the Born approximation, and recognizing the quantum master equation as the Dyson equation of Green's functions theory, we formulate the self-consistent Born approximation to resum the memory-kernel perturbation series in powers of the exciton-environment interaction. Our formulation is in the Liouville space and frequency domain and handles arbitrary exciton-environment spectral densities. In a molecular dimer coupled to an overdamped oscillator environment, we conclude that the self-consistent cycle significantly improves the Born-approximation energy-transfer dynamics. The dynamics in the self-consistent Born approximation agree well with the solutions of hierarchical equations of motion over a wide range of parameters, including the most challenging regimes of strong exciton-environment interactions, slow environments, and low temperatures. This is rationalized by the analytical considerations of coherence-dephasing dynamics in the pure-dephasing model. We find that the self-consistent Born approximation is good (poor) at describing energy transfer modulated by an underdamped vibration resonant (off-resonant) with the exciton energy gap. Nevertheless, it reasonably describes exciton dynamics in the seven-site model of the Fenna-Matthews-Olson complex in a realistic environment comprising both an overdamped continuum and underdamped vibrations.

Canonically Consistent Quantum Master Equation

We put forth a new class of quantum master equations that correctly reproduce the asymptotic state of an open quantum system beyond the infinitesimally weak system-bath coupling limit. Our method is based on incorporating the knowledge of the reduced steady state into its dynamics. The correction not only steers the reduced system toward a correct steady state but also improves the accuracy of the dynamics, thereby refining the archetypal Born-Markov weak-coupling second-order master equations. In case of equilibrium, we use the exact mean-force Gibbs state to correct the Redfield quantum master equation. By benchmarking it with the exact solution of the damped harmonic oscillator, we show that our method also helps correct the long-standing issue of positivity violation, albeit without complete positivity. Our method of a canonically consistent quantum master equation opens a new perspective in the theory of open quantum systems leading to a reduced density matrix accurate beyond the commonly used Redfield and Lindblad equations, while retaining the same conceptual and numerical complexity.

Explaining the Efficiency of Photosynthesis: Quantum Uncertainty or Classical Vibrations?

Photosynthetic organisms are known to use a mechanism of vibrationally assisted exciton energy transfer to efficiently harvest energy from light. The importance of quantum effects in this mechanism is a long-standing topic of debate, which has traditionally focused on the role of excitonic coherences. Here, we address another recent claim: that the efficient energy transfer in the Fenna-Matthews-Olson complex relies on nuclear quantum uncertainty and would not function if the vibrations were classical. We present a counter-example to this claim, showing by trajectory-based simulations that a description in terms of quantum electrons and classical nuclei is indeed sufficient to describe the funneling of energy to the reaction center. We analyze and compare these findings to previous classical-nuclear approximations that predicted the absence of an energy funnel and conclude that the key difference and the reason for the discrepancy is the ability of the trajectories to properly account for Newton's third law.

Weak and Ultrastrong Coupling Limits of the Quantum Mean Force Gibbs State

The Gibbs state is widely taken to be the equilibrium state of a system in contact with an environment at temperature T. However, non-negligible interactions between system and environment can give rise to an altered state. Here, we derive general expressions for this mean force Gibbs state, valid for any system that interacts with a bosonic reservoir. First, we derive the state in the weak coupling limit and find that, in general, it maintains coherences with respect to the bare system Hamiltonian. Second, we develop a new expansion method suited to investigate the ultrastrong coupling regime. This allows us to derive the explicit form for the mean force Gibbs state, and we find that it becomes diagonal in the basis set by the system-reservoir interaction instead of the system Hamiltonian. Several examples are discussed including a single qubit, a three-level V-system, and two coupled qubits all interacting with bosonic reservoirs. The results shed light on the presence of coherences in the strong coupling regime, and provide key tools for nanoscale thermodynamics investigations.

Exploiting the Causal Tensor Network Structure of Quantum Processes to Efficiently Simulate Non-Markovian Path Integrals

In the path integral formulation of the evolution of an open quantum system coupled to a Gaussian, noninteracting environment, the dynamical contribution of the latter is encoded in an object called the influence functional. Here, we relate the influence functional to the process tensor-a more general representation of a quantum stochastic process-describing the evolution. Then, we use this connection to motivate a tensor network algorithm for the simulation of multitime correlations in open systems, building on recent work where the influence functional is represented in terms of time evolving matrix product operators. By exploiting the symmetries of the influence functional, we are able to use our algorithm to achieve orders-of-magnitude improvement in the efficiency of the resulting numerical simulation. Our improved algorithm is then applied to compute exact phonon emission spectra for the spin-boson model with strong coupling, demonstrating a significant divergence from spectra derived under commonly used assumptions of memorylessness.

Modelling the ultra-strongly coupled spin-boson model with unphysical modes

A quantum system weakly coupled to a zero-temperature environment will relax, via spontaneous emission, to its ground-state. However, when the coupling to the environment is ultra-strong the ground-state is expected to become dressed with virtual excitations. This regime is difficult to capture with some traditional methods because of the explosion in the number of Matsubara frequencies, i.e., exponential terms in the free-bath correlation function. To access this regime we generalize both the hierarchical equations of motion and pseudomode methods, taking into account this explosion using only a biexponential fitting function. We compare these methods to the reaction coordinate mapping, which helps show how these sometimes neglected Matsubara terms are important to regulate detailed balance and prevent the unphysical emission of virtual excitations. For the pseudomode method, we present a general proof of validity for the use of superficially unphysical Matsubara-modes, which mirror the mathematical essence of the Matsubara frequencies.

Optimizing co-operative multi-environment dynamics in a dark-state-enhanced photosynthetic heat engine

We analyze the role of coherent, non-perturbative system-bath interactions in a photosynthetic heat engine. Using the reaction-coordinate formalism to describe the vibrational phonon-environment in the engine, we analyze the efficiency around an optimal parameter regime predicted in earlier studies. We show that, in the limit of high-temperature photon irradiation, the phonon-assisted population transfer between bright and dark states is suppressed due to dephasing from the photon environment, even in the Markov limit where we expect the influence of each bath to have an independent and additive effect on the dynamics. Manipulating the phonon bath properties via its spectral density enables us to identify both optimal low- and high-frequency regimes where the suppression can be removed. This suppression of transfer and its removal suggests that it is important to consider carefully the non-perturbative and cooperative effects of system-bath environments in designing artificial photosynthetic systems and also that manipulating inter-environmental interactions could provide a new multidimensional "lever" by which photocells and other types of quantum devices can be optimized.

Plasmonic bio-sensing for the Fenna-Matthews-Olson complex

We study theoretically the bio-sensing capabilities of metal nanowire surface plasmons. As a specific example, we couple the nanowire to specific sites (bacteriochlorophyll) of the Fenna-Matthews-Olson (FMO) photosynthetic pigment protein complex. In this hybrid system, we find that when certain sites of the FMO complex are subject to either the suppression of inter-site transitions or are entirely disconnected from the complex, the resulting variations in the excitation transfer rates through the complex can be monitored through the corresponding changes in the scattering spectra of the incident nanowire surface plasmons. We also find that these changes can be further enhanced by changing the ratio of plasmon-site couplings. The change of the Fano lineshape in the scattering spectra further reveals that “site 5” in the FMO complex plays a distinct role from other sites. Our results provide a feasible way, using single photons, to detect mutation-induced, or bleaching-induced, local defects or modifications of the FMO complex, and allows access to both the local and global properties of the excitation transfer in such systems.