In nonequilibrium thermal baths with T 1 ≠ T 2 , steady-state concurrence as a function of scaled temperature gradient of two thermal baths Δ T for different values of average temperature T_m in the strong coupling regime (a) and weak coupling regime of the IHB case (b), CHB case (c). We set ξ = 20 γ in panel (a), ξ = 2 γ in panels (b) and (c), Ω m = 20 γ , Δ Ω = 2 γ and γ A a = γ A b = γ B a = γ B b = γ .