Strain-Spintronics: Modulating Electronic and Magnetic Properties of

2 days ago - Products. Journals A–Z · eBooks · C&EN · C&EN Archives · ACS Legacy Archives · ACS Mobile · Video. User Resources. About Us · ACS ...
0 downloads 0 Views 2MB Size
Subscriber access provided by UNIV OF LOUISIANA

C: Physical Processes in Nanomaterials and Nanostructures

Strain-Spintronics: Modulating Electronic and Magnetic Properties of HfMnCO MXene by Uniaxial Strain 2

2

2

Edirisuriya M. D. Siriwardane, Pragalv Karki, Yen Lee Loh, and Deniz Cakir J. Phys. Chem. C, Just Accepted Manuscript • DOI: 10.1021/acs.jpcc.9b00594 • Publication Date (Web): 26 Apr 2019 Downloaded from http://pubs.acs.org on April 26, 2019

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 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 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.

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 27 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

The Journal of Physical Chemistry

Strain-Spintronics: Modulating Electronic and Magnetic Properties of Hf2MnC2O2 MXene by Uniaxial Strain Edirisuriya M. D. Siriwardane, Pragalv Karki, Yen Lee Loh, and Deniz Çakir∗ Department of Physics and Astrophysics, University of North Dakota, Grand Forks, North Dakota 58202, USA E-mail: [email protected] Abstract Next-generation spintronic nanoscale devices require two-dimensional (2D) materials with robust ferromagnetism. Among 2D materials, MXenes are favorable for spintronic applications due to their high electron conductivity and mobility. A recently reported MXene, Hf2 MnC2 O2 , possesses a high Curie temperature (greater than 800 K) and a high magnetic moment per formula unit (3 µB ). Since 2D materials have greater elastic strain limits than their bulk counterparts, their properties can be tuned effectively using strain engineering. Here, we investigate modifications in the structural, electronic, and magnetic properties produced by uniaxial strain on a Hf2 MnC2 O2 monolayer. The strain-free Hf2 MnC2 O2 nanosheet is an indirect-band-gap semiconductor. Our calculations predict that an indirect-to-direct band gap transition occurs at about 1%-3% tensile strain applied in the armchair direction. At 7% strain applied in the zigzag direction and 9% strain applied in the armchair direction, this semiconductor material becomes a half-metal which is favorable for spintronic applications. Under 4% compressive strain in either direction, a semiconductor-to-metal transition is predicted.

1

ACS Paragon Plus Environment

The Journal of Physical Chemistry 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

More importantly, the Curie temperature of the material can be enhanced significantly by applying tensile strain. For instance, under 8% strain, the Curie temperature becomes greater than 1200 K.

Introduction Spin-polarized transport occurs in ferromagnetic materials, where the density of states of spin-up and spin-down electrons are shifted relative to each other near the Fermi level. 1 Spintronics, which utilizes spin-polarized transport to manipulate spin degrees of freedom, is being investigated intensively in order to develop current- and next-generation technologies such as computer memory units and quantum computing. 1,2 In ferromagnetic materials, a Curie temperature higher than the room temperature is essential to maintain its intrinsic magnetization above the room temperature. 3 High carrier mobility along with enhanced spin-orbit coupling in 2D materials make them potential candidates for spintronic applications. 4 MXenes are 2D transition metal carbides, carbonitrides, or nitrides with the general formula Mn+1 Xn Tx , where M represents an early transition metal, X stands for carbon and/or nitrogen, and T is surface termination group. 5 MXenes are made by etching out the A layers of MAX phases where A is a group-A element such as Al. Since Ti3 C2 Tx was introduced as the first ever MXene, 6 around 20 different MXenes have been synthesized including Ti2 CTx and Hf3 C2 O2 . 7,8 It has been reported that the double transition metal MXenes Ti2 MnC2 (OH)2 , Hf2 MnC2 O2 , and Hf2 VC2 O2 , and also the nitride MXenes Mn2 NF2 , Mn2 NO2 , Mn2 N(OH)2 , and Cr2 NO2 , are ferromagnetic materials exhibiting Curie temperatures well above room temperature. 5,9 Introducing approaches to control the magnetic states in those robust ferromagnetic materials is crucial in quantum information devices. 9,10 Strain engineering can be effectively applied to 2D materials due to the high elastic strain limits of nanosheets as compared to their bulk counterparts. 10–12 Strain engineering has been exploited to tune numerous properties of 2D materials like graphene, transition metal dichalcogenides (TMDs), MXenes, 2

ACS Paragon Plus Environment

Page 2 of 27

Page 3 of 27 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

The Journal of Physical Chemistry

and phosphorene. 10,12–18 Enhanced magnetic moments have been observed in T-VS2 15 and WTe2 19 TMDs and also in Ti2 C and V2 N MXenes 20 under biaxial tensile strain. It has been reported that ferromagnetism in Li-doped TiO2 nanosheets can be tuned via uniaxial strain. 21 The ferromagnetic half metallic Ti2 NO2 2D materials can be converted into an antiferromagnetic semiconductor by applying tensile strain greater than 6%. 22 It is also reported that the magnetic properties of 2D CrX2 (X=S, Se, Te) and p-doped PtSe2 can be modified by using strain engineering. 23,24 The band gap engineering of 2D materials using uniaxial strain has been extensively studied in recent years. HSE06 studies show that Hf3 C2 O2 is a semiconductor with a 0.155 eV band gap. At around −2.5 % compressive strain, this monolayer shows metallic properties, whereas a band gap of 0.56 eV can be realized at a strain value of 5.0%. 25 Ti2 CO2 monolayer undergoes an indirect-to-direct band gap transition under biaxial strain of 4% and the uniaxial strain of 6%. 26 Similarly, the band gap of black phosphorus, PtSe2 and CrX (X=S, Se, Te) monolayers can be significantly modified by exploiting uniaxial strain. 23,27,28 Ferromagnetic semiconductors are useful as sources of spin-polarized carriers 29 for spintronic applications. Among magnetic monolayer materials with high Curie temperatures, Hf2 MnC2 O2 is a semiconductor, where it has been already proven that biaxial strain can be employed to tune the electronic and magnetic properties significantly. 30 In this research work, we applied uniaxial compressive and tensile strains on a Hf2 MnC2 O2 monolayer in two different directions, namely zigzag and armchair, to manipulate the electronic, and magnetic properties. Our Density Functional Theory (DFT) calculation show that the conduction bands are highly sensitive to strain. Due to the conduction band shifting under uniaxial strain, semiconductor-to-metal phase transition occurs at 4% compression in both directions, at 7% tension in zigzag direction, and at 9% tension in the armchair direction. The dispersion of bands at the lower part of conduction band also changes as a result of band shifting, which alters the effective mass of the electrons. Moreover, we carried out Monte Carlo simulations to investigate the variation of Curie temperature under uniaxial strain.

3

ACS Paragon Plus Environment

The Journal of Physical Chemistry 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

We found that uniaxial tensile strain can be utilized to enhance the Curie temperature in Hf2 MnC2 O2 . For instance, the Curie temperature increases from 851 K for the strain-free material to greater than 1200 K at 8% tensile strain.

Computational Method Our first-principles calculations were based on spin-polarized density functional theory as implemented in the Vienna ab initio Simulation Package (VASP) 31–34 with Perdew-BurkeErnzerhof (PBE) 35,36 pseudopotentials in the Generalized Gradient Approximation (GGA) method. A plane wave basis set was considered with Projected Augmented Wave (PAW) method 37,38 with energy threshold value of 400 eV in the calculations. The criteria for energy and force convergence of the ionic steps were set to 10−6 eV and to 10−2 eV/Å, respectively. A vacuum space of around 20 Å was used in the c-axis direction in order to avoid interaction between two neighboring images. The Brillouin zone was sampled on a Gamma-centered 12×6× 1 Monkhorst-Pack k-mesh. The Density of States (DOS) calculations were performed using a 24 ×12×1 k-mesh. The GGA with on-site Coulomb potential method (GGA+U) 39 was used to treat the strong on-site Coulomb interaction of localized d electrons of the transition metal atoms. The U potentials are 2 eV for Hf and 4 eV for Mn which were taken from a previous. 5 In order to infer the strength of exchange couplings between Mn atoms, different super-cell sizes were considered as described in Section . HSE06 40 calculations for several strain values were done to verify the validity of the GGA+U calculations. The same magnetic ground state can be obtained from both GGA+U and HSE06 methods. Density Functional Perturbation Theory (DFPT) as implemented in VASP 41 was employed for elastic constant (Cij ) calculations. The dynamical stability was investigated by calculating phonon dispersions using PHONOPY code. 42

4

ACS Paragon Plus Environment

Page 4 of 27

Page 5 of 27 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

The Journal of Physical Chemistry

Results and Discussion Strain and Structure

Figure 1: (a)Top and (b)side views of Hf2 MnC2 O2 monolayer. Table 1: The elastic constants C11 , C22 , C12 , monolayer thickness (t) and 3D Young’s modulus (Y3D ) values under -4%, 0%, 4% and 8% uniaxial strains. Strain (%) C11 (GPa) C22 (GPa) C12 (GPa) t (nm) Y3Dx (GPa) Y3Dy (GPa)

-4 429.7 394.3 91.8 0.911 408.3 374.6

zigzag (x) 0 4 8 385.2 328.1 258.5 384.9 373.4 347.8 88.2 82.4 82.5 0.901 0.894 0.888 364.9 309.9 238.9 364.7 352.7 321.4

-4 407.9 386.5 104.7 0.910 379.6 359.6

armchair (y) 0 4 8 385.2 312.9 214.7 384.9 381.7 382.1 88.2 82.4 83.6 0.901 0.894 0.887 364.9 295.1 196.4 364.7 359.9 349.5

In order to apply uniaxial strain, the rectangular cell structure for Hf2 MnC2 O2 shown in Fig. 1 was considered. The GGA+U calculations give the optimized monolayer lattice parameters as a0 = 3.253 Å and b0 = 5.634 Å. As can be seen in Fig. 1 (b), Mn is sandwiched 5

ACS Paragon Plus Environment

The Journal of Physical Chemistry

between two O-Hf-C layers where O is the surface termination. Thus, there are four Hf, four O, four C, and two Mn atoms in the rectangular simulation cell where Hf makes bonds with both C and O atoms while Mn makes bonds only with C atoms. The computed Hf-C, Mn-C and Hf-O bond lengths are 2.319, 2.201, and 2.115 Å, respectively. 2.6

(a)

2.6

Hf1-C1 Hf2-C1 Mn1-C1 Mn2-C1 Hf1-O1 Hf2-O1

2.5

Bond Length (Å)

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 27

2.5

2.4

2.4

2.3

2.3

2.2

2.2

2.1

2.1

2 -4 -2

0

2

4

6

Strain (%)

8

(b)

10

2 -4 -2

0

2

4

6

Strain (%)

8

10

Figure 2: Variation of the bond lengths under uniaxial strain in the (a) x and (b) y directions. Strain is applied in the zigzag (x ) or armchair (y ) directions as x =

a−a0 a0

and y =

b−b0 b0

where a and b are the lattice constants of the strained monolayer in the x and y directions, respectively. Due to the symmetry, studying Hf1-O1, Hf2-O1, Hf1-C1, Hf2-C2, Mn1-C1 and Mn2-C1 is sufficient for investigating the variation of the bond lengths due to the applied uniaxial strain. As can be seen in Fig. 2(a), transition metal (TM)-C bond lengths change differently according to the direction of uniaxial strain. For instance, Hf1-C1 and Mn1C1 increase under tensile strain in the y direction whereas those quantities decrease under tensile strain in the x direction. On the contrary, Hf2-C1 and Mn2-C1 bond lengths rise (decline) as tensile strain is applied in the x (y) direction. These distinct behaviors exist due to the Poisson effect where the structure contracts (expands) in one direction when a tension (compression) is applied in the perpendicular direction. Our calculations show that the Poisson ratio is around 0.22, which indicates that the bond length change in the opposite direction due to the strain is very small. Despite the fact that TM-C bond lengths 6

ACS Paragon Plus Environment

Page 7 of 27

are behaving differently as above, both Hf-O bond lengths increase (decrease) under uniaxial tensile (compressive) strain applied in the zigzag and armchair directions. The strain energy (Es ) under uniaxial strain is calculated using Es = E − E0 where E is the total energy of the strained material and E0 is the total energy of the strain-free material. The stress is calculated as a function of strain using the equation σ =

1 ∂E V0 ∂i

where

V0 is the volume of the strain-free monolayer and i is the strain in the direction i (x or y). Here, the thickness (6.961 Å) was considered as the distance between two oxygen layers. As was illustrated by Fig. 3, the monolayer is able to handle higher stress values when the uniaxial strain is applied in the zigzag direction. This is due to the fact that TM-C bond lengths do not increase much at higher tensile strains in the x direction compared to that in the y direction as illustrate in Fig. 2. Thus, at higher strains in the y direction, bond lengths become weaker. The ultimate stress in the x direction is 46.3 GPa while that in the y direction is 35.7 GPa. Moreover, the corresponding ultimate strain is 18% in both directions. 50

(b)

(a) 40

εx εy

3

Stress (GPa)

4

Strain Energy (eV)

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

The Journal of Physical Chemistry

2

1

0 -4

30

20

10

0

4

8

12

0

16

Strain (%)

0

5

10

15

20

Strain (%)

Figure 3: (a) Strain energy per formula unit and (b) stress as a function of uniaxial strain in the zigzag and armchair directions.

The surface Young’s modulus Y2D in the x and y directions for a rectangular structure can be calculated using Eq. 1 and Eq. 2 where the elastic constants C11 and C22 are the axial compression moduli in the x and y directions, respectively, and C12 is the transverse

7

ACS Paragon Plus Environment

The Journal of Physical Chemistry 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 27

expansion modulus: 11 Y2Dx =

C11 C22 − C212 C22

(1)

Y2Dy =

C11 C22 − C212 C11

(2)

Since the thickness changes under uniaxial strain (see TABLE 1), we normalized Y2D as Y3D = Y2D /t and the elastic constants as Cij = Cij(2D) /t where i and j are 1 and 2, respectively. The thicknesses are calculated by considering multilayer configurations with van der Waals interactions as defined by DFT-D3 with Becke-Jonson damping. 43,44 Thus, TABLE 1 provides the Young’s modulus values and the elastic constants in GPa. At 0% strain, C11 = C22 = 385 GPa and Y3Dx = Y3Dy = 365 GPa due to the hexagonal symmetry of the strain-free Hf2 MnC2 O2 monolayer. When uniaxial strain is applied, those quantities are no longer equal (C11 6= C22 and Y3Dx 6= Y3Dy ), since the symmetry is broken. At 8% strain in either direction, there is a significant difference between x and y components of Y3D , where the latter is higher than the former. Thus, stiffness will be low in the zigzag direction, when a high uniaxial tension is applied (Y3Dx = 238.9 GPa at x = 8% and Y3Dx = 196.4 GPa at y = 8% ). Even though stiffness declines in both directions under tension due to bond weakening, there is still a remarkably high stiffness (Y3Dy > 310 GPa) in the armchair direction at 8% strains compared to that in the zigzag direction. As shown in Fig. S1 ( see the Supplementary Material), we calculated the ab-initio phonon calculations along the high symmetry directions of Brillouin zone for large strain values (10 and 18%) in order to elucidate the dynamical stability Hf2 MnC2 O2 under strain. It is clear that the Hf2 MnC2 O2 monolayer is stable at strain values as large as 10%. Imaginary vibrational frequencies in x = 18% indicates that the nanosheet cannot retain its stability at much higher strains like 18%, even though Fig. 3 indicates that the material’s ultimate strain in both directions is 18%.

8

ACS Paragon Plus Environment

Page 9 of 27 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

The Journal of Physical Chemistry

Uniaxial Strain Effects on Electronic Properties The band structures were obtained by considering the high-symmetry points (see Fig. 4 for an unstrained case) of the rectangular supercells, i.e. Γ (0,0,0), X (0.5,0,0), S (0.5,0.5,0) and Y (0,0.5,0), under uniaxial strains applied in the zigzag (x ) and armchair (y ) directions as shown in Fig. 5. Strain-free Hf2 MnC2 O2 is a semiconductor with an indirect band gap of 0.282 eV. The conduction band minimum (CBM) resides at the S point, while the valence band minimum (VBM) resides at the Γ point. As strain is applied, the lower three bands (see Fig. 4), namely CB1, CB2, and CB3, shift considerably as illustrated in Fig. 5. DOS calculations reveal that TM-d and O-p states hybridize to form CB1, while TM-d orbitals together with p orbitals of Hf, C, and O hybridize to construct CB2 and CB3. CB1 band is occupied by the spin-up electrons while CB2 and CB3 bands are occupied by the spin-down electrons, as shown in Fig. 4.

Figure 4: Spin-resolved electronic band structure and partial density of states (DOS) for strain free rectangular super-cell of Hf2 MnC2 O2 monolayer. Spin-up and spin-down bands are shown with red and blue colors, respectively. The Brillouin Zone for the rectangular super-cell is included with high symmetric points. Fermi energy marks zero energy. TABLE 2 reveals that Hf2 MnC2 O2 monolayer remains an indirect-band-gap semiconductor under tensile strain in the zigzag direction. While tensile strain rises in this direction, CB1 does not move noticeably with respect to the top of the valence band until strain of 4% 9

ACS Paragon Plus Environment

The Journal of Physical Chemistry

Figure 5: The calculated spin-resolved band structure of Hf2 MnC2 O2 for different uniaxial strain values in (a) zigzag and (b) armchair directions. Spin-up and spin-down bands are in red and blue, respectively. Fermi energy marks zero energy.

Band Gap (eV)

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 27

εx εy

0.2

0.1

0 -4

-2

0

2

4

6

8

10

Strain (%)

Figure 6: Variation of band gap as a function of strain. (see Fig. 5). Thus, band gap remains between 0.28 eV and 0.24 within x =0% to x =4% range, as shown in Fig. 6. After x = 4%, the bottom of CB2 becomes CBM, which appears at Y point. Since CB2 moves towards the valence band as tensile strain rises further in the zigzag direction, there is a sudden band gap drop after x = 4%. The CB2 and valence band start crossing the Fermi level at x = 7%. At x = 7%, only the bands with spin-up electrons cross the Fermi level. As a consequence, Hf2 MnC2 O2 becomes a half-metal, where spin-down

10

ACS Paragon Plus Environment

Page 11 of 27

(a)

εx = -4%

2

εx = 0%

εx = 4%

εx = 8%

PDOS

1 0 -1 -2

(b)

Hf-p Mn-p

Hf-d Mn-d -2

0

2

-2

εy = -4%

2

0

C-p O-p 2

-2

εy = 0%

0

2

-2

εy = 4%

0

2

εy = 8%

1

PDOS

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

The Journal of Physical Chemistry

0 -1 -2

-2

0

2

-2

0

2

-2

0

2

-2

0

2

Energy (eV)

Figure 7: The calculated Partial Density of States (PDOS) of Hf2 MnC2 O2 for different uniaxial strain values in (a) zigzag and (b) armchair directions. Fermi energy marks zero energy. electron bands make a band gap of 0.61 eV. The partial density of states (PDOS) studies show that even though the main contributor in conduction bands of spin-up electrons near the Fermi level is Mn-d orbitals at 0% strain, Mn-d, Hf-d and C-p contribute under tensile strains. The contributions from each orbital for the valence bands near the Fermi level barely changes as strain is applied (see Fig. 7). Even though the band gap of Hf2 MnC2 O2 changes in a similar way under the strain in the armchair direction (see Fig. 6), the CBM appears at the different k-points as illustrated in TABLE 2. For the strains greater than 3% in this direction, the CBM resides between Γ and X points since CB3 reaches the Fermi level. Similar to the case of zigzag direction, the spinup electron bands are populated near the Fermi level at higher strains. Thus, Hf2 MnC2 O2 turns into a half-metal at y = 9% with a band gap of 0.82 eV made by the bands with spin-down electrons. It is clear that there is a sharp drop in band gap when compressive strain is applied on

11

ACS Paragon Plus Environment

The Journal of Physical Chemistry 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 27

Table 2: The position of the conduction band minimum (CBM) for each 2D material under uniaxial strain in the zigzag (x) and armchair (y) directions. VBM resides always at the Γ point. Strain (%) -3 -2 -1 0 1 2 3 4 5 6 7 8

CBMx CBMy Γ S Γ S Γ S S S S Γ S Γ S Γ S Γ-X Y Γ-X Y Γ-X Γ-X Γ-X

Table 3: The average Bader charge transfer ∆q for each atom are mentioned for strained and unstrained Hf2 MnC2 O2 . Property ∆qHf (e) ∆qO (e) ∆qC (e) ∆qM n (e)

-4% 3.937 -2.001 -2.687 1.504

-2% 3.937 -2.007 -2.674 1.489

zigzag (x) 0% 2% 4% 3.939 3.938 3.943 -2.011 -2.011 -2.023 -2.669 -2.660 -2.652 1.480 1.466 1.464

8% 3.861 -2.058 -2.635 1.426

-4% 3.939 -2.004 -2.690 1.511

-2% 3.939 -2.007 -2.680 1.496

armchair (y) 0% 2% 3.939 3.936 -2.011 -2.009 -2.669 -2.659 1.480 1.463

4% 3.815 -2.014 -2.521 1.439

8% 3.859 -2.093 -2.474 1.417

the monolayer in either direction. CB1 forms the CBM which appears at the Γ point from x = y = -1% to -3%. CB1 moves downward and crosses the Fermi-level at x = y = -4%. Since a very small portion of the bottom of the conduction band and the top of the valence bands crosses the Fermi level, Hf2 MnC2 O2 is a semi-metal at x = y = -4%. Bottom of CB1 and VBM move towards each other to reduce the separation between the anti-bonding and bonding states. Since Hf-O bond lengths become shorter at negative strains, TM-d and O-p bonding enhances. Materials like Hf2 MnC2 O2 monolayers require special attention due to the presence of partially filled d orbitals. 45 The Hartree electron-electron repulsion potential in DFT implies how an electron interacts with the electron-density of the system. Since the electron-density contains the charge of that electron as well, an unphysical self interaction energy is present 12

ACS Paragon Plus Environment

Page 13 of 27 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

The Journal of Physical Chemistry

in the total energy. 46,47 This self interaction error is significant in the transition metal atoms where the unpaired electrons are delocalized to minimize the self-interaction. Our GGA calculations in Fig. S2 show that all the strained Hf2 MnC2 O2 monolayers are metallic systems (see Supplementary Material) as a result of delocalization of unpaired d-electrons. To overcome this issue, GGA+U method utilizes the on-site Coulomb potential (U), which localizes the transition metal electrons in the partially filled energy states. Furthermore, HSE06 calculations constitute the nonlocal Hartree-Fock exchange potential to minimize the effects of so called self interaction. 48,49 Thus, both GGA+U and HSE06 methods prove that the strain-free Hf2 MnC2 O2 nanosheets are semiconductors, and it retains the semiconductor properties at 4% strains in both zigzag and armchair directions, as illustrated by Fig. S3 (see Supplementary Material). As mentioned in the TABLE S1 and S2, hybrid functional calculations provide higher band gaps. For instance, 0.484 eV and 0.324 eV band gaps are obtained for x = 0% and x = 4%, respectively. On the contrary, GGA+U provides smaller band gap values of 0.282 eV and 0.239 eV at x = 0% and x = 4%, respectively. HSE06 method confirms that bands are shifting in the same way as in GGA+U method. Thus, zero band gaps exist under x = y =-4% and x =8% and also a small band gap (≈ 0.103 eV) exists when y =8% due to bands shifting. Band structure variations are directly related to the charge transfer in each atomic species. Thus, we studied the Bader charge analysis for several strain values in the zigzag and armchair directions, as mentioned in TABLE 3. 50 It is clear from TABLE 3 that the structural variations under compression in either direction considerably affect the charge transfer of Mn and C atoms as compared to those of Hf and O atoms. Small compressive strains alter ∆qMn and ∆qC noticeably as compared to ∆qHf and ∆qO . This is mainly because of considerable shortening of Mn1-C1 and Mn2-C1 bond lengths, while uniaxial compression is applied in the x and y directions respectively (see Fig. 2). However, Hf-O bond length change within 0% to -4% is relatively small. At large tensile strains, less charge transfer for Hf, Mn and C atoms can be seen due to the fact that charges are more localized when the bond strengths

13

ACS Paragon Plus Environment

The Journal of Physical Chemistry

become weaker. (a) 2.5

(b) m h (m0)

1.5 1

*

*

m e(m0)

2

0

1.2

Γ-X Γ-Y

1 0.8 0.6

0.5 -2

0

2

4

0.4

6

-2

0

εx (%)

4

6

4

6

(d) 1 0.8

*

m h (m0)

*

2.5 2 1.5 1 0.5 0

2

εx (%)

(c) 3 m e (m0)

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 27

-2

0

2

4

0.6 0.4

6

-2

0

2

εy (%)

εy (%)

Figure 8: Variation of the effective mass (m∗ ) of (a) electrons and (b) holes under uniaxial strain in the zigzag direction (x ) and m∗ of (c) electrons and (d) holes under uniaxial strain in the armchair direction (y ). Here, m0 is the rest mass of electron.

The effective mass is inversely proportional to the curvature of the band dispersion. We computed the effective masses of electrons and holes as a function of uniaxial strain. The effective mass (m∗ ) can be obtained by using the equation



2

m ~ =h ¯



∂ 2E ∂k 2

−1 ~e~k

(3)

where E and k are the energy and reciprocal lattice vector in a high symmetric direction ~k = k~e~k . 17 The high symmetric directions Γ−X and Γ−Y in reciprocal lattice correspond to the armchair and zigzag directions in real space, respectively.Thus, by calculating the effective mass of electrons and holes in Γ−X and Γ−Y directions, we are able to find a direction with low effective mass, which enhances the transport and thermoelectric properties. 27,51,52 Regardless of the strain direction, Hf2 MnC2 O2 monolayer is a semiconductor in the strain range of [-2%, 6%]. Thus, this strain region was considered for m∗ calculations of this 2D material. For unstrained Hf2 MnC2 O2 , the effective masses of electrons and holes are larger along the Γ−Y direction, since the bottom of CB1 in the Γ−Y direction is flatter than that in the Γ−X direction. When strain is applied in the zigzag direction, m∗e remains 14

ACS Paragon Plus Environment

Page 15 of 27 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

The Journal of Physical Chemistry

around 0.3m0 -0.4m0 in the Γ−X direction, while m∗e is between 2m0 and 2.5m0 in the other direction. Same kind of behavior of m∗e can be seen for y case, except at 6% strain, where m∗e is approximately same in both directions as a result of change in energy order of CB1 and CB2. Therefore, as evidenced by Fig. 8, the armchair direction contains low effective mass electrons.

Uniaxial Strain Effects on Magnetic Properties According to GGA+U calculations, the total magnetic moment (µ) per formula unit is 3.00µB (µB is the Bohr magneton). The main contribution to µ is from the Mn atom, which is in the +4 oxidation state. Our spin-polarized calculations show that contributions from other ions are negligible. The electronic configuration for Mn4+ is [Ar] 3d3 . The uncoupled 3d electrons occupy the t2g states (dxy , dyz , and dxz ), while the eg states (dx2 − dy2 and dz2 ) are empty. 5 The magnetic moment behaves approximately in the same way under uniaxial strain in both directions [see Fig.9(a)]. While tension does not affect µ, compressive strain decreases µ starting at around −2%. This can be attributed to the increased charge transfer from Mn atom to C atom under compression as mentioned in the TABLE 3. Thus, Mn atoms lose electrons, while C atoms gain them under compressive strain. µ becomes 2.986 µB at x = y = −4%. Since the magnetic moment change is negligible, we can assume that magnetic moment is approximately constant within this range of strain. Our HSE06 calculations also show a constant magnetic moment throughout this range of strain, as stated in TABLE S1 and S2 (see Supplementary Materials). In order to find the magnetic ground state of the monolayer at each strain applied, we calculated the total energies for various ferromagnetic (FM), antiferromagnetic (AFM), and non-magnetic (NM) configurations as illustrated in Fig. 10. In the AFM1 configuration [Fig. 10(b)], spin-down Mn atoms form straight lines, whereas in AFM2 [Fig. 10(c)] they form zigzag lines. AFM3 [Fig. 10(d)] has spin-down stripes that are two sites thick. AFM4 [Fig. 10(e)] is the same as AFM1 rotated by 120◦ . As shown in Fig. 11 and S4 (see Supple15

ACS Paragon Plus Environment

The Journal of Physical Chemistry

(a)

3

(b)

εx εy

εx εy

1200

TC (K)

2.995

μ (μΒ)

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 27

2.99

1000

2.985 800

2.98 -4

-2

0

2

4

6

8

10

-4

Strain (%)

-2

0

2

4

6

8

Strain (%)

Figure 9: (a) The total magnetic moment and (b) Curie temperature (TC ) of Hf2 MnC2 O2 as functions of strain.

Figure 10: (a) FM, (b) AFM1, (c) AFM2, (d) AFM3 and (e) AFM4 spin configurations with rectangular super-cells indicated. Only Mn atoms are considered, since contribution for magnetic moment from other atoms are negligible. mentary Materials), we find that EFM has the lowest than energy compared to any of the other configurations. Thus the ground state is ferromagnetic under uniaxial strain within the range −4% to 10%. ENM − EFM is higher than 2.3 eV whereas EAFM − EFM is always very low (less than 0.2 eV). We approximated the magnetic behavior of Hf2 MnC2 O2 by a classical spin-half Ising model on a triangular lattice, with strain-induced anisotropy. The Hamiltonian is 53

H = E0 − J1

X

σi σj − J2

hijix

X

σi σ j ,

(4)

hijidiag

where σi = ±1 are Ising spins, hijix indicates that the sum is over nearest neighbors in the horizontal x direction, hijidiag indicates nearest neighbors in other directions, J1 and J2 are

16

ACS Paragon Plus Environment

Page 17 of 27

0.2

(a)

(b)

EAFM1-EFM EAFM2-EFM EAFM3-EFM EAFM4-EFM

0.15

EAFM - EFM (eV)

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

The Journal of Physical Chemistry

0.1

0.05

0 -4 -2

0

2

4

6

8

10 -4 -2

Strain (%)

0

2

4

6

8

10

Strain (%)

Figure 11: The energy differences per formula unit between each AFM and FM phases (EAFM − EFM ) in the (a) zigzag and (b) armchair directions.

(a)

(b)

Figure 12: Magnetization of Hf2 MnC2 O2 as a function of temperature for strain (ε) directions (a) x and (b) y from anisotropic Ising model simulations of 128×128 triangular lattice. couplings along horizontal and diagonal directions as shown in Fig. 10, and E0 represents the non-magnetic part of the energy. According to Eq. 4, the energies of the spin configurations

17

ACS Paragon Plus Environment

The Journal of Physical Chemistry 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 27

in Fig. 10(a,b,e) are

EFM = E0 − J1 − 2J2

(5)

EAFM1 = E0 − J1 + 2J2

(6)

EAFM4 = E0 + J1 .

(7)

Solving for the couplings yields J1 = 41 (2EAFM4 − EAFM1 − EFM ) and J2 = 14 (EAFM1 − EFM ), which allows us to infer J1 and J2 from the results of DFT energy calculations. Thus, we obtain J1 and J2 as a function of uniaxial strain εx or εy . We used the above Ising model to predict the finite-temperature magnetic behavior of Hf2 MnC2 O2 monolayer. Although Ising models on two-dimensional lattices can be solved using analytic methods, 54 here we used Monte Carlo methods (in particular the Wolff singlecluster algorithm 55 ), which are easily generalized to multilayers or 3D models. Figure 12 shows results of Ising model simulations of 128 × 128 lattices. The magnetization M drops to zero near the Curie temperature Tc , which is strain-dependent. The Binder cumulant crossing-point method 56 was used to determine Tc accurately (see Fig. S4 in Supplementary Materials). Tc increases with increasing strain, as shown in Fig. 9(b). As illustrated by Fig. 9, Tc becomes greater than 1200 K while it converts into a half-metal under the uniaxial tensile strain higher than 8%. Thus, Hf2 MnC2 O2 monolayer can be transformed into an immensely useful material in the spintronic applications.

Conclusion A first principle study based on density functional theory was carried out on Hf2 MnC2 O2 ferromagnetic semiconductors to study the structural, electronic and magnetic property variations. According to the direction of uniaxial strain, different bond lengths are altered differently. It is clear that TM-C bonds are weakened significantly by the tensile strain in the

18

ACS Paragon Plus Environment

Page 19 of 27 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

The Journal of Physical Chemistry

armchair direction compared to the tensile strain in the zigzag direction. Thus, the maximum stress that the system can handle becomes lower when the tensile strain is applied in the armchair direction. The tensile strain (compressive strain) in both directions lowers (enhances) the Young’s modulus. The phonon dispersions indicate that Hf2 MnC2 O2 is dynamically stable under high strains like 10%. As a result of the conduction bands shifting under uniaxial strain, variation of CBM location can be observed while VBM always resides at the Γ point. Consequently, indirect-to-direct band gap transition occurs in 1% to 3% range of tensile strain in the armchair direction. Hf2 MnC2 O2 monolayer becomes a half-metal at x = 7% and y = 9%. Semiconductor-to-metal transition occurs at x = y = −4%. It is worth noticing that Hf2 MnC2 O2 becomes half-metal within the strain range between x = 7% and x = 10% in which Hf2 MnC2 O2 was found to be dynamically stable. Our effective mass calculations show that the armchair direction is preferable for the electron transport due to small effective masses in the Γ−X direction. Moreover, the magnetic moment remains constant (there is only a negligible change) within -4% to 10% strain range. It is found that Curie temperature is highly tunable under uniaxial strain. For instance, TC is around 706 K at -4% in both directions while it is greater than 1200 K at 8% strain.

Supplementary Material • The Band structure calculations, the band gap and the total magnetic moment values using GGA (without U) and HSE06 • Energy differences between nonmagnetic phase and ferromagnetic phase • Binder plot for finding the Curie temperature of strain free monolayer • Cohesive energy calculations and Phonon dispersions

19

ACS Paragon Plus Environment

The Journal of Physical Chemistry 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

Acknowledgement Computer resources used in this work was provided by Computational Research Center (HPC-Linux cluster) at University of North Dakota. A part of this work was supported by University of North Dakota Early Career Award (Grant number: 20622-4000-02624). We also acknowledge financial support from ND EPSCoR through NSF grant OIA-1355466. This work was performed, in part, at the Center for Nanoscale Materials, a U.S. Department of Energy Office of Science User Facility, and supported by the U.S. Department of Energy, Office of Science, under contract no. DE-AC02-06CH11357.

References 1. Prinz, G. A. Magnetoelectronics. Science 1998, 282, 1660–1663. 2. Sarma, S. D.; Fabian, J.; Hu, X.; Zutic, I. Theoretical Perspectives on Spintronics and Spin-polarized Transport. EEE Trans. Magn. 2000, 36, 2821–2826. 3. Kraftmakher, Y. Curie Point of Ferromagnets. Eur. Phys. J 1999, 18, 448. 4. Han, W. Perspectives for Spintronics in 2D Materials. APL Materials 2016, 4, 032401. 5. Dong, L.; Kumar, H.; Anasori, B.; Gogotsi, Y.; Shenoy, V. B. Rational Design of Two-Dimensional Metallic and Semiconducting Spintronic Materials Based on Ordered Double-Transition-Metal MXenes. J. Phys. Chem. Lett. 2017, 8, 422–428. 6. Naguib, M.; Kurtoglu, M.; Presser, V.; Lu, J.; Niu, J.; Heon, M.; Hultman, L.; Gogotsi, Y.; Barsoum, M. W. Two-Dimensional Nanocrystals Produced by Exfoliation of Ti3AlC2. Adv. Mater. 2011, 23, 4248–4253. 7. Anasori, B.; Lukatskaya, M. R.; Gogotsi, Y. 2D Metal Carbides and Nitrides (MXenes) for Energy Storage. Nat Rev Mater. 2017, 2, 16098.

20

ACS Paragon Plus Environment

Page 20 of 27

Page 21 of 27 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

The Journal of Physical Chemistry

8. Zhou, J.; Zha, X.; Zhou, X.; Chen, F.; Gao, G.; Wang, S.; Shen, C.; Chen, T.; Zhi, C.; Eklund, P. et al. Synthesis and Electrochemical Properties of Two-Dimensional Hafnium Carbide. ACS Nano 2017, 11, 3841–3850. 9. Kumar, H.; Frey, N. C.; Dong, L.; Anasori, B.; Gogotsi, Y.; Shenoy, V. B. Tunable Magnetism and Transport Properties in Nitride MXenes. ACS Nano 2017, 11, 7648– 7655. 10. Zhou, Y.; Wang, Z.; Yang, P.; Zu, X.; Yang, L.; Sun, X.; Gao, F. Tensile Strain Switched Ferromagnetism in Layered NbS2 and NbSe2. ACS Nano 2012, 6, 9727–9736. 11. Siriwardane, E. M.; Karki, P.; Sevik, C.; Çakır, D. Electronic and Mechanical Properties of Stiff Rhenium Carbide Monolayers: A First-principles Investigation. Appl. Surf. Sci. 2018, 458, 762 – 768. 12. Akinwande, D.; Brennan, C. J.; Bunch, J. S.; Egberts, P.; Felts, J. R.; Gao, H.; Huang, R.; Kim, J.-S.; Li, T.; Li, Y. et al. A Review on Mechanics and Mechanical Properties of 2D Materials—Graphene and Beyond. EML 2017, 13, 42 – 77. 13. Amorim, B.; Cortijo, A.; de Juan, F.; Grushin, A.; Guinea, F.; Gutiérrez-Rubio, A.; Ochoa, H.; Parente, V.; Roldán, R.; San-Jose, P. et al. Novel Effects of Strains in Graphene and Other Two Dimensional Materials. Phys. Rep. 2016, 617, 1 – 54. 14. Lu, P.; Wu, X.; Guo, W.; Zeng, X. C. Strain-dependent Electronic and Magnetic Properties of MoS2 Monolayer, Bilayer, Nanoribbons and Nanotubes. Phys. Chem. Chem. Phys. 2012, 14, 13035–13040. 15. Kan, M.; Wang, B.; Lee, Y. H.; Sun, Q. A Density Functional Theory Study of the Tunable Structure, Magnetism and Metal-insulator Phase Transition in VS2 Monolayers Induced by In-plane Biaxial Strain. Nano Res. 2015, 8, 1348–1356.

21

ACS Paragon Plus Environment

The Journal of Physical Chemistry 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

16. Khan, S. A.; Amin, B.; Gan, L.-Y.; Ahmad, I. Strain Engineering of Electronic Structures and Photocatalytic Responses of MXenes Functionalized by Oxygen. Phys. Chem. Chem. Phys. 2017, 19, 14738–14744. 17. Deniz Çakir, H. S.; Peeters, F. M. Tuning of the Electronic and Optical Properties of Single-layer Black Phosphorus by Strain. Phys. Rev. B 2014, 90, 205421. 18. Aierken, Y.; Çakir, D.; Peeters, F. M. Strain Enhancement of Acoustic Phonon Limited Mobility in Monolayer TiS3. Phys. Chem. Chem. Phys. 2016, 18, 14434–14441. 19. Amin, B.; Kaloni, T. P.; Schwingenschlögl, U. Strain Engineering of WS2, WSe2, and WTe2. RSC Adv. 2014, 4, 34561–34565. 20. Gao, G.; Ding, G.; Li, J.; Yao, K.; Wu, M.; Qian, M. Monolayer MXenes: Promising Half-metals and Spin Gapless Semiconductors. Nanoscale 2016, 8, 8986–8994. 21. Yuan, J.; Liu, Y.; Wu, P.; Zhou, W. Tuning the Ferromagnetism of a Single Layered Titanium Dioxide Nanosheet with Hole Doping and Uniaxial Strain. J. Phys. Condens. Matter 2018, 30, 305804. 22. Chen, W.; Li, H.-F.; Shi, X.; Pan, H. Tension-Tailored Electronic and Magnetic Switching of 2D Ti2NO2. J. Phys. Chem. C 2017, 121, 25729–25735. 23. Chen, W.; Kawazoe, Y.; Shi, X.; Pan, H. Two-Dimensional Pentagonal CrX (X = S, Se or Te) Monolayers: Antiferromagnetic Semiconductors for Spintronics and Photocatalysts. Phys. Chem. Chem. Phys. 2018, 20, 18348–18354. 24. Zulfiqar, M.; Zhao, Y.; Li, G.; Nazir, S.; Ni, J. Tunable Conductivity and Half Metallic Ferromagnetism in Monolayer Platinum Diselenide: A First-Principles Study. J. Phys. Chem. C 2016, 120, 25030–25036. 25. Zha, X.-H.; Zhou, J.; Luo, K.; Lang, J.; Huang, Q.; Zhou, X.; Francisco, J. S.;

22

ACS Paragon Plus Environment

Page 22 of 27

Page 23 of 27 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

The Journal of Physical Chemistry

He, J.; Du, S. Controllable Magnitude and Anisotropy of the Electrical Conductivity of Hf3C2O2 MXene. J. Phys. Condens. Matter 2017, 29, 165701. 26. Yu, X.-f.; Cheng, J.-b.; Liu, Z.-b.; Li, Q.-z.; Li, W.-z.; Yang, X.; Xiao, B. The Band Gap Modulation of Monolayer Ti2CO2 by Strain. RSC Adv. 2015, 5, 30438–30444. 27. Peng, X.; Wei, Q.; Copple, A. Strain-engineered Direct-Indirect Band Gap Transition and Its Mechanism in Two-Dimensional Phosphorene. Phys. Rev. B 2014, 90, 085402. 28. Li, P.; Li, L.; Zeng, X. C. Tuning the Electronic Properties of Monolayer and Bilayer PtSe2 via Strain Engineering. J. Mater. Chem. C 2016, 4, 3106–3112. 29. Wolf, S. A.; Awschalom, D. D.; Buhrman, R. A.; Daughton, J. M.; von Molnár, S.; Roukes, M. L.; Chtchelkanova, A. Y.; Treger, D. M. Spintronics: A Spin-Based Electronics Vision for the Future. Science 2001, 294, 1488–1495. 30. Siriwardane, E. M. D.; Çakır, D. Strain Engineering of Electronic and Magnetic Properties of Double-transition Metal Ferromagnetic Semiconductor MXenes. J. Appl. Phys. 2019, 125, 082527. 31. Kresse, G.; Hafner, J. Ab initio Molecular Dynamics for Liquid Metals. Phys. Rev. B 1993, 47, 558–561. 32. Kresse, G.; Hafner, J. Ab initio Molecular-Dynamics Simulation of the Liquid-Metal– Amorphous-Semiconductor Transition in Germanium. Phys. Rev. B 1994, 49, 14251– 14269. 33. G. Kresse, J. F. Efficiency of ab-initio Total Energy Calculations for Metals and Semiconductors Using a Plane-wave Basis Set. Comput. Mater. Sci. 1996, 6, 15–50. 34. Kresse, G.; Furthmüller, J. Efficient Iterative Schemes for ab initio Total-energy Calculations Using a Plane-wave Basis Set. Phys. Rev. B 1996, 54, 11169–11186.

23

ACS Paragon Plus Environment

The Journal of Physical Chemistry 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

35. Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865–3868. 36. Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple [Phys. Rev. Lett. 77, 3865 (1996)]. Phys. Rev. Lett. 1997, 78, 1396–1396. 37. Blöchl, P. E. Projector Augmented-wave Method. Phys. Rev. B 1994, 50, 17953–17979. 38. Kresse, G.; Joubert, D. From Ultrasoft Pseudopotentials To the Projector Augmentedwave Method. Phys. Rev. B 1999, 59, 1758–1775. 39. Dudarev, S. L.; Botton, G. A.; Savrasov, S. Y.; Humphreys, C. J.; Sutton, A. P. Electronenergy-loss Spectra and the Structural Stability of Nickel Oxide: An LSDA+U study. Phys. Rev. B 1998, 57, 1505–1509. 40. Krukau, A. V.; Vydrov, O. A.; Izmaylov, A. F.; Scuseria, G. E. Influence of the Exchange Screening Parameter on the Performance of Screened Hybrid Functionals. J. Chem. Phys. 2006, 125, 224106. 41. Baroni, S.; de Gironcoli, S.; Dal Corso, A.; Giannozzi, P. Phonons and Related Crystal Properties from Density-functional Perturbation Theory. Rev. Mod. Phys. 2001, 73, 515–562. 42. Togo, A.; Tanaka, I. First Principles Phonon Calculations in Materials Science. Scr. Mater 2015, 108, 1 – 5. 43. Grimme, S.; Antony, J.; Ehrlich, S.; Krieg, H. A Consistent and Accurate ab initio Parametrization of Density Functional Dispersion Correction (DFT-D) for the 94 Elements H-Pu. J. Chem. Phys. 2010, 132, 154104. 44. Grimme, S.; Ehrlich, S.; Goerigk, L. Effect of the Damping Function in Dispersion Corrected Density Functional Theory. J. Comput. Chem. 2011, 32, 1456–1465.

24

ACS Paragon Plus Environment

Page 24 of 27

Page 25 of 27 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

The Journal of Physical Chemistry

45. van der Marel, D.; Sawatzky, G. A. Electron-electron Interaction and Localization in d and f Transition Metals. Phys. Rev. B 1988, 37, 10674–10684. 46. Bao, J. L.; Gagliardi, L.; Truhlar, D. G. Self-Interaction Error in Density Functional Theory: An Appraisal. J. Phys. Chem. Lett 2018, 9, 2353–2358. 47. Lundberg, M.; Siegbahn, P. E. M. Quantifying the Effects of the Self-interaction Error in DFT: When Do the Delocalized States Appear? J. Chem. Phys. 2005, 122, 224103. 48. Schmidt, T.; Kraisler, E.; Kronik, L.; Kümmel, S. One-Electron Self-Interaction and the Asymptotics of the Kohn–Sham Potential: an Impaired Relation. Phys. Chem. Chem. Phys. 2014, 16, 14357–14367. 49. Arbuznikov, A. V. Hybrid Exchange Correlation Functionals and Potentials: Concept Elaboration. J STRUCT CHEM+ 2007, 48, S1–S31. 50. Yu, M.; Trinkle, D. R. Accurate and Efficient Algorithm for Bader Charge Integration. J. Chem. Phys. 2011, 134, 064111. 51. Gandi, A. N.; Alshareef, H. N.; Schwingenschlögl, U. Thermoelectric Performance of the MXenes M2CO2 (M = Ti, Zr, or Hf). Chem. Mater. 2016, 28, 1647–1652. 52. Pei, Y.; LaLonde, A. D.; Wang, H.; Snyder, G. J. Low Effective Mass Leading to High Thermoelectric Performance. Energy Environ. Sci. 2012, 5, 7963–7969. 53. Farsal, D.; Snina, M.; Badia, M.; Bennai, M. Critical Properties of Two-Dimensional Anisotropic Ising Model on a Square Lattice. J Supercond Nov Magn 2017, 30, 2187– 2195. 54. Husimi, K.; Syozi, I. The Statistics of Honeycomb and Triangular Lattice. Progr. Theor. Exp. Phys 1950, 5, 177–186. 55. Wolff, U. Collective Monte Carlo Updating for Spin Systems. Phys. Rev. Lett. 1989, 62, 361–364. 25

ACS Paragon Plus Environment

The Journal of Physical Chemistry 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

56. Binder, K. Critical Properties from Monte Carlo Coarse Graining and Renormalization. Phys. Rev. Lett. 1981, 47, 693–696.

26

ACS Paragon Plus Environment

Page 26 of 27

Page 27 of 27 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

The Journal of Physical Chemistry

For Table of Contents Only

27

ACS Paragon Plus Environment