Magnetic and Electronic Properties of Heavy Lanthanides (Gd, Tb, Dy, Er, Ho, Tm)

: Rare earth metals (REM) occupy a special and important place in our lives. This became especially noticeable during the rapid development of industry in the industrial era of the twentieth century. The tendency of development of the rare-earth metals market certainly remains in the XXI century. According to experts estimates the industry demand for chemical compounds based on them will tend to grow during the nearest years until it reaches the market balance. At the same time, the practical use of high-purity rare-earth metals requires the most accurate understanding of the physical properties of metals, especially magnetic ones. Despite a certain decline in interest in the study of high-purity REM single crystals during the last decade, a number of scientiﬁc groups (Ames Lab, Lomonosov Moscow State University (MSU), Baikov Institute of Metallurgy and Materials Science Russian Academy of Science (RAS)) are still conducting high-purity studies on high-purity metal samples. The present article is a combination of a review work covering the analysis of the main works devoted to the study of heavy REMs from gadolinium to thulium, as well as original results obtained at MSU. The paper considers the electronic properties of metals in terms of calculating the density of states, analyzes the regularities of the magnetic phase diagrams of metals, gives the original dependences of the Neel temperature and tricritical temperatures for Gd, Tb, Dy, Er, Ho, Tm, and also introduces a phenomenological parameter that would serve as an indicator of the phase transformation in heavy REMs. value linear dependence of the Neel of bulk heavy REMs on de function G 2/3 all ﬁeld, each The ﬁeld.


Introduction
The family of lanthanides (rare earth metals) includes 13 or 15 elements, from cerium to ytterbium or from lanthanum to lutetium. In the Periodic Table of elements, they are placed between barium and hafnium, so REMs usually have a separate row at the bottom of the table. From the very meaning of Mendeleev's table, one would think that all REMs have exactly the same properties. However, this was a mistake because the technology for obtaining these substances in pure form was virtually non-existent. REMs were produced only in the form of complex similar compounds, which created the illusion of their physical and chemical similarity. When REM became available in a sufficiently pure form, it turned out that the differences between them can be as great as those between other chemical elements. For this reason, intensive research on REM began only in the 1950s. Today, these metals are widely used in almost all industries and, most importantly from a scientific point of view, serve as model objects for testing and refining various physical theories.
At present, it cannot be claimed that a complete theory of pure metallic REM atoms or ions, much less their crystals or alloys and compounds, has been created, although many theoretical approaches to explain various characteristics of REM crystals have been proposed. First theoretical works on heavy REMs were published in the 1960s [1][2][3].
where J ( If one assumes that the conduction band is almost empty (almost free electrons approximation), then E k = where Θ is the angle between the vectors → k и → q . From Equations (2) and (3) Using the spherical coordinate system (except for the angle Θ one can introduce α which is the angle between → R and → q ) and replacing the summation by integration, in the low-temperature approximation, when the Fermi distribution can be considered rectangular, one can obtain J e f f ( → R) = mJ 2 8π 3 2 R 2 1 4R 2 (sin 2k F R − 2k F R cos 2k F R) (5) where k F is the value of the wave vector on the Fermi surface. If one introduces the Ruderman-Kittel function in the form: then J e f f ( Formulae (5) and (7) are two variants for RKKY interaction. It is characterized by two features: it is long-range and its amplitude decreases according to the power-law; the interaction is oscillatory, sign-variable, i.e., a particular magnetic ion can be connected with its distant neighbors both ferro-and antiferromagnetically. Determining the specific type of interaction in a given rare-earth metal requires taking into account the interaction of the ion with several neighbors of varying degrees of remoteness. The total electron energy due to the exchange interaction looks more complex than the equation for the RKKY interaction obtained in the second order of perturbation theory [26].
The exchange interaction between conduction electrons and 4f-electrons depends on the orientation of the resultant spin of the unoccupied 4f-layer of the REM ion. The helicoidal and sinusoidal arrangement of the magnetic moments of the ions is characterized by a periodic change in the spin projection on the hexagonal c-axis or an axis perpendicular Crystals 2021, 11, 82 4 of 26 to it. As a consequence, the energy of the conduction electrons is periodic with a period equal to the period of the magnetic superlattice, which is much greater than the period of the crystal lattice. As magnetic order arises and increases, the Fermi surface will continuously deform. The effects of strong changes in the Fermi surface due to the formation of helicoidal or sinusoidal structures at temperatures below the Neel temperature have been considered in detail and quantified by Dzyaloshinskii [27,28]. These effects should be considered when explaining the very appearance of spiral ordering. The specifics of the Fermi surface in heavy REMs determine the character of the magnetic structure immediately below the magnetic ordering temperature. In the course of calculations of the Fermi surface of heavy REMs [29,30] the presence of the so-called "bridges" of the full Fermi surface was discovered. In [31] it was shown that the change in the thickness of this "bridge" from one metal to another provides the correct trend in the behavior of the magnetic properties of REM; in particular, the absence of such a bridge can explain why gadolinium does not have a periodic ordering. The numerical results obtained so far, as well as the absence of a rigorous theory of rare-earth magnetism, allow only to say that the properties of the Fermi surface play a fundamental role in establishing periodic magnetic structures in heavy REMs.
The magnetic structure is significantly transformed when the temperature decreases due to changes in the energy spectrum of conduction electrons and due to the influence of magnetocrystalline and magnetoelastic interactions [26]. Polarization of conduction electrons also strongly changes in REM as a result of magnetic transitions. Thus, magnetic ordering and electronic structure in heavy REM are interrelated.
In this paper, based on the method proposed in [32], densities of states for all heavy REM were calculated ( Figure 1). The Quantum ESPRESSO utility, an integrated set of opensource computer programs for calculating electronic structures and modeling nanoscale materials based on density functional theory (DFT), plane waves (PW), and pseudopotentials, was used to calculate the density of states. A unit cell of two atoms was used to calculate the electronic structures of chemical elements with an hcp crystal lattice. The parameters of the crystal lattice a and c were taken from [33]. Calculations were performed using the projector-augmented waves (PAW) method [34,35], combining the ideas of the pseudopotential method and Linearized Augmented Planewave (LAPW) [36]. The method reproduces the exact wave function of the valence electrons, combining the flexibility of the pseudopotential method with the accuracy of the full potential methods. In the calculations of the exchange-correlation (XC) energy functional the local density approximation (LDA) was used, which depends exclusively on the electron density value at each point in space (and not, for example, on density derivatives or Kohn-Sham orbitals [37]). The specific implementation of LDA is taken from [38]. To calculate the density of states (DOS), it is necessary to calculate the Brillouin zone integrals. For this purpose, it is necessary to replace them with a weighted sum over special points in k-space. The number of points in k-space (i.e., k-grid) in which the Brillouin zone is discretized determines the accuracy of the calculation. The choice of the k-point grid depends on the problem under consideration. More grid points (smaller grid interval) gives better total energy convergence. Here the Brillouin zone is discretized using the scheme proposed in [39], to obtain a grid with equal intervals. The DOS calculation was performed in two stages. First, a self-consistent charge density was calculated with a small number of points in the k-lattice, then nonself-consistent calculations were performed using an already fixed self-consistent charge density. The use of an 8 × 8 × 6 face-centered k-lattice proved to be sufficient to achieve acceptable accuracy. The kinetic energy cutoff parameters for the wave functions were chosen equal to 60 Ry, for charge density and potential-240 Ry.  The plots show that at E=EF the DOS value is maximal or close to maximal. This insignificant change in the Fermi energy is because the concentration of free electrons per unit volume changes weakly in the series of heavy REM from gadolinium to thulium. In [40] it was predicted that the periodic electronic structure caused by the s-f interaction, which determines the magnetic properties of the considered metals, affects the shape of the Fermi surface. Self-intersections of the Fermi surface take place even in the case of infinitesimal perturbation. Thus, the exchange interaction, by removing degeneracy, leads to the splitting of the peaks in the figure n and to the formation of non-intersecting cavities of the Fermi surface. It also follows from the aforementioned work that the half-width of the energy gap between peaks is proportional to the integral of the s-f interaction and to The plots show that at E = E F the DOS value is maximal or close to maximal. This insignificant change in the Fermi energy is because the concentration of free electrons per unit volume changes weakly in the series of heavy REM from gadolinium to thulium. In [40] it was predicted that the periodic electronic structure caused by the s-f interaction, which determines the magnetic properties of the considered metals, affects the shape of the Fermi surface. Self-intersections of the Fermi surface take place even in the case of infinitesimal perturbation. Thus, the exchange interaction, by removing degeneracy, leads to the splitting of the peaks in the figure n and to the formation of non-intersecting cavities of the Fermi surface. It also follows from the aforementioned work that the half-width of the energy gap between peaks is proportional to the integral of the s-f interaction and to the mean value of the magnetic moment of the f-electrons. This fact explains the strong difference of the DOS curve of gadolinium from the other elements. The total orbital momentum of gadolinium is 0, and the spin has the highest value compared to the other elements (7/2). Therefore, the weak spin-orbital interaction in gadolinium leads to the observed absence of split peaks in Figure 1 and to their significant broadening due to the maximum value of the average magnetic moment per electron. The calculated density of states for gadolinium is in good agreement with previously published data [41][42][43]. Gd has strongly correlated electrons, with a split between spin-up and spin-down f-bands. Gd is at half-filling of f-band. Without proper consideration of strongly correlated f-electrons, even magnetic moment per Gd atom is incorrect in DFT. Subtle differences in the wave functions of conduction electrons and their energy bands in different metals lead to a variety of magnetic structures and properties of these metals.
Theoretical investigation of the problem of magnetic ordering in REM is also carried out within the framework of the molecular field model [3,5,[44][45][46]. It correctly describes important properties of ferromagnets, such as the presence of spontaneous magnetization and its temperature dependence, the high-temperature behavior of heat capacity, and magnetic susceptibility [47]. When the nature of the initial Hamiltonian and the interaction are properly accounted for, it is possible to use this theory to obtain magnetic structures that correspond to the experimental data, as well as to determine the boundaries between the phases. At the same time, it does not accurately describe the magnetization at low temperatures and in the vicinity of the Curie temperature. In addition, this theory predicts a reversal to zero of the magnetic part of the heat capacity above the Curie temperature, which does not correspond to reality. At the same time, this model is applicable to describe the magnetic properties of REM because of the strong atomic localization of f-electrons.
The complete spin Hamiltonian can be represented aŝ where the different terms describe the following contributions:Ĥ ex -exchange energy, H c f -crystal field,Ĥ me -magnetoelastic energy. The exchange contribution can be represented aŝ where ( → R i − → R j ) is exchange interaction between spins localized at nodes → R i и → R j , andŝ i , s j are the corresponding total spin operators. Overlaying this interaction on the set of states of ions renumbered using the quantum number of the total momentum, one can obtain: whereĴ i =L i +Ŝ i is is the total-momentum operator of the i-th ion. The microscopic theory gives in the simplest form: where χ( → q ) is the generalized susceptibility of conduction electrons (response function), N is the total number of unit cells in the crystal, U( → q ) is the matrix element of the average energy between the f-and s-electrons. Equation (10) describes isotropic exchange. Anisotropy can result from several mechanisms [48]. Firstly, because of the dependence of the s-f matrix elements on the orbital states of the f-electrons. Estimates show that this effect does not exceed 10% of the average Crystals 2021, 11, 82 7 of 26 exchange energy. Secondly, in ordered phases, the polarization of conduction electrons results in splitting between the longitudinal and transverse components of J ij . This effect is small because the exchange splitting of the conduction zone is~10 −4 of the typical zone width. A similar contribution is made by the spin-orbital interaction of the conduction electrons. Another, very small, contribution may be made by the modification of the s-f exchange matrix element by the spin-orbital interaction of conduction electrons [49] or indirect exchange through spin-orbital interaction of conduction electrons [50]. Although the individual contributions of the different anisotropy mechanisms are small, their sum may be significant enough to result in appreciable anisotropic effects.
The terms describing the interaction with the crystal field can be included by decomposing the electrostatic potential around the lattice node occupied by the REM ion into a series of spherical harmonics and passing from spherical harmonics to equivalent orbital moment operators by the Wigner-Eckart theorem. The number of members of the series is limited by the point group of symmetry of the node and the orbital momentum of the f-electron (l = 3). The general equation for the Hamiltonian describing the crystal field contribution can be written asĤ Here the first summation is performed on nodes i, and θ m l (J i ) are operators equivalent to the corresponding spherical harmonics obtained by replacing the Cartesian components of spherical harmonics with Cartesian components of J, B lm = α l r l A lm (13) where α l is the Stevens factor, which depends on the specific environment of the ion, r l means the probable value according to the radial distribution, A lm are the crystal field parameters, which can be calculated from the charge distribution around the ion. In many cases magnetic ordering can be accompanied by magnetostrictive stresses, which modify the crystal field, lowering the ion energy. The resulting magnetoelastic contributions are usually written in terms of equilibrium stresses and magnetostrictive coefficients [51].
Various magnetic structures observed in REM arise as a result of the combined action of the three types of energy discussed above: exchange energy, crystal field, and magnetoelastic energy. The exchange energy plays a decisive role in the case of noncollinear structures, since the Ruderman-Kittel function is sign-variable, and the force which it determines is long-range. Other energies tend to be localized in the lattice, determining moment orientation and leading to ferromagnetism, or at least commensurability of the crystal and magnetic structures. Magnetic ordering in heavy REMs can be derived from the molecular field model [44,46,51].
In general terms, all structures can be described by the following equation: where the summation is carried out for all periods α (with the corresponding wave vectors → q α ), which exist simultaneously, and the vector amplitudes → J qα are complex. The general relation (14) includes ferromagnetism with moments lying in the basis plane; conical spirals modulated along the c-axis (including the limiting case of pure modulation described by a longitudinal spin density wave); ferromagnetic structures directed along the c-axis. In phenomenological theory, it is more convenient to write these structures as functions of the angle Θ between neighboring moments along the axis: Axes ξ, η, and ζ are chosen in the hexagonal system so that the axis ζ is directed along the axis c of the crystal, the axis ξ is the easy magnetization axis in the basal plane, η is perpendicular to ξ. In this case, relations (15) do not describe the ferrimagnetic structure and spin waves.
In molecular field theory, the exchange energy per ion is obtained from Equation (10) as: where J ( From Equation (16) one can see that the system will choose such a value of → q m as the magnetic period, for which J ( → q ) is minimal. Minimizing energy at a given value → q m considering also the axial anisotropy leads to the following values of the angle Θ: The latter equation describes the cone structure in general; the coefficients describe the axial anisotropy of the crystal and are calculated from the equations for the crystal field. Stability conditions for these structures are discussed in [4]. At high temperatures, the member containing V 2 , is maximal, the magnitude and sign of V 2 determine whether the moments lie in the basal plane (V 2 > 0), whether they form a cone structure or are located along the hexagonal axis, which corresponds to a ferrimagnetic structure or a longitudinal wave of spin density (V 2 0). Regular spiral or cone structures are obtained if one neglects anisotropy in the basal plane and magnetoelastic effects when calculating energy. Accounting for these relatively small influences can lead to helicity modification, the formation of "bunches" of spins near-certain directions. This phenomenon was discovered in neutron diffraction, and later a theory was proposed for the case of a flat helix [52]. A similar consideration, taking into account the magnetoelastic effects, can be carried out for the general case of a conical structure [53]. When a magnetic field is applied to spiral or conical structures, distortion with a small excess moment along the field occurs. In fields exceeding the first critical value (at a given temperature), a fan phase emerges, in which moments are located near the field direction, and in fields exceeding the second critical value, the fan phase breaks down, passing into collinear ferromagnetic [2]. If the anisotropy in the basis plane is large enough, the fan phase may not occur.
In general, the magnetic and crystal structures of rare-earths are incommensurable with each other. However, because the periods of these structures change unequally with temperature, there are situations when the commensurability can take place. The reasons leading to the transition to the commensurable phase are as follows: • axial anisotropy, together with the magnetoelastic forces, favors the maximum momentum at each node for a given configuration of exchange forces, which can only be achieved in a commensurate structure; • exchange interaction also favors the maximum momentum, providing a change with a periodicity not inconsistent with the maximum I(q). At higher temperatures, increasing entropy in the incommensurable phase likely plays a role in maintaining incommensurability → q m and the development of instability at potential transition points, preventing the occurrence of a phase transition.
With a further decrease in temperature, the growth of the anisotropic terms of the Hamiltonian leads to phase transitions of the first order to other structures as in terbium and dysprosium a planar spiral structure undergoes the transition to a ferromagnetic phase with the of easy magnetization direction lying in the basal plane. Magnetostrictive stresses play a decisive role in this transition.
In some cases, a conical structure (ferromagnetic spiral) arises at low temperatures. The conditions for the emergence of such a phase were considered within the framework of the molecular field model [4]. It was shown that the reason for the transition to a conical structure lies in the peculiarities of the temperature dependence of the axial anisotropy coefficients. In the case of a sufficiently slow temperature change of the anisotropy coefficients near the helix-cone transition temperature, the instability can lead to the appearance of spin waves with a wave vector close to → q m , which, in turn, leads to an intermediate phase between the spiral and the cone [53]. As it is shown in [54], this intermediate phase is a canted spiral [55].
In pure REMs, other intermediate phases are observed that are quite simply described, for example, by a superposition of spiral and longitudinal spin density wave structures with equal periods. Such a phase occurs, in particular, in erbium. The calculation of the magnetic structure of erbium was performed within the framework of the molecular field model [56], both the effective crystal field parameters and the interplanar two-ion interaction parameter were taken into account. The parameters were estimated from the experimentally obtained values of transition points and magnetization. The details of the obtained magnetic structure were compared with the results of neutronography [57]. The calculation predicted the real existing structures, although there were quantitative discrepancies with the experiments.
The molecular field method cannot be applied to strictly describe the properties of magnetics in the two limiting regions-near the points of phase transitions and at low temperatures. At the same time, the low-temperature properties of heavy REMs are could be determined by spin waves' theory [58]. When applied to rare-earth metals, it (both for the case of ferro-and antiferromagnetic ordering) should be designed considering anisotropic and magnetoelastic forces. Various aspects of the theory of spin waves in rare-earth metals were considered in [59]. In the case of a ferromagnetic rare-earth metal, the simple theory must be supplemented by considering the energy contributions of the crystal field and magnetoelastic forces. The first of these contributions was taken into account in [1], in which dispersion relation for zero temperature is proposed. Anisotropy leads to the appearance of a gap in the spectrum of spin waves, the predicted value of which, for example, for dysprosium is about 24 K. The temperature and field dependences of the spectra are usually obtained by renormalizing the coefficients and exchange integrals. Magnetoelastic forces can be considered only in two marginal approximations. In one of them, the frozen lattice approximation, it is assumed that in the case of small k, the lattice deformations are frozen so that they have statistically equilibrium values independently of the oscillations of magnetic moments. In another extreme case, it is assumed that for small k the strains can follow the momentum oscillations. Variants of the theory differ in the number and type of terms containing magnetostriction constants. Even greater difficulties arise when building a theory of spin waves in antiferromagnetic structures. For REMs, the dispersion relations are even more complicated [59]. At present, the theory of spin waves in complex structures in the presence of many commensurate contributions to the energy of the crystal is not yet sufficiently developed for rare-earth metals to serve as a convenient model object for their verification.
In [60] it is theoretically shown that REM with different anisotropy configurations at the Neel temperature can undergo both first-and second-order phase transitions. For the hexagonal lattice in the molecular field approximation [61] a wide range of different new types of phase diagrams has been proposed. It was found that the Neel temperature in the presence of a strong exchange within the basal plane can split in the magnetic field into two lines of phase transformations of the first and second order, which was subsequently proved experimentally [62][63][64] The dominant character of a magnetic ordering turned out to be connected with the sign of the coefficient determining the relation of the quadrupole moment of f-electrons with the direction of the total mechanical momentum J. In terbium, dysprosium, and holmium the coefficient is negative, while in erbium and thulium it is positive. This explains the fact that in the last two metals the higher order crystal field components favor the orientation of the moments parallel to the c-axis. At the same time, the particular arrangement of the moments depends on the relative magnitudes and signs of these components. In thulium, all crystal field effects only favor a single direction of spins parallel to the c-axis. In other heavy REMs, there is always a probability of finding a moment with a component outside the plane, and it increases in the series from holmium to terbium.
Spin ordering in heavy REMs is mainly of two types. First, longitudinal oscillations of the spin amplitude along the c-axis (erbium, thulium). Second, a helical structure with moments lying in the basal plane (terbium, dysprosium, holmium). More complex structures are combinations of these two structures. Although the existence of periodic spin structures is supported by long-range exchange forces, their exact form is determined by anisotropy and magnetostriction, each of which favors the ferromagnetic ordering. At lower temperatures, the energies associated with anisotropy and magnetostriction become larger and larger, and their contribution to the decrease in spin turning angles with decreasing temperature increases, hence they should be important in determining the possible transition to a collinear FM state at the lower Curie temperature. In [59,65,66] it is assumed that magnetostriction is the dominant process that stabilizes the FM state. At the AFM-FM transition, the energy barrier separating the magnetically ordered states is compensated mainly by the magnetoelastic energy increasing during cooling [67,68]. Spontaneous magnetostriction decreases the thermodynamic potential of the FM phase more strongly than the AFM, which leads to a decrease in the TD potential difference and phase change. It was shown by Dzyaloshinskii [28], that the magnetostriction along the hexagonal axis c in REM is due to the dependence of the exchange energy on the lattice parameter c. Thus, the AFM-FM transition arises not only due to changes in the energy spectrum of the conduction electrons but also due to an increase in the magnetoelastic energy of the exchange interaction when cooling below the Neel temperature.
Despite all said, the science of rare-earth magnetism is still on the threshold of understanding how the interaction of different groups of electrons leads to manifold magnetic structures and observable excitation spectra. The goal of a fully microscopic theory of lanthanide magnetism is to achieve the best possible understanding of magnetic structures, excitation spectra, and microscopic magnetic properties, expressed in terms of electronic structures of f-and s-electrons, as well as their interactions with each other and the crystal lattice.

Magnetic Phase Diagrams of Heavy Lanthanides
The basic physical and chemical properties of REMs are now sufficiently well studied. At the same time, a great variety of new magnetic structures and features of phase transitions make these metals irreplaceable objects for studying fundamental questions of the physics of magnetic phenomena.
Among heavy REMs, all elements except ytterbium (not a typical representative of this family) have hcp structure at room temperature with a space group P6 3 /mmc with a sequence of layers ABAB. In the vicinity of the melting point (1585, 1629, 1680, and 1734 K for gadolinium, terbium, dysprosium, and holmium, respectively) these 4 metals change the crystal structure to bcc. Similar transitions in erbium and thulium are not observed. It is important to note that the melting temperatures of the metals monotonically increase from gadolinium to erbium with a slight decrease (by 5 K) for thulium compared to erbium. At the same time, the boiling point, on the contrary, monotonically decreases except for erbium, where its boiling point slightly (by 53 K) exceeds the analogous value of holmium [23].
Literature analysis shows that when cooling below room temperatures, all metals become ferromagnetic at some temperature characteristic for each metal (the Curie temperature), and all except for Gd pass into this magnetic state from the antiferromagnetic phase, which is formed directly from the paramagnetic phase. The Curie temperature represents the highest temperature at which ferromagnetism is observed in the absence of an external magnetic field. The temperature dependence of the critical field was previously studied both directly by measuring magnetization [93], and indirectly [94]. In the general case, the critical field increases almost linearly with temperature up to a maximum value in the vicinity of the Neel temperature, and then rapidly decreases when the transition to the paramagnetic state occurs.
Let us consider the behavior of each metal of the heavy REM series separately. Gadolinium is a true ferromagnet over the entire temperature range, but it has [72] changes in the spin orientation, which was found in the magnetocrystalline anisotropy studies: in the temperature range 232-294 K the moments are parallel to the c-axis and below 232 K the angle of rotation of the moment relative to the c-axis increases, reaching a maximum of 75 • at 180 K, then the rotation angle decreases and remains within 32 • at 4.2 K. The question about the magnetic structure of gadolinium in the region of temperatures immediately below the Curie temperature is still open [28,95]. Gadolinium is the only metal among heavy REM that does not have an AFM phase and goes directly to the FM state at about 293 K [96]. Some evidence for the existence of the AFM phase in the very low critical field near 240 K was presented in [79,97], however, the results were not reproduced in other similar studies and were correlated with the presence of impurities and defects in the structure of the investigated samples. It is possible that these observations can also be associated with the anomalous behavior below 240 K along the easy magnetization direction [98][99][100]. According to neutronographic data [99], the c-axis is the easy one in the temperature range from the Curie temperature to the spin-reorientation transition temperature 232 K (Figure 2).
In Tb (Figure 3) spiral AFM structure exists in the range~222-228 K in a magnetic field less than 155 Oe. The critical field values of 300-800 Oe published earlier [77,101], can probably be related to the presence of impurities in the samples, which play the role of capture centers and thus prevent the destruction of the weakly bound AFM phase. Dependence of the critical field on temperature H crit (T) has a symmetric form with an obvious maximum at 226 K. As the magnetic field increases in the same temperature interval in fields up to 5 kOe, a wide region of an intermediate fan phase existence is observed. It exists between the FM and PM order, that is, in the temperature range 222-227 K and the maximum field value of 5 kOe. In the low-field region, it is replaced by helical AFM ordering. The fan phase disappears at~227.3 K in the field of~5 kOe. A characteristic feature of the phase diagram presented in this paper is a rather wide region of the fan phase existence. It is possible to find only one phase diagram containing the fan-type phase [101], other authors did not even include this phase in the diagram, since its existence was probably proved only by only one experimental technique. However, the final solution of the fan existence in terbium remains open, since ab initio calculations do not allow calculations of fan-like structures in heavy REMs. They only show that, if one simulates a collinear AFM structure in the specified range of temperatures and magnetic fields, one can conclude that its existence is more likely along the hard magnetization axis than in the basal plane. No evidence for the existence of any collinear AFM structures could be found in terbium. Thus, detailed studies of elastic thermal neutron scattering on single-crystalline terbium are necessary to more accurately determine the type of magnetic ordering in this temperature and field range. In Tb (Figure 3) spiral AFM structure exists in the range ~222-228 K in a magnetic field less than 155 Oe. The critical field values of 300-800 Oe published earlier [77,101], can probably be related to the presence of impurities in the samples, which play the role of capture centers and thus prevent the destruction of the weakly bound AFM phase. Dependence of the critical field on temperature Hcrit(T) has a symmetric form with an obvious maximum at 226 K. As the magnetic field increases in the same temperature interval in fields up to 5 kOe, a wide region of an intermediate fan phase existence is observed. It exists between the FM and PM order, that is, in the temperature range 222-227 K and the maximum field value of 5 kOe. In the low-field region, it is replaced by helical AFM ordering. The fan phase disappears at ~227.3 K in the field of ~5 kOe. A characteristic feature of the phase diagram presented in this paper is a rather wide region of the fan phase existence. It is possible to find only one phase diagram containing the fan-type phase [101], other authors did not even include this phase in the diagram, since its existence was probably proved only by only one experimental technique. However, the final solution of the fan existence in terbium remains open, since ab initio calculations do not allow calculations of fan-like structures in heavy REMs. They only show that, if one simulates a collinear AFM structure in the specified range of temperatures and magnetic fields, one can conclude that its existence is more likely along the hard magnetization axis than in the basal plane. No evidence for the existence of any collinear AFM structures could be found in terbium. Thus, detailed studies of elastic thermal neutron scattering on single-crystalline terbium are necessary to more accurately determine the type of magnetic ordering in this temperature and field range.
In fields above 5 kOe, there are only two magnetic phases. After the fan phase is completely suppressed, the boundary between the FM and PM phases begins to slowly shift towards higher temperatures-up to 228 K in a 75 kOe-with increasing magnetic field, which is the typical behavior for the Curie temperature [71,102]. In fields above 5 kOe, there are only two magnetic phases. After the fan phase is completely suppressed, the boundary between the FM and PM phases begins to slowly shift towards higher temperatures-up to 228 K in a 75 kOe-with increasing magnetic field, which is the typical behavior for the Curie temperature [71,102]. In dysprosium, the region of helical AFM phase existence (85-178.5 K) is much wider than in terbium (Figure 4). At a high magnetic field in dysprosium, intermediate fan phases exist at temperatures above 125 K. The first intermediate fan phase, the existence of which has been confirmed by different experimental data of magnetization, MCE and ultrasonic measurements [73,103,104] exists in the temperature range from 125 to 180 K in  In dysprosium, the region of helical AFM phase existence (85-178.5 K) is much wider than in terbium (Figure 4). At a high magnetic field in dysprosium, intermediate fan phases exist at temperatures above 125 K. The first intermediate fan phase, the existence of which has been confirmed by different experimental data of magnetization, MCE and ultrasonic measurements [73,103,104] exists in the temperature range from 125 to 180 K in the range of magnetic fields from 5 to 23 kOe. At the same time, recent neutronographic studies and theoretical calculations [46,105] indicate heterogeneity and the presence of a broader region of magnetic ordering of this type. In this case, probably, the vast region of the fan phase existence consists of areas with long and short-range order depending on the chemical purity of the sample, which in some works are interpreted as different types of fan ordering [73,106,107]. An extensive intermediate phase region, probably a fan with a helical structure, exists in the range 125-180 K and 5-25 kOe. A similar phase is observed in holmium [108], for which the existence of a helifan structure was theoretically proved. The specific boundaries of the helifan-type phase existence require clarification by neutronography. In addition, a second intermediate phase in the region of 180 K in the field range of 6-23 kOe is observed in dysprosium. A similar result was previously published in [88], which suggested the existence of a vortex magnetic structure in this region, but it found no further experimental confirmation, although the existence of such unstable structures had previously been predicted theoretically [109] in the case of a magnetic field applied along the hard magnetization c-axis. Because the vortex magnetic ordering is no longer found in any of the heavy REMs, it is reasonable to assume the existence of another intermediate fan phase in this region, which is indirectly confirmed by the results of neutronography [105].  The phase detected in the temperature range from 90 to 130 K is probably intermediate between FM ordering and fan phase with long-range ordering. Likely, it exhibits predominantly FM ordering with portions of fan phase with short-range ordering. [46,105] In the Curie temperature region, according to theoretical calculations [46], there is a canted FM ordering. The existence of such a phase has not yet been established experimentally, and the reason for its occurrence may be a combined magnetic-structural phase transition in dysprosium near 90 K, which causes local distortion of the FM structure. At the same time, the region of such a phase existence is very limited; in lower temperatures, the common collinear FM ordering is observed.
In the region of temperatures above the Neel temperature (179 K) dysprosium is paramagnetic. Some works suggest the presence of antiferromagnetic clusters in the PM region [73,110], but these assumptions require experimental verification.
Holmium is ferromagnetic below the Curie temperature (TC = 20 K) and antiferromagnetic between the Curie and the Neel temperature (TN = 131 K) in the range of magnetic fields from 0 to ~30 kOe ( Figure 5). In fields up to 2 kOe at low temperatures (13÷20 The phase detected in the temperature range from 90 to 130 K is probably intermediate between FM ordering and fan phase with long-range ordering. Likely, it exhibits predominantly FM ordering with portions of fan phase with short-range ordering. [46,105] In the Curie temperature region, according to theoretical calculations [46], there is a canted FM ordering. The existence of such a phase has not yet been established experimentally, and the reason for its occurrence may be a combined magnetic-structural phase transition in dysprosium near 90 K, which causes local distortion of the FM structure. At the same time, the region of such a phase existence is very limited; in lower temperatures, the common collinear FM ordering is observed. In the region of temperatures above the Neel temperature (179 K) dysprosium is paramagnetic. Some works suggest the presence of antiferromagnetic clusters in the PM region [73,110], but these assumptions require experimental verification.
Holmium is ferromagnetic below the Curie temperature (T C = 20 K) and antiferromagnetic between the Curie and the Neel temperature (T N = 131 K) in the range of magnetic fields from 0 to~30 kOe ( Figure 5). In fields up to 2 kOe at low temperatures (13÷20 K) a spiral ordering of FM type (cone) with the cone axis directed along the hard magnetization c-axis is observed. At temperatures below 50 K magnetic satellites of the 5th and 7th orders are observed, which indicates a strong deviation of the magnetic moments in the direction of the easy magnetization axis in the basal plane [111]. Application of the field in the basal plane leads to the transformation of the conic ordering into an FM structure [112]. Anomalies were also found to represent transitions into intermediate phases of the fan type with helicoid structure (helifan) and "simple" fan (as the magnetic field increases) as the AFM phase transforms into ferromagnetic phase in the temperature range from 45 to 95 K. According to theoretical estimations, the existence of helifan type 3/2 is recognized as the most probable in holmium, however, neutronographic studies are required to accurately determine its type. A further field increase causes the fan structure to collapse and the transition to full FM ordering [111,112]. In holmium, there is an intermediate region of coexistence of the two phases between the fan and FM phases. In comparison with the five intermediate phases observed in [113], detected from magnetoresistance measurements, only one phase is observed in the present work, and it can be assumed that, despite the high purity of the crystals used in the present work, impurity effects still play a significant role in those regions of the phase diagram that are most sensitive to the magnetic field (intermediate metastable phases). These impurities can be the cause of mixed-phase states, such as "ferro+fan" which combines features of both fan and FM structures.
Two regions of spin-slip structures existence are observed between 20 and 35 K (I) and between 35 and 42 K (II). The extensive region in the phase diagram between ~95 K and the Neel temperature in the field range from 30 to 80 kOe is associated with the presence of a spin-slip structure (III), an indication of which was observed from magnetic measurements at 98 K in [84]. The existence of such structures has previously been predicted [114,115], but so far there has been no experimental confirmation of these structures' existence using a set of different experimental techniques. The spin-slip model explains the transitions in magnetic helices in REM in terms of simple commensurable structures. The spin-slip theory is based on experimentally observed localized spin defects due to, among other things, the purity of the samples [116]. Therefore, the number of varieties of such spin configurations can be extremely large. In [117] it has been suggested that the appearance of spin-slip type defects in a simple (antiferromagnetic) helix can lead to var- A further field increase causes the fan structure to collapse and the transition to full FM ordering [111,112]. In holmium, there is an intermediate region of coexistence of the two phases between the fan and FM phases. In comparison with the five intermediate phases observed in [113], detected from magnetoresistance measurements, only one phase is observed in the present work, and it can be assumed that, despite the high purity of the crystals used in the present work, impurity effects still play a significant role in those regions of the phase diagram that are most sensitive to the magnetic field (intermediate metastable phases). These impurities can be the cause of mixed-phase states, such as "ferro+fan" which combines features of both fan and FM structures.
Two regions of spin-slip structures existence are observed between 20 and 35 K (I) and between 35 and 42 K (II). The extensive region in the phase diagram between~95 K and the Neel temperature in the field range from 30 to 80 kOe is associated with the presence of a spin-slip structure (III), an indication of which was observed from magnetic measurements at 98 K in [84]. The existence of such structures has previously been predicted [114,115], but so far there has been no experimental confirmation of these structures' existence using a set of different experimental techniques. The spin-slip model explains the transitions in magnetic helices in REM in terms of simple commensurable structures. The spin-slip theory is based on experimentally observed localized spin defects due to, among other things, the purity of the samples [116]. Therefore, the number of varieties of such spin configurations can be extremely large. In [117] it has been suggested that the appearance of spin-slip type defects in a simple (antiferromagnetic) helix can lead to various kinds of spontaneous ferrimagnetic structures (in the limiting case) with macroscopic momentum, which appears to be observed in thulium [106,[118][119][120].
Outside the temperature range of 19-128 K and the field range of 10 to 80 kOe holmium is a typical ferromagnet. At temperatures above 131.7 K holmium is paramagnetic ( Figure 5).
The most complex magnetic structures are observed in erbium ( Figure 6). Below 20 K erbium has a ferromagnetic helix (cone) type structure, which in higher fields consistently transfers to a fan and then to a collinear FM phase. The value of the Curie temperature according to different sources varies in the range from 16 to 20 K, which, apparently, is associated with the purity of the studied samples. [121,122] Erbium, like all heavy REMs, has an extensive AFM ordering region in the range of 20-85 K. The subregions of the long-range magnetic order are also distinguished in this region. In the range of 20-50 K, both transverse and longitudinal components of the magnetic moment oscillate, i.e., a complex spiral is realized. In this temperature range, spin-slip defects of various types and, in the limit, ferrimagnetic ordering are likely to occur. This region is intermediate: a helicoidal component of the magnetic moment in the basal plane exists here. In the interval 52-84 K, a longitudinal spin-wave is observed in erbium, which is characterized by the oscillation of the longitudinal component of the magnetic moment along the caxis [57,123,124]. The spin-wave is a limiting case of the helical structure. Below 52 K, this wave is supplemented by oscillations of the magnetic moment components in the basal plane of the helix type, with the period of the helicoid coinciding with the period of the longitudinal spin-wave. A whole cascade of transitions between spin-slip phases of various configurations exists in erbium in the temperature range 40-70 K [125]. In the same temperature range, an additional transformation of the magnetic structure associated with the transition of erbium to a commensurate phase with wave vector 6/23 is observed. Apparently, the coexistence of phases with different wave vectors of the magnetic structure is observed in this temperature range. The temperature boundary of this effect (in the region of 40 K) is determined by the temperature region of a significant increase in magnetic anisotropy, sufficient to suppress the effect of commensurability [55]. In erbium, the presence of certain phase transitions depends particularly on the purity of the samples: thus, when in [126,127] such commensurate points were found as 2/7, 3/11, 4/15, 5/19, 6/23, 1 4 и5/21, then in [63] 5/19, 1 4 и5/21 are only found. Above the AFM ordering boundary, there is a narrow region of phases with fantype helicoidal ordering (helifan), as well as a number of intermediate phases containing features of cone, fan, and ferromagnetic ordering. The number of commensurate structures in erbium does not exceed 7-8 [126,128]. Impurity atoms can form a so-called impurity magnetic center and distort the magnetic structure in their environment [129,130]. Applying or removing the field leads to irreversible changes in the local magnetic moment on the impurity atom (using tin atoms as an example).
The magnitude of the cone half-dissolution angle Θ in erbium is related to the structural arrangement of the AFM structure. Exposure to a magnetic field can induce a series of consecutive spin-slip structures in the basal plane. The mechanism of this relationship can be explained as follows. Magnetizing of the crystal along the c-axis leads to a magnetostrictive change in the initial sample volume and the lattice parameter c. This violates the commensurability conditions of a given magnetic structure and can lead to a transition to a structure with a different wave vector. The resulting violation of the initial balance of crystal field energy, exchange, and magnetoelasticity causes a jump change in the Θ angle value. An explanation related to the influence of the magnetic field along the c-axis on the magnetic structure in the basal plane, determined by the competition of interplanar exchange and sixth-order anisotropy, is also possible. The transformation of the initial structure when the cone half-solution angle changes are related to the different dependence of the exchange (~sin 2 Θ) and sixth-order anisotropy (~sin 6 Θ) on this angle [131]. This explanation is also true in the case of holmium.
is supplemented by oscillations of the magnetic moment components in the basal plane of the helix type, with the period of the helicoid coinciding with the period of the longitudinal spin-wave. A whole cascade of transitions between spin-slip phases of various configurations exists in erbium in the temperature range 40-70 K [125]. In the same temperature range, an additional transformation of the magnetic structure associated with the transition of erbium to a commensurate phase with wave vector 6/23 is observed. Apparently, the coexistence of phases with different wave vectors of the magnetic structure is observed in this temperature range. The temperature boundary of this effect (in the region of 40 K) is determined by the temperature region of a significant increase in magnetic anisotropy, sufficient to suppress the effect of commensurability [55]. In erbium, the presence of certain phase transitions depends particularly on the purity of the samples: thus, when in [126,127] such commensurate points were found as 2/7, 3/11, 4/15, 5/19, 6/23, ¼ и 5/21, then in [63] 5/19, ¼ и 5/21 are only found. Above the AFM ordering boundary, there is a narrow region of phases with fan-type helicoidal ordering (helifan), as well as a number of intermediate phases containing features of cone, fan, and ferromagnetic ordering. The number of commensurate structures in erbium does not exceed 7-8 [126,128]. Impurity atoms can form a so-called impurity magnetic center and distort the magnetic structure in their environment [129,130]. Applying or removing the field leads to irreversible changes in the local magnetic moment on the impurity atom (using tin atoms as an example).
The magnitude of the cone half-dissolution angle Θ in erbium is related to the structural arrangement of the AFM structure. Exposure to a magnetic field can induce a series of consecutive spin-slip structures in the basal plane. The mechanism of this relationship Thulium (Figure 7) also has an extensive AFM ordering region in the range from zero to the Neel temperature (56 K). In the region of low temperatures and up to 32 K, the metal has an antiphase domain structure in which the magnetic moments of the four atomic layers lying in the hcp basal plane are oriented upward, while the magnetic moments of the three closest layers are oriented downward. Thus, there is no variety of spin-slip phases in Tm: they are all transformed into a ferrimagnetic structure that is commensurate with the crystal lattice. As the temperature increases, the antiphase structure begins to break down due to the appearance of the FM component, and in the range of 32-40 K, the magnetic structure is incommensurate with the crystal structure. This intermediate phase contains features of both low-temperature ferrimagnetic ordering and the static longitudinal spin-wave, which is observed in thulium up to 56 K. A number of authors did not observe this intermediate phase in Tm at all [132], which seems to be related to the purity of the material. The axis of the magnetic moment oscillations in the spin-wave is the hexagonal c-axis.
The peculiarities of the magnetic behavior of erbium and thulium, as compared to terbium, dysprosium, and holmium, is that in these metals the magnetic anisotropy tends to orient the magnetic moments along the c-axis. As a result, in Tm, in the FM region below 32 K, the c-axis is the easy one, while in erbium it is the axis of the FM helix below 20 K. Between the Curie and Neel temperatures, the c-axis in these metals is the spin-wave modulation axis. If the field is applied along the c-axis, there is a gradual collapsing of the FM helix (cone) in erbium, which continues in fields above 15 T. [133] As a rule, the magnetic moments are finally oriented parallel to each other in fields above 15 T: in fields above 17 T at helium temperatures, a ferromagnetic structure is realized in erbium [133].
break down due to the appearance of the FM component, and in the range of 32-40 K, the magnetic structure is incommensurate with the crystal structure. This intermediate phase contains features of both low-temperature ferrimagnetic ordering and the static longitudinal spin-wave, which is observed in thulium up to 56 K. A number of authors did not observe this intermediate phase in Tm at all [132], which seems to be related to the purity of the material. The axis of the magnetic moment oscillations in the spin-wave is the hexagonal c-axis. Figure 7. Magnetic phase diagram of Tm. The experimental data has been taken from [118].
The peculiarities of the magnetic behavior of erbium and thulium, as compared to terbium, dysprosium, and holmium, is that in these metals the magnetic anisotropy tends to orient the magnetic moments along the c-axis. As a result, in Tm, in the FM region below 32 K, the c-axis is the easy one, while in erbium it is the axis of the FM helix below 20 K. The analysis of magnetic phase diagrams allows us, despite the above-mentioned differences between the individual elements, to highlight a number of common properties. First of all, this is the uniformity of the crystal structure (hcp). The interatomic distances in all metals are fairly constant, although the ionic radii noticeably decrease with the order number of the element due to lanthanide compression. In addition, in all heavy REM FM and AFM spiral phases coexist (except for gadolinium). The main types of helical structures are presented in [55]. The period of modulated structures usually changes with temperature, taking incommensurable or commensurable values with respect to the period of the crystal lattice. Anisotropy in the basal plane can cause a significant transformation of a simple helical structure and lead to "locking" of the wave vector during the phase transition from incommensurable to commensurable structure. Increasing anisotropy leads to "binding" of individual spins into "bunches" and, in the limit, to a complete collapse of moments to one of the easy directions. The ideal spiral structure, characterized by a constant angle of moment rotation, is a very crude approximation. A more accurate approximation is a model in which the angle is modulated as it moves along the c-axis, not randomly, but in a well-defined way. The energies of these phases are close, since the critical fields in which the collapse of the helical structure are only a few percent of the effective exchange fields.
In five of the considered heavy REMs, different types of ferromagnetic spin ordering were found. In terbium and dysprosium, there is a parallel orientation of moments in the basal plane in the direction of light magnetization. In holmium and erbium, no complete ordering of moments is observed, and the structure has a ferromagnetic helical arrangement (ferromagnetic cone). One of the momentum components is directed parallel to the c-axis and has a ferromagnetic ordering, while the helical configuration is preserved in the basal plane. For the same reasons, in [134] it was proposed to introduce a weak axial component (along the c-axis) for dysprosium. The appearance of this component was observed at temperatures below 100 K (with poor repeatability on three samples); accordingly, the appearance or disappearance of this component may be related to the presence of impurities.
The character of phase transitions is significantly affected by the magnetic field, pressure, defects, and a number of other factors [135][136][137][138]. First-order magnetic phase transitions (FOPT) take place between phases of any symmetry and, in addition to hysteresis phenomena, are usually accompanied by metastable phenomena and the appearance of various intermediate phases. Neutronographic measurements [139][140][141][142] showed that, for example, in holmium, there is a much stronger tendency to spiral ordering than in terbium. Spin waves in erbium and thulium reveal odd higher harmonics of the wave vector as the temperature decreases, corresponding to a symmetric transformation of the sinusoidal structure into an antiphase one. In the framework of the molecular field model it was shown [45], that during the development of a spin wave-type structure, the amplitude of the sinusoidally varying momentum reaches the value inherent in a free ion at T ≈ 0.5T N . Since this value of momentum is maximum, a further increase in energy can be achieved by the development of higher harmonics (conversion of the sinusoid into a periodic sequence of rectangular form, which is observed experimentally).
In the present work, the values of the tricritical temperatures were theoretically calculated for all heavy REMs for the first time. The magnetization curves of heavy REM usually have a maximum at the para-antiferromagnetic transition (the Neel temperature), which is accompanied at lower temperatures by a very sharp increase in magnetization when the ferromagnetic phase is formed, usually through intermediate metastable phases. Because of this general behavior of magnetization, it is correct to apply the relations for determining the critical temperatures obtained in [2]. The results are presented in Table 1. In the case of REM, the critical temperature is the point at which the difference between the PM and FM ordering disappears, that is, the boundary of the FOPT transforms into the boundary of SOPT. By definition, dH cr dT = 0 at T = T cr . For the existence of a tricritical point, the simultaneous existence of three phases is necessary, and the condition that the derivative of the critical parameter is equal to zero in temperature is not always satisfied (a typical example is the triple point of water). The variety of magnetic phase transitions in heavy REM does not exclude the existence of tricritical points in phase diagrams. Moreover, a case of comparison of experimental and theoretical values of the tricritical temperature for dysprosium was found in [73]. It is pointed out that the theoretical value of 168.4 K obtained by the authors is lower than the experimental value of 172 K obtained in the same paper, but is slightly higher than the 165 K obtained in [143] on samples of the worst quality. At the same time, in [73] the experimental value of the tricritical temperature is determined indirectly, by the type of phase boundaries in the magnetic phase diagram. The accuracy of this method is insufficient, so the experimental determination of tricritical temperature values requires the explicit detection of three phases with a physically distinguishable phase boundary. The theoretical value of the terbium tricritical point does not correspond to the definition of the tricritical point, since at 236 K either the transition from FM to PM phase or a completely paramagnetic state is observed in Tb. It can be noted that in all the rest heavy REM the tricritical point is more likely to be in the region of helical AFM ordering or on its boundary with other phases. As in terbium, this region is the narrowest in the series of heavy REMs, the experimental detection of the tricritical point in terbium is of particular interest and difficulty.
The PM-APM transition in heavy REM has the features of both the first-and secondorder phase transitions. Within the molecular field model, the following equation for the transition temperature from the paramagnetic to the antiferromagnetic state can be derived [45]: This formula is valid for any antiferromagnetic structures. If one assumes that the elementary exchange ( → q m ) is constant for the entire series of heavy REMs, then one can conclude that there is a linear dependence of the Neel temperature T N on de Gennes function G = g j − 1 2 J(J + 1)-is the average value of the spin projection on the total mechanical momentum. However, due to the limitations of the molecular field model, the best agreement with the experimental data is observed for the ratio T N ∼ G 2/3 [144]. In the present study, the validity of this relation was experimentally proved for high-purity single crystals of heavy REM (Figure 8). This dependence was first proposed in [145] but tested only in the case of alloys [144,146,147]. We can conclude that in heavy REM single crystals there is a one-to-one functional dependence between the AFM-PM transition temperature and de Genne's function. This fact can be explained by the character of the temperature dependence of the spin rotation angle (the angle between the magnetic moments of the neighboring magnetic layers) in the helical antiferromagnetic structure on de Genne's function. The value of the rotation angle weakly depends on the temperature at small values of the function, with increasing G the change becomes more noticeable, at a certain value the rotation angle decreases sharply to zero at some finite temperature.
At the same time in [84], it is found that the Neel temperature in holmium is split by the magnetic field (similarly, in erbium). This result confirms the effect [135], observed in REMs with a parameter of the order of 4 (holmium, dysprosium, terbium), occurrence in the field.
The diversity and complexity of magnetic phase transitions in heavy REM require a unified theory of mixed magnetic phase transitions, by analogy with the Bean-Rodbell and Kittel models for FOPT [148][149][150] and Landau theory [151] for second-order phase transitions (SOPT). Thus it is reasonable to find a phenomenological parameter that would serve as an indicator of phase transformation in heavy REM. In the case of REM, taking into account the phase transition cascade, which is observed in most of the considered metals, it is probably not quite correct to speak about a specific phase transition temperature, but rather about the region of material phase transformation depending on the magnetic field.
In this sense, the magnetic contribution to the heat capacity Cmagn (T,H) can be interpreted as the order parameter in the theory of first-and second-order magnetic phase We can conclude that in heavy REM single crystals there is a one-to-one functional dependence between the AFM-PM transition temperature and de Genne's function. This fact can be explained by the character of the temperature dependence of the spin rotation angle (the angle between the magnetic moments of the neighboring magnetic layers) in the helical antiferromagnetic structure on de Genne's function. The value of the rotation angle weakly depends on the temperature at small values of the function, with increasing G the change becomes more noticeable, at a certain value the rotation angle decreases sharply to zero at some finite temperature.
At the same time in [84], it is found that the Neel temperature in holmium is split by the magnetic field (similarly, in erbium). This result confirms the effect [135], observed in REMs with a parameter of the order of 4 (holmium, dysprosium, terbium), occurrence in the field.
The diversity and complexity of magnetic phase transitions in heavy REM require a unified theory of mixed magnetic phase transitions, by analogy with the Bean-Rodbell and Kittel models for FOPT [148][149][150] and Landau theory [151] for second-order phase transitions (SOPT). Thus it is reasonable to find a phenomenological parameter that would serve as an indicator of phase transformation in heavy REM. In the case of REM, taking into account the phase transition cascade, which is observed in most of the considered metals, it is probably not quite correct to speak about a specific phase transition temperature, but rather about the region of material phase transformation depending on the magnetic field.
In this sense, the magnetic contribution to the heat capacity C magn (T,H) can be interpreted as the order parameter in the theory of first-and second-order magnetic phase transformations in non-zero fields. Figure 9 summarizes the displacements of the maximum in the magnetic specific heat capacity of heavy REM, where T/T pk is the ratio of the measured temperature to the temperature in which the maximum of the magnetic specific heat at the certain magnetic field is observed. The experimental data are compared with the theoretical calculations (solid black line) taken from Figure 9 of Ref. [152]. One can see that the common trend is observed for all metals; the critical field value is observed, i.e., the value where the maximum begins to shift towards higher temperatures. The following values are experimentally observed: Gd-1.8 T; Tb-0.3 T; Dy-0.3 T; Ho-2 T; Er-2.3 T, comparing with a theoretical value of 0.4 T from [152].
Thus, it can be apparently seen that the field-induced non-monotonic behavior of the maximum in the magnetic heat capacity is an effect common to all of them. Obviously, the appropriate mathematical apparatus of such a theory needs to be elaborated.  [152]. Thus, it can be apparently seen that the field-induced non-monotonic behavior of the maximum in the magnetic heat capacity is an effect common to all of them. Obviously, the appropriate mathematical apparatus of such a theory needs to be elaborated. Figure 9. Position of the maximum in the magnetic specific heat of heavy REM in dependence on the magnetic field. The experimental data for Gd is taken from [72,153], Tb [71], Dy [46,73], Ho [62], Er [154], Tm [155]. Lines drawn through the data points are guides for the eye. The theoretical curve (solid black line) is from [152].

Conclusions
For the first time, the densities of states (DOS) of heavy REMs were calculated using the projector-augmented waves method (PAW). DOS for pure metals do not differ significantly, except for gadolinium. This fact distinguishes it from the series of heavy REM (absence of antiferromagnetic ordering), while at the same time explaining the general laws of the existence of magnetic phases in heavy lanthanides. From this point of view, for the first time the last published magnetic phase diagrams of metals from gadolinium to thulium, three of which were obtained by the authors of this paper, are analyzed together. In low fields, hexagonal anisotropy provides a ferromagnetic ordering of spins; as the field increases, a stable helical antiferromagnetic ordering is observed. Deformation of the spiral results in a fan phase with a helical structure (helifan). In high magnetic fields bulk heavy lanthanides are predominantly in fan or ferromagnetic states. The spin slip type phases arise due to the dominance of anisotropy energy over exchange energy, the existence of structural and chemical defects. For all elements, there must be a theoretical Figure 9. Position of the maximum in the magnetic specific heat of heavy REM in dependence on the magnetic field. The experimental data for Gd is taken from [72,153], Tb [71], Dy [46,73], Ho [62], Er [154], Tm [155]. Lines drawn through the data points are guides for the eye. The theoretical curve (solid black line) is from [152].

Conclusions
For the first time, the densities of states (DOS) of heavy REMs were calculated using the projector-augmented waves method (PAW). DOS for pure metals do not differ significantly, except for gadolinium. This fact distinguishes it from the series of heavy REM (absence of antiferromagnetic ordering), while at the same time explaining the general laws of the existence of magnetic phases in heavy lanthanides. From this point of view, for the first time the last published magnetic phase diagrams of metals from gadolinium to thulium, three of which were obtained by the authors of this paper, are analyzed together. In low fields, hexagonal anisotropy provides a ferromagnetic ordering of spins; as the field increases, a stable helical antiferromagnetic ordering is observed. Deformation of the spiral results in a fan phase with a helical structure (helifan). In high magnetic fields bulk heavy lanthanides are predominantly in fan or ferromagnetic states. The spin slip type phases arise due to the dominance of anisotropy energy over exchange energy, the existence of structural and chemical defects. For all elements, there must be a theoretical tricritical point (the point of coexistence of the three magnetic phases), the experimental detection of which remains outside the scope of this paper and requires a separate study. In this case, the theoretical value of the tricritical temperature monotonically decreases from terbium to thulium. The linear dependence of the Neel temperature of high-purity bulk heavy REMs on de Gennes function G 2/3 was shown for the first time. In all heavy rare earth metals, there is a common, triggered by an external magnetic field, non-monotonic response of the maximum magnetic heat capacity with a minimum of the maximum temperature shift dependence at some value of the critical field, different for each element. The magnetic heat capacity can be considered as an order parameter in developing a theory of phase transitions in a magnetic field.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.

Conflicts of Interest:
The authors declare no conflict of interest.