“…The low-dimensional Heisenberg spin models, involving quantum fluctuations between spins, play a significant role in this regard because they have been proven to be ideal candidates for a rigorous investigation of the entangled states under the influence of the external stimuli such as magnetic field (homogeneous or inhomogeneous) and/or temperature [ 12 , 13 , 14 , 15 , 16 , 17 , 18 , 19 , 20 , 21 , 22 , 23 , 24 ]. Moreover, many analytical and numerical calculations have been performed to examine the tuning of the quantum and thermal bipartite entanglement by varying the exchange anisotropy parameter [ 19 , 20 , 21 , 22 , 23 , 24 , 25 , 26 , 27 , 28 ], the uniaxial single-ion anisotropy [ 16 , 17 ], the Dzyaloshinskii–Moriya interaction (spin-orbit coupling) [ 18 , 19 , 20 , 26 , 27 ], the next-nearest-neighbour interaction [ 13 , 14 , 29 ], as well as by introducing impurities into the system [ 28 , 30 ].…”