Up to now, there was a tremendous interest for investigating entanglement of various Heisenberg models
[1]}, [2]}, [3]}, [4]}, [5]}, [6]}, [7]}.There have been a great number of studies focusing on the applications of entanglement in quantum teleportation, spintronics and quantum information processing [8]}, [9]}, [10]}, [11]}. Spin models are known as simple popular quantum systems in solid state physics and quantum comuting which exhibits remarkable entanglement features [12]}, [13]}, [14]}, [15]}, [16]}, specially quantum antiferromagnetic spin chains have been the subject of numerous experimental as well as theoretical studies for a long period of time [17]}. A wide variety of the analytical and numerical techniques have been developed in order to study Ising-Heisenberg spin models [1]}, [19]}, [20]}, [21]}, [22]}, [23]} due to enhance the applications of quantum information theory in solid state physics and condensed matter such as investigating thermal entanglement for the mixed spin chains [17]}, [25]}, [26]}, [27]}, [28]}, [29]}.
