Two-Dimensional Manganese Nitride Monolayer ... - ACS Publications

Jun 6, 2018 - 2 > 0, confirming the MnN monolayer is mechanically stable. .... where Jij represents the magnetic exchange coupling constant between th...
10 downloads 0 Views 1MB Size
Subscriber access provided by Kaohsiung Medical University

C: Physical Processes in Nanomaterials and Nanostructures

Two-Dimensional Manganese Nitride Monolayer With Room Temperature Rigid Ferromagnetism Upon Strain Zhenming Xu, and Hong Zhu J. Phys. Chem. C, Just Accepted Manuscript • DOI: 10.1021/acs.jpcc.8b02323 • Publication Date (Web): 06 Jun 2018 Downloaded from http://pubs.acs.org on June 6, 2018

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

1

Two-Dimensional Manganese Nitride Monolayer with Room

2

Temperature Rigid Ferromagnetism upon Strain

3 Zhenming Xu, Hong Zhu*

4 5 6

University of Michigan–Shanghai Jiao Tong University Joint Institute,

7

Shanghai Jiao Tong University, 800, Dongchuan Road, Shanghai, 200240, China

8 9

Abstract: Developing low-dimensional spintronic materials with room temperature

10

magnetic ordering and large spin polarization is the key for the fabrication of practical

11

spintronic devices with high circuit integration density and speed. Here,

12

first-principles calculations were performed to systematically investigate a

13

two-dimensional hexagonal MnN monolayer with room temperature magnetic

14

ordering and 100% spin polarization. MnN monolayer is thermally, dynamically and

15

mechanically stable. MnN monolayer is intrinsically half-metallic with very wide

16

band gap. The Curie temperature of MnN monolayer is estimated to be ~368 K,

17

which is more than the room temperature and insensitive to strain. MnN monolayer

18

shows half-metallic and excellent magnetic stability upon external strains from -10%

19

to 10%. Our calculations indicate that MnN monolayer could be a promising material

20

for the room temperature spintronic devices.

21 22

Introduction

23 24

Graphene, discovered by Andre Geim and co-workers by a simple method of scotch

25

tape stripping in 2004, opens a new era of studying and applying two-dimensional

26

(2D) materials1-4. The reduction of material dimension gives rise to the novel physical

27

and chemical properties, which are quite different from that of the three-dimensional

28

(3D) structure5-6. 2D materials can be widely applied in many fields, such as

*

Corresponding author:

E-mail: [email protected]; University of Michigan–Shanghai Jiao Tong University Joint Institute, Shanghai Jiao Tong University, 800, Dongchuan Road, Shanghai, 200240, China

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

1

nanoelectronic and optical devices7-10, energy storage11-16, catalysis17-20 and chemical

2

sensors21-23.

3 4

Nanoelectronic devices, especially for spintronics, which are urgently needed for

5

developing the high speed and low energy consuming information technology, are

6

one of the most important application areas for 2D materials24. It is due to the fact that

7

fabricating spintronic devices at nanoscale based on the low-dimensional spintronic

8

materials can achieve high integration density and speed24. However, most 2D

9

materials like graphene1, 3, graphdiyne25-26, phosphorene27-28, h-BN29-30, borophene31,

10

etc. intrinsically show non-magnetic or weakly magnetic characteristics, which limit

11

their applications in spintronics. In general, a practical spintronic material should have

12

a sizable room temperature magnetic moment, which means that a large spin

13

polarization of electron state around the Fermi energy level is essential24. Among

14

many reported 2D materials, only a few 2D transition metal carbides, nitrides, sulfide

15

and trihalides show intrinsic ferromagnetisms, such as Ti2C32, Cr2C33, Mn2C34, FeC235,

16

Ti2N32, CrN36, MoN237, MoS238-39, NbS240, and TMX2 transition metal dihalides (X =

17

Cl, Br and I) monolayer41-42, in which the unpaired d-electrons in transition metals

18

contribute to the formation of the spin-polarized electronic states. In addition, earlier

19

reports suggest that transition metal (TM) incorporation into metal-free 2D materials

20

such as graphene43-44 and carbon nitride45-46 could induce local magnetic moment and

21

also lead to half-metallicity.

22 23

Here, we report a room temperature ferromagnetic and half-metallic 2D hexagonal

24

manganese nitride with 100% spin polarization and high Curie temperature, which

25

may be synthesized by a spontaneous graphitic conversion from the ultrathin

26

(111)-oriented cubic MnN in the future experiments47. Similar spontaneous

27

graphitization of ultrathin cubic structures can be found in diamond, boron nitride,

28

and rock salt47. In this work, using first-principles calculations and ab-initio molecular

29

dynamics (AIMD) simulations based on density functional theory (DFT), we have

30

systematically investigated the hexagonal MnN monolayer to ascertain its potential

31

application for spintronic devices. Firstly, the thermal, dynamical, and mechanical

32

stability of MnN monolayer were confirmed by AIMD simulation, phonon dispersion

33

calculation, and mechanical parameter calculation. Then, the spin-polarized electronic 2

ACS Paragon Plus Environment

Page 2 of 24

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

1

structures, magnetic properties and Curie temperature of MnN monolayer with and

2

without strain are discussed in detail.

3 4

Computational methods

5 6

All atomic and electronic structure calculations were carried out by using the

7

projector augmented wave (PAW) method48 in the framework of DFT49, as

8

implemented in Vienna ab-initio Simulation Package (VASP). The generalized

9

gradient

10

approximation

(GGA)50 49

exchange-correlation functional 6

2

5

2

and

Perdew–Burke–Ernzerhof

(PBE)

is used. The electron configurations of Mn and N

3

11

atom are 3p 4s 3d and 2s 3p , respectively. To consider the strong correlation effects

12

of transition metal Mn element, structural optimizations were performed by using the

13

spin-dependent GGA plus Hubbard correction U (GGA+U) method51. The Hubbard U

14

parameter of Mn atom is set to 5.5 eV based on the previous study on manganese (III)

15

compound52. For more accurate descriptions of electronic structures, the hybrid

16

functional HSE0653-54 is employed for the band-structure and density of electron state

17

(DOS) calculations. A vacuum region of 15 Å is applied to avoid interactions between

18

the neighboring configurations. After convergence tests (see details in Supporting

19

Information, Table S1), the plane-wave energy cutoff is set to 520 eV. The

20

Monkhorst–Pack method55 with 11×11×1 and denser 25×25×1 k-point meshes are

21

mployed for the Brillouin zone sampling, respectively for the structural relaxations

22

and electronic structure calculations. The convergence criterions of energy and force

23

calculations are set to 10−5 eV and 0.01 eV Å−1, respectively. The phonon calculations

24

were performed with a 3×3×1 supercell based on the density functional perturbation

25

theory (DFPT)56, as implemented in the PHONOPY code57. In addition, ab-initio

26

molecular dynamics (AIMD) simulations based on the Born-Oppenheimer

27

approximation were performed at the room temperature of 300 K in a statistical

28

ensemble with fixed particle number, volume and temperature (NVT) by using a

29

3×3×1 supercell of MnN monolayer as initial structure. The time step is set to 3 fs,

30

and the supercell systems were simulated for 10000 steps, with a total of 30 ps. A

31

plane-wave energy cutoff of 350 eV, Gamma-point only k-point mesh and periodic

32

boundary condition were also used dor AIMD simulations.

33 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

1 2

Results and discussion

3

Geometry and stability

4 5

Our proposed graphene-like MnN monolayer belongs to the hexagonal honeycomb

6

structure, and the corresponding atomic structure is depicted in Figure 1. The

7

optimized lattice constant a of MnN monolayer is 3.388 Å, and the corresponding

8

Mn-N bond length is 1.956 Å. To confirm the dynamical stability of MnN monolayer,

9

its phonon dispersion was firstly calculated, as depicted in Figure 2a. All the

10

vibrational modes of our predicted MnN monolayer show positive frequencies,

11

indicating MnN monolayer is stable or at least structurally metastable. The vibrational

12

modes with lower frequencies (less than 450 cm−1) are composed of three acoustic

13

branches and one optical branch. Viewed from the total and partial density of phonon

14

states, Mn atom with greater atomic mass mainly devotes to the acoustic vibrational

15

modes with lower frequencies. While the high-frequency optical branches especially

16

for those above 800 cm-1 are mainly contributed by N atom with smaller atomic mass,

17

separated from the lower frequency optical vibrational modes by a phonon gap of 350

18

cm−1.

19 20

Then, AIMD simulations were performed by using a 3×3×1 supercell of MnN

21

monolayer at 300 K for 30 ps to examine its thermal stability. During the whole

22

AIMD simulation process, the total potential energies of MnN monolayer vibrate

23

around a constant of ~ -8.38 eV/atom, and the fluctuations of total potential energy

24

during the last 10 ps are ~20 meV/atom, as shown in Figure 2b. In addition, the

25

supercell structure after 30 ps AIMD simulations nearly maintains the initial

26

hexagonal honeycomb structure with some weak puckers, as illustrated in the inset of

27

Figure 2b, which indicates the thermal stability of MnN monolayer. The elastic

28

constants of MnN monolayer were also calculated to further verify its mechanical

29

stability58-59. According to the Born elastic theory60-61, the elastic constant constitutes a

30

symmetric 6 × 6 tensor matrix in the linear elastic range. Our calculated elastic

31

constants are C11 = 36.663 GPa, C22 = 36.663 GPa, C12 = 32.193 GPa, and C44 =

32

2.235 GPa, respectively (see details in Supporting Information, Table S2). All these

33

elastic constants comply well with the Born criteria for the 2D hexagonal crystal60-61, 4

ACS Paragon Plus Environment

Page 4 of 24

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

1

that is C11 > 0, C22 > 0, C44 > 0 and C11C22 −C122 > 0, confirming MnN monolayer is

2

mechanically stable. All above analyses of the dynamical, thermal, and mechanical

3

properties of MnN monolayer consistently demonstrate its stability under ambient

4

environment.

5 6

Figure 1. Schematic diagrams of the top and side views of 2D MnN 2×2×1 supercell structure, the purple and blue

7

balls represent Mn and N atom, respectively,

8 9

Figure 2. (a) Phonon spectrum, total and partial density of phonon states; (b) total potential energy fluctuations of

10

MnN monolayer during 30 ps AIMD simulation at 300 K, and the final puckered structures after AIMD simulation

11

are also shown.

12 5

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

1

Page 6 of 24

Electronic structure and magnetic properties at equilibrium lattice

2 3

Electronic structure and magnetic properties of the material at equilibrium lattice are

4

important to investigate its potential application for spintronic devices. We firstly

5

found out that the total energy of spin-polarized MnN monolayer is lower than that of

6

the spin-unpolarized MnN monolayer by more than 0.40 eV/atom. Then, we

7

calculated the magnetocrystalline anisotropy energies (MAEs) of MnN unit cell by

8

considering the effect of spin-orbital coupling. Herein, two magnetization directions

9

in plane, namely [100] and [110] direction, and two directions out of plane, namely

10

[001] and [111] directions, are considered. The calculated MAEs by the hybrid

11

functional HSE06 for [100], [110], and [111] directions are 135.59, 134.27, and

12

160.00 μeV/Mn atom relative to the [001] direction, respectively, as shown in Table

13

S3, which indicates the magnetic moment along the [001] direction (out-of-plane) is

14

much more stable than other directions. Thus, the easy axis (EA) of MnN monolayer

15

is along the [001] direction (out-of-plane), which is consistent with 2D Fe2Si62 and

16

FeC235 sheet.

17 18

The 4×4 supercells of MnN monolayer with one ferromagnetic (FM) and six

19

antiferromagnetic (AFM) spin orderings along the [001] direction were constructed

20

based on the previous work of 2D hexagonal ferromagnetic materials42,

21

schematically shown in Figure S1. The exchange energy (Eex) between FM and AFM

22

spin ordering structures is thus calculated by the formulas of Eex = EAFM – EFM based

23

on the GGA+U method, which has been already successfully applied for 2D metal

24

halides42. If the exchange energy is positive, the ground state of 2D MnN monolayer

25

system is ferromagnetic. To evaluate the sensitiveness of the magnetic state to the

26

Hubbard U parameter, Eex of the MnN 4×4 supercell with different U values have

27

been calculated, as shown in Figure S2. It confirms that FM state is more preferable

28

than the AFM state for all considered U values and Eex converges to 0.410 eV/Mn

29

atom around U = 6.0 eV, which validates our choice for Hubbard U parameter for

30

Manganese of 5.5 eV based on the previous study on manganese (III) compound52.

31

According to these calculated positive exchange energies of other antiferromagnetic

32

(AFM) spin orderings (Figure S1), MnN monolayer prefers the FM ground state with

33

magnetic moment along the [001] direction. Therefore, we only focus on the FM 6

ACS Paragon Plus Environment

63

, as

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

1

ground state with magnetic moment along the [001] direction to study the electronic

2

structure and magnetic properties in the following.

3 4

The spin-polarized band structures and total density of electron states (TDOS) were

5

calculated at HSE06 level for better understanding of the electronic structures and

6

magnetic properties of MnN monolayer. Viewed from Figure 3, MnN monolayer

7

exhibits completely different electron conductivities of the spin-up and spin-down

8

channels. The spin-down channel shows a very large direct band gap of 5.249 eV with

9

an insulator feature, while the spin-up channel shows metallic characteristic due to the

10

considerable occupied electron states cross the Fermi energy level. Therefore, MnN

11

monolayer is intrinsically half-metallic, and hence has 100% spin polarization near

12

the Fermi level, which is quite different from its rocksalt-structured phase without any

13

spin polarization (Figure S3). Moerover, the partial density of electron states (PDOS)

14

and partial charge distributions were also calculated to understand the bonding

15

mechanism of Mn-N (see details in Supporting Information, Figure S4).

16 17

We note from Figure S4 that Mn 3 orbital is non-degenerate and does not bond to

18

other orbitals. The double-degenerate Mn 3dxz and 3dyz spin-up and -down states

19

overlap with N pz states and form π-type bonding orbitals parallel to the plane of MnN

20

monolayer. The Mn 3dxy and Mn 3   states are doubly degenerate, hybridizing

21

with a Mn 4s state and forming three Mn sd2 hybrid orbitals. These sd2 hybrid orbitals

22

overlap with three N sp2 hybrid states and form a MnN hexagonal sheet. This Mn

23

3d-orbital split is a typical picture for the trigonal-type complexes, which is quite

24

different from the common 3d orbital splitting of the octahedral symmetry crystal

25

with lower energy orbitals t2g (3dxy, 3dxz and 3dyz) and higher energy orbitals eg

26

(3   and 3  )64-65. This is because that the symmetry and potential field of Mn

27

atom in this 2D MnN monolayer is different from that of 3D octahedral crystal, which

28

is consistent with the previous relevant report of the CrN monolayer36.

29 30

Mn 4s, 3dxy and 3   orbitals strongly overlap with N 2s, 2px and 2py orbitals,

31

respectively, and two electrons from 3dxy and 3   orbital of Mn atom exactly

32

pair with two electrons from 2px and 2py orbital of N atom, respectively, leaving three

33

unpaired 3d electrons of Mn atom and one unpaired 2p electron of N atom. Due to 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

1

100% spin polarization of these four unpaired valence electrons near the Fermi level,

2

MnN monolayer has a net magnetic moment of 4  per unit cell. Here, Mn atom

3

contributes to three quarters of the magnetism with major spin-densities, as the spin

4

charge density plots shown in Figure 3c. MnN monolayer with a wide half metallic

5

gap and large spin polarization can effectively serve in room temperature spintronic

6

devices.

7 8

Figure 3. Electronic structure of MnN monolayer at HSE06 level. (a) Spin-polarized energy band structure of MnN

9

monolayer; (b) Spin-polarized total density of electron state (TDOS) of MnN monolayer; Spin-up and spin-down

10

is depicted in black and red, respectively. The dashed lines in (a) and (b) indicate the Fermi level. (c) Top and side

11

views of the spatial spin density distribution isosurfaces (0.01 e Å−3) of MnN monolayer. The yellow and

12

light-blue isosurfaces correspond to spin-up and spin-down states, respectively.

13

8

ACS Paragon Plus Environment

Page 8 of 24

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

1

Curie temperature

2 3

Curie temperature of the ferromagnetism–paramagnetism transition for ferromagnetic

4

materials near or above 300 K is another critical determining factor for the practical

5

room temperature spintronic devices. Because of the magnetic moment of MnN

6

monolayer perferring the out-of-plane [001] direction (see details in Supporting

7

Information, Table S3), and the calculated MAEs by the HSE06 method being

8

comparable to the reported 2D Ising-like ferromagnetic CrXTe3 and Fe2Si sheets62, 66,

9

the Kosterlitz-Thouless transition is unlikely to occur in MnN monolayer67. Thus,

10

Curie temperature of MnN monolayer can be estimated by using statistical Monte

11

Carlo (MC) simulations based on the Ising model. This calculation method has been

12

widely used to estimate the Curie temperature of 2D materials35, 68. In the 2D Ising

13

model, the nearest neighboring spin interactions can be approximatively equivalent to

14

an exchange constant J, and the spin system energy is described by the classical

15

Hamiltonian of Ising model: H = −     ,

16

where J represents the magnetic exchange coupling constant between two

17

neighboring i and j spins, and S and S are spin parallels or anti-parallels to the

18

vacuum direction of 2D model, respectively for i and j site. The exchange coupling

19

constant J can be evaluated from the exchange energy between FM and AFM spin

20

arrangements in Figure S1 in Supporting Information. Here, we adapted the highest

21

energy AFM spin arrangement for exchange coupling constant J calculation based on

22

the previous work about 2D hexagonal crystal structure42, 63, 69: J =

2  

23

where Eex is 0.408 eV per Mn atom (unit cell), z is the number of the nearest

24

neighbors around one Mn atom (z = 6 for 2D hexagonal crystal structure), and S is the

25

net magnetic moments of 4  per unit cell of MnN monolayer. Therefore, the

26

corresponding estimated exchange coupling constants J is determined here to be 8.499

27

meV.

28

9

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 10 of 24

1

Then, MC simulations for 105 steps by using the Metropolis algorithm and a 2D

2

hexagonal 8080 grid points with periodic boundary conditions were performed to

3

study the ferromagnetic systems. The calculated magnetizations per unit cell (Mn

4

atom) of the unstrained MnN monolayer as a function of different temperatures are

5

depicted in Figure 5a. Below the ice point of 273.15 K, all spins in MnN monolayer

6

are aligned in parallel, providing a FM state. However, when temperature exceeding

7

400 K, spins tend to be randomly aligned due to the stronger temperature perturbation,

8

resulting in zero net magnetization (PM state). Our calculated Curie temperature of

9

the FM-PM phase transition for MnN monolayer is ~368 K, which is higher than the

10

room temperature. This indicates MnN monolayer meets the basic condition of

11

applying ferromagnetic materials for the practical room temperature spintronic

12

devices. In addition, our calculated Curie temperature of MnN monolayer is

13

comparable with some other 2D ferromagnetic materials, such as metal halides MX2

14

(X = Cl, Br and I) (23–138 K)42, CrSnTe3 (170 K)66, and FeC2 (245 K)35.

15

16 17

Figure 5. (a) Variations of the calculated magnetizations (in  per unit cell) of MnN monolayer with respect to

18

different temperatures and equi-biaxial strains. (b) Total magnetic moment (in  per unit cell) at 0 K and Curie

19

temperature of MnN monolayer with respect to different equi-biaxial strains

20 21 22

Strain effect

23 24

Strain effect on the electronic structure and magnetism of 2D materials is another

25

important concern for their applications in spintronic devices37,

26

studied the effect of strain on the electronic structure and magnetic properties of MnN 10

ACS Paragon Plus Environment

40, 70-72

. Here, we

Page 11 of 24 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

1

monolayer, and applied a series of equi-biaxial strains from −10% to 30% to MnN

2

monolayer whose atomic positions were fully relaxed. The strain is defined as ε = (a –

3

a0)/a0 × 100% = (b – b0)/b0× 100%, where a, b and a0, b0 are the lattice constants of

4

the unit cell under strain and at equilibrium, respectively. Here, the positive and

5

negative ε represents the tensile and compressive strain, respectively. The lattice

6

constant a, Mn-Mn and Mn-N bond length of the strained MnN monolayer systems

7

are summarized in Table S4 in Supporting Information. The energy-strain and stress–

8

strain curves of MnN monolayer were calculated, as shown in Figure S5 in

9

Supporting Information. Applying strain on the equilibrium monolayer can increase

10

its total energy and stress, as expected. When applying equi-biaxial strains from 20%

11

to 30%, the stress of MnN monolayer doesn’t further change. We note that an

12

equi-biaxial strain of 30% doesn’t exceed the elastic limit of MnN monolayer,

13

because no rapid attenuation of its stress response upon strain is observed in Figure

14

S5. Considering it is difficult to reach those strains more than 10% in the

15

experimental operations and practical applications, so we only focus our discussions

16

on the electronic and magnetic property changes upon equi-biaxial strains ranging

17

from −10% to 10%. The total magnetic moments ( per unit cell) at HSE06 level of

18

MnN monolayer with respect to the realistic equi-biaxial strains ranging from −10%

19

to 10% were calculated, as shown in Figure 5b and summarized in Table S3 in

20

Supporting Information. It can be seen that both these compressive and tensile strains

21

don't change the total magnetic moment of MnN monolayer at all.

22 23

For better understanding the rigidity of the total magnetic moment under different

24

equi-biaxial strains, the spin-polarized total density of electron states (TDOS) at

25

HSE06 level of MnN monolayer under different equi-biaxial strains ranging from −10%

26

to 10% were calculated, as depicted in Figure 6. The robust band gaps of the

27

spin-down channels and the metal-like spin-up channels of those strained MnN

28

monolayer are preserved, displaying 100% spin-polarization near the Fermi level,

29

which are similar to the unstrained case. Band gaps and band edges of the spin-down

30

electron state at HSE06 level for MnN monolayer under different strains are depicted

31

in Figure 7. A closer look shows the band gaps of the spin-down electron states

32

decrease upon compressive or tensile strains (Figure 7a), especially for the

33

− 8 % and 10 % strained systems, which are similar to those cases of the band gaps 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

1

of 2D MoS273-74 and phosphorene (in the zigzag direction)75 semiconductors

2

decreasing upon the applications of tensile or compressive strains. It can be seen from

3

Figure 7b that the valence band maximum (VBM) of the spin-down electron states

4

with respect to the Fermi level goes up from -3.340 eV (unstrained system) to -2.891

5

eV (-6% compressive strain), then goes down to -3.937 eV (-10% compressive strain)

6

when the compressive strains increasing from 0% to -10%. Meanwhile, the

7

conduction band minimum (CBM) significantly shifts down from 1.909 to 0.625 eV

8

when the compressive strains increasing from 0% to -10%, whereas the shifts of band

9

edge don't display a linear function upon strain. The nonlinear shifts of

10

energy-band-edge upon strain can be explained by the Heitler-London's exchange

11

energy model75, that the energy-band-edge shifts are closely related to the

12

bonding/antibonding nature of the near-band-edge orbitals upon strain76. Even MnN

13

monolayer under moderate strains (-2% and 2%), its VBM and CBM still shift

14

significantly. During the stretching period with increased tensile strains from 0% to

15

10%, VBM markedly goes up from -3.340 to -2.415 eV, and CBM also shifts up from

16

1.909 to 2.194 eV (8% tensile strain), then goes down to 2.177 eV (10% tensile strain).

17

The increasing magnitude of VBM is larger than that of CBM, thus the corresponding

18

band gaps of spin-down states of those MnN monolayers under compressive strains

19

significantly reduce. Moreover, we also calculated the DOS of these strained MnN

20

monolayer by GGA+U method for validating those calculated results by HSE06

21

method, as shown in Figure 7 and Figure S6 in Supporting Information. The

22

calculated CBM edges of spin-down states by HSE06 method are in good accordance

23

with those by GGA+U method, while the VBM edges and band gaps calculated by

24

HSE06 method are larger than those by GGA+U method, which is common for DFT

25

calculation. Although there are absolute deviations between the band gaps calculated

26

by these two methods, the variation trend of band gaps calculated by HSE06 method

27

is basically consistent with that by GGA+U method. In this section, we only focus on

28

the variation of band gap upon different strains rather than the absolute value, so both

29

the HSE06 and GGA+U method are competent for producing the variation of band

30

gap.

31 32

In general, strains from -10% to 10% don't significantly change the electronic

33

structure of MnN monolayer (Figure 6), and 100 % spin-polarizations near the Fermi 12

ACS Paragon Plus Environment

Page 12 of 24

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

1

level of these strained system are fully preserved by the robust band gaps of

2

spin-down channels, so their magnetic properties are completely same as the

3

unstrained system (Figure 5b). MnN monolayer shows excellent FM half metallic and

4

magnetic stability upon external strains from -10% to 10%, so the lattice mismatches

5

between MnN monolayer and substrate materials in the experimental synthetic

6

process don't affect its magnetic property at all. This is different from most reported

7

two-dimensional materials with variational magnetic moments upon strain32,

8

such as half-metallic 2D Fe2Si crystal with variable magnetic moments of 3.02-3.11

9

μ! under moderate strains (-3% to 3%)62. Our predicted MnN monolayer is promising

10

to be used for the spin electronic devices where the stable spin magnetic moment

11

against strain under room environment is needed.

37, 77

,

12 13

In addition, the dependence of Curie temperature of MnN monolayer with respect to

14

strain has been also studied. The energy difference between AFM and FM spin

15

ordering structures under compressive strain sharply goes down from 0.325 to 0.057

16

eV per unit cell, and the tensile strained systems slightly go down to 0.275 eV per unit

17

cell, as summarized in Table S4 in Supporting Information. It can be seen from Figure

18

5b that the response curve of Curie temperature v.s. strain is not monotonous, and it

19

reduces more rapidly under enhanced compressive strain than those under enhanced

20

tensile strain. In other words, Curie temperatures of those strained MnN monolayers

21

decrease with the biaxial tensile or compressive strain, indicating that there is no

22

room to ramp up Curie temperature of MnN monolayer by applying strain. Wu et. al.,

23

hold that the tensile strain reduces Curie temperature of 2D Fe2Si crystal, and the

24

compressive strain enhances Curie temperature, due to the monotonous variation of

25

Fe-Fe bond length resulting in the monotonous variation of the interaction between

26

neighbor Fe atoms62, which is different from our MnN crystal with non-monotonous

27

response. For 2D FeCl2, NbS2 and NbSe2 crystals, their Curie temperatures can be

28

dramatically enhanced by the tensile strain based on DFT calculations42, 69, in contrast

29

to 2D Fe2Si crystal and our MnN monolayer. Experimentally, the enhancement of

30

Curie temperature by tensile stress has been carried out in Fe64Ni36 Invar alloy

31

because of the increase of Fe-Fe interatomic distance78. While for the single layer

32

CrSnTe3, it is found that both compressive and tensile strains can reduce its Curie

33

temperature66, which is consistent with our MnN monolayer. Therefore, the energy 13

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

1

differences between the strained AFM and FM spin ordering structures (effective

2

couple interactions between neighboring magnetic atoms) are not the exclusive

3

responses to structural deformation (bond length of neighbor magnetic atoms), and

4

the variation of electronic structure upon strain may be another important factor.

5

Interestingly, except for the cases of the strained MnN monolayer under larger

6

compressive strain, the variation tendency of Curie temperature for the other strained

7

MnN monolayers is generally in accordance with the variation curve of band gap and

8

VBM edge of spin-down electron state (Figure 7). Since the protective effect of 100%

9

spin polarization near the Fermi level gradually diminished by the reduced band gap

10

under enhanced compressive or tensile strain, the energy differences between FM and

11

AFM structure gets smaller, and FM spin ordering structures are more likely to switch

12

to AFM or PM spin disordered structures. It is worth to point out that some other

13

important factors may also affect Curie temperature of our strained MnN monolayer.

14

More future efforts should be made to thoroughly figure out the general influence

15

mechanism of structural deformation, electronic structure and some other factors on

16

Curie temperature of 2D materials.

14

ACS Paragon Plus Environment

Page 14 of 24

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

1 2

Figure 6. Spin-polarized total density of electron state (TDOS) near the Fermi level at HSE06 level of MnN

3

monolayer with respect to different equi-biaxial strains. The dashed line means the Fermi level is set to zero eV.

15

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

1 2

Figure 7. (a) Band gap (in eV) and (b) band edge (in eV) of spin-down electron state of MnN monolayer under

3

different equal-biaxial strains. All energies are referenced to vacuum level.

4 5

Summary

6 7

A graphene-like two-dimensional hexagonal MnN monolayer with unique electronic

8

structure and stable magnetic properties was studied by first-principles calculations.

9

According to AIMD simulations, phonon spectrum and elastic constants analyses,

10

MnN monolayer is thermally, dynamically and mechanically stable. Electronic

11

structure calculations indicate that MnN monolayer is intrinsically half-metallic with

12

a net magnetic moment of 4  per unit cell, due to the 100% spin polarization and

13

large band gap near the Fermi level. The Curie temperature of ~368 K of MnN

14

monolayer is estimated by MC simulations, which is higher than the room

15

temperature and

16

two-dimensional materials with variational magnetic moments upon strain, MnN

17

monolayer can maintain the FM half metallicity and a constant magnetic moment

not enhanced

by strain.

Different from most reported

16

ACS Paragon Plus Environment

Page 16 of 24

Page 17 of 24 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

1

even under ± 10 % strains, due to the 100% spin polarization electron states near the

2

Fermi level being fully preserved by the robust band gap of spin-down electron states.

3

Our calculations indicate that the MnN monolayer could be a promising material for

4

the room temperature spintronic devices.

5 6 7

Author Information

8 9

Corresponding Authors

10

*

11

The authors declared that they have no conflicts of interest to this work.

E-mail: [email protected]

12 13 14

Supporting Information

15 16

Energy cut-off and k-point convergence. Elastic constants. Magnetic anisotropy

17

energies. Lattice constant, magnetic moment, average energy difference between

18

AFM and FM states, exchange coupling constant, energy band edge and band gap of

19

spin-down electron state respect to different equi-biaxial strains. AFM spin

20

arrangements. Energy difference between AFM and FM spin arrangements with

21

different

22

rocksalt-structured phase calculated by GGA. Energies and stresses response under

23

different equi-biaxial strains. Spin-polarized TDOS at GGA+U level respect to

24

different equi-biaxial strains.

Hubbard

correction

U

values.

Electronic

structures

of

MnN

25 26 27

Acknowledgements

28 29

This work is supported by the National Natural Science Foundation of China

30

(51602196) and the Shanghai Sailing Program (16YF1406000). All simulations were

31

performed at the Shanghai Jiao Tong University High Performance Computing Center.

32

The authors also thank Prof. Weidong Luo for the discussions.

33 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

1

References

2 3

1.

4

20 (12), 2277-2289.

5

2.

6

graphene oxide: synthesis, properties, and applications. Adv. Mater. 2010, 22 (35), 3906-3924.

7

3.

8

two-dimensional nanomaterial. Angew. Chem. Int. Ed. Engl. 2009, 48 (42), 7752-7777.

9

4.

Loh, K. P.; Bao, Q.; Ang, P. K.; Yang, J., The chemistry of graphene. J. Mater. Chem. 2010,

Zhu, Y.; Murali, S.; Cai, W.; Li, X.; Suk, J. W.; Potts, J. R.; Ruoff, R. S., Graphene and

Rao, C. N.; Sood, A. K.; Subrahmanyam, K. S.; Govindaraj, A., Graphene: the new

Novoselov, K. S.; Jiang, D.; Schedin, F.; Booth, T. J.; Khotkevich, V. V.; Morozov, S. V.;

10

Geim, A. K., Two-dimensional atomic crystals. Proc. Natl. Acad. Sci. U. S. A. 2005, 102 (30),

11

10451-10453.

12

5.

13

2013, 113 (5), 3766-3798.

14

6.

15

9451-9469.

16

7.

17

phosphorene. Journal of Materials Chemistry C 2017, 5 (9), 2291-2300.

18

8.

19

der Waals heterostructures. Science 2016, 353 (6298), 9439-9439.

20

9.

21

heterostructures and devices. Nature Reviews Materials 2016, 1, 16042.

22

10. Xu, J.; Shim, J.; Park, J. H.; Lee, S., MXene Electrode for the Integration of WSe2 and

23

MoS2 Field Effect Transistors. Adv. Funct. Mater. 2016, 26 (29), 5328-5334.

24

11. Anasori, B.; Lukatskaya, M. R.; Gogotsi, Y., 2D metal carbides and nitrides (MXenes) for

25

energy storage. Nature Reviews Materials 2017, 2, 16098.

26

12. Naguib, M.; Halim, J.; Lu, J.; Cook, K. M.; Hultman, L.; Gogotsi, Y.; Barsoum, M. W., New

27

two-dimensional niobium and vanadium carbides as promising materials for Li-ion batteries. J.

28

Am. Chem. Soc. 2013, 135 (43), 15966-15969.

29

13. Liu, F.; Zhou, J.; Wang, S.; Wang, B.; Shen, C.; Wang, L.; Hu, Q.; Huang, Q.; Zhou, A.,

30

Preparation of High-Purity V2C MXene and Electrochemical Properties as Li-Ion Batteries. J.

31

Electrochem. Soc. 2017, 164 (4), A709-A713.

32

14. Mashtalir, O.; Lukatskaya, M. R.; Zhao, M. Q.; Barsoum, M. W.; Gogotsi, Y., Amine‐

33

Assisted Delamination of Nb2C MXene for Li‐Ion Energy Storage Devices. Adv. Mater. 2015,

Xu, M.; Liang, T.; Shi, M.; Chen, H., Graphene-like two-dimensional materials. Chem. Rev.

Zhang, H., Ultrathin Two-Dimensional Nanomaterials. ACS Nano 2015, 9 (10),

Lin, H.-F.; Liu, L.-M.; Zhao, J., 2D lateral heterostructures of monolayer and bilayer

Novoselov, K. S.; Mishchenko, A.; Carvalho, A.; Castro Neto, A. H., 2D materials and van

Liu, Y.; Weiss, N. O.; Duan, X.; Cheng, H.-C.; Huang, Y.; Duan, X., Van der Waals

18

ACS Paragon Plus Environment

Page 18 of 24

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

1

27 (23), 3501-3506.

2

15. Liu, Y.; Zhao, Y.; Jiao, L.; Chen, J., A graphene-like MoS2/graphene nanocomposite as a

3

highperformance anode for lithium ion batteries. Journal of Materials Chemistry A 2014, 2 (32),

4

13109-13115.

5

16. Xu, Z.; Lv, X.; Li, J.; Chen, J.; Liu, Q., A promising anode material for sodium-ion battery

6

with high capacity and high diffusion ability: Graphyne and Graphdiyne. RSC Adv. 2016, 6,

7

25594–25600.

8

17. Rahman, M. Z.; Kwong, C. W.; Davey, K.; Qiao, S. Z., 2D phosphorene as a water

9

splitting photocatalyst: fundamentals to applications. Energy & Environmental Science 2016, 9

10

(3), 709-728.

11

18. Deng, D.; Novoselov, K. S.; Fu, Q.; Zheng, N.; Tian, Z.; Bao, X., Catalysis with

12

two-dimensional materials and their heterostructures. Nature nanotechnology 2016, 11 (3),

13

218-230.

14

19. Lukowski, M. A.; Daniel, A. S.; Meng, F.; Forticaux, A.; Li, L.; Jin, S., Enhanced hydrogen

15

evolution catalysis from chemically exfoliated metallic MoS2 nanosheets. J. Am. Chem. Soc.

16

2013, 135 (28), 10274-10277.

17

20. Xue-Rong Shi, S.-G. W., Jia Hu, Hui Wang, Yan-Yan Chen,Zhangfeng Qin, Jianguo Wang,

18

Density functional theory study on water–gas-shift reaction over molybdenum disulfide.

19

Applied Catalysis A: General 2009, 365, 62–70.

20

21. Gonzalez, C.; Biel, B.; Dappe, Y. J., Adsorption of small inorganic molecules on a

21

defective MoS2 monolayer. PCCP. 2017, 19 (14), 9485-9499.

22

22. Guo, H.; Zhang, W.; Lu, N.; Zhuo, Z.; Zeng, X. C.; Wu, X.; Yang, J., CO2 Capture onh-BN

23

Sheet with High Selectivity Controlled by External Electric Field. The Journal of Physical

24

Chemistry C 2015, 119 (12), 6912-6917.

25

23. Hu, W.; Xia, N.; Wu, X.; Li, Z.; Yang, J., Silicene as a highly sensitive molecule sensor for

26

NH3, NO and NO2. Phys. Chem. Chem. Phys. 2014, 16 (15), 6957-6962.

27

24. Li, X.; Yang, J., First-principles design of spintronics materials. National Science Review

28

2016, 3 (3), 365-381.

29

25. Li, Y.; Xu, L.; Liu, H.; Li, Y., Graphdiyne and graphyne: from theoretical predictions to

30

practical construction. Chem. Soc. Rev. 2014, 43 (8), 2572-2586.

31

26. Ivanovskii, A. L., Graphynes and graphdyines. Prog. Solid State Chem. 2013, 41 (1-2),

32

1-19.

33

27. Carvalho, A.; Wang, M.; Zhu, X.; Rodin, A. S.; Su, H.; Castro Neto, A. H., Phosphorene: 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

Page 20 of 24

1

from theory to applications. Nature Reviews Materials 2016, 1, 16061.

2

28. Ding, B.; Chen, W.; Tang, Z.; Zhang, J., Tuning Phosphorene Nanoribbon Electronic

3

Structure through Edge Oxidization. Journal of Physical Chemistry C 2016, 120, 2149-2158.

4

29. Sun, Q.; Dai, Y.; Ma, Y.; Wei, W.; Huang, B., Lateral heterojunctions within monolayer

5

h-BN/graphene: a first-principles study. RSC Advances 2015, 5 (42), 33037-33043.

6

30. Lin, Y.; Connell, J. W., Advances in 2D boron nitride nanostructures: nanosheets,

7

nanoribbons, nanomeshes, and hybrids with graphene. Nanoscale 2012, 4 (22), 6908-6939.

8

31. Andrew J. Mannix, X.-F. Z., Brian Kiraly, Joshua D. Wood, Diego Alducin, Benjamin

9

D. Myers, Xiaolong Liu, Brandon L. Fisher, Ulises Santiago, Jeffrey R. Guest, et al.

10

Synthesis of borophenes: Anisotropic, two-dimensional boron polymorphs. Nature 2015,

11

350 (6267), 1513-1516.

12

32. Gao, G.; Ding, G.; Li, J.; Yao, K.; Wu, M.; Qian, M., Monolayer MXenes: promising

13

half-metals and spin gapless semiconductors. Nanoscale 2016, 8 (16), 8986-8994.

14

33. Si,

15

Functionalization-Induced Metal–Insulator Transition in Graphene-like Two-Dimensional Cr2C

16

Crystals. ACS applied materials & interfaces 2015, 7 (31), 17510-17515.

17

34. He, J.; Lyu, P.; Nachtigall, P., New two-dimensional Mn-based MXenes with

18

room-temperature ferromagnetism and half-metallicity. Journal of Materials Chemistry C 2016,

19

4 (47), 11143-11149.

20

35. Zhao, T.; Zhou, J.; Wang, Q.; Kawazoe, Y.; Jena, P., Ferromagnetic and Half-Metallic

21

FeC2 Monolayer Containing C2 Dimers. ACS applied materials & interfaces 2016, 8 (39),

22

26207-26212.

23

36. Kuklin, A. V.; Kuzubov, A. A.; Kovaleva, E. A.; Mikhaleva, N. S.; Tomilin, F. N.; Lee, H.;

24

Avramov, P. V., Two-dimensional hexagonal CrN with promising magnetic and optical

25

properties: A theoretical prediction. Nanoscale 2017, 9 (2), 621-630.

26

37. Wang, Y.; Wang, S.-S.; Lu, Y.; Jiang, J.; Yang, S. A., Strain-Induced Isostructural and

27

Magnetic Phase Transitions in Monolayer MoN2. Nano Lett. 2016, 16 (7), 4576-4582.

28

38. Luxa, J.; Jankovsky, O.; Sedmidubsky, D.; Medlin, R.; Marysko, M.; Pumera, M.; Sofer, Z.,

29

Origin of exotic ferromagnetic behavior in exfoliated layered transition metal dichalcogenides

30

MoS2 and WS2. Nanoscale 2016, 8 (4), 1960-1967.

31

39. Zhang, X.; Bao, Z.; Ye, X.; Xu, W.; Wang, Q.; Liu, Y., Half-metallic properties of 3d

32

transition metal atom-intercalated graphene@MS2 (M = W, Mo) hybrid structures. Nanoscale

33

2017, 9 (31), 11231-11238.

C.;

Zhou,

J.;

Sun,

Z.,

Half-Metallic

Ferromagnetism

20

ACS Paragon Plus Environment

and

Surface

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

1

40. Xu, Y.; Liu, X.; Guo, W., Tensile strain induced switching of magnetic states in NbSe2 and

2

NbS2 single layers. Nanoscale 2014, 6 (21), 12929-12933.

3

41. He, J.; Ma, S.; Lyu, P.; Nachtigall, P., Unusual Dirac half-metallicity with intrinsic

4

ferromagnetism in vanadium trihalide monolayers. Journal of Materials Chemistry C 2016, 4

5

(13), 2518-2526.

6

42. Kulish, V. V.; Huang, W., Single-layer metal halides MX2 (X = Cl, Br, I): stability and

7

tunable magnetism from first principles and Monte Carlo simulations. Journal of Materials

8

Chemistry C 2017, 5 (34), 8734-8741.

9

43. Luan, H.-X.; Zhang, C.-W.; Li, S.-S.; Zhang, R.-W.; Wang, P.-J., First-principles study on

10

ferromagnetism in W-doped graphene. RSC Advances 2013, 3 (48), 26261-26265.

11

44. Maassen, J.; Ji, W.; Guo, H., Graphene spintronics: the role of ferromagnetic electrodes.

12

Nano Lett. 2011, 11 (1), 151-155.

13

45. Choudhuri, I.; Kumar, S.; Mahata, A.; Rawat, K. S.; Pathak, B., Transition-metal

14

embedded carbon nitride monolayers: high-temperature ferromagnetism and half-metallicity.

15

Nanoscale 2016, 8 (29), 14117-14126.

16

46. Ghosh, D.; Periyasamy, G.; Pandey, B.; Pati, S. K., Computational studies on magnetism

17

and the optical properties of transition metal embedded graphitic carbon nitride sheets. J.

18

Mater. Chem. C 2014, 2 (37), 7943-7951.

19

47. Sorokin, P. B.; Kvashnin, A. G.; Zhu, Z.; Tománek, D., Spontaneous graphitization of

20

ultrathin cubic structures: A computational study. Nano Lett. 2014, 14 (12), 7126-7130.

21

48. Blöchl, P. E., Projector augmented-wave method. Physical Review B 1994, 50 (24),

22

17953-17979.

23

49. Kohn, W.; Sham, L. J., Self-Consistent Equations Including Exchange and Correlation

24

Effects. Phys. Rev. 1965, 140 (4A), A1133-A1138.

25

50. Perdew, J. P.; Burke, K.; Ernzerhof, M., Generalized Gradient Approximation Made

26

Simple. Phys. Rev. Lett. 1996, 77 (18), 3865-3868.

27

51. Dudarev, S. L.; Botton, G. A.; Savrasov, S. Y.; Humphreys, C. J.; Sutton, A. P.,

28

Electron-energy-loss spectra and the structural stability of nickel oxide: An LSDA+U study.

29

Physical Review B 1998, 57 (3), 1505-1509.

30

52. Seo, D.-H.; Gwon, H.; Kim, S.-W.; Kim, J.; Kang, K., Multicomponent Olivine Cathode for

31

Lithium Rechargeable Batteries: A First-Principles Study. Chem. Mater. 2010, 22 (2), 518-523.

32

53. Heyd, J.; Scuseria, G. E.; Ernzerhof, M., Hybrid functionals based on a screened

33

Coulomb potential. The Journal of Chemical Physics 2003, 118 (18), 8207-8215. 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

Page 22 of 24

1

54. Henderson, T. M.; Paier, J.; Scuseria, G. E., Accurate treatment of solids with the HSE

2

screened hybrid. physica status solidi (b) 2011, 248 (4), 767-774.

3

55. Monkhorst, H. J.; Pack, J. D., Special points for Brillouin-zone integrations. Physical

4

Review B 1976, 13 (12), 5188-5192.

5

56. Baroni, S.; de Gironcoli, S.; Dal Corso, A.; Giannozzi, P., Phonons and related crystal

6

properties from density-functional perturbation theory. Rev. Mod. Phys. 2001, 73 (2), 515-562.

7

57. Togo, A.; Tanaka, I., First principles phonon calculations in materials science. Scripta

8

Mater. 2015, 108, 1-5.

9

58. Liu, C. S.; Zhu, H. H.; Ye, X. J.; Yan, X. H., Prediction of a new BeC monolayer with

10

perfectly planar tetracoordinate carbons. Nanoscale 2017, 9 (18), 5854-5858.

11

59. Teng, Z. W.; Liu, C. S.; Yan, X. H., A CO monolayer: first-principles design of a new direct

12

band-gap semiconductor with excellent mechanical properties. Nanoscale 2017, 9 (17),

13

5445-5450.

14

60. Born, M.; Huang, K., Dynamic Theory of Crystal Lattice, Clarendon. Oxford: 1954.

15

61. Wang, J.; Yip, S.; Phillpot, S. R.; Wolf, D., Crystal instabilities at finite strain. Phys. Rev.

16

Lett. 1993, 71 (25), 4182-4185.

17

62. Sun, Y.; Zhuo, Z.; Wu, X.; Yang, J., Room-Temperature Ferromagnetism in

18

Two-Dimensional Fe2Si Nanosheet with Enhanced Spin-Polarization Ratio. Nano Lett. 2017,

19

17 (5), 2771-2777.

20

63. Zhang, C.-w.; Yan, S.-s., First-Principles Study of Ferromagnetism in Two-Dimensional

21

Silicene with Hydrogenation. The Journal of Physical Chemistry C 2012, 116 (6), 4163-4166.

22

64. Srivastava, A. M.; Brik, M. G., Crystal field studies of the Mn

23

perovskite, LaAlO3. Opt. Mater. 2013, 35 (8), 1544-1548.

24

65. Brik, M. G.; Srivastava, A. M., Comparative crystal field analysis of energy level schemes

25

and nephelauxetic effect for Cr4+, Cr3+, and Mn4+ ions in Y2Sn2O7 pyrochlore. Opt. Mater. 2013,

26

35 (6), 1251-1256.

27

66. Zhuang, H. L.; Xie, Y.; Kent, P. R. C.; Ganesh, P., Computational discovery of

28

ferromagnetic semiconducting single-layer CrSnTe3. Physical Review B 2015, 92 (3), 035407.

29

67. Carteaux, V.; Moussa, F.; Spiesser, M., 2D Ising-Like Ferromagnetic Behaviour for the

30

Lamellar Cr2Si2Te6 Compound: A Neutron Scattering Investigation. Europhysics Letters. 1995,

31

29 (3), 251-256.

32

68. Kan, M.; Zhou, J.; Sun, Q.; Kawazoe, Y.; Jena, P., The Intrinsic Ferromagnetism in a

33

MnO2 Monolayer. The journal of physical chemistry letters 2013, 4 (20), 3382-3386.

4+

22

ACS Paragon Plus Environment

energy levels in the

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

1

69. Zhou, Y.; Wang, Z.; Yang, P.; Zu, X.; Yang, L.; Sun, X.; Gao, F., Tensile strain switched

2

ferromagnetism in layered NbS2 and NbSe2. ACS Nano 2012, 6 (11), 9727-9736.

3

70. Liu, C.-S.; Teng, Z.-W.; Ye, X.-J.; Yan, X.-H., Two-dimensional tetragonal AlP monolayer:

4

strain-tunable direct–indirect band-gap and semiconductor–metal transitions. Journal of

5

Materials Chemistry C 2017, 5 (24), 5999-6004.

6

71. Conley, H. J.; Wang, B.; Ziegler, J. I.; Haglund, R. F., Jr.; Pantelides, S. T.; Bolotin, K. I.,

7

Bandgap engineering of strained monolayer and bilayer MoS2. Nano Lett. 2013, 13 (8),

8

3626-3630.

9

72. Zhu, L.; Wang, S.-S.; Guan, S.; Liu, Y.; Zhang, T.; Chen, G.; Yang, S. A., Blue

10

Phosphorene Oxide: Strain-tunable Quantum Phase Transitions and Novel 2D Emergent

11

Fermions. Nano Lett. 2016, 16 (10), 6548-6554.

12

73. He, K.; Poole, C.; Mak, K. F.; Shan, J., Experimental demonstration of continuous

13

electronic structure tuning via strain in atomically thin MoS2. Nano Lett. 2013, 13 (6),

14

2931-2936.

15

74. Scalise, E.; Houssa, M.; Pourtois, G.; Afanas’ev, V.; Stesmans, A., Strain-induced

16

semiconductor to metal transition in the two-dimensional honeycomb structure of MoS2. Nano

17

Research 2011, 5 (1), 43-48.

18

75. Peng, X.; Wei, Q.; Copple, A., Strain-engineered direct-indirect band gap transition and its

19

mechanism in two-dimensional phosphorene. Physical Review B 2014, 90 (8), 085402.

20

76. Peng, X.; Copple, A., Origination of the direct-indirect band gap transition in strained

21

wurtzite and zinc-blende GaAs nanowires: A first principles study. Physical Review B 2013, 87

22

(11), 115308.

23

77. Guo, H.; Lu, N.; Wang, L.; Wu, X.; Zeng, X. C., Tuning Electronic and Magnetic Properties

24

of Early Transition-Metal Dichalcogenides via Tensile Strain. The Journal of Physical

25

Chemistry C 2014, 118 (13), 7242-7249.

26

78. Gorria, P.; Martínez-Blanco, D.; Pérez, M. J.; Blanco, J. A.; Hernando, A.; Laguna-Marco,

27

M. A.; Haskel, D.; Souza-Neto, N.; Smith, R. I.; Marshall, W. G.; et al. Stress-induced large

28

Curie temperature enhancement in Fe64Ni36 Invar alloy. Physical Review B 2009, 80 (6),

29

064421.

30 31

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

1 2

TOC Graphic

24

ACS Paragon Plus Environment

Page 24 of 24