2
studied by the DFT calculations. All Janus 2H-GdXY mono-
layers show FM semiconductor characters with high Tcbe-
yond 260 K. Interestingly, only the 2H-GdBrCl monolayer
possess the PMA behavior, while the 2H-GdICl and 2H-
GdIBr monolayers exhibit the in-plane magnetic anisotropy
(IMA) character. The competition between Gd atom-p/dor-
bitals and halogen atom-porbital can result in a transition of
easy axis direction from the [001] to [100] plane for the Janus
2H-GdXY monolayers. A spontaneous and robust valley po-
larization is observed, which can be effectively adjusted by
the external biaxial strains.
II. COMPUTATIONAL DETAILS
All density functional theory (DFT) calculations are per-
formed using the projected augmented wave (PAW) [24,25]
approach as implemented in the Vienna ab-initio package
(VASP) [26,27]. Considering the exchange and correlation
functional interactions, the Perdew-Burke-Ernzerhof (PBE)
within generalized gradient approximation (GGA) is applied
[28]. The van der Waals (vdW) correction is considered for
the bulk GdXY using the Grimme (DFT-D3) method [29]. The
plane wave cut-off energy is set to 500 eV and a vacuum space
of 18 Å is applied along the z-axis [001] direction to avoid the
interactions between adjacent layers. The crystal structure of
the Janus 2H-GdXY monolayers is completely relaxed until a
force of less than 0.01 eV/Å per atom and an energy differ-
ence of less than 10−6eV between two convergence steps is
observed. The Brillouin zone is sampled using converged Γ-
centered k-meshes with a density of 144 k-points (12×12×1)
for structural relaxation and 576 k-points (24×24×1) for the
electronic calculations [30]. The electron configurations in-
cluding 5s25p64f75d16s2for Gd [31], 4s24p5for Br, 3s23p5
for Cl, and 5s25p5for I atom are considered. The SOC ef-
fect is included in the calculations to investigate electronic,
magnetic and valley-related properties of the Janus 2H-GdXY
monolayers. The rotationally invariant local spin density ap-
proximation (LSDA)+Hubbard (U) method is employed to
treat the strongly correlated corrections to the Gd 4 felec-
trons [32] and the corresponding on-site U/exchange inter-
action Jparameters is set at 9.20 eV/1.20 eV [13,33]. The
Phonon dispersion spectrum of the Janus 2H-GdXY mono-
layers are obtained by the PHONOPY code [34,35] using a
2×2×1 supercell. Ab initio molecular dynamic (AIMD) sim-
ulations [36] adopt the NVT ensemble [37] based on the Nosé-
Hothermostat [38] controlled the temperature of systems at
300 K with a total of 8.0 ps at 2.0 f s per time step. The
VASPKIT code is used to process some of the VASP data [39].
The Tcof the Janus 2H-GdXY monolayers are estimated by
using the Monte Carlo simulation package MCSOLVER [40]
based on the Wolff algorithm. The Berry curvature and opti-
cal properties are calculated based on Fukui’s method [41] by
VASPBERRY code which is developed by Prof. Kim [42].
III. RESULTS AND DISCUSSION
Figure 1(a) displays the crystal structure of the Janus 2H-
GdXY monolayers. Clearly, the loss of reflection symme-
try of Gd atom reduces the symmetry of the systems. The
optimized lattice constants for the 2H-GdBrCl, 2H-GdICl,
and 2H-GdIBr monolayers are 3.835, 3.96, and 4.019 Å, re-
spectively. The FM and anti-ferromagnetic (AFM) ordering
of crystal structures for the Janus 2H-GdXY monolayers is
shown in Fig. 1(b) and 1(c). The ferromagnetic stability ener-
gies (∆E=EAFM−EFM) between the FM and AFM ordering
for the 2H-GdBrCl, 2H-GdICl, and 2H-GdIBr monolayers
are 165.2, 154.9, and 152.2 meV, respectively, strongly sug-
gesting the existence of FM coupling among three systems.
Figure S1 of the Supplemental Material (SM) [43] shows the
spin density images of the Janus 2H-GdXY monolayers. One
can see that the magnetic moments are mainly contributed by
the Gd and Cl (Br) atoms for the Janus 2H-GdICl (2H-GdIBr)
monolayers and by all Gd, Cl and Br atoms for the Janus 2H-
GdBrCl monolayer, suggesting that the valence electrons tend
to gather around the Cl (Br) atom with stronger electroneg-
ativity. In order to prove the stability of Janus 2H-GdXY
monolayers, the calculations of phonon dispersion and AIMD
simulation are carried out. As shown in Fig. 1(d)-1(f), the
phonon dispersions exhibit the positive value in the whole
Brillouin zone, strongly suggesting the dynamical stability for
the Janus 2H-GdXY monolayers. In addition, with time evo-
lution, the small fluctuations (about ±1 eV) of free energy, the
total magnetic moment is kept at about 128.0 µBand the orig-
inal configuration does not show large distortion have implied
a good thermal stability of the Janus 2H-GdXY monolayers
[Fig. S2(a)-S2(c) of the SM [43]]. To estimate the possibility
of mechanical exfoliation, the exfoliation energy of the Janus
2H-GdXY monolayers is calculated in four-layer slab models
with AB-stacking [44] [Fig. S3(a) of the SM [43]] from their
layered bulk crystals. Considering the small value of sepa-
ration distance, the exfoliation process is performed with the
fixed atomic positions. As shown in Fig. S3(b)-S3(d) of the
SM [43], the increase of separation distance (d−d0) leads
to an obvious increase in the energy differences ∆E=Ed−Ed0,
which converge to 0.232, 0.243, and 0.239 J/m2, respectively,
for the Janus 2H-GdBrCl, 2H-GdICl, and 2H-GdIBr mono-
layers, respectively. These cleavage energies are remarkably
lower than the measured value (0.36 J/m2) for the graphite
[45], which can be further confirmed by the variation of cleav-
age strength (the first-order derivative of cleavage energy).
This means that the Janus 2H-GdXY monolayers are easily
to exfoliate experimentally.
The band structures of the Janus 2H-GdXY monolayers
without considering SOC are shown in in Fig. 2(a)-2(c). The
energies at the K and K0valleys are equal, suggesting that
the valley splitting does not appear without considering SOC.
The typical bipolar magnetic semiconductor character [BMS:
the valence and conduction bands possess opposite spin po-
larization when approaching the Fermi level (EF)] [46] with
an indirect band gap can be observed, where the valence band
maximum (VBM) located at the Γpoint and the conductor
band minimum (CBM) located at the M point have the oppo-