Subscriber access provided by Brought to you by ST ANDREWS UNIVERSITY LIBRARY
Article
Spin Dependent Transport and Optical Properties of Transparent Half Metallic g-CN Films 4
3
Arqum Hashmi, M. Umar Farooq, Tao Hu, and Jisang Hong J. Phys. Chem. C, Just Accepted Manuscript • DOI: 10.1021/jp510179p • Publication Date (Web): 05 Dec 2014 Downloaded from http://pubs.acs.org on December 7, 2014
Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.
The Journal of Physical Chemistry C 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 25
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
Spin dependent transport and optical properties of
2
transparent half-metallic g-C4N3 films
3
Arqum Hashmi, M. Umar Farooq, Tao Hu, and Jisang Hong
4
Department of Physics, Pukyong National University, Busan 608-737, Korea
5
ABSTRACT
6
Half-metallic materials play an essential role in the spin dependent transport study and so far it is
7
believed that the 3d transitional metal elements are necessary to obtain ferromagnetic state.
8
Herein, we demonstrate that the metal free two-dimensional graphitic carbon nitride (g-C4N3)
9
film show half metallic state and this is independent of the film thickness and stacking order. We
10
find that the half-metallicity results in a very large magnetoresistance even at finite bias.
11
Moreover, we have also investigated various optical properties. The g-C4N3 thin films have a
12
very high refractive index for parallel polarization and extremely weak reflectivity for
13
perpendicular polarization. Additionally, the g-C4N3 thin films are optically transparent because
14
of very weak absorption coefficients. Overall, our findings suggest that the transparent half-
15
metallic g-C4N3 thin films can be used for both spintronics and optical applications because of
16
potential multifunctionality.
17
ACS Paragon Plus Environment
1
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 2 of 25
18 19
INTRODUCTION
20
Spintronics is one of the most promising fields in condensed matter physics because the spin
21
degree of freedom can provoke various applications of innovative technologies. Regarding the
22
spintronics, half-metallic materials have many advantages because one spin channel is
23
conductive, but the other spin channel is insulating.1,2 This will result in a giant
24
magnetoresistance which is an essential quantity for spintronics. Various types of materials show
25
half-metallic state, but most of them are based on the conventional 3d transition metal
26
elements.3–8
27 28
In general, 2p electron systems show a non-magnetic behavior in a bulk state. Of course, several
29
studies proposed defect induced magnetic state in two-dimensional geometry.9–11 However, those
30
studies were about local magnetic state, not ferromagnetism arising from long range ordering. In
31
the spin dependent transport study, a long spin relaxation time is an important factor for spin
32
information delivery. By virtue of the fact that the spin relaxation time is inversely proportional
33
to the strength of spin-orbit coupling, the half-metallicity in 2p electron system will bring an
34
interesting issue in the spin dependent transport property. Nonetheless, it is rare to find such
35
reports on the half-metallicity from pure 2p electronic materials.12,13
36 37
Interestingly, the graphitic carbon nitride, namely g-C4N3, was synthesized and showed a half-
38
metallic state in free standing monolayer.14,15 After then, many other works were reported on
39
free standing structure.16–20 Recently, the influence of Li adsorption and the mechanism of stable
40
adsorption of hybridized system with BN layer were presented as well.21,22 However,
ACS Paragon Plus Environment
2
Page 3 of 25
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
41
manipulating a single layer structure will require highly precise experimental technique and also
42
complex fabrication process for device applications. It will be more practical, realistic, and easier
43
to treat thicker layers than free standing monolayer. So far, no studies on this issue have been
44
performed. In this respect, our objective is to propose that the g-C4N3 film can maintain half-
45
metallic state even in multilayer structure and this half-metallicity can be applied for spintronics.
46
We will prove this by calculating the spin dependent current and magnetoresistance at finite bias.
47
Along with these magnetic properties, if the multilayer g-C4N3 system can be utilized for optical
48
device applications, it will be a remarkable finding and attract great research interests in diverse
49
fields in condensed matter physics because of potential multifunctionality. Thus, we will study
50
the frequency dependent optical properties such as dielectric functions, refractive index,
51
reflectivity, and absorption coefficients.
52 53
NUMERICAL METHOD
54
We employed the Vienna ab initio simulation package (VASP) using a plane wave basis
55
set.23,24 The core electron interactions were described using the projector augmented wave
56
(PAW) method. For layer structured materials, it is of necessity to include van der Waals
57
interaction. Thus, in this report, we explore the structural, adsorption, electronic and magnetic
58
properties of multilayer g-C4N3 using non-empirical vdW-DF method proposed by Dion et al.25
59
as implemented in the VASP code. We used most recently developed less repulsive optB88-like
60
exchange functional method
61
interlayer distances were achieved.25,27,28 A vacuum distance of 15 Å was imposed and we
62
carried the self-consistent calculations with 9x9x1 k-mesh. All results have been obtained with
63
very high plane-wave energy cut-off of 700 eV. The structure optimization has been performed
26
because significant improvements in binding energies and
ACS Paragon Plus Environment
3
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 4 of 25
64
until the residual forces on individual atoms were smaller than 0.01 eVÅ-1. The convergence
65
criterion for energy was set to 0.1 meV. It is worth noting that the g-C4N3 is known to have a
66
(2x2) reconstructed surface geometry. Thus, the ground state lattice constant of g-C4N3 is 9.68 Å
67
and entire calculations were performed with this lattice constant. In addition, ferromagnetic
68
(FM), antiferromagnetic (AFM) and non-magnetic (NM) spin configurations were considered to
69
find out the magnetic ground state. The thickness varies from bilayer to 4 monolayers (MLs). In
70
this type of layer structured material, it is necessary to find the stacking order. Thus, three
71
different stacking orders, so called, AA, AB and AC for bilayer system and ABA, ABB and
72
ABC stacking orders for the trilayer film were examined. Using the result of trilayer structure,
73
we adopted the stacking order of ABAB for 4 ML film. For the transport calculations, we
74
employ the non-equilibrium Green's function (NEGF) method based on density functional
75
theory. All of the electron transport calculations were performed using the density functional
76
theory code OpenMX, which is based on localized pseudoatomic orbitals (PAOs) and norm-
77
conserving pseudopotentials.29,30 Two primitive orbitals (s and p orbitals) generated by a
78
confinement scheme with cutoff of 5 Bohr radius are employed.31,32 The Perdew-Burke-
79
Ernzerhof (PBE) exchange-correlation functional, which is derived within the generalized
80
gradient approximation (GGA) is used.33 To obtain bias dependent current, we have divided the
81
sheet of g-C4N3 into three parts; two leads and scattering part.34 The central sheet sandwiched
82
between two g-C4N3 leads is used for the NEGF calculations. For spin dependent transport study,
83
we changed the relative spin direction in the two leads and the bias dependent current in both
84
parallel and antiparallel configurations is obtained
85
ACS Paragon Plus Environment
4
Page 5 of 25
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
86
NUMERICAL RESULTS
87
The g-C4N3 bilayer structure can have three different configurations such as AA, AB, and AC
88
stacking as shown in Figure 1 (a) - (c). For clarification, we also present the (1x1) cell of free
89
standing g-C4N3 layer in Figure 1(d). The C and N atoms located on the hexagonal ring are
90
referred as C1 and N1, respectively, whereas C atom located outside the ring is denoted as C2. In
91
AA stacking, carbon and nitrogen atoms of two layers have the same lateral positions. In AB
92
stacking, C1 and N1 atoms are on the hexagonal rings of the other layer while C2 atoms have no
93
direct coupling. In AC stacking, carbon and nitrogen atoms of both layers are on the bridge sites
94
of each other. Figure 1 (e) – (g) shows three different stacking orders for the trilayer film. Here,
95
we considered ABA, ABB, and ABC stacking. Since the AB stacking was the most stable
96
configuration in bilayer system, we constructed trilayer structure by adding one more layer to the
97
AB stacked bilayer g-C4N3. Table 1, the ABA stacking became the most stable structure. Thus,
98
this suggests that the multilayer film may display ABAB stacking with increase of the film
99
thickness. So, we only considered ABAB stacking in 4 ML film.
100 101
Table I presents the total energy differences among FM, AFM and NM spin arrangements and
102
magnetic moment per (1x1) cell. Here, it should be remarked that we have calculated the total
103
energy with (2x2) cell and the same spin configurations are considered as Du et al. did.15 In
104
bilayer and trilayer systems, the total energy of AB and ABA stacking in FM state was set to
105
zero, respectively as a reference energy.
106
ACS Paragon Plus Environment
5
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 6 of 25
107
TABLE 1: Stacking order dependent total energy difference (in meV/cell) and energy
108
difference per atom (in meV) shown in parentheses. The magnetic moment (in µB) per (1x1)
109
unit cell . FM
AFM
NM
µB / cell
AB
0
365 (6.5)
357 (6.3)
1.000
AC
29 (0.5)
340 (6.1)
342 (6.1)
0.968
AA
42 (0.8)
328 (5.8)
328 (5.8)
0.921
ABA
0
518 (6.2)
552 (6.6)
1.000
ABB
62 (0.7)
532 (6.3)
534 (6.3)
0.935
ABC
71 (0.8)
563 (6.7)
568 (6.8)
0.969
0
700 (6.3)
735 (6.6)
1.000
Stacking 2ML
3ML
4ML ABAB 110 111
The most stable bilayer structure is found in the AB stacking. FM ground state was still
112
preserved in the bilayer structure and the magnetic ground state is independent of the stacking
113
order. The calculated total magnetic moment of the bilayer g-C4N3 film in 2x2 supercell is 8.00
114
µB , 7.75 µB and 7.37 µB in AB, AC and AA stacking respectively. Note that the magnetic
115
moment of pristine free standing g-C4N3 monolayer is 1.0 µB per (1x1) unit cell.15,19,22 In bilayer
116
system, we found an integer magnetic moment in AB stacking. This implies that the most stable
117
AB stacking film has a half-metallic state. On the other hand, we find the same behavior even in
118
3 and 4 ML films. This result tells that the g-C4N3 displays half-metallic state in multilayer
119
films. From Table I, we find that the order of spin polarization in bilayer is AB>AC>AA while it
120
becomes ABA>ABC>ABB in trilayer structures. There is an interesting relationship between
ACS Paragon Plus Environment
6
Page 7 of 25
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
121
spin polarization and stacking order. First of all, the AB stacking order produces the largest spin
122
polarization and the smallest spin polarization is obtained in identical stacking (AA order) in
123
bilayer system. If we add one more layer on AB bilayer structure to form a trilayer film, the ABB
124
stacking will have the least spin polarization because of the identical BB stacking. In addition,
125
the BC arrangement in ABC stacking is like an AC bilayer stacking if we start from the layer B.
126
On the other hand, BA arrangement in ABA stacking is equivalent to AB bilayer stacking with
127
respect to the B layer. Thus, we find the order of spin polarization such as ABA>ABC>ABB.
128 129
It is noteworthy that the stability of FM state seems increasing as the film thickness increases.
130
However, in reality, the energy difference per atom remains almost the same as shown in
131
parentheses of Table 1. Hence, the ferromagnetic stability shows no significant thickness
132
dependency. Indeed, we have checked another AFM configuration. In this case, the AFM
133
interaction between different layers of g-C4N3 is considered. The FM state is still stable spin
134
configuration and the calculated energy difference was almost 20 meV as compared with AFM
135
state. In addition, this value is also thickness and stacking order independent. We now show the
136
geometric feature and adsorption energy,
137
Ntot represent the total energy of the fully optimized N layer g-C4N3 film, the energy of free
138
standing single layer g-C4N3, the number of layers and the total number of atoms in the supercell,
139
respectively. The optimized interlayer distances and adsorption energies are presented in Table 2.
140
TABLE 2: Optimized interlayer distances (in Å) and adsorption energies (Ea in meV/atom) 2ML AB stacking AC stacking AA stacking
=
( ( × )
where EN, , N and
Interlayer Distance
3.27
3.30
3.32
Ea
-0.032
-0.031
-0.031
ACS Paragon Plus Environment
7
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 25
3ML
ABA stacking
ABB stacking
ABC stacking
Interlayer Distance
(AB) 3.27
(AB) 3.30
(AB) 3.30
(BA) 3.27
(BB) 3.31
(BC) 3.31
Ea
-0.044
-0.043
-0.043
4ML
ABAB stacking
_
_
Interlayer Distance
3.27
_
_
Ea
-0.051
_
_
141 142
We have observed that the multilayer films still maintain corrugated layer structure as found in
143
free standing single layer g-C4N3 structure.15,21 From the calculated adsorption energy, one can
144
understand that the g-C4N3 multilayer is a weakly bounded system. This suggests that the main
145
interaction among layers should be described by the van der Waals interaction because the
146
adsorption energy per atom is very small as compared with that of bilayer graphene.35
147 148
Figure 2 and 3 shows the band structures and density of states (DOS) for bilayer and trilayer
149
respectively. Indeed, the band structures and DOS in 4ML films were almost the same as those
150
of bilayer and trilayer systems, thus we here only presented the results of bilayer and trilayer
151
films. Band structure and DOS of AB and ABA stacking of bilayer and trilayer show the half
152
metallic behavior as shown in Figure 2(a) and 3(a), respectively. One can clearly find that the
153
majority spin electrons are fully occupied below the Fermi level, whereas the minority spin
154
bands cross the Fermi level. In AA and AC configurations for bilayer and ABB and ABC
155
configurations for trilayer, we find a similar trend although the majority spin bands barely cross
156
the Fermi level. Nonetheless, the minority spin bands maintain the same features as found in AB
157
and ABA stacking. This can be confirmed by the spectral shape of DOS. This half-metallic
ACS Paragon Plus Environment
8
Page 9 of 25
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
158
characteristic in thin film geometry will bring an interesting issue in the spin dependent transport
159
phenomenon because a giant magnetoresistance can be achieved.
160 161
We now discuss the spin dependent transport properties at finite bias. In Figure 4 (a) and (b), the
162
schematic illustration for two leads and scattering part for both parallel and antiparallel
163
configurations is displayed, respectively. We apply the bias voltages along the negative x-axis,
164
so the electrons move to the positive x-axis. The periodic boundary condition is used along the y-
165
axis and the vacuum distance of 12 Å is applied along the perpendicular to the film surface. We
166
performed the NEGF calculation by taking 150 numbers of poles along the x-direction and 9 k-
167
points along the y direction and the convergence criterion is kept 10-8 Hartree. It is obvious that
168
the electronic band structure controls the transport property of a material. As we discussed
169
already, the most stable configuration of multilayer g-C4N3 shows the half-metallic state and no
170
substantial changes in their ground state band structures were found. Thus, we expect that the
171
essential transport properties for multilayer g-C4N3 will be the same. Hence, we only calculated
172
the transport properties of single layer g-C4N3. In Figure 4 (c), we present the spin dependent I-V
173
curve originated from each spin component as a function of external bias. In parallel spin
174
configuration, the majority spin current is completely suppressed and this is a straightforward
175
result because the majority spin band has a gap as shown in Figure. 2 and only the minority spin
176
electrons contribute to the current. In antiparallel spin configuration, the current is significantly
177
suppressed because of the reversed spin direction in the opposite lead. From these I-V curves of
178
each spin component, we calculate bias dependent magnetoresistance (MR) ratio defined by
179
MR(%) = (Ip-Iap)/Ip where the Ip(ap) represents the total current of parallel (antiparallel) spin
180
configuration. Then, the maximum MR ratio is 100 % in the theoretical point of view when the
ACS Paragon Plus Environment
9
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 25
181
total current in an antiparallel configuration is zero at any bias voltage. However, in general, the
182
Iap can be measured at finite bias, so the 100 % MR ratio is hardly achievable at finite bias. In
183
our system, we find that the MR ratio is still large enough even at bias voltage of 0.5 eV and this
184
giant MR ratio suggests that the g-C4N3 film can be an ideal candidate for spintronics.
185 186
To explore the optical properties, the most important quantity is the frequency dependent
187
dielectric function () = () + () because various optical properties such as
188
refractive index, reflectivity, and absorption function can be extracted from this value. Here, we
189
have considered both parallel (E )׀׀and perpendicular (E⊥) electric field polarization to the film
190
surface. From the calculated electronic band structures, one can easily understand that the
191
contribution to the dielectric function due to the intraband transition is negligible if the photon
192
energy is larger than approximately 0.6 eV. Thus, the interband transition plays an essential role
193
in a wide range of photon energy for the frequency dependent dielectric function. Additionally,
194
we expect very weak interband optical transition in the visible range because of the band
195
characteristic. This implies that the g-C4N3 films will be optically transparent. Thus, the g-C4N3
196
films can be utilized for both spintronics and optical application purposes.
197 198
Figure 5 shows the real and imaginary part of frequency dependent dielectric functions for both
199
bilayer and trilayer by changing the stacking order. At a zero frequency, the static dielectric
200
functions (ε1(0)) are subject to the thickness and polarization. In bilayer, the ε1(0) is in the range
201
from 13 to 25 for parallel polarization while it is about 4 for perpendicular polarization. In
202
trilayer, the ε1(0) is in the range from 35 to 50 for parallel polarization while it is about 6 for
203
perpendicular polarization. As the photon energy increases, the thickness and stacking order
ACS Paragon Plus Environment
10
Page 11 of 25
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
204
dependence of the dielectric function disappear, but only the electric field polarization
205
dependence remained. One can see noticeable asymmetric behavior in the imaginary part of the
206
dielectric function. This may indicate an anisotropic optical behavior with respect to the electric
207
polarization. In Figure 6 (a) and (b), the calculated refractive index is presented for bilayer and
208
trilayer films, respectively. The real part of the refractive index of parallel electric polarization
209
was larger than that of perpendicular polarization and the refractive index of perpendicular
210
polarization is frequency dependent, but it is almost constant of 1.2 in the visible range for
211
parallel polarization. Interestingly, we find very large refractive index of 2.6 near the blue light
212
for parallel polarization. This large value of refractive index can have potential application for
213
optical device. The reflectivity value is also an important quantity for reflecting and anti-
214
reflecting applications. For instance, it is highly desirable to obtain zero reflectivity for anti-
215
reflection coating purposes. As shown in Figure 6(c) and (d), the perpendicular polarization
216
component has very weak reflectivity. On the other hand, we can find high reflectivity in the
217
ultraviolet region for parallel polarization. In Figure 6(e) and (f), we present the frequency
218
dependent absorption coefficients. We also find an anisotropic behavior in the frequency
219
dependent functions. As displayed, the absorption coefficient for perpendicular polarization is
220
quite small compared with that of parallel polarization. Since the intensity of propagating light
221
decays exponentially as a function of length, the calculated absorption coefficients indicate that
222
the g-C4N3 films can be optically transparent even in much thicker films.
223 224
CONCLUSION
225
In summary, we investigated the structural, magnetic, optical and spin dependent transport
226
properties of g-C4N3 films. The most stable structure was found in AB type stacking order. The
ACS Paragon Plus Environment
11
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 25
227
adsorption energy calculations indicated that the interlayer interaction should be described by the
228
van der Waals interaction. We found metal free half metallic behavior and the essential
229
electronic structure was insensitive to the film thickness and stacking order. Due to the half-
230
metallic behavior, the minority spin electrons contribute to the spin dependent current at finite
231
bias in parallel spin configuration whereas the current is strongly suppressed in an antiparallel
232
configuration. Consequently, we find very large magnetoresistance up to the bias voltage of 0.5
233
eV. Along with these magnetic properties, we also explored the optical properties. In general,
234
the optical properties are insensitive to the film thickness and stacking order, but they are
235
strongly subjective to the electric polarization. For instance, the refractive index for
236
perpendicular polarization seems frequency independent in the visible range. Interestingly, a
237
large refractive index of 2.6 was found in parallel polarization, whereas an extremely low
238
reflectivity was obtained in perpendicular polarization. In addition, the multilayer g-C4N3 films
239
are optically transparent. Overall, we propose that the metal free transparent half metallic g-C4N3
240
films can be utilized for spintronics and optical device applications.
241
ACS Paragon Plus Environment
12
Page 13 of 25
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
242
AUTHOR INFORMATION
243
Corresponding Author
244
*E-mail:
[email protected] 245
Notes
246
The authors declare no competing financial interest.
247
ACKNOWLEDGMENT
248
This research was supported by Basic Science Research Program through the National
249
Research Foundation of Korea(NRF) funded by the Ministry of Education, Science and
250
Technology (No. 2013R1A1A2006071) and by the Supercomputing Center/Korea Institute of
251
Science and Technology Information with supercomputing resources including technical support
252
(KSC-2014-C3-015)
ACS Paragon Plus Environment
13
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
253 254 255 256
Page 14 of 25
REFERENCES (1)
De Groot, R. A.; Mueller, F. M.; Engen, P. G. van; Buschow, K. H. J. New Class of
Materials: Half-Metallic Ferromagnets. Phys. Rev. Lett. 1983, 50, 2024–2027. (2)
Groot, R. A. de; Mueller, F. M.; Engen, P. G. van; Buschow, K. H. J. Half‐metallic
257
Ferromagnets and Their Magneto‐optical Properties (invited). J. Appl. Phys. 1984, 55, 2151–
258
2154.
259
(3)
Kang, J.-S.; Park, J.-G.; Olson, C. G.; Youn, S. J.; Min, B. I. Valence Band and Sb 4d
260
Core Level Photoemission of the XMnSb-Type Heusler Compounds (X=Pt,Pd,Ni). J. Phys.
261
Condens. Matter 1995, 7, 3789.
262 263 264 265 266
(4)
Hanssen, K. E. H. M.; Mijnarends, P. E. Positron-Annihilation Study of the Half-Metallic
Ferromagnet NiMnSb: Theory. Phys. Rev. B 1986, 34, 5009–5016. (5)
Galanakis, I.; Dederichs, P. H.; Papanikolaou, N. Slater-Pauling Behavior and Origin of
the Half-Metallicity of the Full-Heusler Alloys. Phys. Rev. B 2002, 66, 174429. (6)
Medvedeva, J. E.; Freeman, A. J.; Cui, X. Y.; Stampfl, C.; Newman, N. Half-Metallicity
267
and Efficient Spin Injection in AlN/GaN∶Cr (0001) Heterostructure. Phys. Rev. Lett. 2005, 94,
268
146602.
269
(7)
270 271
Mallajosyula, S. S.; Pati, S. K. Vanadium−Benzimidazole-Modified sDNA: A One-
Dimensional Half-Metallic Ferromagnet. J. Phys. Chem. B 2007, 111, 13877–13880. (8)
Du, C.; Adur, R.; Wang, H.; Hauser, A. J.; Yang, F.; Hammel, P. C. Control of
272
Magnetocrystalline Anisotropy by Epitaxial Strain in Double Perovskite Sr2FeMoO6 Films.
273
Phys. Rev. Lett. 2013, 110, 147204.
ACS Paragon Plus Environment
14
Page 15 of 25
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
274 275 276 277 278 279 280 281
The Journal of Physical Chemistry
(9)
Kim, D.; Hashmi, A.; Hwang, C.; Hong, J. Magnetization Reversal and Spintronics of
Ni/Graphene/Co Induced by Doped Graphene. Appl. Phys. Lett. 2013, 102, 112403. (10) Yang, J.; Kim, D.; Hong, J.; Qian, X. Magnetism in Boron Nitride Monolayer: Adatom and Vacancy Defect. Surf. Sci. 2010, 604, 1603–1607. (11) Si, M. S.; Xue, D. S. Magnetic Properties of Vacancies in a Graphitic Boron Nitride Sheet by First-Principles Pseudopotential Calculations. Phys. Rev. B 2007, 75, 193409. (12) Li, Y.; Zhou, Z.; Shen, P.; Chen, Z. Spin Gapless Semiconductor−Metal−Half-Metal Properties in Nitrogen-Doped Zigzag Graphene Nanoribbons. ACS Nano 2009, 3, 1952–1958.
282
(13) Dutta, S.; Pati, S. K. Half-Metallicity in Undoped and Boron Doped Graphene
283
Nanoribbons in the Presence of Semilocal Exchange-Correlation Interactions. J. Phys. Chem. B
284
2008, 112, 1333–1335.
285
(14) Lee, J. S.; Wang, X.; Luo, H.; Dai, S. Fluidic Carbon Precursors for Formation of
286
Functional Carbon under Ambient Pressure Based on Ionic Liquids. Adv. Mater. 2010, 22, 1004–
287
1007.
288 289 290 291 292 293
(15) Du, A.; Sanvito, S.; Smith, S. C. First-Principles Prediction of Metal-Free Magnetism and Intrinsic Half-Metallicity in Graphitic Carbon Nitride. Phys. Rev. Lett. 2012, 108, 197207. (16) Sun, S.-J. Interior Edges Induced Half-Metallic Ferromagnetism in Graphitic Carbon Nitride Structures. J. Magn. Magn. Mater. 2013, 344, 39–43. (17) Li, Y.; Sanvito, S.; Hou, S. Origin of the Half-Metallic Properties of Graphitic Carbon Nitride in Bulk and Confined Forms. J. Mater. Chem. C 2013, 1, 3655–3660.
ACS Paragon Plus Environment
15
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
294 295 296 297 298 299 300 301 302 303 304 305 306 307 308 309 310 311
Page 16 of 25
(18) Li, X.; Zhou, J.; Wang, Q.; Kawazoe, Y.; Jena, P. Patterning Graphitic C–N Sheets into a Kagome Lattice for Magnetic Materials. J. Phys. Chem. Lett. 2013, 4, 259–263. (19) Li, X.; Zhang, S.; Wang, Q. Stability and Physical Properties of a Tri-Ring Based Porous G-C4N3 Sheet. Phys. Chem. Chem. Phys. 2013, 15, 7142–7146. (20) Wu, M.; Wang, Q.; Sun, Q.; Jena, P. Functionalized Graphitic Carbon Nitride for Efficient Energy Storage. J. Phys. Chem. C 2013, 117, 6055–6059. (21) Hashmi, A.; Hu, T.; Hong, J. Transition from Half Metal to Semiconductor in Li Doped G-C4N3. J. Appl. Phys. 2014, 115, 124312. (22) Hashmi, A.; Hong, J. Metal Free Half Metallicity in 2D System: Structural and Magnetic Properties of G-C4N3 on BN. Sci. Rep. 2014, 4. (23) Kresse, G.; Furthmüller, J. Efficiency of Ab-Initio Total Energy Calculations for Metals and Semiconductors Using a Plane-Wave Basis Set. Comput. Mater. Sci. 1996, 6, 15–50. (24) 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. (25) Dion, M.; Rydberg, H.; Schröder, E.; Langreth, D. C.; Lundqvist, B. I. Van Der Waals Density Functional for General Geometries. Phys. Rev. Lett. 2004, 92, 246401. (26) Klimeš, J.; Bowler, D. R.; Michaelides, A. Chemical Accuracy for the van Der Waals Density Functional. J. Phys. Condens. Matter 2010, 22, 022201.
ACS Paragon Plus Environment
16
Page 17 of 25
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
312
(27) Graziano, G.; Klimeš, J.; Fernandez-Alonso, F.; Michaelides, A. Improved Description of
313
Soft Layered Materials with van Der Waals Density Functional Theory. J. Phys. Condens.
314
Matter 2012, 24, 424216.
315
(28) Liu, W.; Carrasco, J.; Santra, B.; Michaelides, A.; Scheffler, M.; Tkatchenko, A. Benzene
316
Adsorbed on Metals: Concerted Effect of Covalency and van Der Waals Bonding. Phys. Rev. B
317
2012, 86, 245405.
318 319 320 321 322 323 324 325 326 327 328 329 330 331
(29) The Code, OpenMX, Pseudo-Atomic Basis Functions and Pseudopotentials Are Available on Web Site (http://www.openmx-Square). (30) Troullier, N.; Martins, J. L. Efficient Pseudopotentials for Plane-Wave Calculations. Phys. Rev. B 1991, 43, 1993–2006. (31) Ozaki, T. Variationally Optimized Atomic Orbitals for Large-Scale Electronic Structures. Phys. Rev. B 2003, 67, 155108. (32) Ozaki, T.; Kino, H. Numerical Atomic Basis Orbitals from H to Kr. Phys. Rev. B 2004, 69, 195113. (33) Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865–3868. (34) Ozaki, T.; Nishio, K.; Kino, H. Efficient Implementation of the Nonequilibrium Green Function Method for Electronic Transport Calculations. Phys. Rev. B 2010, 81, 035116. (35) Birowska, M.; Milowska, K.; Majewski, J. A. Van Der Waals Density Functionals for Graphene Layers and Graphite Acta Physica Polonica A 2011, 845.
ACS Paragon Plus Environment
17
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 25
332
Figure 1. Schematic illustrations of (a) AA stacking (b) AB stacking (c) AC stacking in bilayer
333
film (d) 1x1 unit cell of g-C4N3 (e) ABA stacking (f) ABB stacking (g) ABC stacking in trilayer
334
film.
335
Figure 2. Calculated Band structures of (a) AB stacking (b) AC stacking (c) AA stacking and
336
DOS of (d) AB stacking (e) AC stacking (f) AA stacking.
337
Figure 3. Calculated Band structures of (a) ABA stacking (b) ABB stacking (c) ABCstacking
338
and DOS of (d) ABA stacking (e) ABB stacking (f) ABC stacking.
339
Figure 4. Schematic illustrations of device components for transport calculations with two leads
340
and scattering part (central region). The arrow indicates the direction of electron flow in (a)
341
parallel spin configuration (b) antiparallel spin configurations. (c) I-V curve and
342
magnetoresistance (MR) with the applied voltage. The left and right scales are for current and
343
MR ratio, respectively.
344
Figure 5. Real part of the frequency dependent dielectric function of (a) bilayer (b) trialyer
345
and imaginary part of the frequency dependent dielectric function of (c) bilayer (d) trilayer.
346
Figure 6. Frequency and polarization dependent (a) - (b) refractive index (c) - (d) reflectivity (e)
347
–(f) absorption coefficient for bilayer and trilayer, respectively.
348 349 350 351
ACS Paragon Plus Environment
18
Page 19 of 25
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
352
The Journal of Physical Chemistry
Table of Contents (TOC)
353
ACS Paragon Plus Environment
19
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
ACS Paragon Plus Environment
Page 20 of 25
Page 21 of 25 3
3
2
2
2
E n e r g y (e V )
(a )
(b )
(c )
1
1
1
0
0
0
-1
-1
-1
-2
-2
-2
G
-3
K
G
M
G
-3
4 0
M
G
-2 0
-2
0
E n e r g y (e V )
2
G
M
(f)
2 0
2 0
0
0
-2 0
-2 0
-4 0 4
K
4 0
0
-4
G
(e )
2 0
-4 0
K
-3
4 0
(d )
D O S ( S t a t e s /e V .a t o m .s p in )
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
The Journal of Physical Chemistry 3
-4
-2 0 2 ACS ParagonE Plus Environment n e r g y (e V )
-4 0 4
-4
-2
0
E n e r g y (e V )
2
4
The Journal of Physical Chemistry 3
3
3
2
2
2
E n e r g y (e V )
(a )
(b )
(c )
1
1
1
0
0
0
-1
-1
-1
-2
-2
-2
G
-3
4 0
D O S ( S t a t e s /e V .a t o m .s p in )
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
Page 22 of 25
K
G
M
G
-3
4 0
(d )
K
M
G
G
-3
4 0
(e )
2 0
2 0
0
0
0
-2 0
-2 0
-2 0
-4 0
-4 0
-4 0
-2
0
E n e r g y (e V )
2
4
-4
-2 0 2 E n Plus e r g y Environment (e V ) ACS Paragon
4
G
M
(f)
2 0
-4
K
-4
-2
0
E n e r g y (e V )
2
4
Page 23 of 25
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
The Journal of Physical Chemistry
ACS Paragon Plus Environment
The Journal of Physical Chemistry
5 0
E
(a )
5 0
//
A A -s ta c k in g A B -s ta c k in g A C -s ta c k in g
4 0
E
//
A B A -s ta c k in g A B B -s ta c k in g A B C -s ta c k in g E
⊥
3 0
A A -s ta c k in g A B -s ta c k in g A C -s ta c k in g
2 0
E
(b )
4 0
⊥
3 0
ε1(ω)
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
Page 24 of 25
A B A -s ta c k in g A B B -s ta c k in g A B C -s ta c k in g
1 0
2 0
1 0
0
0
-1 0
-1 0 0
1
2
3
4
5
3 0
0
1
2
3
4
5
3 0
E
//
(c )
A A -s ta c k in g A B -s ta c k in g A C -s ta c k in g E
2 0
⊥
E
//
(d )
A B A -s ta c k in g A B B -s ta c k in g A B C -s ta c k in g E
2 0
⊥
ε2(ω)
A A -s ta c k in g A B -s ta c k in g A C -s ta c k in g
1 0
A B A -s ta c k in g A B B -s ta c k in g A B C -s ta c k in g
1 0
0
0 0
1
2
3
P h o to n E n e rg y (e V )
4
5 ACS Paragon Plus 0 Environment 1
2
3
P h o to n E n e rg y (e V )
4
5
Page 25 of 25
The Journal of Physical Chemistry
8
R e f r e a c t i v e I n d e x N ( ω)
(a ) 6
Im g
4
E
E
//
A A -s ta c k in g A B -s ta c k in g A C -s ta c k in g E //
A A -s ta c k in g A B -s ta c k in g A C -s ta c k in g E ⊥
6
A A -s ta c k in g A B -s ta c k in g A C -s ta c k in g
A A -s ta c k in g A B -s ta c k in g A C -s ta c k in g
4
Im g
2
E
//
⊥
A B A -s ta c k in g A B B -s ta c k in g A B C -s ta c k in g E ⊥
A B A -s ta c k in g A B B -s ta c k in g A B C -s ta c k in g
A B A -s ta c k in g A B B -s ta c k in g A B C -s ta c k in g
0 1
2
0 .6
3
4
5
0
1
2
0 .6 E
A A - s t a c k in g A B - s t a c k in g A C - s t a c k in g
0 .4 E
E
(d )
//
(c ) R ( ω)
E
A B A -s ta c k in g A B B -s ta c k in g A B C -s ta c k in g E //
2
0
R e fle c tiv ity
R e a l
(b )
⊥
0
3
4
5
3
4
5
3
4
5
//
A B A - s t a c k in g A B B - s t a c k in g A B C - s t a c k in g
0 .4
⊥
E
⊥
A A - s t a c k in g A B - s t a c k in g A C - s t a c k in g
A B A - s t a c k in g A B B - s t a c k in g A B C - s t a c k in g
0 .2
0 .2
0 .0
0 .0 0
1
2
3
4
5
0
8
A b s o rp tio n c o e ffic ie n t / 1 0 5 (c m -1 )
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
8 R e a l
1
2
8 E
A A - s t a c k in g A B - s t a c k in g A C - s t a c k in g
6 E
//
A B A - s t a c k in g A B B - s t a c k in g A B C - s t a c k in g
6 E
⊥
A A - s t a c k in g A B - s t a c k in g A C - s t a c k in g 4
E
(f)
//
(e )
⊥
A B A - s t a c k in g A B B - s t a c k in g A B C - s t a c k in g 4
2
2
0
0 0
1
2
3
P h o to n E n e rg y (e V )
4
5
0
ACS Paragon Plus Environment
1
2
P h o to n E n e rg y (e V )