Qubit-qubit entanglement in the Tavis-Cummings model with two independent resonators

In this paper, a three-qubit Tavis-Cummings model with two independent lossless single-mode resonators is considered. It is assumed that the initial states of the resonator fields are thermal fields and the qubits are in genuine entangled W- and GHZ-type states. For the model under study and the specified initial states of the resonator fields and qubits, we have exactly solved the quantum Liouville equation for the full density matrix. The full density matrix was used to calculate the entanglement parameters — negativity and fidelity. The computer simulation results showed that in the investigated model, the entanglement for all initial qubit states breaks down rapidly with increasing intensity of the thermal fields of the resonators compared to the previously investigated three-qubit models. Moreover, a sudden death of entanglement is observed even for vacuum resonator fields.