Room-Temperature Ferromagnetism in Two-Dimensional Fe2Si

Apr 25, 2017 - †CAS Key Laboratory of Materials for Energy Conversion, School of Chemistry and Materials Sciences, and CAS Center for Excellence in ...
0 downloads 11 Views 1MB Size
Subscriber access provided by UB + Fachbibliothek Chemie | (FU-Bibliothekssystem)

Communication

Room-Temperature Ferromagnetism in 2D Fe2Si Nanosheet with Enhanced Spin-Polarization Ratio Yingjie Sun, Zhiwen Zhuo, Xiaojun Wu, and Jinlong Yang Nano Lett., Just Accepted Manuscript • Publication Date (Web): 25 Apr 2017 Downloaded from http://pubs.acs.org on April 25, 2017

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

Nano Letters is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 23 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Nano Letters

Room-Temperature Ferromagnetism in 2D Fe2Si Nanosheet with Enhanced Spin-Polarization Ratio Yingjie Sun,#,† Zhiwen Zhuo,#,† Xiaojun Wu,*,†,‡,§ and Jinlong Yang ‡, § †

CAS Key Laboratory of Materials for Energy Conversion, School of Chemistry and Materials

Sciences, and CAS Center for Excellence in Nanoscience, University of Science and Technology of China, Hefei, Anhui 230026, China. ‡

Hefei National Laboratory for Physical Science at the Microscale, University of Science and

Technology of China, Hefei, Anhui 230026, China. §

Synergetic Innovation Center of Quantum Information & Quantum Physics, University of

Science and Technology of China, Hefei, Anhui 230026, China KEYWORDS. First-principles calculations, two-dimensional Fe2Si, ferromagnetism

ABSTRACT. Searching experimental feasible two-dimensional (2D) ferromagnetic crystals with large spin-polarization ratio, high Curie temperature and large magnetic anisotropic energy is one key to develop next-generation spintronic nano-devices. Here, 2D Fe2Si nanosheet, one counterpart of Hapkeite mineral discovered in meteorite, with novel magnetism is proposed on the basis of first-principles calculations. 2D Fe2Si crystal has a slightly buckled triangular lattice

ACS Paragon Plus Environment

1

Nano Letters 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 2 of 23

with planar hexacoordinated Si and Fe atoms. The spin-polarized calculations with hybrid HSE06 function method indicate that 2D Fe2Si is a ferromagnetic half metal at its ground state with 100% spin-polarization ratio at Fermi energy level. The phonon spectrum calculation and ab-initio molecular dynamic simulation shows that 2D Fe2Si crystal has a high thermodynamic stability and its 2D lattice can be retained at the temperature up to 1200K. Monte Carlo simulations based on the Ising model predict a Curie temperature over 780 K in 2D Fe2Si crystal, which can be further tuned by applying a biaxial strain. Moreover, 2D structure and strong inplane Fe-Fe interaction endow Fe2Si nanosheet sizable magnetocrystalline anisotropy energy with the magnitude of at least two orders larger than those of Fe, Co and Ni bulks. These novel magnetic properties render the 2D Fe2Si crystal a very promising material for developing practical spintronic nano-devices

INTRODUCTION Spintronics, which uses the spin of electrons for information storage, transportation and processing, has attracted intensive interests from both science and industry in the past decades.1 In tradition, ferromagnetic crystals with large spin-polarization ratio, high Curie temperature, and large magnetocrystalline anisotropic energy (MAE) are fundamental to build practical spintronic devices that could work at room temperature2-6. To realize this concept at nanoscale, it is one key issue to develop low-dimensional ferromagnetic crystals that possess the above qualities. However, searching experimentally feasible low-dimensional ferromagnetic materials still remains a big challenge for the low yield, small domains and rich boundaries/defects of prepared samples, or the difficulty in obtaining crystalline structures.7-9

ACS Paragon Plus Environment

2

Page 3 of 23 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Nano Letters

Since the discovery of graphene, intensive research efforts have been devoted to develop spintronics based on two-dimensional (2D) crystals. 10-14 Graphene have attracted many interests for its gate tunable spin transport,15 long spin diffusion lengths of about 4 micros at room temperature.16 However, significant magnetism is rarely observed in experiment due to the delocalization of s and p orbitals in carbon materials. Alternatively, 2D inorganic graphene-like crystals, such as transition metal dichalcogenides17-19, oxides20,22, nitrides23,24, carbide25-27, halide28-29, silicide30-33, and complex ternary compounds34-39, present great opportunities in spintronics applications for the possible existence of d orbital itinerant magnetism transition metal atoms and coupling between them. Here, we propose a new 2D Fe2Si nanosheet, one counterpart of Fe2Si alloy, with roomtemperature ferromagnetism, enhanced spin-polarization ratio and sizeable MAE, suitable for spintronics application in the nanoscale. Fe2Si, also named as hapkeite mineral, is one unique member of Fe-Si alloys, which have many technologically important compounds with excellent soft magnetic properties, e.g. Fe3Si, Fe2Si, Fe5Si3, FeSi, β-FeSi2, and α-FeSi2.40-48 Fe2Si is the compound has to be formed under the space weathering condition on airless or meteorite impacts. Therefore, it was predicted that this compound should exist on the lunar surface, and was indeed found in a lunar meteorite named as Dhofar 280 in 2003.49 Experimentally, single phase Fe2Si alloy can be obtained by a rapid quenching technique.50 With novel high-saturation magnettization, low coercivity, high permeability, low magnetostrictive coefficient, low magnetocrystalline anisotropy, and possible half-metallic like band structure, Fe2Si has great potential applications in high-power transformers, microwaves, motors, and high-frequency applications etc. 51,52

ACS Paragon Plus Environment

3

Nano Letters 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 23

Based on first-principles calculations and ab-initio molecular dynamic (AIMD) simulation, we show that graphene-like 2D Fe2Si is a ferromagnetic crystal with high structural stability up to the temperature of 1200 K. The Curie temperature estimated with Monte Carlo simulation based on Ising model reaches 590 K, which can be further enhanced by applying biaxial strain. In particular, 2D Fe2Si crystal is half metal with 100% spin polarization ration at Fermi energy and its MAE is two orders higher than those of Fe, Co, and Ni, enabling 2D Fe2Si crystal a very promising candidate for spintronics at nanoscale.

COMPUTATIONAL METHODS Our first-principles calculations are performed on the basis of density functional theory (DFT) by using Vienna ab initio Simulation Package (VASP).53 Projector augmented-wave (PAW) pseudopotential54 is used to account electron–ion interactions. The exchange and correlation terms are described using general gradient approximation (GGA) in the scheme of Perdew-Burke-Ernzerhof (PBE) functional.55 Test calculations with PBEsol functional present similar results as those of PBE functional. To deal with the strong-correlated effect on d-orbital electrons of Fe, the screened hybrid Heyd-Scuseria-Ernzerhof (HSE06) functional56,57 is also used. A vacuum separation larger than 15 Å along the z direction is employed in order to avoid interactions between two Fe2Si images in nearest-neighbor unit cell. A cutoff energy of 500 eV is used to expand the wave functions in plane waves. Both the lattice constant and positions of all atoms are relaxed until the force is less than 10-2 eV/Å. The convergence criterion of the total energy is 10−5eV. A Γ-centered Monkhorst-Pack scheme58 with k-point mesh of 15 × 15 × 1 is used. The phonon dispersion spectrum is calculated with finite displacement method method59 on 4×4×1 supercell with 48 atoms by using PHONOPY package.60 The AIMD simulation61with a

ACS Paragon Plus Environment

4

Page 5 of 23 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Nano Letters

canonical ensemble is performed by using the algorithm of Nóse and the temperature of the system is maintained at 300, 600, 900, 1200, or 1500 K, respectively. The time step is 1.0 fs and the total simulation time is 10.0 ps.

RESULTS AND DISCUSSION 2D Fe2Si sheet can be obtained by extracted one stoichiometric layer from hexagonal Fe2Si crystal, as illustrated in Figure 1. The hexagonal Fe2Si crystal has layered structure with P3-M1 symmetric group (Group 146), composing uneven graphene-like honeycomb sublattice of FeSi intercalated with hexagonal sublattice of Fe atoms, as shown in Figure 1a and 1b. The optimized lattice constants with PBE calculations are a=b=3.93 Å, and c=4.84 Å, agreeing well with previous theoretical results.48,52 However, these values are significantly different from experimental values of a=b=4.05 Å , and c=5.09 Å.48 Note that GGA fails to correctly describe the structural and electronic properties of the strong-correlated system, such as underestimation of the lattice constants, local magnetic moments and band gaps, which can be improved by using HSE06 method.62 The optimized lattice constants of Fe2Si bulk with HSE06 method are a=b=4.07 and c=5.11 Å, agreeing well with experimental values.52 The spinpolarized electronic band structures of Fe2Si calculated with PBE method exhibit half-metalliclike property, as shown in Figure S1 (see Supporting Information), which is consistent with previous theoretical study at the same calculation level.52 Considering the strong-correlated effect, however, the electronic band structure calculated with the screened HSE06 hybrid function indicates that the hexagonal Fe2Si is a ferromagnetic metal, as shown in Figure 1c. Based on the calculated density of states (DOS) projected on atomic orbitals in Figure 1c, the magnetism is mainly contributed by d orbital of Fe atoms, and the local magnetic moments are

ACS Paragon Plus Environment

5

Nano Letters 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 23

2.92, 2.35, and 2.27 µB for Fe atoms in two hexagonal sublattices and honeycomb sublattice in the unit cell by using HSE06 method, respectively. These values are significantly larger than those obtained with PBE method (2.57, 0.10, and 0.70 µB).52 Note that The spin polarization ratio is defined as δ=|ρ↑-ρ↓|/(ρ↑+ρ↓), where ρ↑ and ρ↓ are the absolute values of DOS at the Fermi energy level for spin-up and down channels, respectively. The calculated spin-polarization ratio is about 19.5% at the Fermi energy level by using HSE06 method.

Figure 1. (a) and (b) are the top and side views of the optimized structure of hexagonal Fe2Si bulk, respectively. (c) The spin-polarized band structures and density of states of hexagonal Fe2Si bulk calculated with HSE06 functional. (d) The top and side views of the optimized structure of 2D Fe2Si crystal. The DOS of major and minor spin polarized channel are labelled with black up and down arrows, respectively. The fermi energy level is set as zero. Γ (0.0, 0.0, 0.0), A (0.0, 0.0, 1/2), H (-1/3, 2/3, 1/2), K (-1/3, 2/3, 0.0), M (0.0, 1/2, 0.0), L (0, 1/2, 1/2) refer to the high-symmetric k-points in the first Brillouin zone.

ACS Paragon Plus Environment

6

Page 7 of 23 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Nano Letters

Extracting one stoichiometric layer from hexagonal Fe2Si bulk, the optimized structure of 2D Fe2Si crystal is displayed in Figure 1d. Distinctly different from the structure of stoichiometric layer in bulk, 2D Fe2Si crystal has a slightly buckled structure with same symmetry group at its bulk counterpart (P3-M1). Quasi-planar hexa-coordinated Fe and Si atoms are packed in a triangular lattice where each Si atom is surrounded by six neighboring Fe atoms and Fe atoms form a graphene-like sublattice, sharing the similar structure as those of metal carbide and silides.63-66 Using HSE06 method, the optimized lattice constants increase to a=b=4.24 Å, where the Fe-Si and Fe-Fe bond lengths are 2.46 and 2.51 Å. Note that the optimized lattice constants are a=b=3.97 Å, and the Fe-Si and Fe-Fe bond lengths are 2.33 and 2.43 Å, respectively by using PBE method, which are significantly different from HSE06 method due to the strong-correlated effect of d electrons in system. Test calculations with PBE+U method with U=3.5 eV present similar results as HSE06 method, as summarized in Table S1. The thickness of 2D Fe2Si is about 0.56 Å at HSE06 calculation levels. To ascertain the structural stability of 2D Fe2Si crystal, the cohesive energy is calculated with the definition of Ecoh = (ESi + 2×EFe – EFe2Si)/3, where ESi, EFe, and EFe2Si are the total energies of Si, Fe atom, and 2D Fe2Si unitcell, respectively. The calculated cohesive energy is 4.10 eV per atom. This value is smaller than that of graphene (7.85 eV/atom), but is significantly larger than those of silicene (3.98 eV/atom), germane (3.26 eV/atom),65 and planar hexa-coordinated Cu2Si (3.46 eV/atom)64 and Cu2Ge (3.17 eV/atom)66, respectively, implying the possible stability of 2D Fe2Si crystal. The thermodynamic stability of 2D Fe2Si is evaluated with phonon spectrum calculation and AIMD simulation. Figure 2a displays the calculated phonon dispersion of 2D Fe2Si crystal without any presence of imaginary-vibration mode. This indicates that 2D Fe2Si crystal with

ACS Paragon Plus Environment

7

Nano Letters 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 23

buckled structure is dynamically stable. Test calculation is performed on 2D Fe2Si crystal with perfect planar structure and significant imaginary-vibration mode can be found. The AIMD simulation is further performed on 2D Fe2Si at the temperatures of 300, 600, 900, 1200, and 1500 K. Figure 2b to 2d display the snapshot of structures at 10.0 ps at the temperatures of 900, 1200, and 1500 K, respectively. And, the final structures of 300 and 600K are provided in Figure S2. It is clear that 2D Fe2Si retains its structures at the temperature up to 1200 K, and the structure is partially destroyed at the temperature of 1500 K in the simulating period, indicating that the structure of 2D Fe2Si sheet is thermally stable at high temperature at least up to 1200 K.

Figure 2. (a) The calculated phonon spectrum of 2D Fe2Si crystal. (b) to (d) are the structural snapshots of 2D Fe2Si crystal at time of 10 ps during AIMD simulation under the temperatures of 900, 1200, and 1500 K, respectively. Γ (0.0, 0.0, 0.0), M (0.0, 1/2, 0.0), K (−1/3, 2/3, 0.0) refer to the high-symmetric k-points in the first Brillouin zone. In the following, the electronic properties of 2D Fe2Si crystal are investigated. Two magnetic orders, i.e. ferromagnetic (FM) and antiferromagnetic (AFM) coupling orders, are considered to

ACS Paragon Plus Environment

8

Page 9 of 23 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Nano Letters

determine the magnetic ground state of 2D Fe2Si. Our calculations show that 2D Fe2Si has a FM ground state. The AFM and FM energy difference is 0.41 eV per cell. Figure 3a displays the calculated electronic band structures and DOS of 2D Fe2Si crystal at FM ground state by using HSE06 method. 2D Fe2Si crystal is half metal with 100% spin-polarization ratio at the Fermi energy level, which is over 5 times larger than that of hexagonal Fe2Si bulk (19.5%) at the HSE06 calculation level. The enhanced spin-polarization ratio endows 2D Fe2Si crystal great advantages in spintronics application at nanoscale.

Figure 3. (a) The band structures, total and partial DOS projected on atomic orbitals of 2D Fe2Si crystal are calculated by using HSE06 method. The fermi energy level is set to zero. “up” and “down” arrows denote the spin-up and down polarization, respectively. (b) The deformation charge density distribution with the isosurface value of 0.007 a.u. The yellow and green denote the charge density increase and decrease, respectively. (c) The spin charge density distribution with the isosurface value of 0.03 a.u. (d) The simulated magnetic moment (M) and specific heat (Cv) with respect to temperature for 2D Fe2Si crystal.

ACS Paragon Plus Environment

9

Nano Letters 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 23

The major magnetism in 2D Fe2Si crystal is contributed by Fe atoms, as illustrated by the spin charge density distribution in Figure 3c. The local magnetic moment on Fe atom is about 3.03 µB, which is larger than that in its bulk. The partial DOS projected on atomic orbitals indicate that the spin-polarization at the Fermi energy level is mainly contributed by dxz and dyz orbitals of Fe atoms. To further understand this novel magnetism behavior, crystal field theory is used to analysis the electronic occupation of d orbital in 2D Fe2Si crystal. Note that Fe atoms are hexcoordinated, the energies of five d orbitals in Fe atom are splitted into a(d2), e1(dxz, dyz), and e2(dxy, dx2-y2), as schemed in the insert of Figure 3a. A freestanding Fe atom has a valence electronic configuration of 3d64s2. The deformation charge density, shown in Figure 3b, and Bader charge analysis result indicate that each Fe atom losses about 0.28 electron to Si atom. Consequently, Fe atom is approximately +1 charged, and the left seven electrons occupy a, e1, and e2 orbitals with a and e1 orbital half-filled in a high spin state, resulting enhanced local magnetic moment of 3.07 µB and spin-polarization ratio. Table 1. Magnetic anisotropy energies in µeV/Fe and the easy axis (EA) for 2D Fe2Si crystal are summarized. E(100)-E(001) E(010)-E(001) E(110)-E(001) E(111)-E(001) EA 560

574

550

325

(001)

MAE is an important parameter of magnetic materials to determine the low-temperature magnetic orientation with respect to the lattice structure, which is directly related to the thermal stability of magnetic data storage. Large values of MAE per atom imply the possibility to reduce the grain size per bit of information. Note that reduced dimensionality and symmetry would significantly enhance the MAE value of material, a sizeable MAE is expected in 2D Fe2Si crystal

ACS Paragon Plus Environment

10

Page 11 of 23 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Nano Letters

when comparing with its bulk counterpart. In addition, spin-orbital coupling (SOC) calculations are performed on 2D Fe2Si crystal to obtain the relative stability along three magnetization directions in plane, i.e. (100), (010), and (111) direction, and two directions out of plane, i.e. (001) and (111) directions, as summarized in Table 1. The easy axis (EA) of 2D Fe2Si is along the out-of-plane (001) direction. The calculated MAEs for (100), (010), (110), and (111) directions are 0.560, 0.574, 0.550, and 0.325 meV per Fe atom based on PBE calculation level, respectively, which are about one to two orders larger than those of cubic Fe (1.4 µeV per atom), Co( 65 µeV per atom), and Ni (2.7 µeV per atom).67,68 Moreover, the magnitudes are comparable to those of Fe, Co monolayer deposited on Rh(111) and Pt(111) substrates (0.08 to 0.37 meV per metal atom).69 The MAE values are also calculated with HSE06 method, which is more expensive in computation. The MAE values with HSE06 method are smaller than that obtained with PBE method, which are 172, 198, 331 and 225µeV per Fe atom for (100), (010), (110) and (111) directions, respectively. The large MAE values make 2D Fe2Si suitable for magneto electronics applications. To develop a practical spintronics working at ambient environment, the Curie temperature (TC) of materials should be comparable to or above the room temperature. Since the easy axis of 2D Fe2Si is out-of-plane (001) direction and the calculated MAE is comparable to previous reported 2D CrXTe3,70 which behaviors as a 2D Ising-like ferromagnet,71 the Kosterlize-Thouless transition is unlikely to take place in Fe2Si sheet.71 Here, Monte Carlo (MC) simulation based on Ising model is used to estimate the TC of 2D Fe2Si crystal with FM state, which has been used to predict TC of 2D magnetic crystal in previous simulation.6, 37, 70 The Hamiltonian of Ising model is defined as H = − ∑ 

(1)

ACS Paragon Plus Environment

11

Nano Letters 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 23

, where Jij is the nearest-neighboring exchange parameter and M=3.0 is the spin magnetic moment on Fe atom. Using the energy difference between AFM and FM ground states based on HSE06 method, the estimated exchange parameter J is about 7.7 meV. To obtain Tc, we calculate the specific heat Cv at first after the system reaches equilibrium at a given temperature. A supercell of 60 × 60 × 1 with periodic boundary condition is used in MC simulation. In the simulations, the spins on all magnetic sites flip randomly. The Tc is obtained by locating the peak position in Cv(T) plot. From the simulated Cv(T) curve in Figure 3d, the Curie temperature of 2D Fe2Si is estimated to be 780 K. Based on PBE results of energy difference between AFM an FM states, the exchange parameter J is about 8.5 meV and the estimated TC reaches 860 K. It should be pointed out that the MAE in 2D Fe2Si is large but not infinite as assumed in Ising model. As a result, the calculated temperature with Ising model is likely overestimated. Generally, to fabricate 2D crystal based electronics, a high in-plane stiffness is required to avoid the curling and obtain free-standing membranes. Here, 2D Young’s modulus is calculated to estimate the in-plane stiffness of 2D Fe2Si. -The calculated in-plane Young’s modulus is 71 Nm-1, which is about 24% of in-plane Young’s modulus of graphene (295 Nm-1). Note this value is larger than that of silicene (61 Nm-1), germanene (42 Nm-1), and comparable to those of Cu2X and Ni2X (X=Si, and Ge) sheet (60 to 93 Nm-1).64-66 Then, the gravity induced out-of-plane deformation can be estimated by the equation of 

h/L ≈ (ρgL/ )

(2)

, where ρ = 1.7 ×10-6 kg/m2 is the density of 2D Fe2Si crystal, and L is the size the nanosheet. Taking L≈100 µm, we obtain h/L≈2.84×10-4, comparable to that of graphene.73 Therefore, it is possible that 2D Fe2Si crystal keeps free-standing planar structure by withstanding its weight.

ACS Paragon Plus Environment

12

Page 13 of 23 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Nano Letters

Figure 4. (a) and (c) are the band structures and total DOS of 2D Fe2Si crystal under the biaxial strain of -3% and 3% by using HSE06 method, respectively. The Fermi energy level is set to zero. (b) and (d) are the simulated magnetic moment (M) and specific heat (Cv) with respect to temperature for 2D Fe2Si crystal under the biaxial strain of -3% and 3%, respectively. At last, the strain effect on the electronic and magnetic properties of 2D Fe2Si crystal is considered. Figure 4 displays the calculated band structures and DOS of 2D Fe2Si crystal under biaxial strain of -3% and 3%. Other results are displayed in Figure S3 (see Supporting Information). The biaxial deformation within the range from -3% to 3% won’t change the FM half metallicity of 2D Fe2Si crystal. Meanwhile, the energy difference between AFM and FM states varies from 0.501 eV to 0.289 eV, as summarized in Table S2. Consequently, the Curie temperature is enhanced to 820 K when the biaxial strain is -3%, but decreases to 680 K when the biaxial strain reaches 3%. Such behavior change with the biaxial strain can be understood with the structural deformation of 2D Fe2Si crystal. As summarized in Table S1, when the biaxial strain changes from -3% to 3%, the Fe-Si and Fe-Fe bond lengths increases

ACS Paragon Plus Environment

13

Nano Letters 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 23

monotonously from 2.39 to 2.53 and 2.44 to 2.57 Å, respectively, which result the couple interaction between neighbor Fe atoms. Experimentally, it is still a big challenge to obtain 2D Fe2Si nanosheet directly by exfoliating Fe2Si bulk. Note that Fe2Si ultrathin films have been obtained in experiment,51 and Fe2Si ultrathin film containing one stoichiometric layer spontaneously transforms into 2D Fe2Si nanosheet, it is possible to obtain 2D Fe2Si sheet by reducing the thickness of Fe2Si ultrathin film to the atomic scale, which has been used to obtain graphene-like ZnO nanosheet.74 Alternatively, it is still possible to obtain 2D Fe2Si nanosheet by using atomic layer deposition (ALD) and molecular beam epitaxy (MBE) techniques, which have been used to prepare 2D inorganic materials, such as transition metal disulfides, on selected substrates.75-77

CONCLUSIONS In summary, we presented a new 2D Fe2Si crystal with room-temperature ferromagnetism and enhanced spin-polarized ratio by using first-principles calculations. 2D Fe2Si crystal can be looked as a stoichiometric layer of its bulk counterpart. Distinctly different from Fe2Si bulk, quasi-planar hexa-coordinated Si and Fe atoms are packed in a slightly buckled triangular lattice in 2D Fe2Si crystal. Each Si atom is surrounded by six neighboring Fe atoms, and Fe atoms form a graphene-like structure. Both phonon spectrum calculation and AIMD simulation indicate that 2D Fe2Si crystal has a high thermodynamic stability, and can retain its lattice structure at the temperature up to 1200 K. Spin-polarized calculations with hybrid HSE06 function show that 2D Fe2Si crystal is FM half metal, and the spin-polarization ratio is over 5 times larger than that of its bulk counterpart. 2D Fe2Si crystal possesses both high Curie temperature of about 780 K, and

ACS Paragon Plus Environment

14

Page 15 of 23 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Nano Letters

a sizeable MAE, which is one to two orders larger than those of Co, Fe, and Ni bulk. These novel magnetic properties provide great opportunities for spintronics and electronics devices.

ASSOCIATED CONTENT Supporting Information. Supporting Information Available: The spin-polarized band structures of hexagonal Fe2Si calculated with PBE functional; the structural snapshots of 2D Fe2Si crystal of AIMD simulation under the temperatures of 300 and 600 K, respectively; the band structures of 2D Fe2Si crystal under a biaxial strain; the bond length, local magnetic moment, and AFM-FM energy differences of deformed 2D Fe2Si crystal by applying biaxial strain. This material is available free of charge via the Internet at http://pubs.acs.org. AUTHOR INFORMATION Corresponding Author * [email protected] Author Contributions #

These authors contributed equally.

Notes The authors declare no competing financial interest. ACKNOWLEDGMENT

ACS Paragon Plus Environment

15

Nano Letters 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 16 of 23

This work is supported by the NSFC (21573204, 21421063, 51172223), MOST (2016YFA0200602), Strategic Priority Research Program of CAS (XDB01020300), the Fundamental Research Funds for the Central Universities (WK2060140014), National Key Basic Research Program (2012CB922001), National Program for Support of Top-notch Young Professional, External Cooperation Program of BIC CAS (211134KYSB20130017), and by USTCSCC, SCCAS, Tianjin, and Shanghai Supercomputer Centers

REFERENCES 1.

Fert, A. Rev. Mod. Phys. 2008, 80, 1517-1530.

2.

Meservey, R.; Tedrow, P. M. Solid State Commun. 1972, 11, 333-336.

3.

Chen, J.-R.; Odenthal, P. M.; Swartz, A. G.; Floyd, G. C.; Wen, H.; Luo, K. Y.; Kawakami, R. K. Nano Lett. 2013, 13, 3106−3110.

4.

Wang, D.; Chen, Q.; Xing, G.; Yi, J.; Bakaul, S. R.; Ding, J.; Wang, J.; Wu, T. Nano Lett. 2012, 12, 3994-4000.

5.

Seixas, L.; Carvalho, A.; Neto, A. H. C. Phys. Rev. B 2015, 91. 155138.

6.

Li, X. X.; Wu, X. J. ; Yang, J. L. J. Am. Chem. Soc. 2014, 136, 5664−5669.

7.

Abel, M.; Clair, S.; Ourdjini, O.; Mossoyan, M.; Porte, L. J. Am. Chem. Soc. 2011, 133, 1203-1205.

8.

Bieri, M.; Blankenburg, S.; Kivala, M.; Pignedoli, C. A.; Ruffieux, P.; Muellen, K.; Fasel, R. Chem. Commun. 2011, 47, 10239-10241.

ACS Paragon Plus Environment

16

Page 17 of 23 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Nano Letters

9.

Lee, J. S.; Wang, X.; Luo, H.; Dai, S. Adv. Mater. 2010, 22, 1004-1007.

10. Mas-Balleste, R.; Gomez-Navarro, C.;Gomez-Herrero,J.; Zamora, F. Nanoscale, 2011, 3, 20. 11. Tang, Q.; Zhou, Z,; Chen, Z. F. WIREs Comput Mol Sci. 2015, 5,360–379. 12. Nicolosi, V.; Chhowalla, M.; Kanatzidis, M. G.; Strano, M. S.; Coleman, J. N. Science 2013, 340, 1226419. 13. Xu, M. S.; Liang, T.; Shi, M. M.; Chen, H. Z. Chem. Rev. 2013, 113, 3766-3798. 14. Li, X. L.; Wu, X. J. WIREs Comput. Mol. Sci. 2016, 6, 441-455. 15. Cho, S.; Chen, Y. F.; Fuhrer, M. S. Appl. Phys. Lett. 2007, 91, 123105. 16. Tombros, N.; Jozsa, C.; Popinciuc, M.; Jonkman, H. T.; van Wees, B. J. Nature 2007, 448, 571-574. 17. Ataca, C.; Sahin, H.; Ciraci, S. J. Phys. Chem. C 2012, 116, 8983-8999. 18. Kan, M.; Adhikari, S.; Sun, Q. Phys. Chem. Chem. Phys. 2014,16, 4990-4994. 19. Ma, Y. D.; Ying Dai, Y.; Guo, M.; Niu, C. W.; Zhu, Y. T.; Huang, B. B. ACS Nano. 2012, 6, 1695-1701. 20. Kan, M.; Zhou, J.; Sun, Q.; Kawazoe, Y.; Jena, P. J. Phys. Chem. Lett. 2013, 4, 33823386. 21. Kan, E. J.; Li, M.; Hu, S. L.; Xiao, C. Y.; Xiang, H. J.; Deng, K. M. J. Phys. Chem. Lett. 2013, 4, 1120−1125.

ACS Paragon Plus Environment

17

Nano Letters 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 18 of 23

22. Tang, Z. K.; Li, X. B. ;Zhang, D. Y.; Zhang, Y. N.; Liu, L. M. J. Mater. Chem. C 2015, 3, 3189-3197. 23. Zhang, X. M.; Yu, Z. M.; Wang, S. S.; Guan, S.; Yang, H. Y.; Yao, Y. G.; Yang, S. Y. A. J. Mater. Chem. A 2016, 4, 15224-15231. 24. Zhang, S. H.; Li, Y. W.; Zhao, T. S.; Wang, Q. Sci. Rep. 2014, 4, 5241. 25. Ji, X.; Xu, K.; Chen, C.; Zhang, B.; Wan, H. Z.; Ruan, Y. J.; Miao, L.; Jiang, J. J. J. Mater. Chem. A 2015, 3, 9909-9914. 26. Si, C.; Zhou, J.; Sun, Z. M. ACS Appl. Mater. Interfaces 2015, 7, 17510-17515. 27. Khazaei, M.; Arai, M.; Sasaki, T.; Chung, C. Y.; Venkataramanan, N. S.; Estili, M.; Sakka, Y.; Kawazoe, Y. Adv. Funct. Mater. 2013, 23, 2185-2192. 28. Torun, E.; Sahin, H.; Singh, S. K.; Peeters, F. M. Appl. Phys. Lett. 2015,106, 192404. 29. Lin, S. H.; Kuo, J. L. Phys. Chem. Chem. Phys. 2014, 16, 20763-20771. 30. Bui, V. Q.; Pham, T. T.; S. Nguyen, H. V. S.; Hung M. Le, H. M. J. Phys. Chem. C 2013, 117, 23364-23371. 31. Dzade, N. Y.; Obodo, K. O.; Adjokatse, S. K.; Ashu, A. C.; Amankwah, E.; Atiso,C. D.; Bello, A. A.; Igumbor, E.; Nzabarinda, S. B.; Obodo, J. T.; Ogbuu, A. O.; Femi, O. E.; Udeigwe, J. O.; Waghmare, U. V. J. Phys.: Condens. Matter 2010, 22,375502 . 32. Liu, Z. F.; Wang, X. Q.; Cai, J. T.; Zhu, H. J. J. Phys. Chem. C 2015, 119, 1517-152. 33. Han, N. N.; Liu, H. S.; Jijun Zhao, J. J. J. Supercond. Nov. Magn. 2015, 28, 1755-1758.

ACS Paragon Plus Environment

18

Page 19 of 23 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Nano Letters

34. Sivadas, N.; Daniels, M. W.; Swendsen, R. H.; Okamoto, S.; Xiao, D. Phys. Rev. B 2015, 91, 235425. 35. Chen, X. F.; Qi, J. S.; Shi, D. N. Phys. Lett. A 2015, 379, 60-63. 36. Zhuang, H. L. L.; Xie, Y.; Kent, P. R. C.; Ganesh, P. Phys. Rev. B 2015, 92, 035407. 37. Li, X. X.; Wu, X. J.; Yang, J. L. J. Am. Chem. Soc. 2014, 136, 11065-11069. 38. Sachs, B.; Wehling, T. O.; Novoselov, K. S.; Lichtenstein, A. I.; Katsnelson, M. I. Phys. Rev. B 2013, 88, 201402. 39. Tang, Z. K.; Liu, W. W.; Zhang, D. Y.; Lau, W. M.; Liu, L. M. RSC Adv. 2015, 5, 7715477158. 40. Bauer, F.; Hiltl, M.; Rappengluck, M. A.; Neumair, A.; Ernstson, K. Meteorit. Planet. Sci. 2013, 48, A52-A52. 41. Khalaff, K.; Schubert, K. J. Less-Common Metals 1974, 35, 341-345. 42. Shinjo, T.; Nakamura, Y.; Shikazono, N. J Phys. Soc. Japan 1963, 18, 797-801. 43. Hines, W. A.; Menotti, A. H.; Budnick, J. I.; Burch, T. J.; Litrenta, T.; Niculescu, V.; Raj, K. Phys. Rev. B 1976, 13, 4060-4068. 44. Christensen, N. E. Phys. Rev. B 1990, 42, 7148-7153. 45. Takarabe, K.; Teranishi, R.; Oinuma, J.; Mori, Y.; Suemasu, T.; Chichibu, S.; Hasegawa, F. Phys. Rev. B 2002, 65, 165215. 46. Girlanda, R.; Piparo, E.; Balzarotti, A. J. Appl. Phys. 1994, 76, 2837-2840.

ACS Paragon Plus Environment

19

Nano Letters 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 20 of 23

47. Sandalov, I.; Zamkova, N.; Zhandun, V.; Tarasov, I.; Varnakov, S.; Yakovlev, I.; Solovyov, L.; Ovchinnikov, S. Phys. Rev. B 2015, 92. 205129. 48. Kudielka, H. Z. Kristallogr. -Cryst. Mater. 1977, 145, 177–189. 49. Anand, M.; Taylor, L. A.; Nazarov, M. A.; Shu, J.; Mao, H. K.; Hemley, R. J. Proc. Natl. Acad. Sci. USA 2004, 101, 6847-6851. 50. Varga, L. K.; Mazaleyrat, F.; Kovac, J.; Greneche, J. M. J. Phys.: Condens. Matter 2002, 14, 1985-2000. 51. Chen, Y. T.; Tan, Y. C. J. Alloys Comp. 2014, 615, 946-949. 52. Tang, C. P.; Tam, K. V.; Xiong, S. J.; Cao, J.; Zhang, X. AIP Adv. 2016, 6, 065317. 53. Kresse, G.; Hafner, J. Phys. Rev. B 1993, 48, 13115-13118. 54. Blöchl, P. E. Phys. Rev. B 1994, 50, 17953-17979. 55. Perdew, J. P.; Burke, K.; Ernzerhof, M. Phys. Rev. Lett. 1996, 77, 3865-3868. 56. Heyd, J.; Scuseria, G. E.; Ernzerhof, M. J. Chem. Phys. 2003, 118, 8207-8215. 57. Heyd, J.; Scuseria, G. E.; Ernzerhof, M. J. Chem. Phys. 2006, 124. 219906. 58. Monkhorst, H. J.; Pack, J. D. Phys. Rev. B 1976, 13, 5188-5192. 59. Ackland, G. J.; Warren, M. C.; Clark, S. J. J. Phys.: Condens. Matter 1997, 9, 7861–7872. 60. Togo, A.; Oba, F.; Tanaka, I. Phys. Rev. B 2008, 78. 134106.

ACS Paragon Plus Environment

20

Page 21 of 23 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Nano Letters

61. Bucher, D.; Pierce, L. C. T.; McCammon, J. A.; Markwick, P. R. L. J. Chem. Theory Comput. 2011, 7, 890–897. 62. Marsman, M.; Paier, J.; Stroppa, A.; Kresse, G. J. Phys.: Condens. Matter 2008, 20, 064201. 63. Li, Y. F.; Liao, Y. L.; Chen, Z. F. Angew. Chem. Int. Ed. 2014, 53, 7248-7252. 64. Yang, L. M.; Bacic, V.; Popov, I. A.; Boldyrev, A. I.; Heine, T.; Frauenheim, T.; Ganz, E. J. Am. Chem. Soc. 2015, 137, 2757-2762. 65. Yang, L. M.; Popov, I. A.; Frauenheim, T.; Boldyrev, A. I.; Heine, T.; Baciac, V.; Ganz, E. Phys. Chem. Chem. Phys. 2015,17, 26043-26048. 66. Yang, L. M.; Popov, I. A.; Boldyrev, A. I.; Heine, T.; Frauenheim, T.; Ganz, E. Phys. Chem. Chem. Phys. 2015, 17, 17545-17551. 67. Daalderop, G. H. O.; Kelly, P. J.; Schuurmans, M. F. H.; Jansen, H. J. F. J. Phys. Colloques 1988, 49, C8-93-C8-94. 68. Daalderop, G. H. O.; Kelly, P. J.; Schuurmans, M. F. H. Phys. Rev. B 1990, 41, 1191911937. 69. Lehnert, A.; Dennler, S.; Blonski, P.; Rusponi, S.; Etzkorn, M.; Moulas, G.; Bencok, P.; Gambardella, P.; Brune, H.; Hafner, J. Phys. Rev. B 2010, 82, 094409. 70. Zhuang, H. L.; Xie, Y.; Kent, P. R. C.; Ganesh, P. Phys. Rev. B 2015, 92, 035407. 71. Carteaux, V.; Moussa, F.; Spiesser, M. Europhys. Lett. 1995, 29, 251-256.

ACS Paragon Plus Environment

21

Nano Letters 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 22 of 23

72. Kosterlitz, J. M.; Thouless, D. J. J. Phys. C: Solid State Phys. 1973, 6, 1181. 73. Booth, T. J.; Blake, P.; Nair, R. R.; Jiang, D.; Hill, E. W.;Bangert, U.; Bleloch, A.; Gass, M.; Novoselov, K. S.; Katsnelson, M. I.; Geim, A. K. Nano Lett. 2008, 8, 2442-2446. 74. Tusche, C.; Meyerheim, H. L.; Kirschner, J. Phys. Rev. Lett. 2007, 99, 026102. 75. Johnson, R. W.; Hultqvist, A.; Bent, S. F. Mater. Today 2014, 17, 236-246. 76. Dasgupt, N. P.; Lee, H. B. R.; Bent, S. F.; Weiss, P. S. Chem. Mater. 2016, 28, 1943-1947. 77. Ugeda, M. M.; Bradley, A. J.; Shi, S. F.; Jornada, F. H. da; Zhang, Y.; Qiu, D. Y.; Ruan, W.; Mo, S. K.; Hussain, Z.; Shen, Z. X.; Wang, F.; Louie, S. G.; Crommie, M. F. Nat. Mater. 2014, 13, 1091-1095.

ACS Paragon Plus Environment

22

Page 23 of 23 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Nano Letters

SYNOPSIS Table of Contents (TOC)

ACS Paragon Plus Environment

23