Weifu Cen, Lin Lyu, Yinye Yang1, * and Meihui Fan
(1. Guidance Center of Undergraduate Innovation and Entrepreneurship, Guizhou Minzu University, Guiyang 550025, China;2. Special and Key Laboratory of Guizhou Provincial Higher Education for Photoelectronic Information Processing and Analysis,Guizhou Minzu University, Guiyang 550025, China)
Abstract: A system study of the three-dimensional normal stress for regulating electronic structure and magnetic property of Fe2Ge is studied. The density states of Fe more than 92% contribution come from Fe 3d, the density states of Ge mainly contributed from Ge 4p and Ge 4s, and the Fe 3d spin induces the Ge 4p electron transfer. The inductive effect increases germanium electron energy, weakens the Fe spin density of states, opposes the stability of the ferromagnetic state. The magnetic moment varies from 5 to 3 μB with the stress charges from -30 to 30 GPa. The charge of Fe is negative whereas the Ge atom is positively charged, the Fe atom loses charge, the charge transfers to the Ge atom. The unevenly distributed charge forms the new occupy state and spin polarization state in the Fe2Ge electron structure system. The Fe is the electron donor, the total electron is transferred to Ge, but the total numbers of gain electron and total numbers of lost electron are not equal, so the Fe2Ge electron system may have hybridization between the Fe 3d state and Ge 4p state. The magnetic of Fe2Ge mainly comes from the unoccupied Fe 3d orbital, the Fe 3d is positive spin-polarization state and the spin-polarization strength is decreased, the Ge 4p is negative spin-polarization state and the spin-polarization strength are increased. Moreover, electrons-spin polarization is relevant to the structure parameters of the Fe2Ge system, and controls spin-polarized electronic behavior by means of adjusting ferromagnetic.
Keywords: Fe2Ge; electronic structure; magnetic property; stress controlling; first-principles
The CASTEP is used for structure optimization, where the electron-core interactions are described by the projector augmented wave method for the potentials, and the exchange correlation energy was calculated within the generalized gradient approximation (GGA) of the Perdew-Burke-Ernzerhof form[23]. The cutoff energy for the plane wave basis set used to expand the Kohn-Sham was 350 eV for all calculations. Structural relaxations and total energy calculations were performed ensuring that the Hellmann-Feynman forces acting on ions were less than 10-2eV/? for all studied structures.
Since iron is a transitional metals element with atomic number 26 and its outer shell (3d64s2) contains 3d electrons, the GGA approaches fail to describe the strongly correlated localized 3d electrons and predicting a metallic ground state for the iron germanide. Thereby, to account for the strong on-site coulomb repulsion among the localized 3d electron in Fe2Ge, we used a Hubbarad-U constant as described by Obukuro[23]. The GGA + U cases, represented by the Hubbard-like term U and the exchange term J led to an improvement of the ground state properties such as the band gap and the magnetic moments. In all cases investigated, the density of states of bulk Fe2Ge was calculated and a good agreement found the optimum U constant was 3. Here, we research the magnetism of Fe2Ge using the method of GGA+U with the three-dimensional normal stress regulating the magnetism.
The investigated Fe2Ge compound is a hexagonal structure with lattice constants ofa=b=0.369 1 nm,c= 0.540 1 nm, the space group is P63/mmc (No.194)[21], the structure presenting the two-dimensional hexagonal structure with the Fe atom at their centers, the molecular structure diagram was shown in Fig.1.
Fig.1 The molecular structure of Fe2Ge
The mechanical behavior of solids is usually defined by constitutive stress-strain relations, the stress and strain follow Hook’s law if the stress is a linear function of strain. Strainεis characterized the lattice constant along with the change of stressσ, the deformation energy[24]uas follows:
The relationship between total energy and stress (E-S) of Fe2Ge was shown in Fig.2. The total energy increased with the tensile stress and compressive stress increasing, the rock bottom of the E-S curve appeared at 0 GPa, the E-S curve was not symmetric to the point of 0 GPa, those varies were similar with the change of elastic potential energy, the elastic potential energy was minimum, the Fe2Ge magnetic system was the most stable under the stress of 0 GPa. The elastic potential energy is the macroscopic representation of molecular potential in the perspective of physics, the E-S curve indicates the molecules potential energy is minimum under 0 GPa, the change rate of E-S curve with compressive is less than that with tensile stress.
Fig.2TherelationshipbetweentotalenergyandstressofFe2Ge
The relationship between strain and lattice constants (L-S) of Fe2Ge was shown in Fig.3. The strain of lattice constants ofa,bandcpresented the quasi linear change increased with stress increasing in the stress range of -30-0 GPa and 0-30 GPa. However, the strain change trend of lattice constantsaandbare the same, differ fromc, the change rate of lattice constantscis greater than lattice constants ofaandb. The most significant changes occurred in the stress range of 0 and 10 GPa, the lattice constantsaandbare rapidly decreased and have a turning point at 6 GPa, the lattice constantcis decreased in the stress range of 0-6 GPa and is slightly increased in the stress range of 6-10 GPa. The Poisson ratio of Fe2Ge is 0.914 22, so Fe2Ge is an elastic material. In this paper, the strain of lattice constant is changed in the range of 0- 4.382 84%, the strain is almost linearly with stress change, the stress range is in the elastic range.
The integrated spin is defined as the difference between spin up DOS and spin down DOS, the|integrated spin| is defined as totals of spin up DOS and spin down DOS, the relationship between strain and integrated spin density (I-S) of Fe2Ge was shown in Fig.4. |integrated spin| and integrated spin are almost linearly decreased with the stress increasing from -30 GPa to 30 GPa, this change is similar to the change of lattice constant, the results mean that the changing of lattice constant is one of the important reasons of |integrated spin| and integrated spin change. But the maximum change between the curve of I-S and the curve of L-S appeared at the stress range of 0-10 GPa, in this stress range the curve of I-S is decreased rapidly. The cause of change was a pressing need to better understand the density states changed with stressing.
Fig.3TherelationshipbetweenstrainandstressofFe2Ge
Fig.4 The integration spin density of Fe2Ge
The partial density of states (PDOS) of the atoms is calculated to further investigate the interaction between the Fe atom and Ge atom. The PDOS of Fe2Ge deduced from band structure calculations and PDOS of states for the Fe and Ge are presented in Fig.5, under the stress of 0 GPa. Here, the Fermi level is taken as a reference. The Fe PDOS is dominated by the Fe 3d orbital, whereas the PDOS on Ge atoms show only small contributions from Ge(4s) and Ge(4p) orbitals. The magnetic moment is mainly dominated by the 3d orbitals of Fe atom. The structure asymmetry exists in the maximum PDOS of Fe 3d orbitals of spin-up and spin-down near the Fermi level, and the Fe 3d spin up DOS was different with the Fe-d spin down DOS, this asymmetrical spin density states represented the magnetic of Fe2Ge. In the magnetic system of Fe2Ge, the increasing energy of electron weakens the spin density of states for Fe and opposed the stability of the ferromagnetic state. Importantly, the significant difference existed in the spin-up and spin-down configuration of Ge-p, the spin states estimated came from the Ge-p electron transfer by Fe-d spin induced. The magnetic moment of Fe2Ge is computed as well and the magnetic moment about 4 μB per Fe2Ge molecule (the studied system has 2 Fe2Ge molecule). Closing to the value of 2 μB per Fe atom quoted for bulk Fe2Ge, the results correspond to Ref.[19]. The spin orbit, the charge transfer, and the orbital electric are influenced by the stress. For Fe 3d orbital, the spin up density is chiefly concentrated in -8-0 eV, the spin down density is chiefly concentrated in -5-3 eV. The Fe-d density state change on the stress range of -30 - 30 GPa was present in Fig. 6. The highest unoccupied state of spin up moves to the high energy direction, and the occupied state of spin down moves to the low energy direction in the process of the stress increasing from -30 GPa to 30 GPa, this variation depends mainly on atom characteristic in the system of Fe2Ge. When the pressure is 2 GPa and 6 GPa, the electron density varies with energy increasing which is different with the variation under other stress. This may be the main causes of the integrated spin and the lattice constant occurring significant change. For further understanding the variation characteristics of the magnetic moments with the stress changing, we analyzed the electronic characteristic and the interaction between charge and electrons from the atoms population, in Fe2Ge electronic structure system.
(a) Spin partial density of states for Fe
(b) Spin partial density of states for Ge
Fig.5ThespinpartialdensityofstatesfortheFeandGe
(a) The Fe-d density state change on the stress range of -30~-2 GPa
(b) The Fe-d density state change on the stress range of 2~-30 GPa
Fig.6TheFe-ddensitystatechangeonthestress
Mulliken atomic charge calculation plays an important role in the molecular system because atomic charges affect spin electronic structure and much more electronic properties. The total atomic charges of Fe2Ge obtained by atoms population analysis, the atomic charge of Fe is increased whereas the Ge atoms is decreased, the Ge atom occupies the higher positive value in the process of the stress increasing from -30 GPa to 30 GPa. In Fig.7, the Fe atom lost the charge, the charge transferred to the Ge atom. The amount of transferred charge for Fe is decreased,the amount of gained charge for Ge is increased. This results mean that the unevenly distributed charge formed the new occupy state and spin polarization state.
Fig.7 The charge population
In the Fe2Ge electron system, the total electron of Fe more than 92% of the contribution comes from Fe 3d, the electron of Ge mainly contributes from Ge 4p and Ge 4s. Fig.8 shows the total electron values for the Fe is decrease, for the Ge is increased in the process of the stress increasing from -30 to 30 GPa. The Fe is the electron donor, the total electron is transferred to Ge, but the total numbers of gain and lose electron is not equal. So the Fe2Ge electron system exists hybridization between the Fe 3d state and Ge 4p state, suggesting a strong s-p coupling interaction between them. Such a strong s-p interaction makes the energy levels spin split with the compressive stress increasing, this results led to the system energy of the compressive stress lower than the energy of the tensile stress (as shown in Fig.1).
The Ge 4p orbital is an occupied, in nonmagnetic, the magnetic comes from the unoccupied Fe 3d orbital, but the increasing energy of electron weakens the spin density of states for Fe and opposed the ferromagnetic state, strengthens the spin density of states for Ge and applauded the ferromagnetic state. In Fig.9, the Fe 3d is positive spin-polarization state and the spin-polarization strength is decreased, the Ge 4p is negative spin-polarization state and the spin-polarization strength is increased in the process of the stress increasing from -30 to 30 GPa. The total magnetic moment for the ferromagnetic states is very close to the expected value of 2.0 μB per Fe since the Fe 3d minority states are unoccupied. The partially filled spin-down band split off from the fully occupied spin-up band via Hund exchange. So the stress is one of the important factors for affecting the electrons-spin polarization. Moreover, electrons-spin polarization is relevant to the structure constants of the Fe2Ge system, and control spin-polarized electronic behavior by means of adjusting ferromagnetic.
Fig.8 The total electron
Fig.9 The spin
In summary, using first-principle calculations we investigated the effect of electronic structures and magnetic property of Fe2Ge by stress regulating. In this paper we consider the stress range is in the elastic range from -30 to 30 GPa, because Iron is a transitional metal, the outer shell (3d64s2) contains 3d electrons, using the method of GGA+U calculated magnetism, using a Hubbarad-U constant described the strongly correlated. The density states of Fe more than 92% contribution come from Fe 3d, the density states of Ge mainly contribute from Ge 4p and Ge 4s, and the Fe 3d spin is induced Ge 4p electron transfer by analysis of their electronic structure. The magnetic moment of Fe2Ge is computed as well and the magnetism moment is almost linearly decrease from 5 to 3 μB per Fe2Ge molecule with the stress increasing from -30 to 30 GPa. The charge of Fe is negative whereas the Ge atoms are positively charged, the Fe atom lost charge, the charge transferred to the Ge atom, the unevenly distribute charge formed the new occupy state and spin polarization state in the Fe2Ge electron structure system. But the total numbers of gain and lose electron is not equal, so the Fe2Ge electron system is exist hybridization between the Fe 3d state and Ge 4p state, the magnetic come from the unoccupied Fe 3d orbital, the Fe 3d is positive spin-polarization state and the spin-polarization strength is decrease, the Ge 4p are negative spin-polarization state and the spin-polarization strength are increase in the process of the stress increasing from -30 to 30 GPa.
Journal of Harbin Institute of Technology(New Series)2018年6期