The classical and quantum synchronization between two scattering modes in Bose-Einstein condensates generated by the standing-wave laser beams

The classical and quantum synchronization between two scattering modes in Bose-Einstein condensates generated by the standing-wave laser beams

Lin Zhang    Xiaoting Xu    Xing Liu    Weiping Zhang School of physics and Information Technology, Shaanxi Normal University, Xi’an 710119, China Department of Physics and Astronomy, Shanghai Jiao Tong University, Shanghai 200240, P. R. China Quantum Institute for Light and Atoms, Department of Physics, East China Normal University, No.500, Shanghai 200241, P. R. China Collaborative Innovation Center of Extreme Optics, Shanxi University, Taiyuan, Shanxi 030006, P. R. China

Both classical and quantum dynamics of the synchronization between two nonlinear mechanical modes scattered from Bose-Einstein condensates (BECs) by the standing-wave laser beam are comparatively investigated. As the ultra-low dissipations of the momentum modes in the atomic BECs, the synchronized dynamics are studied in a framework of closed-system theory in order to track down both the classical and the quantum synchronizations from an angle of quantum control. The classical synchronization and the relevant dynamics of measure synchronization, the quantum synchronization and two different types of measures proposed by Mari and estimated by mutual information based on -function are studied respectively in order to reveal both the macroscopic and the microscopic signatures of synchronized behaviors in a closed quantum system. The results demonstrate that the “revival and collapse” of the quantum fluctuations beyond the classical mean-value dynamics due to long-lasting mode coherence discriminates the quantum synchronization from the classical one, which not only excludes the possibilities of an exact synchronization and a perfect density overlap in phase space, but also leads to upper limitations to Mari measure and large unceasing fluctuations to mutual information between two scattering modes. We reveal a close dynamic connection between Mari measure and the mutual information of two nonlinear momentum modes in closed BEC systems by demonstrating an opposite mean-value behavior but a similar fluctuation variation with respect to their respective evolutionary scales.

42.50.Pq, 37.10.Vz, 37.30.+i, 42.65.Pc

I Introduction

In recent years, many attentions are paid to the quantum control problems for exploring classically driven devices which are actually working in a quantum framework Brif (); Wiseman (). For the quantum control on a collection of nearly identical microscopic quantum units from a macroscopic level down to atomic/molecular scale, only the synchronized behavior emerged from the microscopic scale can be identified by the classical devices as a macroscopic signature, and thus the synchronization behaviors at the macroscopic scale, that is, the classical synchronization (CS), play a critical role to trace down (classically detect or control) the coherent quantum dynamics among the microscopic quantum units. However, the synchronized dynamics between different quantum units, which induce dynamical transitions of macroscopic temporal-spacial orders, is intrinsically determined by the microscopic synchronization within a decoherence time for a quantum evolution. Therefore, the quantum synchronization (QS) as another important signature for the quantum coherent dynamics have been extensively studied with the nanomechanical resonators Bagheri (); Shim (); Ying (), the spins and the atomic qubits Orth (); Hush (). The dynamics of the coupled oscillators leading to QS have also been extensively considered in the optomechanical systems Holmes (); Shlomi (); Heinrich (); Zhang ().

In this paper, we want to push this topic by exploring both classical and quantum natures of the synchronization between two nonlinear collective modes in a macroscopic quantum matter: Bose-Einstein condensate (BEC). In BEC, the synchronized behaviors between a large number of atoms are complicated because the dimension of BEC is huge, but the main features of synchronization can also be revealed by only analyzing the dynamics between several collective modes coherently excited from it, especially by the simplest case of a two-mode excitation. In a two-mode quantum system, one mode can be treated as a passive oscillator which is easily initialized to a ground state and the other is the active one which produces quantum control actions via their mutual couplings. In this driven-response closed loops, the dynamical measures of two quantum oscillators to be synchronized within a decoherence time are the useful signatures to guide efficient controls on the quantum evolution of the passive oscillator starting from a specific initial state. Recently, many measures are presented for the judgement of QS Mari (); Ameri () and most of the works considered two coupled nonlinear quantum oscillators exposed to different dissipative baths Goy (); Zhi (); Agrawal (); Lee (); Math (); Walter (); Gie (). These studies are mainly concerned about QS phenomena in the open quantum systems, and the results demonstrated that QS of an open quantum system is eventually determined by the statistical properties of the external baths Breuer (). Further studies indicate that two oscillators coupled to a common reservoir are more likely to be synchronized than that they separately couple to two different reservoirs Giorgi (). These dissipative synchronizations are all related to the stable steady states or the dynamical attractors to which the coupled systems can eventually settle down due to the dissipation and the decoherence induced by the baths Xue (). However, for a control problem, a full real-time evolution of the correlated states between coupled units within a decoherence time is more important than the longtime asymptotic behavior. Therefore the damping and the couplings with the baths can be neglected in a short time interval which enables an assumption of a closed quantum system. As the conservation of the phase volume for a closed quantum system, the synchronization without dissipations will be very different from that of an open quantum system. In this paper, we will consider the dynamical synchronization between two autonomous nonlinear modes generated in the atomic BECs Weiping () by investigating the full synchronized behaviors in both the classical and the quantum regimes.

As discussed above, a complete synchronization for a closed quantum system must include both the macroscopic and the microscopic signatures of the synchronized dynamics, that is, both the CS and QS. In the classical regime, the synchronization phenomena for a dissipative system have been studied for nearly half a centaury from the coupled Huygens clock (in 1665) to chemical reactions, biological and social behaviors Janson (); Pikovsky (). Many CS measures have been proposed in differen ways, such as complete synchronization, generalized synchronization, phase synchronization, lag synchronization, projective synchronization and so on Pikovsky (); Zhou (). However, from a microscopic point of view, how to estimate QS of quantum systems is still a problem. A natural way to define a measure for QS is to generalize the classical counterparts directly to the quantum regime. However, the generalization of the classical measures to quantum case is not always straightforward because not all the classical counterparts in the quantum regime can be uniquely found, for example, the phase synchronization is related to the Hermite phase operator which, obviously, is difficult to be defined in the quantum regime Barnett (). While for the coupled conservative Hamiltonian systems, an interesting CS called measure synchronization (MS) was found Zanette () and, recently, its QS case was considered in the BEC system Qiu (). Up to now, two types of measure for QS are proposed in the literatures, one is based on the distance between two quantum states in the Hilbert space and the other relies on the quantum state correlations. However the distance of the quantum states in Hilbert space is not unique and different norms or inner products can be adopted for different QS definitions along this line. Moreover the QS measures defined by using state correlations are even more complicated because there exist many quantum correlations for two quantum states, such as quantum discord, quantum fidelity, quantum coherence, quantum entropy, and even different quantum entanglements Manzano (). Therefore, the measure to judge QS between two coupled quantum systems might not be theoretically unique and which one we should pick is only determined by what microscopic property we want to resolve. In this paper, we will investigate the dynamics of QS between two atomic scattering modes in BEC by using -function as a tool. As -function is positive and bounded, we can easily combine two types of QS measure together (the distance and the correlation) to explore the dynamical features of QS in a closed quantum system by investigating either the “distance” (overlap) from a classical aspect or the “correlation” based on the mutual information theory.

Ii The Theoretical model

Theoretically, we consider only two low-energy scattering modes and stimulated in BEC by a standing-wave laser beam. The coherent atomic BEC can be prepared either in a number state or in a coherent state (see Ref.Weiping () for the state preparation of BEC). The Hamiltonian describing two scattering modes within a BEC beam in a number state can be described by Weiping ()


where and are the mechanical frequencies of two momentum modes induced by the collective atomic recoil Lin (), is the linear coupling rate with being the Rabi frequency of the atom interacting with the standing-wave field and being the detuning of the laser field from the atomic transition frequency, and is the nonlinear coupling rate which is due to the interatomic dipole-dipole interaction (photon exchange) Weiping (); Smerzi (). For the scattering modes generated in a BEC of number state, the total atomic number is conserved in a sense of .

System (1) can be treated as two coupled nonlinear oscillators in Hilbert space of , where denotes two free Kerr-like nonlinear oscillators with free energy of in a number-state representation of . The first term in describes the mutual coupling of two nonlinear Boson oscillators which induces mode transitions and the second term of describes intensity interaction analogous to two light modes in Kerr medium [20]. Surely, the above model can also be obtained in the BEC system with a double-well potential with two-mode approximation (left-well and right-well modes) Smerzi (). Usually, for a two-mode system, an interaction picture with respect to the middle energy of can be introduced as


where the free half-level gap between the two modes is defined by .

Surely, we can also investigate the scattering modes generated from BEC in a coherent state that is prepared by coherently mixing of several BEC beams with different atomic number of Weiping (). In this sense, the above quantum system (2) can be generalized into a two-mode system in a Hilbert space of . In the following sections we will consider both cases in the interaction picture by omitting the prime of for simplicity.

Iii Nonlinear dynamics of CS: measure synchronization

Firstly, we consider the macroscopic synchronized behavior of two scattering modes in a classical point of view and the Hamiltonian system (2) admits a well-known CS: measure synchronization (MS) Zanette (). The operator equations for the two momentum modes in the BEC beam with a total atomic number of will be determined by


where . Equivalently, we can write it by two parts, the linear and the nonlinear parts, as

If the nonlinear interatomic scattering coefficient (due to dipole-dipole interaction) , the Hamiltonian will reduce to a typical JC model with an effective Rabi frequency of . The nonlinear dipole-dipole coefficient will modify and by introducing a number-dependent Kerr interaction between two momentum modes, leading to an approximate Rabi frequency of .



Figure 1: The dynamics of the quadratures (solid lines), (dashed lines) and their corresponding errors ( and ) in the phase space under different nonlinear coupling rates (upper panel), (middle panel) and (lower panel). The mode detuning is and the initial conditions are and .

The classical dynamics of the two-coupled modes can be investigated by merely putting them into a coherent state of and just replacing with in Eq.(3) and Eq.(4), and then we have


where the normalized variables for and the parameters are rescaled by the total atomic number of . For a classical dynamics, the MS between two scattering modes beyond dissipation can be found in Fig.1 when the coupling rate enters into a certain parametric region. The right-hand column of Fig.1 shows the quadrature dynamics of , and the corresponding quadrature errors of


are shown in the left column in phase space. Fig.1 illustrates that the quadrature errors between modes and generally conduct a quasi-periodic motion and never can reach to a zero-measure orbit (a periodic circle or a constant point) in the phase space when the two modes start from different initial states, which indicates only a partial dynamical synchronization between two coupled modes without damping. As Eq.(5) and Eq.(6) are invariant for exchange of and , we can see a complete classical synchronization (when and ) only happens for starting from exactly the same initial states of and (symmetric case). The dynamics of the two quadratures of modes and perform a partial synchronization from local anti-phase motion to in-phase motion when the nonlinear coupling rate changes from to as shown in the right column of Fig.1. For some specific parameters and initial states, the error dynamics will follow strict periodic orbits in the phase space with a fixed relative phase between two mechanical modes as shown in the middle panel of Fig.1 or in Fig.2(g). We notice that the local anti-phase motion of two modes covers a larger ring area of error orbit than that for a in-phase motion in the phase space, which indicates a weaker energy exchanging between two modes of synchronized motions than that of the anti-synchronized motions. As the conservation of phase volume, the separated trajectories of two modes can not run into a same orbit due to the energy exchange between them through the mode coupling. Therefore, when two modes synchronized, two interwinding trajectories are formed and they will cover the same area in the phase space since the energy exchange meets the detailed equilibrium.

The impossibility of a complete synchronization in a conservative Hamiltonian system was first investigated by Hampton and Zanette Zanette () and they found the interesting MS of two coupled Hamiltonian systems. Since then MS in many non-dissipative systems have been extensively studied by many authors Wangwy (); Tianjing (); Vincent (). As the conservation of an overall phase volume, MS exhibits a typical trajectory overlap of two coupled oscillators in the phase space.

Figure 2: (a) The temporal dynamics of and and the corresponding trajectories in the phase space (inset) with the parameters of ; (b) The covering areas of modes (black dots) and (green circles) in phase space versus the coupling rate of under mode detuning of . Inset: A sample phase picture of the trajectories in a limited time interval with ; (c) The trajectory distance between modes and modified by the coupling rate of under mode detuning of ; (d)-(f) The long-time orbital overlap in the phase space modified by the nonlinear couplings of (d) , (e) and (f) with ; (g) An accidental periodic orbit of two modes with very specific parameters of . All the parameters are scaled by and the initial conditions are the same as in Fig.1.

By using the polar form of a complex variables as


we can easily check the trajectory overlap of two scattering modes by modifying the coupling rate as well as the mode detuning as shown in Fig.2. A typical dynamics of , and their longtime trajectories (inset) are given in Fig.2(a). The covering areas of the quasi-periodic trajectories of modes and in the phase space can be estimated by

where () means the maximum (minimum) value of during a time interval of . As shown in Fig.2(b), the two modes keep nearly a same phase area ( with differences below ) in a large parametric region of , displaying an invariant measure of two modes. We can also define a trajectory distance between two modes by calculating to see the trajectory overlap in the phase space as shown in Fig.2(c). If , the ring radius of trajectory is larger than that of mode and there is no overlap between them. When , the overlap begins and it reaches its maximum value at a critical point C shown in Fig.2(c). Fig.2(c) reveals that the overlap of the quasi-periodic trajectories in the phase space happens only in a very specific parametric window of and , indicating the resonant exchange of energy between two scattering modes. Fig.2(d)-(f) demonstrate the trajectory overlap with the increment of coupling rate between two symmetric modes of . We can find a clear evidence of orbital attraction in MS as shown in Fig.2(c) for that the orbital merging is more quickly than orbital separation with respect to the coupling rate of (the slope of dramatically decreases at the point of C in Fig.2(c)). For some specific initial states and parameters, the trajectories of modes and will run into the interesting synchronized periodic orbits (the correlated Lissajous curves with zero measures) such as shown by Fig.2(g). We can numerically verify that there are many correlated Lissajous orbits in the phase space determined by the initial conditions and the parameters of and . Therefore, according to the classical dynamics of Eq.(5) and Eq.(6), the two nonlinear scattering modes can obtain a partial synchronized dynamics in certain parametric windows in a sense of covering the same phase space of their trajectories with invariant measures or conducting correlated Lissajous orbits with zero measures.

Figure 3: The amplitudes of and of the two mechanical modes modified by and . (a) The density plot of the amplitude of in the plane of and ; The amplitude of changes with (b) or with (c) under the specific coupling or detuning indicated by the dashed line 1 and 2 in (a); (d) The density plot of amplitude of phase difference in the plane of and ; (e)(f) The amplitude of versus or for fixed or . All the frequencies are scalded by and the initial conditions are set to be the same as those in the previous figures.

As for a non-dissipative Hamiltonian system, the trajectory should be very sensitive to the initial conditions because the total phase volume of two coupled systems is determined by the initial states. The synchronized behavior of MS can exhibit a different characteristic if we consider the dynamics via intensity imbalance and relative phase between two modes starting from different initial states. Then the trajectory overlap means that two scattering modes can possess a correlated intensity relation or lock to a constant relative phase as implied by Fig.1 and Fig.2. The instaneous intensity imbalance and phase difference between two mechanical modes can be defined by


As a general quasi-periodic behavior of (see the inset of Fig.3(b)), the oscillating amplitude of changing with the detuning and is calculated in Fig.3(a). Fig.3(b) and Fig.3(c) are the oscillating amplitudes of with respect to and for the specific and , respectively, indicated by the dashed lines of 1 and 2 in Fig.3(a). Fig.3(a) reveals that the amplitude of the intensity imbalance can reach a constant value (AC behavior) and it will disappear for or at a large value of (DC behavior). Similarly, the amplitudes of the relative phase are also shown in the lower frames of Fig.3. Fig.3(d) presents the density plot of oscillating amplitude of in the parametric plane of and . Fig.3(e) and (f) display the amplitudes of changing with respect to parameters and indicated by two dashed lines of 1 and 2 in Fig.3(d), respectively. We can see that, in some parametric windows of and , the oscillating amplitude of is randomly dependent on the parameters. Fig.3(d)-(f) demonstrate that only in specific parametric windows of and can the phase amplitude obtain a stable value for a partial phase synchronization.

Figure 4: The evolutions of the intensity imbalance and phase difference in the phase space with (a) , (b) and (c) , starting from the random initial conditions. The mode coupling rate is and all the frequencies are scalded by . The thick black lines are the separatrices of the orbits and the black dots are the stable steady states.

Alternatively, the intensity imbalance or the phase error being an indicator of the MS Zanette () can be revealed by the equations of


where we introduce the amplitude ratio of as a convenient variable to describe the amplitude balance between two mechanical modes with a number conservation of . We can easily verify that and a fixed amplitude synchronization for only happens when or . In the case of , the solution for the relative phase will be , , where the initial vales and . Specifically, is a steady state of Eq.(10) and (11) under for . For , a quartic equation of should be solved to get the steady states in this case (the complete synchronization is related to the stable steady states corresponding the fixed points shown in Fig.4). Fig.4 displays the orbits in the phase space of with random initial and when the frequency detuning changes from negative to positive. The orbital flow in the phase space indicates that the amplitude synchronization or the phase synchronization only happens under specific initial conditions, such as or as shown in Fig.4 by the separatrices (thick black lines). The phase transitions from self-trapping case on mode in Fig.4(a) to the tunneling case in Fig.4(b) and to the self-trapping case on mode in Fig.4(c) are all connected by the asymptotic synchronization lines passing through the same dynamical frustrated points of and (). The trajectory flow shown in Fig.4 reveals that only a partial phase synchronization can easily be obtained for (Fig.4(b)) and the exact synchronized motion between two conservative modes is indeed an accidental case corresponding only to the stable steady states which can be shared by both modes.

As described in Ref.Zanette (), a complete CS that two different orbits of the system collapse into the same one is impossible in the coupled Hamiltonian systems because the phase volume must be preserved (no damping). However, a weak or partial synchronization can be established between two coupled conservative modes with a signature of two attractive interwinding orbits formed in phase space due to mutual exchange of energy. According to the above calculations on mean-value dynamics, we see that the partial two-mode synchronized motion can safely be acquired through the mode coupling under specific parametric conditions in a MS point of view. If the initial conditions are specially chosen, a full synchronization to reach the stable steady states can still be realized in this system classically. However, in order to consider the synchronization from a quantum aspect, we should give up the trajectory picture to investigate the probability density distribution in the phase space instead. Recently, a work studied the MS in a quantum regime to reveal some quantum properties of two coupled modes in BEC system Qiu (). As the fluctuations are always involved in the quantum dynamics, the CS measure by just investigating the overlap of the mean-value trajectories in the phase space is obviously too rough for a quantum dynamics.

Iv The coherent dynamics of quantum synchronization

iv.1 QS of scattering modes from BEC in a number state

Because of the intrinsic fluctuations in a quantum system, the QS is essentially different from the classical one and a complete synchronization of two coupled quantum modes is definitely avoided due to the uncertainty principle. In order to show the dynamics of QS in a closed quantum system, the total state of two excited modes in BEC should follow the Schrödinger equation of


where the wave function in the Hilbert space of can be expanded by


whose form is due to the conservation of ( denotes the total atomic number of BEC in number state and thus the dimension of is ) and the unitary evolution keeps . Substituting Eq.(13) into Eq.(12), we have the equation of motion for the coefficients as


Clearly, the above coupled equations produce a symmetric tridiagonal matrix of Hamiltonian in the Fock representation and the algorithm to calculate its eigenvalues are well-known (e.g., bisection algorithm bisect ()). Fig.5 gives a simple case of energy levels for BEC with only atoms () modified by the nonlinear couplings of . There are four energy levels in this case whose eigenstates are the superposition of four bare number states of , , and .

Figure 5: The energy levels of the BEC modes of atoms changing with different nonlinear couplings of (a) , (b) , (c) , and with (d) .

Clearly, the population distributions of modes and are dynamically determined by the avoided level crossings of the states which can easily be modified by and . A clear degenerate case of shown in Fig.5(b) indicates a transitionless evolution between four bare states without energy avoided crossings for


This corresponds to the classical case of for a synchronized intensity balance (see Eq.(10) for ) that only the relative phase is linearly increasing just like that in Eq.(15).

Actually, an overlap of the quantum dynamics of modes and in the phase-space should be investigated by using the probability distributions of the total quantum sates. The -function based on the coherent state representation possesses a good behavior to provide positive and bounded probability distributions in the phase space. The coherent properties of a quantum state can be revealed by a continuous function which is a projection of the quantum state onto a reference coherent state by , and its probability density is called -function as

If the -function goes to 1 in the phase space, then the state will be completely described by a coherent state with an overall amplitude and phase. Therefore, the quantum MS of two quantum states can then be traced by the overlap of probability distributions in the phase space revealed by -function. The -function for the two-mode state of BEC in a number state is


where the density operator . In terms of , the above -function can be written as


By using Cauchy’s inequality for a certain , the above -function has an upper limit function of


where the normalization condition in the Hilbert space is used. Clearly, the above bounded -function gives a central symmetric distribution for modes and and no phase information is included.

Figure 6: The temporal density distributions of the -functions for two modes with specific (solid lines) or with specific (dashed lines) at different evolution times . The parameters are and . The time unit of the evolution is .

The detailed evolution of the -function can be tracked down by exactly solving coefficients of s of Eq.(14) under certain initial conditions. As the total atomic number is conserved in the space of , the initial condition for mode is assumed to be in the ground number state and mode is in a vacuum state, that is, for and . Then the initial two-mode -function is


which gives two well-separated probability distributions of two modes with no overlap in the phase space as shown in Fig.6 (). The evolution of the probability distributions in the phase space with for modes and are subsequently shown in Fig.6. As the high dimension of , only the joint distributions of one-mode with a specific value of the other can be depicted in Fig.6. The main property shown in Fig.6 is that the local overlaps of the phase distributions can be obtained during the state evolution, and a transient complete overlap can reach at a certain time (e.g. in Fig.6). The gradually merging of distributions of modes and in the phase space indicates that only a partial QS can be maintained between two coupling modes in a Hilbert space of after a long time evolution. Surely, the unperfect overlaps of the phase distributions can be enhanced by properly adjusting the coupling rate of and the frequency detuning as well as the initial states of two modes. As Fig.6 only shows temporal joint distributions for one mode when the other one is in a specific value, so a marginal single-mode quasi-probability function of or can be achieved by integrating all the values of the other


The above results clearly show that the marginal -functions for and modes produce symmetric circular distributions which conduct breathing motions in the phase space as shown in Fig.7. The breathing modes of the ringlike distribution mask all the irregular details of distribution patterns induced by the nonlinear dynamics of two coupling modes. Although no stable distribution can be reached beyond dissipations, a correlated out-of-phase breathing motion of the ring distributions of and is revealed and a transient perfect overlap of the marginal -functions at certain time () is also found. The dynamics of the probability distributions demonstrated by the marginal -functions in Fig.7 also verifies the partial QS behavior between two scattering momentum modes in BEC.


Figure 7: The reversed dynamical phase distributions of two modes revealed by marginal -functions (upper panel) and (lower panel) with the same parameters, time scale and initial conditions as that in Fig.6.

In order to reveal the avoidance of a perfect QS between two scattering modes generated in the BEC, we can investigate the error (or deviation) operators between and modes defined by

where the zero-point-of-fluctuation units of and for quadratures of modes ( and ) are used respectively. Then the error operators obey and which leads to a fluctuation uncertainty for any quantum state in the Hilbert space as


where the deviation of an arbitrary operator is defined by , its variance is and the fluctuation above mean value is characterized by . The inequality of Eq.(22) indicates that if one quadrature () of the two modes is synchronized (), then the other () will diverge, which means that a complete QS between two coupled modes is impossible in quantum regime. As the quantum states in a closed Hilbert space of are always on a pure states following Schrödinger equation, the uncertainty of fluctuation will always plays the role to destroy the dynamical synchronization between two modes in quantum regime.

In order to reveal why two nonlinear classical oscillators can synchronize completely via their mutual couplings but, in the quantum regime, they never do, we now consider the dynamics of the error operators with both the mean values and their variances, simultaneously. The classical dynamics of the two modes are always characterized by the mean values of their positions and momenta , and, often, the CS is estimated only by


because the higher-order correlations decays very quickly with respect to the classical time scale and the variances are always much smaller than the mean values in a real system. Therefore, a complete CS is only estimated by a vanishing of the mean errors after a longtime evolution. A partial synchronization, such as the phase or amplitude synchronization, is achieved when the differences between amplitudes of or the phases of are finally locked to a constant value, i.e., the deviation

asymptotically converges to a constant of or (see Fig.4).

Figure 8: The dynamics of mean values and their relevant uncertainties for (a) (b) (c) in the symmetric case of and , and (d) (e) (f) for with and . The initial condition is determined by Eq.(19).

However, for a quantum dynamics in a closed system, the high-order coherence of the error operators should be included for a full synchronization judgement, which will lead non-negligible modulations to the mean-value dynamics between two modes. The fluctuations of the error operators beyond mean values can be estimated by


where and are position variances for modes and , and


where and are their momentum variances. Fig.8 shows the mean-value dynamics of the error operators and their corresponding fluctuations in a parametric region where the partial CS in Fig.3 occurs. Fig.8(a)(b) indicate that the varying fluctuations over the mean-value dynamics sometimes plays a dominant role in the dynamics. We can see that the error fluctuations become very large at the “revival times” that introduce large deviations of and to destroy the mean-value synchronization between two modes. Both error quadratures conduct revival and collapse dynamics due to coherent superpositions of different mode frequencies of , . The mode number of also exhibits a collapse and revival dynamics and the number fluctuations increase during the collapse region and depresses during the revival period of time (see Fig.8(c)).

We can find a quantum trapping case of this model that, if the coupling parameter satisfies , all the non-diagonal terms in Eq.(14) vanish and the quantum state will confine to the initial Fock state because it is the eigenstate of the total Hamiltonian. Surely, this case is demanding and accidental because any fluctuations on the parameters will destroy it and the transition from initial state to other states will be activated. However, if , the oscillating amplitudes and frequencies of will be dramatically decreased and the fluctuations of all the variables will be suppressed. Fig.8(d)-(f) demonstrate this case for giving and we can see a clear suppression of varying amplitudes of the mean values and their variances, both of them losing the manifest “collapse and revival” behaviors compared with that in Fig.8(a)-(c).

Because of the unceasing coherent revival of the quantum fluctuations, QS is anyway different from its classical counterpart. In order to compare CS with QS, Mari et. al. Mari () introduced a generalized criterion which is based on the fluctuations of the error operators as

This quantity gives the orbital deviation between two modes in the classical field and can measure the level of QS with a quantum limit of


In the classical dynamics, is unbounded for a regular classical orbit in the phase space (no deviation fluctuations) and will easily break the upper limit of 1/2 such as in an irregular orbit of chaotic state.

Figure 9: The dynamical measure of synchronization for different subspace of (a) and (b) with the same parameters as in Fig.8. Insets: the zoomed-in images of from the dashed rectangles.

Fig.9 demonstrates with a limits of Eq.(26) in the model of two scattering modes of BEC in different number states of . Clearly, the upper limit of holds safely during the quantum dynamics and the measure is closely below the quantum uncertainty limitation: the middle term of Eq.(26) (the dashed blue curves in Fig.9). Fig.9(b) is for the trapping case of . The calculation verifies that the quantum fluctuation is a permanent property of quantum dynamics and always exclude a full QS between the coupled quantum systems.

Recently, Ameri et al. Ameri () demonstrated that the fluctuation measure of has a similar behavior to the coherent measure of the mutual information for modes and , which, in a steady state case, can be defined by

where the Von Neumann entropy and . In our present model, the density operator is time-dependent and can then be given by

where is the transient eigenvalue of the density operator. As the present quantum system will keep a unitary evolution in the absence of dissipation, the entropy reduces to , , and the mutual information becomes .

Figure 10: The dynamical measure of mutual information for different subspace of (a) and (b) with the same parameters as in Fig.8. Insets: the zoomed-in images of from the dashed rectangles.

Now we introduce the mutual information expressed by -function to investigate the dynamical synchronization of two quantum modes as

where the joint -function and the marginal -functions , have been given by Eq.(17) and Eq.(20)(21), respectively. The mutual information concerns about the information sharing between two modes with a limit of

where is the dimension of the Hilbert space of . If the -function can separate at a transient time of by satisfying , then and the two modes will be desynchronized at time without sharing any information. This measure implies that QS is somehow related to the quantum entanglement between two modes. For a pure quantum state under unitary evolution, the mutual information can be expressed by


By using the marginal -function defined by Eq.(20), above mutual information can be calculated by


Fig.10 calculates the dynamics of the mutual information under the same parameters as that in Fig.9 just for a dynamical comparison. We can see, strictly, the mutual information , exhibits a different dynamics from that of measure . However, both measures indicate correlated dynamics of synchronization with a similar decreasing or increasing fluctuation at the same time as that shown in Fig.8. Although Fig.9 and Fig.10 seem to have opposite mean-value tendencies, their relevant “collapse and revival” behaviors keep close connections. Fig.9 and Fig.10 also show that the fluctuation amplitudes of two different measures depends heavily on the initial states of the whole system (see the trapping case of is very different from ). Furthermore, based on -functions, we can see that QS is somehow related to the quantum correlations (e.g., entanglement, discord, mutual entropy) and connected to a mathematical problem of dynamical variable separation for a joint probability function, such as .

iv.2 QS of scattering modes from BEC in a coherent state

As the importance of the initial state of BEC, we now consider QS again for two scattering modes generated from a BEC in a coherent state. In this case, the whole Hilbert space of the system is a direct sum of the subspaces of with different atomic number of , and the wave function in the whole Hilbert space should be expanded by


where the number states of . Therefore the dynamical coefficients meet


where the atomic number of the two modes have no limit in the whole Hilbert space. But within a certain subspace of specific , the above Eq.(29) returns back to Eq.(14). By using the coefficients