Dynamics of a driven spin coupled to an antiferromagnetic spin bath

We study the behavior of the Rabi oscillations of a driven central spin (qubit) coupled to an antiferromagnetic spin bath (environment). It is found that the decoherence behavior of the central spin depends on the detuning, driving strength, the qubit-bath coupling and an important factor, associated with the number of the coupled atoms, the detailed lattice structure, and the temperature of the environment. If the detuning exists, the Rabi oscillations may show the behavior of collapses and revivals; however, if the detuning is zero, such a behavior will not appear. We investigate the weighted frequency distribution of the time evolution of the central spin inversion and give this phenomenon of collapses and revivals a reasonable explanation. We also discuss the decoherence and the pointer states of the qubit from the perspectives of the von Neumann entropy. It is found that the eigenstates of the qubit self-Hamiltonian emerge as the pointer states in the weak system-environment coupling limit.


Introduction
One of the most promising candidates for quantum computation is the implementation using spin systems as quantum bits (qubits) [1,2,3,4,5,6,7]. Combining with nanostructure technology, they have the potential advantage of being scalable to a large system. In particular, the spin in a quantum dot exhibits a relatively long coherence time compared with fast gate operation times. This makes it a good candidate for a quantum information carrier [8,9,10,11]. However the influence of the environment, especially the spin environment, on a spin, which usually causes the decoherence of the spin, is inevitable. Typically, the material environment must be present in order to host the spin or to locally control the electric or magnetic fields experienced by the spin. Therefore the decoherence behavior of a central spin or several spins interacting with a spin bath has attracted much attention in recent years [12,13,14,15,16,17,18,19,20,21,22,23,24,25].
The problem of a central spin coupled to an antiferromagnetic (AF) spin environment was investigated in Refs. [26,27,28]. In these investigations, a pure dephasing model in which the self-Hamiltonian of the undriven central spin commutes with the interaction Hamiltonian to the environment was considered. A similar problem of a spin- 1 2 impurity embedded in an AF environment was studied in the context of quantum frustration of decoherence in Ref. [29]. There the impurity spin is coupled locally in real space to just one spin of the AF environment, in contrast to the central spin model [30,31,32,33,34,35,36,37,38,39,40,41] where the central spin is coupled isotropically with equal strength to all the spins of the environment. In this paper, we consider a more general case of a central spin (qubit) driven by an external field and coupled to an AF spin bath (environment). To control the quantum states of a qubit, some types of time-dependent controllable manipulations , which can be electrical [42], optical [43], or magnetic [44], are desirable. Thus it is important and necessary to study the behavior of a qubit (central spin) or several qubits (spins) interacting with an environment (a spin bath) in the presence of a driving field [45,46,47,48]. Recently, the coherent dynamics of a single central spin (a nitrogen-vacancy center) coupled to a bath of spins (nitrogen impurities) in diamond was studied experimentally [47]. To realize the manipulation of the central spin, the pulsed radio-frequency radiation was used. Under the control of a time-dependent magnetic field, Ref. [48] considered a central spin system coupled to a spin bath. The decay of Rabi oscillations and the loss of entanglement were discussed. However, the correlations between spins in the spin bath were neglected. Recent investigations indicated that the internal dynamics of the spin bath could be crucial to the decoherence of the central spin [12,13,14,15,16,17,18,19,20,21,22]. In our driven spin model, the interactions between constituent spins of the AF environment are taken into account.
The self-Hamiltonian of the central spin, after a transformation to a frame rotating with the frequency of the driving field, can be written in a form of H S = εS z 0 +gS x 0 in the rotating wave approximation. Due to the additional driving term of gS x 0 that provides energy into the system and does not commute with the interaction Hamiltonian that couples to the AF environment in our model, the dynamics of the central spin in this case is dramatically different from the undriven pure dephasing behaviors investigated in Refs. [26,27]. It was shown in Ref. [17] that the form and the rate of Rabi oscillation decay are useful in experimentally determining the intrabath coupling strength for a broad class of solid-state systems [17]. This is also the case in our problem. After the use of the spin-wave approximation to deal with the AF environment, we employed an elegant mathematic technique to trace over the AF environmental degrees of freedom exactly to obtain the reduced density matrix of the driven spin. This enables us to study and describe the decay behaviors of the Rabi oscillations for different initial states and parameters of the central spin-AF environment model beyond the Markovian approximation and Born approximation (perturbation). With the reduced density matrix obtained nonperturbatively, we also investigate and discuss the decoherence and the pointer states of the central spin from the perspective of the von Neumann entropy. We find that the decoherence behavior of the central spin depends on the detuning, driving strength, the coupling between the central spin and the spin environment, and an important factor Ω, associated with the number of the spins, the lattice structure and the temperature of the environment. If the detuning exists, the Rabi oscillations may show a behavior of collapses and revivals; however, if the detuning is zero, such a behavior will not appear. We investigate the weighted frequency distribution of the time evolution of the central spin inversion and give this phenomenon a reasonable explanation. The form and the rate of Rabi oscillation decay is useful in determining the intrabath coupling strength and other related properties of the qubit-environment system [17]. Although we concentrate on the central spin model, our study is applicable to similar models of a pseudo spin or a qubit coupled to an environment. For example, in the study of the decoherence behavior of a flux qubit interacting with a spin environment, the same self-Hamiltonian of the qubit, H S = εS z 0 +gS x 0 , can be used [25]. In such case, ε is the bias energy and g is tunnel splitting of the flux qubit, and the two eigenstates of S z 0 correspond to macroscopically distinct states that have a clockwise or an anticlockwise circulating current [25] which can be denoted as |1 or |0 . To make contact of the driven central spin model with experiments, one may envisage a setup of a small ringshaped flux qubit located at a distance above or below an also ring-shaped AF material that has a common symmetric z-axis with the flux qubit. In this case, the coupling strength between the flux qubit and each of the constituent spins of the AF material (environment) may be regarded to be almost the same.
The paper is organized as follows. In Sec. 2, the model Hamiltonian is introduced and the spin wave approximation is applied to map the spin operators of the AF environment to bosonic operators. After tracing over the environmental modes, the reduced density matrix is obtained and the dynamics of the central spin is calculated. We also investigate the decoherence and the pointer states of the central spin for the cases of zero detuning and nonzero detuning, and calculate the von Neumann entropy which is a measure of the purity of the mixed state. Numerical results and discussions are presented in Sec. 3. Conclusions are given in Sec. 4.

Model and transformed Hamiltonian
We consider a central spin driven by an external microwave magnetic field and embedded in an AF material. To detect the central spin and control its states, the frequency ω c of the microwave magnetic field is tuned to be resonant or near resonant with the central spin. Furthermore, the central spin and the AF environment are made of spin-1 2 atoms. The total Hamiltonian of our model can be written as where H S , H B are the Hamiltonians of the central spin and the AF environment respectively, and H SB is the interaction between them [26,27,49,50]. They can be written as ( = 1) where µ 0 is the Larmor frequency and represents the coupling constant with a local magnetic field in the z direction. The magnetic field creates a local Zeemann splitting, which can then be accessed by the driving field with frequency ω c and (real) coupling strength g. The coupling strength g is proportional to the amplitude of the driving field. We assume that the spin structure of the AF environment may be divided into two interpenetrating sublattices a and b with the property that all nearest neighbors of an atom on a lie on b, and vice versa [51]. S a,i (S b,j ) represents the spin operator of the ith (jth) atom on sublattice a (b). The indices i and j label the N atoms in each sublattice, whereas the vectors δ connect atom i or j with its nearest neighbors. J is the exchange interaction and is positive for AF environment. The effects of the next nearest-neighbor interactions are neglected. For simplicity, significant interaction between the central spin and the environment is assumed to be of the Ising type. This type of interaction has gained additional importance because of its relevance to quantum information processing [24,46,48,52]. The coupling constant between the central spin and AF environment is scaled as J ′ 0 / √ N such that nontrivial finite limit of N → ∞ can exist [26,27,49,52,53]. For convenience, we denote in the following the scaled interaction between the central spin and the AF environment as J 0 = J ′ 0 / √ N . In a frame rotating with the frequency of the driving field, the Hamiltonian of the central spin becomes where the detuning Using the Holstein-Primakoff transformation, we map spin operators of the AF environment onto bosonic operators. We will consider the situation that the environment is in the low-temperature and low-excitation limit such that the spin operators in Eqs. (7) and (8) can be approximated as S + a,i ≈ a i , and S + b,j ≈ b + j . This can be justified because in this limit, the number of excitation is small, and the thermal averages < a + i a i > and < b + i b i > are expected to be of the order O(1/N) and can be safely neglected when N is very large. The Hamiltonians H SB and H B can then be written in the spin-wave approximation [51] as where M is the number of the nearest neighbors of an atom. We note here that in obtaining Hamiltonian (10) in line with the approximations of S + a,i ≈ a i , and S + b,j ≈ b + j in the low excitation limit, we have neglected terms that contain products of four operators. The low excitations correspond to low temperatures, T ≪ T N , where T N is the Néel temperature [54]. Then transforming Eqs. (9)and (10) to the momentum space, we have where γ k = M −1 δ e ik· δ . Furthermore, by using the Bogoliubov transformation, where k , and the Hamiltonians H SB and H B can be diagonalized ( = 1) and be written as where α + k (α k ) and β + k (β k ) are the creation (annihilation) operators of the two different magnons with wavevector k and frequency ω k respectively. For a cubic crystal system in the small k approximation, where l is the side length of cubic primitive cell of the sublattice. The constants − 3 2 NMJ and k ω k can be neglected. Finally, the transformed Hamiltonians become Similar to the famous spin-boson model [55,56], the interaction Hamiltonian does not commute with the self-Hamiltonian of the spin. However, a significant difference from the spin-boson model [55,56] is that the interaction Hamiltonian commutes with the bath Hamiltonian, i.e., [H SB , H B ] = 0. So the problem of the total transformed Hamiltonian, Eqs. (18)- (20), can be solved exactly, even in the case of multi-environment modes and finite environment temperatures.

Reduced density matrix
We assume the initial total density matrix of the composed system is separable, i.e., The density matrix of the AF bath satisfies the Boltzmann distribution ρ B = e −H B /T /Z, where Z is the partition function and the Boltzmann constant has been set to one. If the initial state of the qubit (central spin) is taken as then the reduced density matrix operator of the qubit can be written as where tr B denotes the partial trace taken over the Hilbert space H B of the environment. The partition function Z can be evaluated as where V is the volume of the environment. At a low temperature such that ω max >> T , we may extend the upper limit of the integration to infinity. With x = (2M) 1/2 Jkl/T and N = V /l 3 , we obtain where To obtain the exact density matrix operator, Eq. (23), we need to evaluate e −i(H S +H SB )t |1 and e −i(H S +H SB )t |0 . To proceed, we adopt a special operator technique presented in Ref. [35]. We can see that the Hamiltonian H S + H SB contains operators where S − 0 and S + 0 change the system state from |1 to |0 , and vice versa. It is then obvious that we can write where A and B are functions of operators α + k , α k , β + k , β k , and time t. Using the Schrödinger equation identity and Eq. (27), we obtain with the initial conditions from Eq. (27) given by We note that the coefficients of Eqs. (29) and (30) involve only the operators α + k α k and β + k β k . As a result, we know that A and B are functions of α + k α k , β + k β k , and t. They therefore commute with each other. Consequently, we can treat Eqs. (29) and (30) as coupled complex-number differential equations and solve them in a usual way. This operator approach allows us to solve Eq. (27) and obtain where Following the similar calculations above, we can evaluate the time evolution for the initial spin state of |0 . Let In a similar way, we obtain With Eqs. (27), (33), (34), (36), (37) and (38), the reduced density matrix operator, Eq. (23), can be obtained analytically using a particular mathematical method to deal with the trace over the degrees of freedom of the thermal AF environment. This particular mathematical method will be described in subsection 2.3. Note that our approach also applies to the case in which the qubit (central spin) is initially in a mixed state. For example, if the initial state for the qubit is ρ S (0) = |δ| 2 |1 1| + |γ| 2 |0 0|, the corresponding reduced density matrix is Eq. (23) provided that the second and third terms that contain, respectively, δγ * and δ * γ on its right-hand side are removed.

Dynamics of the central spin
With the time evolution of the density matrix operator, we can investigate the dynamical behavior of the central spin, which is of particular interest from the perspective of practical application. Here we discuss the time dependence of the expectation value of S z 0 (t). It can be written as As demonstrated in Eq. Here we introduce a particular mathematical method to perform the trace over the environment degrees of freedom. Our procedure takes two steps. First, the states with definite eigenvalues, say P 1 and P 2 , of the respective operators k α + k α k and k β + k β k are traced over. Then we sum over all possible eigenvalues P 1 and P 2 of the operators k α + k α k and k β + k β k . In this way, we obtain the final expression of S z 0 (t) , i.e., where Here the operator k α + k α k ( k β + k β k ) in A and B has been replaced by integer P 1 (P 2 ) of its corresponding eigenvalues. The two conditional partition functions are defined as where only the states with eigenvalue P 1 (P 2 ) of the operator k α + k α k ( k β + k β k ) are traced over. Because of the restriction in the trace, the evaluation of Λ(P ) is a little bit involved. To proceed, we define a generating function G(λ) for Λ(P ) in the following manner [57]. For any real number λ (|λ| ≤ 1), we define The Taylor polynomial expansion of function G(λ) at λ = 0 can be written as where G (P ) (0) represents G (P ) (λ)| λ=0 = d P G/dλ P | λ=0 , i.e., the derivative of order P of the the function G(λ) at λ = 0. From Eqs. (45) and (46), the coefficient of λ P in the expansion of G(λ) is Λ(P ) = tr B(P ) e − k ω k α + k α k /T . Therefore, we have It is easy to directly evaluate Eq. (45) to obtain Using the expansion of we obtain It can be shown that . . .

Spin inversion and Rabi oscillation decay
In this section, we present and discuss the results we obtain. We first study the dynamics of the central spin inversion. Figures 1(a) and 1(b) show the time evolution of S z 0 for different values of Ω in the case of resonance (i.e., detuning ε = 0). Figure 2(a) shows the time evolution of S z 0 also in resonance (i.e., detuning ε = 0) but with a different value of the system-environment coupling strength J 0 from that of Fig. 1(a). For the parameters chosen in Figs. 1(a), 1(b) and 2(a), the driving strength g is much larger than the coupling strength, i.e., g ≫ J 0 . As a result, the self-Hamiltonian is dominant over the interaction with the environment. The eigenstates of the self-Hamiltonian are |+ = (|1 + |0 )/ the environment on the central spin as shown in Fig. 1(b). We can see from Figs. 1(a), 1(b) and 2(a) that increasing the factor Ω and increasing the coupling constant J 0 have similar effects. One can observe that the larger the value of Ω is, the stronger the decay of the amplitude of the Rabi oscillations will be. As in the case of the spin-boson model discussed in Ref. [56], the central spin inversion does not oscillate around the value of S z 0 = 0. Its oscillations are, however, biased (shifted) a little bit toward the positive value of the initial S z 0 (0) . With the increase of the value of Ω or J 0 , this effect is enhanced.
It was shown in Ref. [17] that the form and the rate of Rabi oscillation decay are useful in experimentally determining the intrabath coupling strength for a broad class of solid-state systems [17]. This is also the case in our problem. The central spin is a two-level system. However, due to its interaction with the environment, the time evolution of the central spin inversion consists of different frequencies involved in the time series of Eqs. (40)-(42) as well as Eqs. (33)- (35). The frequencies from Eq. (35) can be written as for a pair of integers P 1 and P 2 . The probability distribution of the frequencies is then where n = |P 2 − P 1 |. It is possible that other pairs of integers P ′ 1 and P ′ 2 may exist which correspond to the same frequency κ but with m = |P ′ 2 − P ′ 1 | = n. In such case, we should add an additional probability For example, the frequencies of Eq. (60) can be rewritten as If n = |P 2 − P 1 | is chosen to be 3, then other pairs of integers P ′ 1 and P ′ 2 with m = |P ′ 2 − P ′ 1 | = 17 correspond to the same frequency κ for the parameters ε = 0.5g and J 0 = 0.05g used in Fig. 2(b). Thus, summing the probability distributions with all possible combinations of P 1 and P 2 that correspond to the same frequency κ leads to the final frequency probability distribution. Figure 3(a) shows the frequency probability distribution for the case of ε = 0, J 0 = 0.05g, and Ω = 5. It is obvious that the main frequencies are located near g which is approximately the Rabi frequency. This is also the case for Figs. 1(a), 1(b) and 2(a) where the detuning ε = 0. The contribution of many different frequencies, a consequence of interacting with the AF bath, results in the Rabi oscillation decay of the central spin. For the case of non-zero detuning ε = 0.5g, Fig. 2(b), with J 0 = 0.05g and Ω = 1, shows a decay of oscillations with an approximated Rabi frequency of ε 2 + g 2 ≈ 1.12g. The oscillation residual amplitude approaches almost zero for a time period of about 65 ≤ gt ≤ 80 [see also the inset of Fig. 2(b)]. In such case, different frequencies interfere with each other and produce the zero amplitude variation period, a behavior that is called collapse in quantum optics [58]. At longer times, the oscillation amplitude revives. This is shown clearly in the inset of Fig. 2(b). Reference [47] reported an experimental observation of a similar collapse and revival phenomenon for driven spin oscillations, though the spin bath and the corresponding interactions are different from those discussed here.
The environmental conditions affect the dynamics of the central spin and the completeness of the collapse. Figure 3(b) shows the frequency probability distribution  Fig. 2(b). The inset shows the comparison between the probability distribution of frequencies and that of the case of increasing N 10 times (i.e., Ω = 1×10, J 0 = 0.05g/ √ 10) with other parameters unchanged.
for the central spin inversion with parameters used in Fig. 2(b). The distribution has a left-hand-side cutoff at κ min = g and the center of the distribution is located at κ c = ε 2 + g 2 . One can easily understand this from the frequency relation Eq. (60). It is also obvious from Eq. (60) that the center of the distribution shifts toward larger frequencies with the increase of the detuning [this can also be seen from the comparison between Figs. 3(a) and 3(b)]. If the detuning ε exists, it is possible that the shape of the frequency probability distribution becomes approximately Gaussian, which then results in well-defined collapse and revival behavior regions [56]. If the detuning is zero, only a half side of the distribution exists [see Fig. 3(a)] and the central spin inversion will never show the (complete) collapse and revival behaviors. The frequency probability distribution is determined also by the coupling strength J 0 and the important factor Ω. With the increase of J 0 or Ω, the width of frequency distribution increases and the probability decreases. As a result, the decay of the Rabi oscillations is enhanced. We note that the results presented here depend, of course, on the number of environment atoms N in each sublattice. The influence of the number of atoms N in each sublattice on the dynamics of the central spin comes from two respects. One is the scaled coupling constant J 0 = J ′ 0 // √ N which is proportional to 1/ √ N . The other is the factor Ω, Eq. (26), which is proportional to N. We plot in Fig. 2(a) and Fig. 2 converge toward the same evolution. In other words, a nontrivial finite limit exists in the thermodynamic limit of N → ∞. This has been shown analytically in Refs. [26,27,28] for a spin coupled to an AF environment without an external ac driving field. We show numerically here that this N → ∞ limit also exists in the driven spin case. One can also notice that the difference between the original curve (in solid line) and the two curves (in dashed and dotted lines, respectively) in Fig. 2 is very small in the all time interval shown. This difference in Fig. 2(b) with ε = 0 is also small but is slightly larger at large times. These indicate that the results presented here are very close to the results in the thermodynamic limit. If we also increase the AF atom number N in J 0 and Ω simultaneously but keep other parameters unchanged, then we can find that the number of the contributed frequencies in the probability distribution of frequencies becomes larger (i.e., the separation between two adjacent distributed frequencies becomes smaller or denser) and the value of the distribution becomes smaller [see the insets in Figs. 3(a) and 3(b)]. But the width and the shape of the frequency probability distribution remain almost the same as before [see the insets in Figs. 3(a) and 3(b)]. The fact that the frequency distribution does not change much results in the very similar dynamics of the driven spin (see Fig. 2) when we increase N. The existence of nontrivial finite limits as N → ∞ and the small difference between the original dynamical evolutions with their counterparts in the thermodynamic limit are true for all the dynamic variables ( S z 0 and von Neumann entropy) presented in this article.

von Neumann entropy
The purity of a mixed state can be measured by the von Neumann entropy After diagonalizing the reduced density matrix ρ S (t), we obtain the von Neumann entropy where p k (t) are the eigenvalues of the reduced density matrix ρ S (t) and can be expressed as We investigate the von Neumann entropy of the qubit system in Fig. 4. In Fig. 4(a), for the initial state |ψ(0) = |1 , the entropy E(t) oscillatorily approaches its maximal value of 0.4228 (in blue dashed curve). If the initial state is one of the eigenstates of the qubit self-Hamiltonian, i.e., |ϕ 1 or |ϕ 2 , the von Neumann entropy E(t) remains small with time and approaches a very small value in the long-time limit (in red solid line). In the opposite, if we choose the superposition state |ψ(0) = (|ϕ 1 + |ϕ 2 )/ √ 2 as the initial state, results show that the von Neumann entropy E(t) increases monotonously to the maximal value of ln 2 = 0.6971 (in pink dot-dashed line), representing a completely mixed state. That is to say, this initial state gets maximally entangled with the environment in the course of the evolution. This initial superposition state is thus frangible to the influence of the environment and should be avoided being used in the quantum information processing in the presence of the AF environment. In all the cases, the entropy E(t) increases with the increase of the factor Ω and the coupling constant strength J 0 .
We discuss next the pointer states of the qubit from the perspective of the von Neumann entropy [59,60,61]. For g = 0, i.e., no driving field, it is obvious that the eigenstates of the self-Hamiltonian, µ 0 S z 0 , are the pointer states as the self-Hamiltonian commutes with the interaction Hamiltonian with the environment. In this case, a generic quantum state decays, after the decoherence time, into a mixture of pointer states. Thus the decoherence behavior can be described effectively by the off-diagonal elements (decoherence factor) of the reduced density matrix in the pointer state basis of the system. This is exactly the cases studied in Refs. [26,27], where the off-diagonal elements (decoherence factor) vanish in the long-time limit in the eigenstate basis of the self-Hamiltonian or of the interaction Hamiltonian with the environment regardless of the different initial states. If the driving field strength g = 0, it is difficult to identify the exact pointer states of the system (especially in the non-perturbative and non-Markovian regime) as the self-Hamiltonian does not commute with the interaction Hamiltonian with the bath. Pointer states may be defined as the states which become minimally entangled with the environment in the course of their evolution [59,60]. An operational definition in terms of the von Neumann entropy, introduced in Refs. [59,60,61], is that pointer states are obtained by minimizing the von Neumann entropy over the initial state |ψ(0) and requiring that the answer be robust when varying the time t. In general situations, pointer states result from the interplay between self-evolution and interaction with the environment, and thus their dynamical selection by the environment are complicated. For g >> J 0 or/and ε >> J 0 , the pointer states turn out (approximately) to be the eigenstates of the self-Hamiltonian [24,59]. Here the eigenstates |ϕ 1 and |ϕ 2 of the self-Hamiltonian H S = εS z 0 + gS x 0 defined in Eq. (5) can be written as where Thus for the Hamiltonian parameters chosen in Fig. 4(a) which is in a regime where the self-Hamiltonian of the system dominates, it is expected that the pointer states are very close to the eigenstates |ϕ 1 and |ϕ 2 [24,59]. By changing different initial states near the eigenstates |ϕ 1 or |ϕ 2 with other parameters fixed, it is found that the eigenstate |ϕ 1 or |ϕ 2 produces the minimum entropy increase with time (the red solid line). This suggests that the eigenstates |ϕ 1 and |ϕ 2 are dynamically selected by the environment as the preferred pointer states in the self-Hamiltonian dominant regime. Moreover, when the energy parameters g and ε of the self-Hamiltonian of the central spin system are increased further with respect to the system-bath coupling strength J 0 (i.e., the self-Hamiltonian dominates even more), the values of entropy evolution are lower (the green dashed curve is lower than the pink dot-dashed curve in the zerodetuning case, and the blue dotted curve that is close to zero is lower than the red solid curve in the detuning case) as shown in Fig. 4(b). We have also checked that the four curves with initial state |ϕ 1 or |ϕ 2 in Fig. 4(b) are, respectively, the minimum entropy curves in the course of evolution among those entropy evolution curves with initial states varied near the eigenstate of |ϕ 1 or |ϕ 2 but with other parameters fixed. This confirms furthermore again that in the regime where the self-Hamiltonian of the system dominates (with or without detuning), the eigenstates of the system self-Hamiltonian emerge as the preferred pointer states.

Conclusions
We investigate the decoherence of a qubit in an AF environment, in the presence of a driving field. The difficulty of our problem lies in the fact that the self-Hamiltonian does not commute with the interaction Hamiltonian and the internal dynamics (coupling) of the spin bath are taken into consideration at the same time. The spin-wave approximation is used to map the spin operators of the AF environment onto bosonic operators in the low-temperature and low-excitation limit. Then the resultant model is solved exactly, even in the case of multi-environment modes and finite environment temperatures. Our approach includes the environment dynamics, the qubit dynamics, and the quantum correlations between them. The influence of the AF environment on the qubit is found to depend on two factors, i.e., the coupling constant J 0 and a dimensionless factor Ω = N T 3 4 √ 2π 2 M 3/2 J 3 . Increasing the two factors, the decay time of the Rabi oscillations becomes shorter, the decoherence of the qubit is enhanced, and the qubit state becomes more mixed with time. The time evolution of the Rabi oscillations also depends on the detuning between the driving frequency and the qubit Larmor frequency. If the detuning exists, the Rabi oscillations may show a behavior of collapses and revivals; however, if the detuning is zero, such a behavior will not appears. This can be understood in terms of the weighted frequency distribution investigated here. Also the decoherence and the pointer states of the qubit are discussed from the perspective of the von Neumann entropy. It is found that the eigenstates of the qubit self-Hamiltonian are dynamically selected by the environment as the preferred pointer states in the weak system-environment coupling limit (or in the self-Hamiltonian dominant regime).