In this paper, we prove that the solution of the Landau-Lifshitz flow u(t, x) from H 2 to H 2 converges to some harmonic map as t → ∞. The essential observation is that although there exist infinite numbers of harmonic maps from H 2 to H 2 , the heat flow initiated from u(t, x) for any given t > 0 converges to the same harmonic map as the heat flow initiated from u(0, x). This observation enables us to construct a variant of Tao's caloric gauge to reduce the convergence to harmonic maps for the Landau-Lifshitz flow to the decay of the corresponding heat tension field. The advantage of the strategy used in this paper is that we can see the limit harmonic map directly by evolving u(0, x) along a heat flow without evolving the Landau-Lifshitz flow to the infinite time.flow. When β = 0, α > 0, it reduces to the heat flows of harmonic maps. In this paper, we consider the case when M = H 2 and N = H 2 .Besides the physical motivation, such as the continuous isotropic Heisenberg spin model, the gauge theory, the dynamics of the magnetization field inside ferromagnetic material (Landau-Lifshitz [26]), the Landau-Lifshitz flow (LL) is also a typical model in the differential geometry.We recall the following non-exhaustive list of works. The heat flow (HF) case (α > 0, β = 0) has been intensively studied in the past 60 years, for instance Eells and Sampson [11] for HF from closed manifolds to closed manifolds, Hamilton [16] for HF with Dirichlet boundary condition, Li and Tam [29] for HF from complete manifolds to complete manifolds. When α > 0, β ∈ R, the local well-posedness and partial regularity for weak solutions were considered by many authors for instance [12,35,45]. The local well-posedness of Schrödinger flow was studied by Sulem, Sulem and Bardos [40] for S 2 targets, Ding and Wang [10], McGahagan [34] for general Kähler manifolds. Chang, Shatah and Uhlenbeck [7] obtained the global well-posedness of maps into closed Riemann surfaces for small initial data. The global well-posedness for maps from R d into S 2 with small critical Sobolev norms was proved by Bejenaru, Ionescu, Kenig and Tataru [1, 2]. The one dimensional case was studied by Rodnianski, Rubinstein and Staffilani [38]. The dynamic behavior of LL is known in some cases. For the equivariant Schrödinger flows from R 2 into S 2 with energy below the ground state and equivariant flows from R 2 into H 2 with initial data of finite energy, the global well-posedness and scattering in the gauge sense were proved by Bejenaru, Ionescu, Kenig and Tataru [3, 4]. For the m-equivariant LL from R 2 into S 2 with initial data near the harmonic map, Gustafson, Kang, Tsai [13, 14] proved asymptotic stability for m ≥ 4 and later the case m = 3 was proved by Gustafson, Nakanishi, Tsai [15]. Moreover, there exist blow up solutions near the harmonic maps, see Merle, Raphael, Rodnianski [36] and Perelman [37] for 1-equivariant 2D Schrödinger maps, and Chang, Ding, Ye [6] for 2D symmetric heat flows. Usually the dynamics for flows defined on the Euclidean space and curved ...