{'fraction_non_alphanumeric': 0.04764705882352941, 'fraction_numerical': 0.01626470588235294, 'mean_word_length': 4.034947430771509, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 4, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 4, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Dynamics and Quantum correlations in Two independently driven Rydberg atoms with distinct laser fields

We study the population dynamics in a two-atom setup in which each atom is driven independently by different light fields, but coupling the same Rydberg state. In particular, we look at how an offset in the Rabi frequencies between two atoms influences the dynamics. We find novel features such as amplifying the Rabi frequency of one atom, together with strong Rydberg-Rydberg interactions freezes the dynamics in the second atom. We characterize the Rydberg-biased freezing phenomenon in detail, with effective Hamiltonians obtained for various limits of the system parameters. In the absence of Rabi-offset, the doubly excited state population exhibits a Lorentzian profile as a function of interaction, whereas for very small offsets it shows splitting and thus peaks. Using an effective Hamiltonian as well as the perturbation theory for weak interactions, we show that the peak arises from a competition between Rabi-offset and Rydberg-Rydberg interactions when both are sufficiently small, together with the Rydberg blockade at large interactions. The effective Hamiltonians provide us with analytical results which are in an excellent agreement with full numerical solutions. Also, we analyze the growth and the dynamics of quantum correlations such as entanglement entropy and quantum discord for the coherent dynamics. We extend our studies to the dissipative case in which the spontaneous emission from the Rydberg state is taken into account and in particular, we look at the purity and quantum discord of the steady states. To conclude, our studies reveal that the local manipulation of an atom using Rabi-offset can be an ideal tool to control the quantum correlations and in general, quantum states of the composite two-qubit systems.

A typical Rydberg setup has been modeled as a gas of interacting two-level atoms (qubits), with interactions of either dipolar or van der Waals type. The minimal setup constitutes of either two atoms or two excitations and has been a usual scenario in many of the Rydberg-based experimental studies [5,7,8,19,20,[22][23][24][25][26][27][28][29][30]. On an equivalent note, the atom-based technologies have progressed in such a way that it is possible to probe and manipulate at a single particle level [17,31,32]. Recently, involving Rydberg excitations, single atom addressing has been used to quantify the imperfections in Rabi oscillations [33], proposed to engineer phonon modes in ion crystals [34,35], to achieve controlled local operations to manipulate two atom quantum state [24], controlling resonant dipole-dipole interaction between Rydberg atoms [3], as well as freezing spin excitation dynamics [27] and controlled quantum gates [30,36].
On the other side, to characterize quantum correlations in composite systems is a nontrivial and an important task for developing scalable quantum technologies. In particular, entanglement entropy and quantum discord have gained vital importance because of their applications in condensed matter physics, especially in characterizing various quantum phases [37][38][39][40][41][42][43][44], including topologically ordered ones [42,[45][46][47] and spin liquids [48,49]. While entanglement entropy measures the entanglement between two subsystems [50,51] in pure (b) (a) Figure 1. (a) Schematic picture of two atoms trapped using microtraps with a separation of R. Each of the atoms is assumed to be driven independently by different laser fields coupling the ground state |g to the same Rydberg state |r . Γ is the decay rate from the Rydberg state |r . (b) shows the level scheme for two atoms and the corresponding laser parameters associated with each state including the Rabi couplings. The two primary trajectories by which the state |rr is populated from |gg , via |gr and |rg are shown by curved arrows.
states, discord is treated as a more general non-classical correlation, may possess a nonzero value even for separable mixed states [52,53] and is highly relevant for deterministic quantum computation with a single qubit [44,[54][55][56]. Associated to entanglement entropy in pure states, an alternate measure for entanglement in mixed states is the entanglement of formation [50]. So far, using ultracold atomic setups, the dynamics of entanglement (Rényi) entropy has been measured in quantum quench experiments [57,58].
In this paper, motivated by above developments on single particle control and manipulation, we study a minimal setup of two two-level Rydberg atoms driven continuously and independently by two distinct laser fields. In particular, we look at the effect of an offset in Rabi frequencies between the fields on the population dynamics, in the presence of Rydberg-Rydberg interactions. Interestingly, the population dynamics reveals us arXiv:1901.09265v3 [cond-mat.quant-gas] 2 Apr 2019 qualitatively novel features. An interesting scenario emerges when amplifying the Rabi coupling in one atom freezes the dynamics in the second atom, in the presence of strong atomatom interactions. This phenomenon we term it as Rydbergbiased freezing and is also well captured by effective Hamiltonians at different limits of system parameters. Without the Rabi-offset, the time-averaged doubly excited state population exhibits a Lorentzian profile as a function of interaction strength. Whereas for small offset, the Lorentzian profile gets deformed, showing a non-monotonous behaviour with a peak at small interaction strengths. Anew, obtaining an effective Hamiltonian as well as using the second order perturbation theory, we show that the peak arises from a competition between Rabi-offset and Rydberg-Rydberg interaction when both are sufficiently small, together with Rydberg blockade at large interactions. The striking quality of effective Hamiltonians obtained is that we can attain analytical solutions for population dynamics at various limits. Further, we analyze the growth and temporal evolution of quantum correlations such as the entanglement entropy and the quantum discord for the coherent dynamics. While entanglement entropy serves as a good measure for quantum correlations only in pure states, quantum discord is used for both pure and mixed states. The correlation dynamics also reveals the competition between the inter-atomic interactions and Rabi-offset when both are sufficiently small. We extend our studies to the dissipative case in which the spontaneous emission from the Rydberg state is taken into account and in particular, we look at the purity and quantum discord in the steady states [59,60]. Finally, we conclude and demonstrate that the Rabi-offset can constitute an ideal tool to control quantum correlations between the two qubits.
The paper is structured as follows. In Sec. II we discuss the schematic setup, the model Hamiltonian, the master equation and define the quantum correlations that we analyze. In Sec. III we analyze the coherent dynamics of the system as a function of Rabi-offset and interaction strength. In Sec. IV we obtain effective Hamiltonians describing the features of the dynamics in various limits of system parameters. In Sec. V we study the dynamics of quantum correlations such as entanglement entropy and quantum discord for the coherent dynamics discussed in Sec. III. Finally, in Sec. VI we consider the spontaneous emission from the Rydberg state and look at the dynamics and steady-state quantum correlations. In the appendix, we provide detailed calculations on the perturbation theory in the weak interactions limit as well as analytical results for the steady state density matrices and purity for the system and subsystems.

II. SETUP, MODEL AND QUANTUM CORRELATIONS
We consider two two-level atoms or qubits, each of them being strongly confined in two independent micro traps [see Fig. 1(a)], and are driven by distinct laser fields with Rabi frequencies Ω i and detunings ∆ i , coupling the electronic ground state |g to a Rydberg state |r . Let V 0 = C 6 /R 6 provides us the strength of the van der Waals interaction between the two Rydberg excited atoms. In the frozen gas limit [61,62], ignoring the motional dynamics of the atoms, the internal state dynamics of the setup is governed by the Hamiltonian ( = 1): whereσ ab = |a b| with a, b ∈ {r, g},σ x =σ rg +σ gr . Introducing ∆ 1 = ∆, ∆ 2 = ∆ + δ, Ω 1 = Ω and Ω 2 = Ω + ω with ω > 0, we rewrite the Hamiltonian in Eq. (1) as, In the absence of Rabi couplings (Ω = ω = 0), the eigenstates of the Hamiltonian in Eq. (2) are |gg , |rg , |gr and |rr with eigenvalues E gg = 0, E rg = −∆, E gr = −∆−δ and E rr = −2∆− δ + V 0 . The coherent dynamics of the system is obtained by numerically solving the Schrödinger equation: i∂/∂t|ψ(t) = H|ψ(t) and we use the two-atom basis {|gg , |gr , |rg , |rr }. Throughout we take |ψ(t = 0) = |gg and assume both the fields are at resonance (δ = ∆ = 0) with |g -|r transition and we focus on the effect of Rabi-offset ω on the population dynamics. With above assumption, the Hamiltonian reduces toĤ The level scheme of our two atom setup with laser parameters is shown in Fig. 1(b) and the two primary trajectories by which the state |rr is populated from |gg , via |gr and |rg are shown by curved arrows. When ω vanishes, the system is identical to the scenario in which both the atoms are driven by a global field. In that case, for sufficiently large V 0 (> Ω) starting from |gg , the system undergoes coherent Rabi oscillations between |gg and |+ = (|gr +|rg )/ √ 2 with an enhanced Rabi frequency of √ 2Ω. These dynamics is attributed to the Rydberg blockade. A non zero ω breaks the symmetry between the two atoms and hence, the state |+ lost its significance in the dynamics.
Apart from the state dynamics, we also look at how the quantum correlations develop and evolve in our setup in particular, entanglement entropy and quantum discord. Since the initial state |gg is not an eigenstate of the Hamiltonian in Eq. (3), our scenario is identical to that of a quantum quench problem in which the Rabi frequencies are instantaneously quenched from zero to a finite value at t = 0. Let us label the first atom as subsystem A and the strongly driven second atom as subsystem B. The entanglement entropy of the subsystems are, S A = − Tr(ρ A log 2ρ A ) and S B = − Tr(ρ B log 2ρ B ) whereρ A (ρ B ) is the reduced density matrix for A (B), which is obtained by the partial trace of the total density matrix, ρ = |ψ(t) ψ(t)|. The partial trace over ρ leaves the subsystems A and B in a classical mixture of pure quantum states. In terms of the eigenvalues (λ i ) ofρ A , we have S A = − 2 i=1 λ i log 2 λ i , and for a pure state S A = S B . To define the quantum discord, we briefly sketch the mutual information in the classical information theory. The classical mutual information between two subsystems A and B is de- is the Shannon entropy of the subsystem A (B), and H AB is the joint entropy of A and B. An equivalent expression for mutual information is J = H B − H B|A , where H B|A is the conditional entropy, the information needed to describe B when A is known. Though, in the classical theory I = J, in the quantum version, in which the Shannon entropy is replaced by the von Neumann entropy, there exists a discrepancy between I and J which is quantified by the quantum discord.
In the quantum theory, we have I = S A + S B − S AB and J(B : A) = S B − S B|A , where S AB = − Tr(ρ log 2ρ ) is the von Neumann entropy for the stateρ, and S AB = 0 for a pure state. Given a complete set of von Neuman projective measurements {Π i A } on the subsystem A with probabilities {p i }, the conditional entropy of the subsystem B is S B|A = i p i S B|i , where S B|i is the von Neumann entropy for the reduced den- † and I B is the identity operator. It has been shown that the total classical correlation can be obtained as,J(B : [52]. The maximization (max Note that, the quantum conditional entropy depends on the choice of the observables being measured on the other subsystem, and this results in a discrepancy between I and J(B : A) or J(A : B), which is quantified as the quantum discord. For a bipartite pure state |ψ(t) the quantum discord coincides with the entanglement entropy, i.e., D(A : B) = D(B : A) = S A = S B [63].
Once the spontaneous emission from the Rydberg state is taken into account, the dissipative dynamics and the steady state correlations are analyzed using the master equation for the two-particle density matrixρ, with the Lindblad operator given by where the operator,Ĉ i = √ Γσ i ge with Γ as the spontaneous decay rate of the Rydberg state |r . At the steady state, ∂ tρ = 0 and the dissipative mechanism drives the system eventually into a mixed state even though the system is initially prepared in a pure state. Note that, for a mixed state the discord is no longer the same as the entanglement entropy and infact, S A,B has been ruled out from being a good measure of quantum correlations or entanglement since it fails to distinguish between classical and quantum correlations [50,64,65], whereas the quantum discord remains a good measure.

III. STATE POPULATION DYNAMICS
Here, we analyze the effect of Rabi-offset ω on the Rydberg excitation dynamics, in particular, how it affects the Rydberg blockade. We look at the time averaged populations (see Figs. 2 and 3):P αβ = 1/τ τ 0 P αβ (t)dt with α, β ∈ {r, g} of the states {|αβ }, as a function of ω and V 0 . For ω = 0 and sufficiently large V 0 /Ω, the interaction induced level shift in |rr state results in the well known Rydberg blockade with average populations approachingP gg → 0.5,P rr → 0,P + → 0.5, where P + = P rg + P gr . For ω Ω, we haveP gr ≈P rg independently of V 0 . When ω is significantly large, the results shown in Fig.  2 reveal interesting features. The first thing to notice from Fig.  2(b) is that the blockade condition is merely affected by ω for sufficiently large ω. This is understood as follows: as ω becomes large, the second atom is driven strongly compared to the first one, and that results in the augmentation of P gr at the cost of P rg . Hence, to attain blockade the Rydberg-Rydberg interactions just have to dominate the Rabi coupling of the weakly driven atom, i.e., V 0 > Ω, leaving the blockade condition almost independent of ω, which become more apparent in Sec. IV. It also implies that, the blockade dynamics not necessarily always result in the generation of the symmetric entangled state |+ as in the case for ω = 0. Summing up, we have an interesting scenario: for sufficiently large ω/Ω and V 0 /Ω, the dynamics of the first atom nearly freezes and the system exhibits coherent Rabi oscillations between |gg and |gr with P rr ≈ 0 and P rg ≈ 0 as shown in Figs. 3. Note that, the freezing of the first atom emerges as a combined effect of both the Rydberg blockade from large interactions and the strong driving in the second atom, and we term this as Rydberg-biased freezing.
To gain more insights into the population dynamics, in Fig.  4 we show explicitly the average populations as a function of V 0 /Ω for different values of ω/Ω. When V 0 = 0, the two atom states are simply the product of single atom states and we have P gg (t) = cos 2 Ωt 2 cos 2 (Ω+ω)t 2 , P gr (t) = cos 2 Ωt 2 sin 2 (Ω+ω)t 2 , P rg (t) = sin 2 Ωt 2 cos 2 (Ω+ω)t 2 and P rr (t) = sin 2 Ωt 2 sin 2 (Ω+ω)t 2 for the initial state |ψ(t = 0) = |gg . The corresponding time average values areP gg =P rr = 0.375,P gr =P rg = 0.125 for ω = 0 andP gg =P rr =P gr =P rg = 0.25 for ω 0. The effect of a nonzero ω on the time average populations at V 0 = 0 has an influence as well as when V 0 0, and qualitatively noticeable especially, at small ω and V 0 as seen in Fig. 4.
As previously mentioned, when ω = 0 and V 0 0 , the population is transferred from |gg to |rr state via the entangled state |+ , thus preserving the exchange symmetry between the two atoms at any instant, andP rr decreases monotonously with increasing |V 0 | whereasP gg andP + increases and saturates to 0.5 at large |V 0 |. As soon as ω 0, the symmetry is broken, resulting inP gr P rg and as expected,P gr gets larger with larger ω for any V 0 0. As seen in Fig. 4(b) [also in Fig.  2(b)], for small values of ω,P rr exhibits a non-monotonous behaviour as a function of |V 0 |. The single Lorentzian profile centered at V 0 = 0 ofP rr for ω = 0 exhibits a partial splitting, exhibiting two peaks at ±V p 0 around V 0 = 0. This arises from the competition between the terms associated with V 0 and ω in Eq. (3) for small values of V 0 and ω. This is also evident from the plots ofP gr ,P rg , andP gg shown in Fig. 4. Focusing onP rr [ Fig. 2(b)], for very small ω such thatP gr ≈P rg , increasing |V 0 | from zero suppresses the effect of the offset ω, leads to the recovery ofP rr to the value obtained for ω = 0. Once |V 0 | dominatesP rr starts to decrease as expected from the blockade effect. As a result, the peaks separation 2V p 0 increases with increasing ω until it reaches a value for ω such that V 0 can no longer nullify the effect of ω, see Fig. 5(a). At that point, there is a substantial difference between the magnitudes ofP gr andP rg , and V p 0 decreases with further increase in ω, reaches zero. Thus, when ω becomes greater than a particular value, the populationP rr becomes a Lorentzian function of V 0 . The peak separation, V p 0 ∝ √ ω for small values of ω.
The width of the Lorentzian ν rr as a function of ω is shown in Fig. 5(b). The shaded region for small ω in Fig. 5(b) indicates the existence of double peak structure, and ν rr saturates to a constant value at large values of ω. The latter also indicates that the blockade condition is not affected by ω at large values of ω, as we have discussed above. In Sec. IV, using an effective Hamiltonian obtained in the limit {ω, V 0 } Ω, we explain the initial increment ofP rr in V 0 at small ω and also the same is done in Appendix A using the second order perturbation theory obtained in the weak interaction limit. Alternatively, one could also think, the peaks inP rr vs V 0 emerge as a consequence of quantum interference since both ω and V 0 introduce additional phase shifts in the amplitudes of quantum paths populating the state |rr . Note that, the emergent splitting ofP rr in V 0 axis has a resemblance to the Autler-Townes effect, but here using two particle states and also consequently the antiblockade effect [66]. shows the Lorentzian width (ν rr ) ofP rr as a function of ω/Ω. The width is obtained at the half-maximum of P rr distribution in V 0 /Ω. The shaded region for small ω indicates the existence of double peak structure inP rr shown in the inset of (a).
In the second step, we obtain a period average Hamiltonian, The second step is identical to removing the fast oscillations emerging from large V 0 . Even after integrating out the interaction dependent terms, still the effective HamiltonianĤ (V 0 ) e f f cannot be written as a sum of two single particle terms, indicating the existence of quantum correlations between the two atoms, which we quantify later. The first term in Eq. (9) leaves the first atom in the ground state while driving the second atom with a Rabi coupling ω + Ω, and is vice versa for the second term but with a Rabi coupling Ω for the first atom. The two terms being the correlated Rabi couplings [67], also identical to the density assisted interband tunneling for atoms in optical lattices [68]. When ω Ω, the second term in Eq. (9) is merely a perturbation to the first term, leading to the scenario of Rydberg-biased freezing. Truncating the basis to {|gg , |gr , |rg } and using the effective Hamiltonian in the Schrödinger equation, we can explicitly obtain the time dependent populations as, P gr (t) = Ω + ω 2β 2 sin 2 βt (11) where β 2 = (Ω + ω) 2 + Ω 2 /4. Note that, at t = 0,Û(0) = I, the identity operator which leaves the initial state unchanged, and is also same for the different transformations we consider below. The time average populations become:P gg = 1/2, P gr = (Ω + ω)/2β 2 /2 andP rg = (Ω/2β) 2 /2. From Eqs. (11) and (12) we can see that for ω Ω, we haveP gr ≈P rg , which is consistent with numerical results shown in Fig. 2.
For ω Ω we have P gr (t) ≈ sin 2 βt with β ≈ ω/2, indicating the Rabi oscillations between the states |gg and |gr with a Rabi frequency approximately ω. We compare these results with the numerical solutions obtained by solving the full Hamiltonian in Eq. (3), and is found to be in an excellent agreement when V 0 Ω [see Fig. 3(b)]. Now we derive the effective Hamiltonian in the limit ω {Ω, V 0 } by doing a similar procedure as above. For that we introduce a local unitary operatorÛ = exp(iωσ i=2 x t/2) acting only on the second atom which gives us the new Hamiltonian, and then averaging over a time period of T = 2π/ω, we get The effective HamiltonianĤ (ω) e f f can be written as a sum of two single particle Hamiltonians, with an effective detuning V 0 /2 for the first atom. Note that, in the new rotating frame, the Rabi coupling of the second atom is reduced to Ω. FromĤ (ω) e f f it is clear that, for V 0 Ω, the first atom is merely excited resulting in the Rydberg-biased freezing and also indicates that the two atom correlations get suppressed in the large ω limit (see Sec. V).  Finally, we consider the limit {V 0 , ω} Ω, and introduing the unitary operatorÛ = e i Ω H = V 0 4 (cos Ωt − 1) 2 |gg gg| + (cos Ωt + 1) 2 |rr rr| + sin 2 Ωt ((|gr gr| + |rg rg|) + (|gr rg| + H.c.)) + (cos 2 Ωt − 1)(|gg rr| + H.c.) + ω 2σ 2 x + iV 0 sin Ωt
Then, the effective Hamiltonian after averaging over T = 2π/Ω, in the basis {|gg , |gr , |rg , |rr }: and together with the basis vectors transformations: we can estimate the time average populations at small ω and V 0 . For instance, when V 0 = ω = 0 and with the initial state in the rotating frame, |ψ R (t = 0) =Û|ψ(t = 0) = |gg , the time average populations can be estimated directly fromÛ|gg , which gives usP gg =P rr = 0.375 andP gr =P rg = 0.125, as expected. But, when either ω 0 or V 0 0, we need to first use the Hamiltonian evolution ofĤ (Ω) e f f . Taking V 0 = 0 and ω 0, in the rotating frame, the initial |gg state undergoes coherent Rabi oscillations with |gr according toĤ (Ω) e f f , i.e., |ψ R (t) = cos ωt 2 |gg − i sin ωt 2 |gr . Then, the population in |rr state is obtained by projecting |ψ R (t) alongÛ|rr . The latter provides us the same result discussed in Sec. III. Now, we take ω = 0 and V 0 0, but V 0 Ω. TheĤ (Ω) e f f results in |ψ R (t) = cos V 0 t 8 |gg + i sin V 0 t 8 |rr exp(−i3V 0 t/8) and then projecting on to the rotated basis states given in Eqs. (16), we get for ω = 0 and V 0 Ω: These and then projecting to the rotated basis states given in Eqs.
where η = V 2 0 + 4ω 2 . The time average values becomē andP gr =P rg = 1/8 + ω 2 /(2η 2 ). When ω = 0, we retrieve the old results. For ω 0,P rr is an increasing function of V 0 which explains whyP rr increases initially with |V 0 | at small values of ω. The saturation point V P 0 shown in Fig. 5(a) for small ω can be obtained by equating theP rr given in Eq. (22) to the Lorentzian profile, f (V 0 /Ω) = 3Ω 2 /8[Ω 2 + (V 0 /ν rr ) 2 ], obtained forP rr with ω = 0, in which ν rr is obtained by fitting to the exact numerical results. Doing so, we get which is shown as a dashed line in Fig. 5(a). A comparison of time average populations from exact numerics, perturbation theory (Appendix A) and using the effective Hamiltonian is shown in Fig. 6 for small values of ω/Ω and V 0 /Ω, and they are in good agreement with each other.

V. QUANTUM CORRELATIONS
In this section, we analyze the growth and subsequent evolution of the entanglement entropy S A for the coherent dynamics discussed in Sec. III. At t = 0, we have S A (0) = 0 since the initial state |gg is a separable state. Fig. 7 shows S A (t) for different values of V 0 /Ω and ω/Ω and, is an oscillating function of time. Each row of four plots is for a fixed ω/Ω and different V 0 /Ω, and vice versa for the columns. For V 0 = 0, there is no correlation (S A = 0) between the atoms. When ω = 0 and for any interaction strengths (|V 0 | 0), S A (t) oscillates be-tween 0 and its maximum possible value of log 2 2 = 1. Larger the value of V 0 , the maximum correlation is attained between shorter intervals of time. For sufficiently large V 0 , S A (t) exhibits clean periodic oscillation between 0 and 1, indicating a complete blockade in which the system exhibits coherent Rabi oscillations between the separable |gg state and a maximally entangled |+ state. Making ω 0 significantly changes the growth and dynamics of the quantum correlations depending on the value of V 0 . Not only ω slows down the correlation growth but also lower the maximum correlation that can be attained. When both ω and V 0 are very large compared to Ω, the correlations are very well suppressed due to the Rydbergbiased freezing (see Fig. 7). For sufficiently small V 0 , increasing ω slowed down the growth of S A (t) but did not affect the maximum value of S A attained over time (see along the first column in Fig. 7 and also Fig 8(a) in which the maximum of S A (t) is shown). It indicates that ω effectively reduces the effect of interaction strengths between the two atoms for small V 0 or in other words, there exists a competition between V 0 and ω as we pointed out earlier in Secs. III and IV. The maximum value of S A (t) as a function of V 0 /Ω and ω/Ω is shown in Fig. 8(a). We introduced a lower cutoff for V 0 in the vertical axis in Fig. 8(a) since the time taken to attain the maximum correlation becomes extremely large for such small values of V 0 with large ω.
The above results opens up the possibility that the quantum correlations between the two atoms can be easily controlled by means of the Rabi-offset. To demonstrate that, we consider a time dependent ω [see Fig. 8(b)] as follows: 0, 0 ≤ αt ≤ π/2 cos 2 (αt), π/2 ≤ αt ≤ π 1, π ≤ αt ≤ 2π cos 2 (αt), 2π ≤ αt ≤ 5π/2 0, 5π/2 ≤ αt ≤ 3π (24) where α determines the rate at which ω is varied. We take sufficiently large V 0 such that the two-atom setup is in the fully blockade region thus, S A (t) exhibits periodic oscillation between 0 and 1 in the absence of any Rabi-offset. Starting from the initial state |gg , first we slowly ramp ω(t) to ω max . The ω(t) suppresses the population in |rg and consequently in |+ state [see Fig. 8(c)]. The latter results in a significant loss of quantum correlations between the atoms. As we reduce ω(t) back to zero, the correlations are again build up in the system and completely retrieve its maximum value of 1 as ω vanishes. In the example shown in Figs. 8(b) and (c), the value of α is taken such that not only the correlations are rebuild, but also the initial blockade dynamics is retrieved completely. The latter is verified by calculating the overlap functions | X|ψ(t) | 2 where |X ∈ {|gg , |+ }, see Fig. 8(c).

VI. DISSIPATIVE DYNAMICS
At this point, we discuss the effect of spontaneous emission rate from the Rydberg state |r on the dynamics as well as the quantum correlations in the steady states of the master equation in Eq. (6). Since the dissipation drives the system into a mixed state, entanglement entropy is no longer a good measure for quantum correlations [50,64,65], and we restrict ourselves to quantum discords. As stated before, for a mixed state, D(A : B) may not be always equal to D(B : A). In our two atoms setup, D(A : B) = D(B : A) only when ω = 0, due to the exchange symmetry between the atoms.
First, we look at the dynamics in the Rydberg-biased freezing regime by taking a large value for ω/Ω and V 0 /Ω, see Fig. 9. In Fig 9(a) we show the time evolution of the populations [ρ αβ (t) with α, β ∈ {g, r}] for ω/Ω = V 0 /Ω = 5 and Γ/Ω = 0.1. At shorter times we see the damped Rabi oscillations between the states |gg and |gr [inset of Fig 9(a)] whereas the states |rg and |rr are almost suppressed. Eventually, the system reaches the steady state with almost equal populations between |gg and |gr . Fig 9(b) shows both the correlations and the trace of the square of the density matrices of both the total system and the subsystems as a function of time, for the dynamics shown in Fig 9(a) Tr ρ 2 A remains close to unity during the dissipative evolution with small fluctuations initially, and becomes steady at unity as the system converges to the steady state. This makes sense, because the first atom remains frozen in the ground state due to the Rydberg-biased freezing and thus, in a pure state. Whereas, for the second atom (subsystem B), Tr ρ 2 B decreases and eventually converges to 1/2 indicating that it is in a completely mixed state i.e., a mixture of |g and |e with equal populations. As a consequence, the density matrix of the whole system is not completely mixed, and Tr ρ 2 AB converges to 1/2. The evolution of the corresponding quantum discords D(A : B) and D(B : A) are shown in Fig 9(b). Since the initial state is ρ gg = 1, a pure product state, we have vanishing discords at t = 0. In the initial period of time, they exhibit non-periodic oscillations and converge to a small value due to the Rydberg-biased freezing, as the system approaches steady state. We have D(B : A) > D(A : B) since the purity ofρ A remains larger thanρ B at any instant.
Further, we extend our calculations to wider range of V 0 /Ω and ω/Ω. In Fig. 10 we show the steady state populations: ρ gg , ρ gr , ρ rg , and ρ rr as a function of V 0 /Ω and ω/Ω with Γ/Ω = 0.1. Two features are evident from the Fig. 10 Fig. 10(c)], as discussed above. More insights into the steady states are attained from the purity of the system and the subsystems [ Fig. 11]. For small values of V 0 ( Ω), independently the value of ω, at the steady state, the total system [Tr ρ 2 AB ∼ 0.25] as well as the subsystems [Tr ρ 2 A ∼ 0.5 and Tr ρ 2 B ∼ 0.5] are completely mixed with no quantum correlations between the subsystems, see Fig. 12 for the corresponding discords. Thus, the steady state of the system is a product state, ρ AB = ρ A ⊗ ρ B and consequently Tr ρ 2 AB = Tr ρ 2 A × Tr ρ 2 B . This also indicates that for sufficiently small V 0 , the correlations which initially build up in the  For ω = 0, the correlations maximally saturate at large V 0 due to Blockade, and they start to diminish as ω increases. Discords vanish in the Rydberg-biased freezing regime where the system is described by a product state.
system from the interactions have been washed out eventually by the dissipation. As V 0 increases, for sufficiently small ω, the correlations survive in the steady state as shown in Fig. 12, and their magnitude increases with increase in V 0 , and eventually saturates to a constant value (Γ-dependent) at large values of V 0 (∼ 10Ω). The strong correlations at large V 0 for small ω is attributed to the Rydberg blockade [60]. The maximally saturated correlation attained at large V 0 for a given Rabi-offset decreases with increase in ω as shown in Fig. 11(d).
Concerning the purity, as V 0 increases Tr ρ 2 AB increases and saturates to a value of 0.5. Interestingly, the latter happens independent of the value of ω [ Fig. 11(c)], but a complete picture is accessible only through Tr ρ 2 A and Tr ρ 2 B . For instance, the purity as a function of ω/Ω in the blockade region (V 0 /Ω = 10) is shown in Fig. 11(d). Though, Tr ρ 2 AB ∼ 0.5 is independent of ω at large V 0 , the purity of the subsystems depends strongly on ω/Ω. The purity in the strongly driven atom (subsystem B) decreases as a function of ω/Ω and becomes maximally mixed at sufficiently large values, whereas that of first atom (subsystem A) increases with ω and eventually becomes a pure state at very large values of ω/Ω. Thus, the total system is in a product state for large V 0 and ω, and the purity of the becomes Tr ρ 2 AB = Tr ρ 2 A × Tr ρ 2 B as shown in Fig. 11(d). That means, the quantum correlations decrease with increase in ω for large V 0 , see Fig. 11(d) for the corresponding quantum discords. Note that, we also obtained the analytical results for steady state density matrices, the purity of the system and subsystems, see Appendix B, and are in excellent agreement with the exact numerical calculations shown in Fig. 11.

VII. CONCLUSIONS
In conclusion, we studied the dynamics and the quantum correlations in a minimal setup of two two-level Rydberg atoms driven continuously and independently by two distinct laser fields. In particular, we analyzed the effect of an offset in Rabi frequencies between the fields on the Rydberg excitation dynamics, in the presence of Rydberg-Rydberg interactions. Interestingly, we identified novel features in the system, amplifying the driving in one atoms freezes the dynamics in the second atom, in the blockade regime. Nonmonotonous behaviour of doubly excited state population as a function interaction strength for small Rabi-offsets is observed. The effective Hamiltonians obtained via unitary transformations at various limits of system parameters provide us analytical solutions for the dynamics. They are found to be in excellent agreement with the complete numerical calculations. The quantum correlations for both the coherent and dissipative dynamics are studied. We also demonstrate that the quantum correlations can be controlled using time dependent Rabi-offset. Our studies open up a new possibility to engineer quantum states in Rydberg atom setup, and the immediate question would be to extend the studies to more than two atoms, and also at different geometries.  Figure 13. The time average population in |rr state obtained using the perturbation theory (solid lines) in the weak interaction limit and the exact numerical solution (dotted lines) as a function of V 0 for different ω. The initial increment ofP rr in V 0 /Ω for small ω/Ω, discussed in Sec. III is well captured by the second order perturbation theory.