
2
FIG. 2. a) DFT band structure calculations for C3N along
the Γ −K−M−Γ line in the BZ. b) orientation of the BZ
and the high symmetry points Γ, K,M.
with a short recap of the band structure obtained with
the help of the density functional theory calculations. In
Sec. III, effective k·pHamiltonians at Γand Mpoints are
obtained, using symmetry groups and perturbation the-
ory. Certain optical properties of this material are dis-
cussed in Sec. IV. In Sec. V the spectra of Landau levels
for this material are calculated at the Γ and Mpoints.
Finally, our main results are summarized in Sec. VI.
II. Band structure calculations
The band structure of monolayer C3N has been calcu-
lated before at the DFT level of theory [25, 31, 32] and
also using the GW approach[1, 2, 33]. The main effect of
the GW approach is to enhance the band gap and this
does not affect our main conclusions below. To be self-
contained, we repeat the band structure calculations at
the DFT level. The schematics of the crystal lattice of
single-layer C3N is shown in the Fig. 1. The lattice of
C3N possess P6/mmm space group with a planar hexag-
onal lattice and the unit cell contains six carbon and
two nitrogen atoms. We used the Wien2K package[34]
to perform first-principles calculations based on density
functional theory (DFT). For the exchange-correlation
potential we used the generalized gradient approxima-
tion (GGA) [35]. The optimized input parameters such
as RKmax, lmax, and k-point were selected to be 8.5, 10,
and 14 ×14 ×3, respectively. The convergence accuracy
of self-consistent calculations for the electron charge up
to 0.0001 was chosen and the forces acting on the atoms
were optimized to 0.1dyn/a.u. The optimized lattice con-
stant is a0= 4.86˚
A, in good agreement with previous
studies [31, 36].
The calculated band structure is shown in Fig. 2. The
conduction band minimum is located at the Γ point,
while the valence band maximum can be found at the M
of the BZ. Thus, at the DFT level C3N is an indirect band
gap semiconductor with a band gap of Ebg = 0.48 eV
which is in good agreement with previous works [21, 37].
We have checked that the magnitude of the spin-orbit
coupling is small at the band-edge points of interest and
therefore in the following we will neglect it. The main
effect of spin-orbit coupling is to lift degeneracies at cer-
tain high-symmetry points and lines, e.g., the four-fold
degeneracy of the conduction band at the Γ point would
be split into two, two-fold degenerate bands.
III. Effective k ·p Hamiltonians
We now introduce the k·pfor the Γ point, where the
band edge of the CB is located, and for the Mpoint,
where the band edge of the VB can be found.
A. Γpoint
The pertinent point group at the Γ point of the BZ is
D6h. We obtained the corresponding irreducible repre-
sentations of the nine bands around the Fermi level at
the Γ point with the help of the Wien2k package. Us-
ing this information one can then set up a nine bands
k·pmodel along the lines of Ref. [30], see Appendix A
for details. Here we only mention that there is no k·p
matrix element between the VB and the degenerate CB,
CB+1 which means that direct optical transitions are not
allowed between these two bands. Since it is usually dif-
ficult to work with a nine-band Hamiltonian, we derived
an effective low-energy Hamiltonian which describes the
two (degenerate) conduction bands and the valence band.
Using the L¨owdin partitioning technique [38, 39] we find
that
HΓ
eff =HΓ
0+HΓ
k·p,(1a)
HΓ
0=
εvb 0 0
0εcb 0
0 0 εcb+1
(1b)
HΓ
k·p=
α1q20 0
0 (α2+α3)q2−α3(q+)2
0−α3(q−)2(α2+α3)q2
.(1c)
Here εcb =εcb+1 = 0.386 eV and εvb =−1.50 eV are
band edge energies of the degenerate CB minimum and
VB maximum. The wavenumbers qx,qyare measured
from the Γ point, q±=qx±iqyand q2=q2
x+q2
yand in
α2we took into account the free electron term[29].
Note, that there are no linear-in-qmatrix elements be-
tween the VB and the degenerate CB, CB+1 bands. In
higher order of qthese bands do couple, but this is ne-
glected in the minimal model given in Eq. (1c). The mini-
mal model given in Eq. (1) already captures an important
property of the degenerate CB and CB+1 bands from the
DFT calculations, which is that their effective masses