
2
(a)
(b) (c)
FIG. 2: The Ti2C 2D-threelayer system: (a) view from
above; and side views of (b) ferromagnetic and (c)
antiferromagnetic configurations. Red arrows indicate
the internal magnetic moment of Ti atoms.
works well in determining miRNAs and diagnosis of can-
cer biomarkers.[12] Computational studies of MXenes
have gone further than experimental ones. Quantum-
chemical methods make it possible to investigate com-
pounds that are still difficult to obtain experimentally.
Theoretical studies of organometallic molecules on
MXene surfaces were not found in the literature, but
there are several theoretical studies devoted to the ad-
sorption of atoms on the surface of Ti2C and Ti3C22D
MXene layers. It was found [13] that the adsorption en-
ergies of 3d, 4d, and 5d transition metal atoms on Ti3C2
are in the range of -7.98 to -1.05 eV. 3d-transition met-
als on M2C layers (M = Ti, V, Cr, Zr, Nb, Mo, Hf,
Ta, and W) were studied [14] as single-atom catalysts to
find an alternative to Pt-based catalysts. Also, there are
studies about metal atoms on functionalised layers for
single-atom catalysts [15, 16] and Li-ion storage.[17] It
should be noted that these studies did not take into ac-
count the spin polarisation of surfaces. The results of the
studies are presented in Table I. The calculation param-
eters given in parentheses will be explained in the next
chapter.
Determination of the magnetic moment of each atom
in a system is possible using spin-polarized scanning tun-
nelling microscopy.[18–20] In these works, the Fe atom
was studied in the interaction with InSb(110) surface. In
this study, the iron atom falls into the surface due to the
large lattice constant. In the case of the iron atom on top
of the Cu(001) surface, it was shown [21] that the elec-
tronic and magnetic properties of adatoms are strongly
affected by the tip-surface distance.
In the purely theoretical study of iron chains on
Cu(001) and Cu(111) surfaces [22] a single iron atom acts
as a donor of spin momentum. It retains most of the spin
momentum when it is in the so-called ”ontop” position
and loses it when the atom is in or inside the surface.
The Bader charge transfer analysis for 3d TM atoms on
graphene and graphene/Ni(111) surfaces [23] indicates
that for all elements in the series electrons are trans-
ferred from the adatom to the graphene layer, leaving
the net charge on the adatom positive. Magnetic atoms
on the surface were also studied. Cobalt and iron atoms
were placed onto the Pt(111) surface [24] and also on
Pt(111) and Ir(111) surfaces.[25] There was shown that
adatoms induce polarisation on nearby surface atoms.
The magnetic anisotropy parameters also were calculated
and they are in agreement with the experimentally deter-
mined ones, where inelastic tunnelling spectroscopy was
implemented.[26]
The ability of an atom on a surface to have two sta-
ble states was studied in several articles. TM atoms on
the graphene/Ni(111) surface can be ferromagnetic and
antiferromagnetic toward the surface magnetisation.[23]
It was found that for Ti, V and Cr antiferromagnetic
alignment is preferable while for Mn, Fe and Co it is fer-
romagnetic (exchange energy for Fe, in this case, is 10
meV). By utilizing a combination of scanning tunnelling
spectroscopy and DFT methods, it was shown that a Co
atom on semiconducting black phosphorus [27] has two
states: low-spin and high-spin. It was shown experimen-
tally that the state of the atom can be switched elec-
trically. A holmium atom on the MgO surface exhibits
bistability property with up and down spin states.[28] It
was shown that it is possible to read the states using a
tunnel magnetoresistance and induce particular magnetic
state of the iron atom (just to ”write” it) with current
pulses using a scanning tunnelling microscope.
In this paper, we present the study of the FePc
molecule on the Ti2C MXene layer. Here, we mostly
focus on the magnetic properties of the formed system.
Different magnetic configurations of the FePc/Ti2C hy-
brid system are modelled and compared. Additionally, a
single Fe atom and the H2Pc molecule on the top of Ti2C
are considered. It allows one to understand the role of
the phthalocyanine ligand in the complex.
II. COMPUTATIONAL DETAILS
The DFT calculations for the studied hybrid
FePc/Ti2C system with periodic boundary conditions
were performed employing the Quantum Espresso
6.5 numerical package [29] with the generalised gra-
dient approximation (GGA) [30] realised through
the Perdew-Burke-Ernzerhof (PBE) exchange-correlation
functional.[31] The van der Waals interaction between
the molecule and the Ti2C layer was accounted for within
the Grimme DFT-D3 ad-hoc scheme.[32] Rappe-Rabe-
Kaxiras-Joannopoulos (RRKJ) ultrasoft pseudopoten-
tials from pslibrary [33] were implemented. To treat the
strong on-site Coulomb interaction of TM d-electrons,
we used the DFT+U approach within the Hubbard
model.[34] The U parameter value for the Fe atom (U