Abstract We are motivated by factorization of bosonic quantum dynamics and we study the corresponding Lie algebras, which can potentially be infinite dimensional. To characterize such factorization, we identify conditions for these Lie algebras to be finite dimensional. We consider cases where each free Hamiltonian term is itself an element of the generated Lie algebra. In our approach, we develop new tools to systematically divide skew-hermitian bosonic operators into appropriate subspaces, and construct specific sequences of skew-hermitian operators that are used to gauge the dimensionality of the Lie algebras themselves. The significance of our result relies on conditions that constrain only the independently controlled generators in a particular Hamiltonian, thereby providing an effective algorithm for verifying the finiteness of the generated Lie algebra. In addition, our results are tightly connected to mathematical work where the polynomials of creation and annihilation operators are known as the Weyl algebra. Our work paves the way for better understanding factorization of bosonic dynamics relevant to quantum control and quantum technology.
We study, classify, and explore the mathematical properties of finite-dimensional Lie algebras occurring in the quantum dynamics of single-mode and self-interacting bosonic systems. These Lie algebras are contained in the real skew-hermitian Weyl algebra
A^1, defined as the real subalgebra of the Weyl algebra
A1 consisting of all skew-hermitian polynomials. A central aspect of our analysis is the choice of basis for
A^1, which is composed of skew-symmetric combinations of two elements of the Weyl algebra called monomials, namely strings of creation and annihilation operators combined with their hermitian conjugate. Motivated by the quest for analytical solutions in quantum optimal control and dynamics, we aim at answering the following three fundamental questions: (i) What are the finite-dimensional Lie subalgebras in
A^1 generated by monomials alone? (ii)~What are the finite-dimensional Lie subalgebras in
A^1 that contain the free Hamiltonian? (iii) What are the non-abelian and finite-dimensional Lie subalgebras that can be faithfully realized in
A^1? We answer the first question by providing all possible realizations of all finite-dimensional non-abelian Lie algebras that are generated by monomials alone. We answer the second question by proving that any non-abelian and finite-dimensional subalgebra of
A^1 that contains a free Hamiltonian term must be a subalgebra of the Schrödinger algebra. We partially answer the third question by classifying all nilpotent and non-solvable Lie algebras that can be realized in
A^1, and comment on the remaining cases. Finally, we also discuss the implications of our results for quantum control theory. Our work constitutes an important stepping stone to understanding quantum dynamics of bosonic systems in full generality.
Classical simulability is ultimately determined by both the dynamics of a quantum system and the observables being evaluated. Lie-algebraic simulation exploits the latter to make exact polynomial-time classical simulations by propagating observables through low-dimensional invariant operator spaces. However, its conventional formulation in terms of polynomial-dimensional dynamical Lie algebras does not directly accommodate bosonic systems as their algebras are neither compact nor semisimple. In this contribution, we overcome this limitation, making bosonic systems accessible to the Lie-algebraic programme of exact polynomial-time classical simulation. We prove that expectation values, fixed-order correlation functions, including multi-time correlators and out-of-time-ordered correlators, and gradients are efficiently computable whenever their operator modules have polynomial dimension. This recovers Gaussian quantum optics and extends it to non-Gaussian input states, while identifying exact polynomial regimes of interacting non-Gaussian dynamics including bounded-photon Kerr and pair-hopping Hamiltonians and nilpotent polynomial phase dynamics. We show that unlike in the finite-dimensional spin and fermionic setting treated previously, a finite-dimensional bosonic generator algebra alone does not guarantee finite observable dynamics. We further derive a controlled perturbative hierarchy for squeezing beyond exact sector confinement and confirm the predicted error orders numerically. We also evaluate operator spreading on interacting chains of up to
400 modes and connect a topological doublon band with flux-reversed edge motion. These results provide a unified formalism for classifying, discovering, and systematically approximating tractable bosonic quantum dynamics with classical polynomial-time simulation.
Determining the unitary dynamics accessible from finite Hamiltonian resources is a central problem in Hamiltonian engineering and quantum control. Dynamical Lie algebras (DLAs) connect available control Hamiltonians with the reachable dynamics, but their use as a design tool for modifying Hamiltonian generator sets remains less developed. In this work, we develop a finite-dimensional DLA framework for three generator-set operations: composition, invariance, and reduction. For composition, we construct direct sums of component DLAs using spectral projectors on an auxiliary register. For invariance, we analyze when modifications of Pauli-string generating sets preserve the generated Lie algebra, and introduce algebraic diagnostics for added generators. For reduction, we consider compact reductive DLAs and show how projection onto selected simple ideals gives reduced generating sets whose Lie closures are the corresponding ideal sums. We illustrate these results with finite-dimensional examples and numerical checks, including direct-sum dimension addition, central-spin invariance diagnostics, and DLA-based ansatz reduction for block-local Hamiltonians. The results show how DLA structure can be used to diagnose controllability and guide Hamiltonian generator design under explicit algebraic assumptions.
Hermitian polynomials in the quadrature operators
x and
p are the basic ingredients of continuous-variable quantum information, and increasingly exotic ones are now routinely engineered in superconducting circuits and trapped ions. However, such a polynomial need not define unique quantum dynamics: it may define exactly one, none, or an infinite family of physically distinct ones. Such obstructions are invisible to algebraic arguments and can compromise finite-dimensional truncations, which may fail to converge or silently converge to the wrong dynamics. We develop a method, implemented as an algorithm, which decides which case occurs, and how large the family is, from the coefficients alone, with no functional-analytic input from the user. For Gaussian operators the answer is always one; for non-Gaussian ones, a regime continuous-variable quantum computing cannot do without, the other outcomes occur already for simple and physically natural Hamiltonians. We apply the algorithm to operators drawn from the recent literature.
The ability to navigate and control the state space of a quantum system is fundamental to quantum computing, quantum information processing, and advanced spectroscopy. This paper explores a rigorous mathematical framework for quantifying state space accessibility by leveraging the theory of operator commutation relations, particularly within the context of Lie algebras. We demonstrate that the span of iterated commutators of system Hamiltonians and control operators forms a dynamical Lie algebra, whose dimension and structure directly dictate the reachable states and unitary transformations. For closed quantum systems, the Lie algebra rank condition provides a powerful criterion for complete controllability, asserting that if the dynamical Lie algebra spans the entire space of anti-Hermitian operators, then any state can be reached from any other, and any unitary operation can be implemented. For open quantum systems, while environmental interactions introduce decoherence and dissipation, the Lie algebraic approach can still offer insights into the attainable subspaces and the limits of control, often involving Lie semigroups and Lie wedges. We present methods to analyze the accessibility of both pure and mixed states, illustrating how the algebraic properties encoded in commutation relations provide a precise measure of the manipulability of quantum systems. The implications extend to the design of robust quantum control protocols, optimized state preparation, and the development of error-resilient quantum technologies.
We study quantum dynamics generated by time-dependent Hamiltonians in Krylov space, the minimal subspace in which the evolution takes place. We establish a direct link between dynamics in the time-dependent Krylov subspace and the underlying Lie-algebraic structure of the Hamiltonian. We develop a general framework in which the dynamics in the time-dependent Krylov subspace is generated by ladder operators of the associated Lie algebra. In particular, we identify the minimal conditions under which the exact time-dependent Krylov dynamics is naturally determined by the interaction-picture Hamiltonian and governed by an embedded
sl(2,C) subalgebra. We further show that an exact single-exponential representation of the time-evolution operator gives rise to a distinct time-independent Krylov dynamics in a unitarily related basis, from which the exact time-dependent Krylov dynamics can nevertheless be recovered. We also extend the framework to the oscillator algebra as the simplest extension of the nilpotent Heisenberg--Weyl algebra, and provide further examples, including the translated and dilated harmonic oscillator, systems governed by closed Virasoro subalgebras, a spin in a rotating magnetic field, and higher-dimensional generalizations for multi-level systems. In addition, we introduce a new quantum speed limit to the complexity growth rate generated by a time-dependent generator and show that, for evolutions governed by a Lie algebra, it retains the same functional form as in the time-independent case. Remarkably, saturation of this bound is strongly affected by temporal driving and persists only when the Hamiltonian commutes with itself at different times. These results establish a unified framework for characterizing operator growth and Krylov complexity in time-dependent quantum systems with underlying Lie-algebraic structures.