Higher-order solutions to non-Markovian quantum dynamics via hierarchical functional derivative

Solving realistic quantum systems coupled to an environment is a challenging task. Here we develop a hierarchical functional derivative (HFD) approach for efficiently solving the non-Markovian quantum trajectories of an open quantum system embedded in a bosonic bath. An explicit expression for arbitrary order HFD equation is derived systematically. Moreover, it is found that for an analytically solvable model, this hierarchical equation naturally terminates at a given order and thus becomes exactly solvable. This HFD approach provides a systematic method to study the non-Markovian quantum dynamics of an open system coupled to a bosonic environment.