
hydrogen or halogen passivation [20,21] Doping a GeC
monolayer could be considering as another technique to
modify the band gap. For instance, F and C dopant atoms
in a GeC monolayer disrupt the planar structure and a
surface-functionalized GeC monolayer with low-buckling
results [22]. With this type of doping the band gap is seen
to vary from 2.8 eV to 3.2 eV in calculations using HSE.
In this work, we perform DFT calculations based on
the Kohn-Sham formalism implemented in the Quantum
espresso software package [23,24]. In the calculations, we
tune the buckling parameter to study the electronic, the
thermal, and the optical properties of a GeC monolayer.
The results show that the buckling effects can be consid-
ered as an alternative way for controlling it’s physical prop-
erties, such as the band gap, the thermal conductivity and
the heat capacity.
The structure of the paper is as follows: Sec. 2includes
details of the computational methods, and Sec. 3demon-
strates the calculated electrical, the thermal, and the op-
tical properties for a GeC monolayer with different degree
of buckling. The last section, Sec. 4, is the conclusion.
2. Methodology
A 2×2 supercell of a GeC monolayer with equal number
of Ge and C atoms is considered. The GeC structure is
fully relaxed with high values of cutoffs for the plane-waves
kinetic energy and the charge densities fixed at 1088.5 eV,
and 1.088 ×104eV, respectively [25]. In the relaxation
process, the forces on the atoms are less than 10−5eV/˚
A,
where a dense Monkhorst-Pack grid with 18×18×1 is used.
The distance between GeC monolayers is assumed to be
20 ˚
A in the z-direction, which is long enough to cancel out
interlayer interactions. The generalized gradient approxi-
mation (GGA) is used with the Perdew-Burke-Ernzerhof
(PBE) functionals approximating the exchange and the
correlation terms implemented in QE software [26]. In
the calculations of the band structure and the density of
states (DOS), Self-Consistent Field (SCF) and non-self-
consistent field (NSCF) calculations are performed, respec-
tively. In these calculations, we use a Monkhorst-Pack grid
of 18 ×18 ×1 for the SCF and 100 ×100 ×1 for the NSCF
[27]. The optical properties of a GeC monolayer are ob-
tained using QE with the optical broadening of 0.1 eV.
An ab initio molecular dynamics, AIMD, calculations
are utilized to check the thermodynamic stability. The
calculations, done in the NVT ensemble, are performed
for 10 ps with a time step of 1.0 fs using the heat bath
approach described by Nos´e-Hoover [28].
The optical characteristics of the GeC monolayer can
be calculated using a large number of empty bands which
is taken into account to evaluate the dielectric properties,
ε(ω) = ε1(ω) + iε2(ω), where ε1and ε2are the real and
the imaginary parts of the dielectric function. In the long
wavelength limit q→0, ε2(ω) is given in Refs. [29,30]
ε2(ω) = 2e2π
ωε0X
K,c,v |hΨv
K|~u·~r |Ψc
Ki|2δ(Ec
K−Ev
K−ω).(1)
Herein, ωindicates the frequency of the electromagnetic
waves, ε0refers to the free space permittivity. The labels
vand cindicate the valence and conduction bands, re-
spectively, and ~u and ~r demonstrate the polarization and
the position vectors of the electromagnetic field, respec-
tively. The real part and the imaginary part of the com-
plex dielectric functions are connected to each other by the
Kramers-Kronig relation [31,32]. Both ε1(ω) and ε2(ω)
are obtained from the QE package. Once the dielectric
functions are obtained, the real part of the refractive in-
dex is calculated as [33]
n(ω) = 1
√2 hε2
1(ω) + ε2
2(ω)i
1
2+ε1(ω)!
1
2
.(2)
The optical conductivity is then computed from
σoptical =−i ω
4πhε(ω)−1i.(3)
3. Results
In this section, we show the obtained results for the
electronic, the thermal and the optical properties of a GeC
monolayer with different values for the planar buckling
parameter, ∆. In addition, for the sake of comparison, we
recalculate the physical properties of a flat GeC monolayer,
∆=0.0, and use them as reference points to compare to.
3.1. Electronic states
In a flat, or planar, GeC monolayer (∆ = 0.0), all the
Ge and the C atoms are located in the same xy-plane as is
presented in Fig. 1for both a side view (I), and a top view
(II). If a finite planar buckling is considered (∆ 6= 0.0),
the Ge atoms are located in the same plane and all the C
atoms are situated in another plane. The planar buckling
indicates the vertical distance, ∆, between the Ge and
the C planes. There are several proposed techniques to
control the planar buckling in monolayers experimentally.
One of them is using an applied external electric field on
the monolayer [34].
The degree of the hybridization of the s- and the p-
orbitals can be found using a simple equation, cos(θ) =
s/(s−1) = (p−1)/p, where θis the angle between the
equivalent orbitals, sand p[35]. The considered values
of planar buckling affect the orbital hybridization. Most
flat monolayers have an sp2hybridization, and the orbital
hybridization is approaching an sp3configuration when the
planar buckling is increased [36,37]. A flat GeC monolayer
(∆ = 0.0) has an sp2hybridization as is presented in Fig.
2(a), where the degree of the orbital hybridization is 2
indicating an sp2-hybridization at ∆ = 0.0, and the degree
2