Fractional revivals, multiple-Schrödinger-cat states and quantum carpets in the interaction of a qubit with N qubits

Fractional revivals, multiple-Schrödinger-cat states and quantum carpets in the interaction of a qubit with qubits

Shane Dooley dooleysh@gmail.com Quantum Information Science, School of Physics and Astronomy, University of Leeds, Leeds LS2 9JT, United Kingdom.    Timothy P. Spiller York Centre for Quantum Technologies, Department of Physics, University of York, York YO10 5DD, U.K. Quantum Information Science, School of Physics and Astronomy, University of Leeds, Leeds LS2 9JT, United Kingdom.
July 16, 2019
Abstract

We study the dynamics of a system comprised of a single qubit interacting equally with qubits (a “spin star” system). Although this model can be solved exactly, the exact solution does not give much intuition for the dynamics of the model. Here, we find an approximation that gives some insight into the dynamics for a particular class of initial spin-coherent states of the qubits. We find an effective Hamiltonian for the system that is a finite Kerr (one-axis twisting) Hamiltonian for the qubits. The initial spin coherent state evolves to spin-squeezed states on short time scales, and to “multiple-Schrödinger-cat” states (superpositions of many spin-coherent states) on longer time scales, a manifestation of the phenomenon of fractional revivals of the initial state. The evolution of the system is visualized with phase-space plots ( functions) that, when plotted against time, reveal a “quantum carpet” pattern. Of particular interest is the fact that our approximation captures the qualitative features of the model even for small values of . This suggests the possibility of observing the phenomenon of fractional revival in this model for systems of few qubits.

pacs:
03.67.-a, 05.50.+q, 42.50.Dv, 42.50.Md

I Introduction

Collapse and revival is a well known feature of quantum systems whereby the initial wave packet of the system collapses as it evolves, but at a later time returns (either exactly or approximately) to the initial state: the revival. An example is the evolution of a field mode in a Kerr medium with Hamiltonian and initially in a coherent state. After the revival time , the field mode is again in the initial state Yurke and Stoler (1986). Fractional revival is an additional effect where at a rational fraction of the revival time the state of the system is made up of a number of superposed, displaced copies of the initial coherent state. For an initial coherent state evolving in a Kerr medium, for example, the field mode will be in a superposition of coherent states (a “multiple-Schrödinger-cat state”) at a time where is an even number Gantsog and Tanas (1991); Tara et al. (1993).

Another example of collapse and revival is in the Jaynes-Cummings model for the resonant interaction of a field mode and a two-level atom Gerry and Knight (2005). The phrase “collapse and revival” in this context refers to the collapse and revival of the atomic inversion, but for any initial atom state, the field (initially in a coherent state) returns (approximately) to that coherent state after the revival time. This model also seems to have a limited form of fractional revival: For a judiciously chosen initial atom state, at a quarter of the revival time the field mode is (approximately) in a superposition of two coherent states Gea-Banacloche (1991); Bužek et al. (1992). However, this Schrödinger-cat state generation is really due to the conditional evolution of the field rather than fractional revival. This is because there are two orthogonal initial atom states that result in two different effective Hamiltonians, , for the field evolution Gea-Banacloche (1991). Starting in a superposition of these two atom states leads to Schrödinger-cat states of the field. At no time is the field composed of more than two distinct macroscopic states.

To see fractional revival and multiple-Schrödinger-cat states (with more than two components) in the Jaynes-Cummings model requires sub-Poissonian number statistics for the initial state Averbukh (1992); Góra and Jedrzejek (1993) (i.e., a non-classical, number-squeezed initial state). In this case one finds the usual collapse and revival, but also “super-revivals” at longer times, and multiple Schrödinger-cat states at rational fractions of the “super-revival” time Góra and Jedrzejek (1993). These kinds of fractional revival are again different to the Kerr-type fractional revivals since they come after the first revival. Kerr-type fractional revivals, on the other hand, appear before the first revival.

The phenomenon of fractional revival is well known and well investigated, both theoretically Averbukh and Perelman (1989); Robinett (2004) and experimentally Greiner et al. (2002), for systems with infinite-dimensional Hilbert space, including a recent experimental demonstration with a Kerr non-linearity Kirchmair et al. (2013). There has, however, been less discussion of fractional revival in finite dimensional systems (e.g., a system of qubits), although Agarwal and co-workers Agarwal et al. (1997) and Chumakov and co-workers Chumakov et al. (1999) have studied the evolution of an qubit system in a “finite Kerr medium” analogous to the Kerr Hamiltonian above for the field mode. This type of evolution is also the basis of several proposals to generate Schrödinger-cat states of various finite-dimensional systems Gerry (1998); Ferrini et al. (2008).

In this paper we study the interaction of qubits with a single qubit by an analogue of the Jaynes-Cummings Hamiltonian, assuming throughout an initial spin-coherent state (a “classical” state) of the qubits. In Dooley et al. (2013) it was shown that this system exhibits Jaynes-Cummings-like collapse and revival for a particular class of initial spin coherent states of the qubits (even for moderate values of ). Here we show that the same system with a different class of initial spin coherent states also exhibits fractional revivals. These fractional revivals are of the Kerr-type, and are comparable to these found by Agarwal and coworkers and Chumakov and coworkers rather than the Jaynes-Cummings type that requires a sub-Poissonian initial state. The transition from one regime of collapse and revival (Jaynes-Cummings like) to the other (with Kerr-like fractional revival) is made by changing the initial spin-coherent state parameter, something that is in principle very straightforward in the state preparation. We suggest that fractional revivals and the associated multiple-Schrödinger-cat states could be observed in this model even for few-qubit systems.

Although we discuss the system of qubits without reference to a particular physical system, there are various candidates for the implementation of the model. For example, as we discuss in the conclusion, the qubits may be superconducting qubits, nuclear spins in symmetric molecules, or spins associated with Nitrogen vacancy (NV) centers in diamond.

Ii Fractional Revivals in Systems of Qubits

We begin by defining the operators for an qubit system as where is a Pauli -operator for the ’th qubit and . We also define and . Dicke states are simultaneous eigenstates of and with eigenvalues and , respectively. In what follows we restrict to the eigenspace of the -qubit system. This is the symmetric subspace of the -qubit system: states in this subspace are invariant under permutation of qubits. Alternatively, the system restricted to this subspace can be thought of as a spin- particle. It is an -dimensional subspace for which the Dicke states () form a basis. Spin-coherent states Arecchi et al. (1972) in this subspace are separable states of the qubits with each qubit in the same pure state:

(1)

where is a complex number. In the Dicke-state basis the spin-coherent state can be written as

(2)

where

(3)

By making the transformation we can also express (1) in terms of the usual polar and azimuthal Bloch sphere angles and as follows:

(4)

We note that the spin-coherent states are widely regarded as “classical” states of an spin system Markham and Vedral (2003); Giraud et al. (2008) and that the preparation of any spin coherent state from an initial spin coherent state (e.g., the ground state ) is, in principle, straightforward: the arbitrary spin-coherent state can be generated by an appropriate interaction with an external classical field Arecchi et al. (1972).

ii.1 Interaction of qubits with one qubit

We say that a system of qubits interacts with a single qubit via Hamiltonian

(5)

This Hamiltonian was studied in Hutton and Bose (2004); Breuer et al. (2004); El-Orany and Abdalla (2011); Wu et al. (2014); Dooley et al. (2013) and has been referred to as a “spin star” system.

Given any initial state of this system we can find an exact expression for the state at any later time . This expression is, however, cumbersome and difficult to interpret in its exact form. Here, for convenience, we focus on the resonance () and initial states that are of the form

(6)

where the -qubit system is in a spin-coherent state and the single qubit is initially in an arbitrary pure state, written here in terms of the orthonormal basis states that depend on the phase of the spin-coherent state parameter . We consider separately the evolution of the two orthonormal states by our Hamiltonian since the evolution of an arbitrary initial state (6) is just a superposition of these two solutions.

On resonance it is convenient to first rotate Hamiltonian (5) to an interaction picture where the Hamiltonian is (the time independence of this interaction Hamiltonian depends on the resonance condition). The unitary time evolution operator is thus . Expanding the exponential as a Taylor series, and using the identities , and , gives

(7)

Depending on whether the initial state of the qubit is or we write this unitary operator as or , i.e., with . From (7) we find that

(8)

This is still an exact expression that does not give much insight into the features of the system as it evolves in time. However, we can use two approximations that allow us to write the time-evolved state in a useful way.

First approximation. We say that the spin-coherent state is approximately an eigenstate of the operator with complex eigenvalue and an eigenstate of with eigenvalue .

This is a good approximation when (details given in Appendix Appendix A):

(9)

For large this restriction is not at all severe since a very broad range of values of will satisfy condition (9).

We plot in Fig. 1 the quantities and where . The fact that both of these quantities are small when justifies our approximation.

We note that a similar approximation can be made for a field mode where for a coherent state with Loudon (1973).

Figure 1: (Color online). Top: plotted against and for . Bottom: plotted against and for . The black lines show for (or ) and for (or ).

Using the identities and then allows us to write (8) as

(10)

In other words, the effective Hamiltonian given the initial qubit state is

(11)

Unlike Hamiltonian this effective Hamiltonian is diagonal in the Dicke-state basis of the qubits, in the sense that .

The fact that we can find an effective Hamiltonian given the initial qubit state or is reminiscent of the effective Hamiltonian that can be derived in the strong field approximation of the resonant Jaynes-Cummings model: when the atom is in the state , a “semi-classical eigenstate,” the field mode evolves by the effective Hamiltonian Gea-Banacloche (1991); Chumakov et al. (1994). Indeed, in the parameter regime the above model exhibits Jaynes-Cummings-like collapse and revival, as discussed in Dooley et al. (2013). Making the transformations and and taking , Hamiltonian (5) reduces exactly to the Jaynes-Cummings model Holstein and Primakoff (1940) and the initial spin-coherent state (2) becomes a coherent state Markham and Vedral (2003); Dooley et al. (2013).

Second approximation. Expanding the square root in (11) in powers of the operator we obtain

(12)

where , and for are coefficients whose absolute value is always less than unity. [Here is a double factorial, i.e., the product of all odd positive integers less than or equal to .] The second approximation is to truncate this expression to the first two terms:

(13)

This approximation is valid when (or ) and (details given in Appendix Appendix B). In Fig. 2 we plot , the fidelity of the exact state to the approximation against time for various initial values of and . The time axis is scaled by for reasons that should become clear in the next section.

Figure 2: (colour online). Top: Fidelity of approximation to the exact state plotted against and with and . Bottom: Fidelity plotted against and with and . For (or ) and (or ) the fidelity is high (red).

The truncated effective Hamiltonian (13) can be conveniently rewritten as

(14)

If the initial state is either or (rather than some superposition of the two) then (since both and are eigenstates of ) the first two terms of (14) just give a global phase factor that can be ignored. If the initial state is a superposition of and then the first two terms cannot be ignored since they give a relative phase factor. The last term of (14) is proportional to . For an -qubit system the Hamiltonian is known as a one-axis twisting Hamiltonian Kitagawa and Ueda (1993); Ma et al. (2011) or, because it is analogous to the Kerr Hamiltonian for a field mode – a finite Kerr Hamiltonian Chumakov et al. (1999). Our term proportional to in the effective Hamiltonian (14) is therefore a one-axis twisting term for the -qubit system.

Previous examples of finite-dimensional systems whose Hamiltonians include one-axis twisting terms are a collection of two-level atoms interacting with a far detuned field mode Agarwal et al. (1997); Klimov and Saavedra (1998); Bose-Einstein condensates in a double-well potential Milburn et al. (1997); molecular nano-magnets Wernsdorfer (2008); a collection of NV centers coupled to the vibrational mode of a diamond resonator Bennett et al. (2013). Here we have shown that the one-axis-twisting Hamiltonian can also be an effective Hamiltonian in the interaction of qubits with a single qubit. We note, however, that the coupling parameter is weakened by a factor of for the one-axis-twisting term in .

ii.2 Multiple Cat States and Quantum Carpets

In this subsection (following the method of analysis of Averbukh and Perelman (1989); Robinett (2004)) we investigate the evolution of the initial state under the approximations made in the previous section. We note that the initial spin coherent state has (or, equivalently, ) so that all of our approximation conditions are satisfied so long as and .

The system evolves by the effective truncated Hamiltonian so that the evolved state is

(15)

where the tilde above indicates approximation and where we define

(16)
(17)

So that we can eventually arrive at a more useful expression for we now consider some properties of these functions and . First, when is an even number both and are periodic in time with period . To see this note that when is even both and are even integers so that and are exponentials whose phases are integer multiples of . Since and are periodic in time we have : the system returns to its initial state after period . Similarly, when is an odd number we have and the state has a revival time of .

Focussing on the case when is even, we now consider times where and are coprime integers, i.e., rational fractions of the revival time. Then and are both either periodic functions of the discrete variable with period if is odd, or anti-periodic functions of with (anti-) period if is even:

(18)
(19)

Either way, this means that we can write and in terms of their discrete Fourier transforms,

(20)

where we define

(21)

The inverse transform is

(22)

The advantage of writing the functions in terms of their discrete Fourier transforms is that the phases in the exponentials and are now linear in rather than quadratic in . Substituting (16) and (17) into (20) and using the fact that and are periodic in , it is not difficult to show that

(23)

Substituting (22) and (23) into our expression (15) allows us to write the state appealingly as:

(24)

This is a superposition of terms involving spin-coherent states of the -qubit system distributed uniformly in the azimuthal Bloch sphere angle . Expression (24) indicates that the system undergoes fractional revivals at times since it shows that the state at that time is a superposition of displaced copies of the initial wave packet. This is seen most clearly by rewriting (24) as where

(25)

is a unitary operator that is a sum of unitary displacements .

The explicit value of is

(26)

Sums of the kind in (26) are known as generalized Gauss sums Bernt and Evans (1981) and can be calculated for various values of and . For simplicity we take . In this case Bernt and Evans (1981)

(27)
(28)

The initial state is a spin-coherent state of the combined -qubit system:

(29)

In this case (ignoring global phase factors) the evolved state takes a particularly straightforward form:

(30)
(31)

This is a superposition of spin coherent states of the qubit system uniformly spaced around the equator of the Bloch sphere. This is in agreement with the results of the authors of Agarwal et al. (1997) and Chumakov et al. (1999) for the evolution of a spin-coherent state by a finite Kerr Hamiltonian. Taking , for example, is a Greenberger-Horne-Zeilinger (GHZ) state of the qubits. To visualize such states we plot in Fig. 3 the Husimi function

(32)

of the exact state at various times for and for the initial state . For very short times the initial spin-coherent state evolves to a spin-squeezed state Kitagawa and Ueda (1993); Ma et al. (2011), as shown in Fig. 3(b). At later times we see multiple-Schrödinger-cat states [Figs. 3(c),(d),(e)]. We note that Schrödinger-cat states can also be generated in a different parameter regime () for the initial state, as discussed in Dooley et al. (2013).

Figure 3: (Color online). functions [Eq. (32)] and their stereographic projections at various times for the (exact) state with . (a) The initial state , a spin-coherent state. (b) The state . (c) The state . (d) The state , a GHZ state. (e) The state . (f) The state at the revival time.

Since the operator commutes with our Hamiltonian (5), we know that it (and its powers) are conserved quantities of the system. This means that if the function is initially a narrow distribution at the equator of the sphere, the function is constrained to the equator at all times, as seen in Fig. 3. This suggests a concise visualization of the system dynamics by plotting as a function of time and of (ignoring the variation in the polar angle that plays a less interesting role). The resulting pattern is shown in Fig. 4. Such patterns are know as “quantum carpets” Kaplan et al. (2000); Berry et al. (2001). At times that are rational fractions of the period we see bright spots at the values of where there are spin coherent states.

Figure 4: (Color online). Plotted is , the -function slice at . Here . We see a “quantum carpet.”

ii.3 Small

Our approximations require that (or ), and (or, ). As illustrated in Fig. 2, the fidelity of the exact state to the approximate state gets worse as gets small or as gets close to . Our -function plots show, however, that the qualitative features of the approximation are valid well outside of these parameter regimes.

In Fig. 4, for example, it is clear that the multiple-Schrödinger-cat states persist well beyond . In Figs. 3(e) and 3(f) we plot the functions for at and , respectively. In Fig. 3(e) we see something qualitatively like a superposition of spin-coherent states, although the coherent states are distorted [the likely cause for the decrease in fidelity against the ideal superposition of spin-coherent states [Fig. 2)].

Similarly, although our approximation required that we assume , the functions in Fig. 5 show that our approximation captures the qualitative features of the exact evolution of the system for moderately small values of . It is clear from these plots that, although they are not superpositions of perfect spin-coherent states, they are superpositions of distorted spin-coherent states, still highly non-classical states. In Fig. 6 we plot the “slice” as a function of and of for . The “carpet” pattern, although not as sharp as in Fig. 4, is clearly recognizable. The states at and are plotted in Fig. 5(c) and 5(d), respectively.

Figure 5: (Color online). functions (and their stereographic projections) at various times and for small values of . Each time the initial state is . (a) for . (b) for . (c) for ; (d) for .
Figure 6: (Color online). Plotted is for . Even for this small value of the “carpet” pattern is conspicuous.

Iii Conclusion

In Dooley et al. (2013) is was shown that a system of qubits interacting with a single qubit via Hamiltonian (5) leads to collapse and revival analogous to that in the resonant Jaynes-Cummings model when the initial state of the qubits is a spin-coherent state with . Here we have shown that for an initial state in a different parameter regime, , the whole -qubit system undergoes a different type of collapse and revival that includes fractional revivals: at rational fractions of the revival time the system is (approximately) in a multiple-Schrödinger-cat state. The transition between the two different kinds of collapse and revival is made by changing the spin-coherent state parameter , or, equivalently, by changing , the Bloch sphere polar angle for the initial spin-coherent state.

The fractional revivals in the present case are similar to those found for the evolution of an -qubit system by a finite Kerr Hamiltonian, the analogue of a Kerr Hamiltonian for an optical system Agarwal et al. (1997); Chumakov et al. (1999).

We emphasize that the initial spin-coherent state is, in principle, easily prepared since it consists of all qubits in the same pure state. By comparison, fractional revivals in the Jaynes-Cummings model requires an initial state with a sub-Poissonian photon number distribution Averbukh (1992); Góra and Jedrzejek (1993).

We also point out that the qualitative features of this approximation are valid even for moderately small values of . An interesting application of this result might be the observation of fractional revivals in systems of few-qubits. Physical systems where this might be realized are interactions between superconducting qubits McDermott et al. (2005); Niskanen et al. (2007); Neeley et al. (2010) for which the coupling between qubits is natural. Our interaction Hamiltonian (5) is composed of such equal interactions with the central qubit. Alternatively, a superconducting phase qudit Neeley et al. (2009) (which emulates our dimensional subspace) might be coupled to a single superconducting qubit.

Another candidate system is any highly symmetric molecule that consist of spins equally coupled to a central spin. The trimethyl phosphite molecule, for example, has 9 spins, all equally coupled to a single spin Jones et al. (2009). The tetramethylsilane molecule has 12 spins equally coupled to a single spin Simmons et al. (2010). These values of may be large enough to allow for the observation of fractional revivals.

An example of a system where Hamiltonian (5) may be realized for larger values of is the interaction of a superconducting flux qubit with a thin layer of NV centers Marcos et al. (2010). Rabi oscillations for such an interaction have been observed experimentally for Zhu et al. (2011). Squeezed states or Schrödinger-cat states of the NV centers generated in this system may have applications in magnetic field sensing.

Future work will extend the model to qubits interacting with two qubits (an exactly solvable model), both from the perspective of the evolution of the qubits, and from the perspective of the decoherence and disentanglement of the two qubit system as it interacts with the qubits. Future work will also include a study of various forms of realistic decoherence to assess the possibility of experimentally observing fractional revivals, especially in systems of few qubits.

References

Appendix A

Here we justify the approximations

(33)
(34)

when

(35)

Writing the spin-coherent state in its Dicke basis gives

(36)

where

(37)

The operators and act on the Dicke state as follows:

(38)
(39)

so that we can write

(40)
(41)

We expand the expressions in the square brackets in (40) around the average values of :