Engineering SU(1, 1) ⊗ SU(1, 1) vibrational states

We propose an ideal scheme for preparing vibrational SU(1, 1) ⊗ SU(1, 1) states in a two-dimensional ion trap using red and blue second sideband resolved driving of two orthogonal vibrational modes. Symmetric and asymmetric driving provide two regimes to realize quantum state engineering of the vibrational modes. In one regime, we show that time evolution synthesizes so-called SU(1, 1) Perelomov coherent states, that is separable squeezed states and their superposition too. The other regime allows engineering of lossless 50/50 SU(2) beam splitter states that are entangled states. These ideal dynamics are reversible, thus, the non-classical and entangled states produced by our schemes might be used as resources for interferometry.

Quantum state engineering studies the preparation, manipulation, and characterization of arbitrary quantum states. Technological advances allow the coherent control of dynamics in an increasing collection of physical systems; e.g. trapped ions, superconducting circuits, quantum gases, mechanical oscillators. In particular, trapped ions show high addressability, long coherence times and high fidelity readout necessary for quantum state preparation and manipulation within its own experimental issues 1 . Vibrational state characterization is available in this platform through tomographic reconstruction, experimentally demonstrated for Wigner 2 and Hussimi 3 quasi-probability distributions.
Trapped ions have proved a fertile ground for fundamental research and quantum technologies development 4,5 . Single-mode vibrational number, coherent, and squeezed states have been engineered experimentally 6 and theoretical proposals for the synthesis of arbitrary one 7 -and two-dimensional vibrational states 8 has been produced. In particular, trapped ions might act as vibrational beam splitters 9 producing states identical to the photon states on the output ports of a lossless interferometer with number-state inputs 10 . Vibrational interferometry can be used to either explore fundamental quantum mechanics, e.g. quantum decoherence 11,12 , or produce new quantum technologies, e.g. vibrational thermometers 13 or quantum gyroscopes 14 .
Squeezing and entanglement improve phase sensitivity in interferometry in a manner proportional to the inverse of the excitation number of quanta entering an interferometer [15][16][17] . Here, we are interested in the quantum state engineering of orthogonal vibrational modes that show squeezing and entanglement with an underlying SU(1, 1) ⊗ SU(1, 1) symmetry. We will use blue and red resolved second sideband driving 18,19 to this end. In the following, we will present an effective Hamiltonian describing our proposal in the Lamb-Dicke regime. Then, we will show that the asymmetric coupling model produces the superposition of separable squeezed vibrational mode states where the inner state is intrinsically entangled to the vibrational modes. Afterward, we will show that red driving with symmetric coupling is able to produce lossless 50/50 SU(2) beam splitter vibrational states that are factorized from the internal state of the ion. The ideal dynamics producing these states are reversible and seem to suggest their use as interferometers to characterize different aspects of real-world experiments.

Results
We suggest driving two normal phonon modes of the center of mass motion of a trapped ion 20,21 , such that the effective Hamiltonian, www.nature.com/scientificreports www.nature.com/scientificreports/ describes the interaction of the j-th vibrational mode, with frequency ν j and represented by the annihilation (creation) operator â j ˆ † a ( ) j , with two internal levels of the trapped ion, with energy gap ω 0 and represented by Pauli matrices σ j , through a set of external driving fields of frequency ω j , Lamb-Dicke parameter η j , phase φ j , and Rabi coupling strength Ω j . Moving into the reference frame defined by the uncoupled Hamiltonian, ω σ ν ν = + +ˆˆˆˆ † † H a a a a /2 0 0 3 11 1 2 2 2 , driving the k-th vibrational sideband, ω j = ω 0 − kν j where k > 0 and k < 0 define the so-called red and blue driving, in the Lamb-Dicke regime, η 〈 〉  † a a 1 j j j , and under an optical and mechanical rotating wave approximation, we can approximate 22 for red, and blue sideband driving, with effective couplings and phases, , in that order. For reasons that will become obvious, we choose the second sideband driving 19 , k = ±2, and draw from the idea of simultaneous blue and red driving in the simulation of the quantum Rabi model 23,24 to reach the model Hamiltonian at the core of our proposal, This model is the two-phonon interaction analog of the so-called cross-cavity quantum Rabi model 25 that has been used to propose the quantum simulation of para-oscillators in trapped ions 26,27 . Single-mode, two-phonon interactions in trapped ions have been recently proposed to simulate interaction-induced spectral collapse 28 . It may be possible to explore alternative schemes in the trapped-ion platform; for example, continuous dynamical decoupling schemes that have been proposed as alternatives to produce robust realizations of the two-phonon interaction in the ultrastrong regime 29 or dynamics far away from the Lamb-Dicke regime 30 . It might be even possible to explore realizations with superconducting qubits, where a two-photon quantum Rabi model has been proposed 31 , or single trapped cold atoms, where a proposal to realize the quantum Rabi model has arisen 32 .

Discussion
Asymmetric squeezing. Let us simplify and find uses for our general Hamiltonian. First, we propose to work with a system where the amplitudes and relative phases for blue and red driving in each mode are chosen to provide similar coupling parameters, g j,k ≡ g j , and phases, φ j,k ≡ π/2. Under this assumption, the Hamiltonian describing the system, yields an evolution operator, controlled by the internal state of the trapped ion. For example, choosing an initial state with arbitrary phonon fields and the internal state an eigenstate of the σ 1 operator, |ψ(0)〉 = |ξ 1 , ξ 2 , x ± 〉, the ideal evolution, provides a separable state with different squeezing in each mode. These are the so-called two-mode SU(1, 1) Perelomov coherent states, whose properties have been discussed in detail by Gerry and Benmoussa 33 . Figure 1(a) shows the fidelity defined as the trace distance, where γ q and γ j are the effective atomic and vibrational decay rates. The parameters in the simulations are g 2 = 0.75 g 1 and, for the sake of the example, homogeneous decays γ q = γ j = γ = 0.1 g 1 . While we select these parameter values for the sake of showing the dynamics in a lossy system, they are not far from those of trapped Barium 34 or Cadmium 35 where Rabi frequencies are of the order of 15-100 kHz and the coherence time of the qubits are in the 80-120 μs range. Note that Fig. 1(b-d) show the joint phonon number probability, n m , www.nature.com/scientificreports www.nature.com/scientificreports/ for ideal and Fig. 1(e-g) for lossy evolution at different times. We colored the probability bars and show only a small section of the plot to make comparison simpler.
Obviously, choosing an arbitrary internal initial state, |ψ(0)〉 = |ξ 1 , ξ 2 〉(α|x + 〉 + β|x − 〉), provides a superposition of the form, where the state of the whole system is entangled thanks to the internal state of the ion. Figure 2(a) shows the fidelity between the ideal and lossy evolution. Figure 2(b-d) show the joint phonon number probability, for ideal and Fig. 2(e-g) for lossy evolution. Parameters are identical to those in Fig. 1. The fact that a π-phase change in the driving laser phases ideally reverts the effective model dynamics might suggest its use as a squeezed state interferometer. In a real-world experiment, the system will not return to the original initial state due to all the real-world experimental subtleties. For example, these states might help in the characterization of the environment effect on the trapped center of mass motion and internal states of the ion might mention the obvious. Symmetric squeezing. Now, let us simplify our model one step more and assume a system where the coupling parameters, g j,k ≡ g, and the phases are chosen to be equal. The dynamics are given by a simpler Hamiltonian that produces symmetric squeezing if we follow the procedure introduced before. We can take a second step and suppress blue sideband driving, then the dynamics are described by a Hamiltonian, Again, it is straightforward to construct an evolution operator,  www.nature.com/scientificreports www.nature.com/scientificreports/ 2 using a representation of SU(1, 1), 3 3 provided by the two-mode operators 33 , We will use a Hilbert space partition defined by the raising operator, m acting on four states that we will call vacuum states related to a Bargmann index k. This produces four phonon subspaces labelled by a Bargmann index and a vacuum state: k = 1/2 and |1/2; 0〉 = |0, 0〉, k = 1 and |1; 0〉 01 = |0, 1〉, k = 1 and |1; 0〉 10 = |1, 0〉, k = 3/2 and |1, 0〉 = |1, 1〉. Figure 3 shows a pictorial representation of these phonon subspaces. The whole Hilbert state for the quantized center of mass motion is covered once with the four orthogonal subspaces defined by the bases,  Here, we want to make a stop and remember the action of a lossless SU(2) beam splitter 10,36 , such that for a 50/50 beam splitter, θ = π/2, we can rewrite the SU(1, 1) bases above, and so on. Now, this procedure generates entangled orthogonal vibrational states that are factorized from the internal structure of the ion. Again, this result suggests the use of this system as an interferometer that might provide information about the characteristics of an experiment that might be, in principle, different from those available through the superposition of squeezed states proposed before. Figure 4(a) shows the fidelity, t ( )  , between the ideal evolution of the driving algorithm in Eq. (20) and lossy evolution under homogeneous decay for all components. The initial state of the algorithm is the vacuum state for the subspace defined by the Bargmann parameter k = 1/2. That is, the initial state is the ion in the excited state and both vibrational modes cooled down to the vacuum state. Figure 4(b-d) show the joint phonon number probability, P n,m , at times where the 50/50 SU(2) beam splitter states are expected under ideal time evolution. Figure 4(e-g) show the same probabilities for lossy evolution. www.nature.com/scientificreports www.nature.com/scientificreports/

Conclusion
We proposed a trapped ion model under second sideband resolved blue and red driving of two orthogonal modes of the center mass motion in the Lamb-Dicke regime. For parameter regimes providing an effective model with asymmetric coupling of the vibration modes with the internal state of the ion, we showed that time evolution of arbitrary vibrational states with balanced superposition of the internal states of the ion generates a so-called two-mode SU(1, 1) Perelomov coherent state where the vibrational modes are separable. Uneven superposition of the internal states of the ion produces a superposition where entanglement of the whole state of the ion arises. We also showed that under red driving only, the effective model is able to generate 50/50 SU(2) beam splitter vibrational states factorized from the internal state of the ion. The fact that it is ideally possible to engineer these states with reversible dynamics suggests the use of these models as interferometers to characterize real-world experiments.