Tavis–Cummings model
In quantum optics, the Tavis–Cummings model is a theoretical model to describe an ensemble of identical two-level atoms coupled symmetrically to a single-mode quantized bosonic field. The model extends the Jaynes–Cummings model to larger spin numbers that represent collections of multiple atoms. It differs from the Dicke model in its use of the rotating-wave approximation to conserve the number of excitations of the system.
Originally introduced by Michael Tavis and Fred Cummings in 1968 to unify representations of atomic gases in electromagnetic fields under a single fully quantum Hamiltonian — as Robert Dicke had done previously using perturbation theory — the Tavis–Cummings model's restriction to a single field-mode with negligible counterrotating interactions simplifies the system's mathematics while preserving the breadth of its dynamics.
The model demonstrates superradiance, bright and dark states, Rabi oscillations and spontaneous emission, and other features of interest in quantum electrodynamics, quantum control and computation, atomic and molecular physics, and many-body physics. The model has been experimentally tested to determine the conditions of its viability, and realized in semiconducting and superconducting qubits.
Hamiltonian
The Tavis–Cummings model assumes that for the purposes of electromagnetic interactions, atomic structures are dominated by their dipole, as they are for distant neutral atoms in the weak-field limit. Thus the only atomic quantity under consideration is its angular momentum, not its position nor fine electronic structure. Furthermore, the model asserts the atoms to be sufficiently distant that they don't interact with each-other, only with the electromagnetic field, modeled as a bosonic field.Formal derivation
For two atomic-electronic states separated by a Bohr frequency, then transitions between the ground- and excited-states and are mediated by Pauli operators:,, and, and the Hamiltonian separating these energy states in the th atom is. With independent atoms each subject to this energy gap, the total atomic Hamiltonian is thus with total spin operators.Similarly, in a free field with no modal restrictions, creation and annihilation operators dictate the presence of photons in each mode: with wave number, polarization, and frequency. If the dynamics occur within a sufficiently small cavity, only one mode will couple to the atom, thus the field Hamiltonian simplifies to, just as in the Jaynes-Cummings and Dicke models.[|thumb|Schematic diagram of the Tavis–Cummings physical model, representing an ensemble of two-level atoms (four dots on red and green levels) interacting symmetrically with a single mode photonic field (blue standing wave), isolated within a cavity. The atomic level separation is, the cavity's resonant frequency is, and the coupling strength of atom-field interactions is.|400x400px]
Finally, interactions between atoms and the field is determined by the atomic dipole, rendered quantumly as an operator, and the similarly expressed electric field at the atoms' centers, thus which acts on both the qubit and bosonic degrees of freedom. The dipole operator couples the excited and ground states of each atom, while the electric free field solution is:
, which at a static point evaluates as:
, thus the interaction Hamiltonian expands as
Here, specifies the coupling strength of the total dipole to each electric field mode, and functioning as a Rabi frequency that scales with ensemble size due to the Pythagorean addition of single-atom dipoles. Then, in the rotating frame,, which results in corotating terms,, and counterrotating terms and . When and, the corotating terms accumulate phase very slowly, while the counterrotating terms accumulate phase too fast to significantly affect time-ordered integrals, thus the rotating wave approximation allows counterrotating terms to drop in the rotating frame. The cavity permits only one field mode with energy sufficiently close to the Bohr energy,, so the final form of the interaction Hamiltonian is for dephasing.
In total, the Tavis–Cummings Hamiltonian includes the atomic and photonic self-energies and the atom-field interaction:
Symmetries
The Tavis–Cummings model as described above exhibits two symmetries arising from the Hamiltonian's commutation with excitation number and angular momentum magnitude. Since, it is possible to find simultaneous eigenstates such that:The quantum number is bounded by, and, but due to the infinity of Fock space, excitation number is unbounded above, unlike angular momentum projection quantum numbers. Just as the Jaynes-Cummings Hamiltonian block-diagonalizes into infinite blocks of constant excitation number, the Tavis–Cummings Hamiltonian block-diagonlizes into infinite blocks of size up to with constant, and within these larger blocks, further block-diagonalizes into blocks of size where, with constant cooperation number. The size of each of these smallest blocks determine the bounds of the final quantum number that specifies the eigenenergy:, with signifying the ground state of each irrep, and the maximally excited state.
Dynamics
Under the simplifications of real and quasistatically small, the Hamiltonian becomes, whose matrix elements one can express in a joint Dicke and Fock basis such that and. Necessarily,, so the matrix elements are as follows:From these elements, one can express Schrödinger equations of motion to demonstrate the photon field's ability to mediate entanglement formation between atoms without atom-atom interactions:
for which the fine-tuned, multivariate dependence on quantum numbers demonstrates the difficulty of solving the Tavis-Cumming model's eigensystem. Here, a few approximate methods, and an exact solution involving Stark shifts and Kerr nonlinearities follow.
Spectrum approximations
In 1969, Tavis and Cummings found approximate eigenenergies and eigenstates for a nondimensionalized Hamiltonian in three different regimes of approximation: first, for near the ground-state of each irrep; second, for and when each atom sees a highly saturated "averaged" field; third, for with sparse excitations. In all solutions, eigenstates are related to Dicke-Fock joint states by, for coefficients that are solved from the Hamiltonian spectrum.For close to zero, a differential approach provides the eigenvalues: with, average photon number, and differential coefficient.
For an averaged photon field when in a large photon-rich system, the off-diagonal matrix elements in replace with, and each two-level atom interacts independently with a photon field that conveys no information about the other atoms. For each of these atoms, there are two dressed eigenstates coupled to "pseudophoton" number states:, with single-atom eigenenergies. Superpositions of these single-atom dressed states construct the full eigenstates according to addition of angular momenta, mediated by Clebsch-Gordan coefficients. The full eigenvalues are approximately.
For when the atomic ensemble, and not the photon field, is averaged, then the off-diagonal elements in approximate as, smearing over the atomic degrees of freedom in. Eigenstates are constructed as weighted superpositions of photon number states coupled to a spin state, but the eigenvalues are now.
Bethe ansatz
In 1996, Nikolay Bogoliubov, Robin Bullough, and Jussi Timonen found that adding quadratic excitation-dependent terms to the Tavis–Cummings Hamiltonian allowed for an exact analytic eigensystem. In the limit where these Kerr and Stark shifts vanish, this solution can recover the eigensystem of the unmodified Tavis–Cummings system.Including a Kerr term of the form, or a Stark term equivalently results in a new Hamiltonian:, which obeys the same operator symmetries as does the unmodified Tavis–Cummings Hamiltonian, and which reduces to Tavis–Cummings in the limit. Thus the transformation preserves the dynamics and shares joint eigenvectors with the untransformed Hamiltonian. The transformed Hamiltonian, explicitly, is for new parameters and, which is integrable using quantum inverse methods. Separating the dynamics into two complex-parametrized operator matrices, one acting on the bosonic degrees of freedom and the other on the spin degrees of freedom produces a monodromy matrix whose determinant is directly proportional to, whose trace is proportional to, and whose trace's parametric derivative is proportional to. Manipulating the monodromy matrix allows its spectral parameter to determine the Hamiltonian eigenstates and eigenenergies as the complex roots of a Bethe ansatz. Every for must satisfy the following Bethe equations:
then the Hamiltonian eigenvalues arise from the roots as:
In the limit where, the above Bethe equations simplify to, and the eigenenergies to. Eigenstates follow similarly.
Experiments
The Tavis–Cummings model has seen numerous experimental implementations verifying its phenomena, including several since 2009 virtually realizing the model on quantum computational platforms like superconducting qubits and circuit QED. Such experiments utilize the Tavis–Cummings Hamiltonian's ability to generate superradiance wherein the artificial atoms emit and absorb light from the field coherently, as though they were a single atom with a large total angular momentum. Superradiance, scaling dipole-interaction strength, and other features allow Tavis–Cummings-type dynamics to manifest quantum computationally and metrologically desirable states, such as Dicke states through global interactions, as was explored in the 2003 paper by Tessier et al.One realization by Tuchman et al., in 2006, used a stream of ultracold Rubidium-87 atoms, and observed cooperation number, or 12% its maximum possible value, indicating very high interatomic coherence relative to experimental capabilities of the time. This experiment also confirmed the scaling of the dipole-interaction; at the level of single atoms, dipole interactions are much weaker than monopolar interactions, so the ability of Tavis–Cummings dynamics to counteract the weakness of with quadrature addition of dipoles makes neutral atom control more feasible.
Circuit QED
A seminal result from Fink et al. in 2009 involved 3 transmons as virtual "atoms" with qubit-dependent Bohr frequencies for controllable Josephson energy and experimentally determined single electron charging energy, inside a microwave waveguide resonator which supplies a standing electric field at GHz. To ensure symmetric coupling of the qubits to the field, each transmon was placed at an antinode of the standing wave, and to best conserve excitations by minimizing photon leakage, the resonator was kept ultracold which ensured a high quality factor. Manipulating each qubit's Bohr frequency so only one qubit resonated with the field, the team measured each single-qubit coupling strength, then reintroduced the other qubits to compare the total coupling strength with the average strength of the resonating qubits,, confirming that. In addition, the team observed bright and dark states characterized by high emission rates and zero emissions respectively, for 2 and 3 active qubits, with the 3-qubit bright and dark states each being degenerate.In addition to superconducting qubits, semiconducting qubits have also been a platform for Tavis–Cummings dynamics, such as in a 2018 investigation by van Woerkom et al. at ETH Zürich, in which two qubits constructed of double quantum dots coupled to a SQUID resonator, with the two DQDs separated by a distance of 42μm. The micrometer regime is a far greater distance than that over which semiconducting qubits had previously achieved entanglement, and the difficulty of long-range interactions in semiconducting qubits was at the time a major weakness compared to other quantum computing platforms, for which the Tavis–Cummings model's ability to form entanglement through global atom-field interactions is one solution. By observing the reflection amplitude of field waves between the SQUID array and the DQDs, the team isolated the photon number states as they smoothly coupled to the first qubit to form superpositional Jaynes-Cummings eigenstates when the first qubit tuned to the resonator. Similarly, they observed these hybrid states shift into a pair of bright states and a dark state when the second qubit was tuned to resonance. In addition to physical photons mediating the long-range entanglement at, the team found similar energy shifts at signalling qubits interacting with "virtual" photons, measured by the phase shift of the field rather than the reflection amplitude.