Solving quantum dynamics with a Lie algebra decoupling method 3
effects of interest are weak, such perturbative solutions can be used to accurately model the behaviour of
experiments in the laboratory. Many mathematical methods have been developed to express, manipulate
and truncate exponential operators [7–12]. In order to control the full quantum behaviour of individual
quantum systems, however, such as in quantum optics, quantum information and for quantum technologies,
it can be necessary to go beyond these approximations.
One mathematical technique that can be used to solve quantum dynamics exactly makes use of the
Lie algebra induced by the Hamiltonian. The method is based on a Lie algebra decoupling theorem, which
was put forwards by Wei and Norman in 1963 [13]. At its core is the observation that a finite Lie algebra
generated by a Hamiltonian can be used as a basis that allows for a set of scalar differential equations to
be derived. Should the resulting differential equations have analytic solutions, the system dynamics can be
solved exactly.
In this work, we provide a pedagogical introduction to the Lie algebra decoupling theorem and its
application. Specifically, we apply the theorem towards solving the dynamics of three Hamiltonians that
frequently appear in quantum optics and adjacent fields. The Hamiltonians we consider have both linear
and quadratic interaction terms, where each interaction term depends on time in an arbitrary fashion. As
such, the theorem allows us to provide a treatment of the most general case of Gaussian dynamics.
The work is structured as follows. Section 2provides a mathematical introduction to Lie groups and Lie
algebras. Section 3contains a proof of the decoupling theorem, as well as a straight-forward recipe for how
it can be applied. Then, in Sections 4and 5, we apply the theorem to solve the dynamics of a Hamiltonian
with time-dependent linear and quadratic interaction terms, respectively. There are also ways in which
these solutions can be combined to represent the most general Gaussian dynamics, which we demonstrate in
Section 6. The work is concluded with a discussion of symmetries and extensions to open-system dynamics in
Section 7, as well as some final remarks in Section 8. In the appendices that follow, we provide an extension
of the decoupling theorem to phase space, as well as detailed calculations for each section in the main text.
2. Introduction to Lie algebras and Lie groups
The Lie algebra decoupling theorem can be understood and applied without deep knowledge of Lie algebras,
but for the interested reader, we here provide a basic mathematical introduction to Lie groups, Lie algebras,
and the connection between them. In each case, we start the discussion with a formal definition, then provide
a few examples that commonly appear in physics.
2.1. Lie groups
To understand Lie algebras, it helps to first introduce the concept of a Lie group. We start by recounting
the formal definition of a group.
Definition (group) A group (G, ∗)is a set Gwith group elements gand a binary operation ∗such that
G×G→G, which satisfies three conditions:
(i) Associativity: For any three elements x, y, z ∈G, we have (x∗y)∗z=x∗(y∗z).
(ii) Identity: The group must include an identity element ∈G. It must hold that ∀g∈G, multiplying gby
the identity element leaves ginvariant, such that ∗g=g∗=g.
(iii) Inverse: For each element in the group, there must be an inverse element. That is, for each element
g∈Gthere is some inverse element ¯g∈Gsuch that g∗¯g=¯g∗g=.
Groups arise in many separate context in physics. Perhaps one of the simplest examples of a group is the
set of discrete rotations that leave a square invariant. That is, the group elements correspond to rotation
actions (acting on a square) and the group operation corresponds to the composition of actions. This group
is known as Z4and contains rotations of zero, 90○, 180○, and 270○degrees. This set of rotations satisfies the
group axioms, since any combination of the group elements satisfies the same symmetry. That is, we can
rotate the square by first 90○, then 180○, and it remains invariant. The associativity criterion is satisfied by
a proper representation of the group elements, such as a matrix. Lastly, the identity is equivalent to not
rotating the square, and the inverse can easily be constructed by rotating the square in the reverse direction
by the same amount of degrees.