In thermal equilibrium baths with T 1 = T 2 = T , steady-state concurrence as a function of coupling strength ξ for various temperatures of thermal baths T in panel (a). And steady-state concurrence as a function of temperature of thermal baths T for different values of coupling strength ξ in panel (c) and energy detuning Δ Ω in panels (b) and (d). Panels (a) and (b) are for strong coupling regime; Panels (c) and (d) are for weak coupling regime. We set Δ Ω = 2 γ in panels (a) and (c); ξ = 20 γ in panel (b); ξ = 2 γ in panel (d). Other parameters are set to be γ A a = γ A b = γ B a = γ B b = γ and Ω m = 20 γ .