We present a recursive formula for the computation of the static effective Hamiltonian of a system under a fast-oscillating drive. Our analytical result is well-suited to symbolic calculations performed by a computer and can be implemented to arbitrary order, thus overcoming limitations of existing time-dependent perturbation methods and allowing computations that were impossible before. We also provide a simple diagrammatic tool for calculation and treat illustrative examples. By construction, our method applies directly to both quantum and classical systems; the difference is left to a low-level subroutine. This sheds light on the relationship between seemingly disconnected independently developed methods in the literature and has direct applications in quantum engineering.
We present a recursive formula for the computation of the static effective Hamiltonian of a system under a fast-oscillating drive. Our analytical result is well-suited to symbolic calculations performed by a computer and can be implemented to arbitrary order, thus overcoming limitations of existing time-dependent perturbation methods and allowing computations that were impossible before. We also provide a simple diagrammatic tool for calculation and treat illustrative examples. By construction, our method applies directly to both quantum and classical systems; the difference is left to a low-level subroutine. This sheds light on the relationship between seemingly disconnected independently developed methods in the literature and has direct applications in quantum engineering.
In this work, we present a new diagrammatic method for computing the effective Hamiltonian of driven nonlinear oscillators. At the heart of our method is a self-consistent perturbation expansion developed in phase space, which establishes a direct correspondence between the diagram and algebra. Each diagram corresponds to a Hamiltonian term, the prefactor of which, like those in Feynman diagrams, involves a simple counting of topologically equivalent diagrams. Leveraging the algorithmic simplicity of our diagrammatic method, we provide a readily available computer program that generates the effective Hamiltonian to arbitrary order. We show the consistency of our schemes with existing perturbation methods such as the Schrieffer-Wolff method. Furthermore, we recover the classical harmonic balance scheme from our result in the limit of ℏ→0. Our method contributes to the understanding of dynamic control within quantum systems and achieves precision essential for advancing future quantum information processors. To demonstrate its value and versatility, we analyze five examples from the field of superconducting circuits. These include an experimental proposal for the Hamiltonian stabilization of a three-legged Schrödinger cat, modeling of energy renormalization phenomena in superconducting circuits experiments, a comprehensive characterization of multiphoton resonances in a driven transmon, a proposal for an inductively shunted transmon circuit, and a characterization of classical ultra-subharmonic bifurcation in driven oscillators. Lastly, we benchmark the performance of our method by comparing it with experimental data and exact Floquet numerical diagonalization.
The dynamics of classical and quantum systems, which are driven by a high-frequency $(\ensuremath{\omega})$ field, is investigated. For classical systems, the motion is separated into a slow part and a fast part. The motion for the slow part is computed perturbatively in powers of ${\ensuremath{\omega}}^{\ensuremath{-}1}$ to the order ${\ensuremath{\omega}}^{\ensuremath{-}4},$ and the corresponding time independent Hamiltonian is calculated. Such an effective Hamiltonian for the corresponding quantum problem is computed to the order ${\ensuremath{\omega}}^{\ensuremath{-}4}$ in a high-frequency expansion. Its spectrum is the quasienergy spectrum of the time dependent quantum system. The classical limit of this effective Hamiltonian is the classical effective time independent Hamiltonian. It is demonstrated that this effective Hamiltonian gives the exact quasienergies and quasienergy states of some simple examples, as well as the lowest resonance of a nontrivial model for an atom trap. The theory that is developed in this paper is useful for the analysis of atomic motion in atom traps of various shapes.
Abstract Time‐driven quantum systems are important in many different fields of physics as cold atoms, solid state, optics, etc. Many of their properties are encoded in the time evolution operator or the effective Hamiltonian. Finding these operators usually requires very complex calculations that often involve some approximations. To perform this task, a systematic scheme that can be cast in the form of a symbolic computational algorithm is presented. It is suitable for periodic and non‐periodic potentials and, for convoluted systems, can also be adapted to yield numerical solutions. The method exploits the structure of the associated Lie group and a decomposition of the evolution operator on each group generator. To illustrate the use of the method, five examples are provided: harmonic oscillator with time‐dependent frequency (Paul trap), modulated optical lattice, time‐driven quantum oscillator, a step‐wise driving of a free particle, and the non‐periodic Caldirola‐Kanai Hamiltonian. To the extent of the authors' knowledge, whereas the exact form of Paul trap's evolution operator is well known, its effective Hamiltonian was until now unknown. The remaining four examples accurately reproduce previous results.
Parametric gates and processes engineered from the perspective of the static effective Hamiltonian of a driven system are central to quantum technology. However, the perturbative expansions used to derive static effective models may not be able to efficiently capture all the relevant physics of the original system. In this work, we investigate the conditions for the validity of the usual low-order static effective Hamiltonian used to describe a Kerr oscillator under a squeezing drive. This system is of fundamental and technological interest. In particular, it has been used to stabilize Schrödinger cat states, which have applications for quantum computing. We compare the states and energies of the effective static Hamiltonian with the exact Floquet states and quasi-energies of the driven system and determine the parameter regime where the two descriptions agree. Our work brings to light the physics that is left out by ordinary static effective treatments and that can be explored by state-of-the-art experiments.
Based on a previously developed recursive approach for calculating the short-time expansion of the propagator for systems with time-independent potentials and its time-dependent generalization for simple single-particle systems, in this paper we present a full extension of this formalism to a general quantum system with many degrees of freedom in a time-dependent potential. Furthermore, we also present a recursive approach for the velocity-independent part of the effective potential, which is necessary for calculating diagonal amplitudes and partition functions, as well as an extension from the imaginary-time formalism to the real-time one, which enables to study the dynamical properties of quantum systems. The recursive approach developed here allows an analytic derivation of the short-time expansion to orders that have not been accessible before, using the implemented SPEEDUP symbolic calculation code. The analytically derived results are extensively numerically verified by treating several models in both imaginary and real time.
Deriving effective Hamiltonian models plays an essential role in quantum theory, with particular emphasis in recent years on control and engineering problems. In this work, we present two symbolic methods for computing effective Hamiltonian models: the Non-perturbative Analytical Diagonalization (NPAD) and the Recursive Schrieffer-Wolff Transformation (RSWT). NPAD makes use of the Jacobi iteration and works without the assumptions of perturbation theory while retaining convergence, allowing to treat a very wide range of models. In the perturbation regime, it reduces to RSWT, which takes advantage of an in-built recursive structure where remarkably the number of terms increases only linearly with perturbation order, exponentially decreasing the number of terms compared to the ubiquitous Schrieffer-Wolff method. Both methods consist of elementary algebra and can be easily automated to obtain closed-form expressions. To demonstrate the application of the methods, we study the ZZ interaction of superconducting qubits systems in two different parameter regimes. In the near-resonant regime, the method produces an analytical expression for the effective ZZ interaction strength, with an improvement of at least one order of magnitude compared to neglecting the comparably further detuned state. In the dispersive regime, we calculate perturbative corrections up to the 8th order and accurately identify a regime where ZZ interaction is suppressed in a simple model with only resonator-mediated coupling.