## Confinement of quasi-atomic structures in Ti N and Ti N 2 3 2 MXene Electrides

Chandra M. Adhikari, 1 Dinesh Thapa, 2 Talon D. Alexander, 1 Christopher K. Addaman, 1 Shubo Han, 1 Bishnu P. Bastakoti, 3 Daniel E. Autrey, 1 Svetlana Kilina, 2 Binod K. Rai, 4 and Bhoj R. Gautam 1

1 Department of Chemistry, Physics and Materials Science, Fayetteville State University, Fayetteville, NC 28301, USA 2 Department of Chemistry and Biochemistry, North Dakota State University, Fargo, ND 58108, USA 3 Department of Chemistry, North Carolina A &amp; T State University, Greensboro, NC, 27411, USA 4 Savannah River National Laboratory, Aiken, SC 29808, USA

Metal carbides, nitrides, or carbonitrides of early transition metals, better known as MXenes, possess notable structural, electrical, and magnetic properties. Analyzing electronic structures by calculating structural stability, band structure, density of states, Bader charge transfer, and work functions utilizing first principle calculations, we revealed that titanium nitride Mxenes, namely Ti N and Ti N , have excess anionic electrons in their pseudo-atomic structure inside the crystal 2 3 2 lattice, making them MXene electrides. Bulk Ti N 3 2 has competing antiferromagnetic (AFM) and ferromagnetic(FM) configurations with slightly more stable AFM configurations, while the Ti N 2 MXene is nonmagnetic. Although Ti N 3 2 favors AFM configurations with hexagonal crystal systems having 6 /mmm point group symmetry, Ti N 3 2 does not support altermagnetism. The monolayer of the Ti N 3 2 MXene is a ferromagnetic electride. These unique properties of having non-nuclear interstitial anionic electrons in the electronic structure of titanium nitride MXene have not yet been reported in the literature. Density functional theory calculations show TiN is neither an electride, MXene, or magnetic.

## I. INTRODUCTION

Newly discovered two-dimensional layered inorganic early transition metal carbides, nitrides, and carbonitrides compounds with the general formula of M n +1 X , n where M is a transition metal found in groups IIIB-VIB and periods 4-6 and X being carbon (C) and/or nitrogen (N), have received significant attention in the materials science community because of their very useful properties required in various applications within the electronic, magnetic and optical industries [1]. The name MXenes is adopted for M n +1 X n compounds to emphasize their composition and similarities with the well-known twodimensional (2D) material graphene [2].

from electronic wastes [9]. Tunable bandgap with proper selection of surface termination group, superior surface properties, and high mechanical strength make MXenes excellent sensing materials for strain/stress, gas, electrochemical, optical, and humidity sensors [10]. MXenes can also be used in catalysts for hydrogen evolution reactions (HER), oxygen evolution reactions (OER), and bifunctional with HER at the cathode and OER at the anode due to their metallic and intrinsic hydrophilic nature [11-13].

Since they are layered 2D structures, MXenes have a larger surface-to-volume ratio, similar to graphene. Mxenes have high electrical conductivity and a larger oxidation-reduction reaction active surface, which improves electrode-electrolyte ion exchange and heat dissipation, making them highly efficient materials in the industry of rechargeable batteries and supercapacitors [3, 4]. MXenes have a tunable dielectric response, surface-enhanced Raman scattering, high light-to-heat conversion efficiency, and plasmonic broadband absorption, which make them a good candidate for optoplasmonics [5, 6]. MXenes have biomedical applications as antibacterial agents, biosensing, therapeutic diagnosis, implants, and bioimaging [7]. MXenes provide electromagnetic interference shielding, even at a minimal thickness, due to their high electrical conductivity and multiple internal reflections from their flakes [8]. MXeneenabled technologies work as an environmental cleaning agent, removing organic and inorganic pollutants present in the environment. They also help to recycle and recover precious rare earth elements, eliminating radionuclides and heavy ions such as Cr 6+ , Cd 2+ , and Pb 2+

To understand the MXenes' usability in building magnetic sensors, magnetic storage devices, and in spintronics, one must investigate the magnetism of various MXenes. The magnetic properties prediction of MXenes is quite complicated as they depend not only on the parent early transition metals (M) and carbon and/or nitrogen (X) of MXenes but also on many other things such as route of synthesis, strain level, purity of the sample, defect, doping, crystal symmetry, oxidation state of M, and the presence of surface termination groups [14-17].

Interestingly, since the discovery of MXenes in 2011 at Drexel University [2], most studies on MXenes are concentrated more towards carbides than nitrides. With the currently available technology, carbide MXenes are comparatively easier to synthesize and characterize than nitride MXenes. As far as technological applications are concerned, a study showed that nitride-based MXenes have a comparatively larger active surface, a stronger preference to adhere surface termination groups, and higher electrical conductivity than carbide-based Mxenes [18]. To further deepen the knowledge and unearth the fundamental understanding of nitride-based MXenes, we devote this work to titanium nitride MXenes, specifically Ti N and Ti N . 2 3 2 A study on TiN is also presented for comparison purposes. Here and onwards, TiN, Ti N, and 2 Ti N 3 2 will be collectively referred to as Ti-N compounds.

To better analyze the electronic structure and magnetism of MXenes, a discussion on transition metal to carbon or nitrogen bonding is necessary. Transition metals have comparatively lower electronegativity than carbon and nitrogen. For example, the electronegativity of Titanium (Ti) by the Pauling scale is 1.54, while the same for carbon (C) and nitrogen (N) are 2.55 and 3.04, respectively [19, 20]. The absolute electronegativity differences of (Ti, C) and (Ti, N) pairs are less than 1.7 but more than 0.4. From a simple rule of thumb of bond classification [21], both carbide and nitride Mxenes of titanium have neither a pure ionic bonding nor a covalent bonding but prefer to have polar covalent bonds, while the covalent bonding percentage is higher in the Ti to C bond than the Ti to N bond. The electronic charge distribution in Ti 3 N 2 and Ti N Mxenes is com2 plex and does not follow the bond classification rules as they have unique electronic structures. Our first principle DFTcalculations demonstrate that both Ti N 3 2 and Ti N 2 have excess electrons in their pseudo-atomic structure inside their crystal lattice, known as the quasi-Farbe color centers, making them the first discovered titanium-based MXene electrides. In contrast, no quasi-Farbe color centers are present in TiN crystals, dismissing the possibility of TiN being electride.

An introduction to electride materials, a relatively under-explored material class, is in order. Electrides have fully trapped isolated electrons occupying interstitial space between atoms with an anisotropic electronic distribution, which are (i) neither fully transferred as in ionic, (ii) nor shared equally between elements with almost the same electronegativity as in covalent, (iii) nor shared unequally between element with slightly different electronegativity as in polar covalent compounds, but instead, they have pseudo-atomic structures also called empty spheres of interstitial ionic electron (IAE) sites. As these empty spheres have no nucleus, the trapped anionic electrons in IAEs can easily diffuse, resulting in low work functions and high electrical conductivity [22], which makes them a potential candidate for field emission display in cold cathodes, thermionic emitters, battery anode in high-performance energy storage systems, supercapacitors, high-speed transistors and so on [2326]. High electron density, high electron mobility, and the ability to activate nitrogen molecules make electrides excellent catalysts for organic molecule synthesis [27]. Exploiting the electrons' spin and charge distribution, magnetic electrides can be used in spintronic devices for information processing and storage [28].

Carrying out high-throughput material screening in conjunction with a geometry-based algorithm on a large set of materials from the Materials Project database [29], Zhou et al. reconfirmed already reported 12 electrides, namely, ScCl [30], YCl [30], ZrCl [31], ZrBr [32], Ca N [33, 34], Sr 2 2 N [34], Ba 2 N [34], Sc 2 C [35], Y 2 C [36], Tb C [37], Dy C [37], and Ho C [37] and uncovered 12 2 2 2 new materials as layered electrides, namely, TbCl, TbBr, LaBr , ThI , La Br , Y Cl , Gd Br , Sc Cl , Tb Br , 2 2 2 5 2 3 2 3 5 8 5 8

Er I , Sc 6 7 7 Cl 10 , and Ba LiN [38]. Utilizing high through2 put computational screening followed by the first principle calculations, Thapa et al. recently reported that electride material need not necessarily be carbide, nitride, or halide, but they can be a transition metal oxide or chalcogenides such as Zr 2 X, where X= O, Se, and Te [39].

Considering both rare earth scandium (Sc) and yttrium (Y) as early transition metal candidates for MXenes, only two carbide MXenes, Sc C and Y C, are re2 2 ported as electrides in the literature [35, 36]. A layered rare earth lanthanide gadolinium carbide, Gd C, is also 2 reported to have a two-dimensional electron gas, whose excess anionic electron counts increase with the square of the electric field intensity [40]. The electronegativity of metal cations contributes to determining the electronic structure of electrides. For example, less electronegative yttrium can not open a band gap in Y C and remains a 2 semi-metallic electride. In contrast, high electronegative scandium and aluminum in Sc C and Al C with the same 2 2 structure as Y C open up indirect band gaps. 2 The gap is larger for Al 2 C than Sc C, implying that the band gap 2 increases with an increase in electronegativivity [35]. We here make an extensive examination via the first principle DFT calculations and report a new electride material family: Ti n +1 N n MXene. The first principle DFT calculations confirm the presence of IAEs in Ti 3 N 2 and Ti N, 2 but not in TiN. A few additional benefits of the MXene electride family that stand out from other existing electrides are their thickness tunability due to being a stacked layered structure with c/a &gt; 1, higher charge density, and enhanced charge carrier mobility.

TABLE I. Lattice parameters for computationally optimized structures of TiN, Ti N, Ti N 2 3 2 FM and Ti N 3 2 AFM. The values in parenthesis are experimental values presented in Ref. [41] for TiN and Ref. [42] for Ti 2 N. To the best of our knowledge, no experimental values for lattice parameters are yet available for Ti 3 N . 2

| compounds    | a (A) ˚       | b (A) ˚       | c (A) ˚      |
|--------------|---------------|---------------|--------------|
| TiN          | 4.241 (4.244) | 4.241(4.244)  | 4.241(4.244) |
| Ti 2 N       | 4.149 (4.140) | 4.149 (4.140) | 8.873(8.805) |
| Ti 3 N 2 FM  | 3.032         | 3.032         | 14.630       |
| Ti 3 N 2 AFM | 3.032         | 3.032         | 14.629       |

## II. STRUCTURAL ANALYSIS

Ti N 3 2 forms a hexagonal crystal (p6 /mmc 3 space group) with a c/a ratio of 4.8244 and a c value of 14.6285 ˚ . A Ti N crystallizes into a tetragonal structure 2 (I4 1 /amd space group) with a c/a ratio of 2.1387 and a c value of 8.8725 ˚ . Ti A 3 N 2 and Ti N have layered struc2 tures and experimentally feasible c-values to synthesize 2D-MXenes. The titanium in Ti N 3 2 at the 2 a Wyckoff position cages itself within a TiN 6 octahedron by bonding to six equivalent nitrogen atoms, while the titanium at the 4 f Wyckoff position forms a triangular noncoplanar-

FIG. 1. Structure models of Ti 2 N (A), and Ti N 3 2 (B) showing unit cells along the rotated (100) axis to a side-view without empty spheres (a) and with empty spheres (b) and a topview along the (001) crystal axis (c). Empty spheres in two different interlayer regions in Ti 3 N 2 receive different environments. Thus, we denote them as A and B empty spheres for A and B empty sphere layers. Empty spheres in Ti N receive 2 an indifferent crystallographic environment.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00897v1-figure-78be4371-1290-4549-b52a-49fd03b8fce6.jpg)

environment (not shown in the figure). The contentious symmetry measures show that titanium in the 2 a Wyckoff position has less distortion in the coordination environment than titanium in the 4 f Wyckoff position. Nitrogens in Ti 3 N are at the 4 2 f Wyckoff position and confined at the center of NTi 6 octahedron. In Ti N, each Ti 2 resides in the 8 e Wyckoff position and is bonded to three N atoms in a plane. The Nitrogen in Ti N is at the 4 2 a Wyckoff position and bonds to six Ti, forming an NTi 6 octahedron. We have also analyzed TiN for comparison purposes, although it is not an MXene. TiN crystallizes into cubic (Fm3m space group) with vanishing continuous symmetry measures indicating no distortion in its coordination environment. In TiN, titanium resides in the 4 a Wyckoff position and nitrogen at the 4 b Wyckoff position. In TiN, each titanium and nitrogen are caged in a TiN 6 and NTi 6 octahedron. Lattice parameters for computationally optimized structures of TiN, Ti N, Ti N 2 3 2 FM, and Ti N 3 2 AFM and available relevant experimental values are presented in Table I. As per the need for comparison, one can convert the lattice parameter in a face-centered cubic (FCC) to that of corresponding bodycentered cubic (BCC) structures using ℓ FCC = √ 2 ℓ BCC , where ℓ denotes any of the a , b or c lattice parameters.

To measure the stability of any structure, one may need to analyze the formation enthalpy ( E F ) and cohesive energy ( E C ) of the system, where the former is the energy required to form a crystal from its constituent elements in their most stable bulk state while the latter represents the energy to assemble isolated constituent atoms to a crystal or the same to disassemble a crystal into its constituent isolated atoms. More implicitly, E F is the total energy of the optimized crystal structure minus the bulk energies of the crystal structures of all constituent atoms with their respective weight, and E C is the total energy of the optimized crystal structure minus the sum of the energies of all the isolated atoms. For example, for Ti n +1 N n MXenes, the enthalpy of formation E F (Ti n +1 N ) and the cohesive energy E n C (Ti n +1 N n ) read as follows:

$$E _ { F / C } ( \text{Ti} _ { n + 1 } N _ { n } ) & = E _ { T } ( \text{Ti} _ { n + 1 } N _ { n } ) \\ & - ( n + 1 ) \, E _ { B / I } ( \text{Ti} ) - n \, E _ { \text{B/I} } ( N ) \,, \quad ( 1 )$$

where E T is the total energy of specified MXene, E B (Ti), and E (Ti) represent I the energy of the bulk titanium metal and the energy of an isolated titanium metal atom. A similar description holds for nitrogen. One may want to divide the right side of the expression by 2 n +1 to obtain the corresponding energy per atom. For TiN, Eq. (1) reads

$$\text{ve} \quad E _ { F / C } ( \text{TiN} ) = E _ { T } ( \text{TiN} ) - \text{E} _ { \text{B/I} } ( \text{Ti} ) - \text{E} _ { \text{B/I} } ( \text{N} ) \,. \quad ( 2 )$$

Formation enthalpy for Ti N is 2 -1 3165 eV/atom. . The formation enthalpy for antiferromagnetic (AFM) and ferromagnetic (FM) configurations of Ti N 3 2 are very close to each other, which differs only in the fifth significant figure when expressed in the units of eV/atom. More explicitly, E F, AFM (Ti N ) = 3 2 -1 3589 eV/atom and . E F, FM (Ti N ) = 3 2 -1 3583 eV/atom, making AFM config-. uration slightly more favorable. The formation energy of TiN crystal is -1 7386 eV/atom, making it the most sta-. ble among the three titanium nitrides discussed in this article. The formation enthalpy is negative for each titanium nitride of interest in this paper, signifying that the formation of titanium nitride crystals is exothermic and energetically favorable. The cohesive energy per atom for TiN, Ti N, Ti N 2 3 2 AFM, and Ti N 3 2 FM configurations are -4 4561, . -4 9430, . -4 6218, and . -4 6212 eV respec-. tively.

To facilitate efficient and scalable production of MXenes, in addition to the formation of enthalpy and cohesive energy, one needs to investigate the energy required to exfoliate MXenes from their respective MAX phase, called exfoliation energy. Pristine MXenes are usually exfoliated from Ti n +1 AX , where the A element is usun ally aluminum or gallium. The exfoliation energy for constructing Ti n +1 N n by removing the A layer is calculated

FIG. 2. Atom resolved projected density of states of AFM Ti N 3 2 in bulk phase (a), Ti 2 N (b), and TiN (c) per their formula unit (f.u.).E A and E B represent DOS of A and B empty spheres per f.u. whereas (E=E A +E B ) is the total DOS of empty spheres in the respective crystal structure. Calculation shows a vanishing empty sphere DOS for TiN, indicating that TiN is a nonelectride.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00897v1-figure-adee41f5-36ff-4b8a-915e-f72a23217a0a.jpg)

FIG. 3. Bandstructure of AFM Ti N 3 2 (a), Ti 2 N (b), and TiN (c) along their respective high symmetry points. Fermi level (E F ) is shifted to zero energy level. Each of the titanium nitrides is gapless, indicating that they are metallic.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00897v1-figure-75920577-65ae-4e13-b199-463aabf8e659.jpg)

using the following formula [43, 44]:

$$E _ { \text{ext} } = \frac { 1 } { 2 \mathcal { A } } \left [ 2 E ( A ) + 2 E ( T i _ { n + 1 } N _ { n } ) - E ( T i _ { n + 1 } A N _ { n } ) \right ] \,, \quad \begin{pmatrix} \text{The } \, \mathbb { i } \\ \text{ing } t \mathbb { l } \\ \text{excha} \\ \text{the } P \end{pmatrix}$$

where E(A) is the total energy of an A atom in its preferred bulk configuration. For the aluminum (Al) atom, cubic Fm m 3 (space group # 225) is the preferred crystal structure with a total energy of -3.7478 eV per atom. Similarly, E(Ti n +1 N ) and E(Ti n n +1 AN n ) are total energy values for Ti n +1 X n MXene and Ti n +1 AX n MAX phase respectively. A is the area of the exfoliated layer. The calculated exfoliation energy for Ti 2 N and Ti N 3 2 are 0.0224 and 0.0769 eV/ ˚ A 2 indicating that Ti N can be 2 exfoliated more easily than Ti 3 N 2 from their respective Al-based MAX phases.

## III. COMPUTATIONAL METHODS

The first-principles calculations were performed based on Density Functional Theory (DFT) using the Vi- enna Ab initio Simulation Package (VASP) [45, 46]. The ion-electron interactions were accounted for using the projector augmented wave method [47], and exchange-correlation potentials were approximated by the Perdew-Burke-Ernzerhof functional [48]. Atomic coordinates were fully relaxed with the force and energy convergence criteria of 0.01 eV/ ˚ A and 10 -6 eV, respectively. The plane-wave cutoff energy of 520 eV was used. The Monkhorst-pack Γ-centered k -point mesh of 21 × 21 × 5 and 13 × 13 × 7 were used for Brillouin zone sampling for single unit cells for Ti 3 N and Ti N, respec2 2 tively. The same for TiN was taken as 11 × 11 × 11. A Gaussian smearing of 0.05 eV width was used. The empty spheres of 1.117 ˚ A and 1.552 ˚ A Wigner-Seitz radii ( R WS ) were taken for Ti N 3 2 and Ti N, respectively, and posi2 tioned by examining their electron localization function and matching their centers with localized IAEs. R WS for Ti and N are 1.217 ˚ A and 0.741 ˚ A respectively in both TiN MXenes. The initial crystal structures were adopted from the Materials Project [29]. Crystallographic presentations were made using the Vesta program [49].

The dynamics of collective excitations in a crystal

called phonon are due to the vibratory movements of atoms around their atomic equilibrium positions in crystalline lattices, which help understand the local stability of the crystal. To better understand stability of various configurations of the crystal, we performed the density functional perturbation theory (DFPT) calculations [50, 51] using the Phonopy program [52, 53] and analyzed the phonon dispersion and phonon density of states in the chosen titanium nitride compounds. In DFPT, the system's ground state is lightly perturbed by slightly displacing atoms, and the system's responses to such perturbation are analyzed. The DFPT in the Phonopy program constructs a Hessian matrix or a 'dynamical matrix D ij ( ⃗ q ),' from the second-order derivative of the system's total energy ( E ) with respect to atomic displacements from their equilibrium positions. Mathematically,

$$D _ { i j } ( \vec { q } ) = \frac { \partial ^ { 2 } E } { \partial x _ { i } \partial x _ { j } } \quad \quad ( 4 ) \quad \begin{array} { c c } \text{FI} \\ \text{Til} \\ \ell _ { \text{Fe} } \end{array}$$

where ⃗ q is the momentum transfer of phonon in reciprocal spaces, and x i and x j are atomic displacements of i th and j th atoms from their respective equilibrium positions. Diagonalization of D ij ( ⃗ q ) matrix provides all the phonon modes of the crystal lattice as eigenvectors and phonon frequencies corresponding to all phonon modes of the crystal lattice as eigenvalues. We used 2 × 2 × 2 supercell consisting of 96 atoms for Ti 2 N and 4 × 4 × 4 supercell containing 128 for TiN to study phonon dynamics. The used supercells achieve phonon accuracy with vanishing imaginary vibration modes in both Ti N and 2 TiN. The imaginary frequency does not vanish in Ti N 3 2 crystals for as big as 4 × 4 × 2 supercell of 320 atoms. The force and energy convergence criteria of 0.01 eV/ ˚ A and 10 -6 eV, respectively, were used for ionic relaxation during the phonon frequency calculations. In the case of Ti N 3 2 monolayer, however, the imaginary vibration mode hardly disappears only when an as big as 6 × 6 × 1 supercell is constructed and used.

## IV. ELECRONIC STRUCTURE AND PROPERTIES

Studying the distribution of electrons and analyzing their behaviors in a material is crucial to understanding its electronic structure. To explore the relationship of electronic structure with their electronic properties for Ti-N crystalline systems, we theoretically calculate electronic band structure, density of states (DOS), electron localization function (ELF), work function, and charge transfer.

## A. Band structure and Density of states

Electrons' energy as a function of their wave vectors called band structure, the number of electronic states

FIG. 4. ELF maps for a body-centered cubic unit cell of TiN crystal keeping N at the center of the body. Note that ℓ FCC = √ 2 ℓ BCC . Figures show the TiN unit cell (a), 3 D isosurface plot for TiN with an iso-surface value of 0.6 (b), and 2 D contour plot of ELF along (001) and (110) planes (c, d). The color scale bar shows the ELF value measured in color.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00897v1-figure-362926b0-820b-4c69-9b13-db8771334688.jpg)

FIG. 5. Rotated (100) crystal direction view of 3 D isosurface plot of ELF for Ti 2 Nwith an isosurface value of 0.6 (a) and 2 D contour plot of ELF along (100), (110) and (001) view for Ti 2 N with the same isosurface value. The color scale bar shows the ELF value measured in color. Red-colored centers in the interstitial region confirm the presence of IAEs in Ti 2 N crystal.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00897v1-figure-5f8b7558-4ca3-4e70-a216-34ceec62ecd5.jpg)

per unit electron volt of energy interval as a function of energy called DOS, and the IAEs' contributions to the band structures and DOS at the Fermi level (E F ) are calculated and presented in Figs. 2 and 3 for all three Ti-N crystalline systems. As expected, the band structures

overlap for two different spins in the AFM configuration of the bulk Ti 3 N . Each of Ti N 2 3 2 and Ti N has four IAEs 2 in their crystal structures. Each IAE in Ti N receives an 2 identical crystal lattice environment. However, IAEs in Ti N 3 2 present in different interlayer spaces and have nonidentical contributions to the electronic structure. The layer of IAEs in these two interlayer spaces are denoted as A-type (strongly localized) and B-type (less localized) IAEs as shown in Fig 1(B)(b). Near and at the Fermi level of Ti 3 N , both 2 E A and E B empty spheres have significant contributions to the density of states. The DOS of empty spheres near and at the Fermi level is even larger than the nitrogen for both Ti N 3 2 and Ti N. The major 2 contributions to the DOS near and above the Fermi level come from Ti in each of the titanium nitrides discussed in this paper. Orbital resolved partial density of states shows that the significant contributions to the total DOS come from the 3 d orbital of Ti, IAEs, and 2 p orbital of N (not shown in the figure here). At the valence band (energy band below Fermi level), the hybridization of Ti -3 d with empty spheres, Ti -3 d with N -2 p , and N -2 p with empty spheres indicate the complex chemical bonding feature in Ti 3 N 2 and Ti N. Since empty spheres do 2 not contain nuclei and electrons like real atoms, to keep it simple, we did not resolve the empty sphere's basis functions into their orbitals.

(A) TisN? AFM spin up

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00897v1-figure-4657d1b0-58cb-4824-9700-7aa254698f02.jpg)

(B) TizN2 AFM spin down

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00897v1-figure-7ed69aab-a0e6-4d8b-9f0e-f10982ddebe3.jpg)

FIG. 6. ELF plots for Ti N 3 2 crystal in its AFM configurations for up spin (A) and down spin (B) with an iso-surface value of 0.6. Each (A) and (B) part of the figure has a 3 D ELF map showing in a rotated (110) view (a), a side view of 2 D contour map of ELFs along (110) plane (b), the same along (001) plane slicing through the plane containing empty spheres E B (c) and the same for a plane containing the empty sphere E A (d). The color scale bar shows the ELF value measured in color.

FIG. 7. ELF plots for Ti 3 N crystal in its FM configurations 2 for up spin (A) and down spin (B) with an iso-surface value of 0.6. Each (A) and (B) part of the figure has a 3 D ELF map showing in a rotated (110) view (a), a side view of 2 D contour map of ELFs along (110) plane (b), the same along (001) plane slicing through the plane containing empty spheres E B (c) and the same for a plane containing the empty sphere E A (d). The color scale bar shows the ELF value measured in color.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00897v1-figure-601bb397-5266-4a04-9da6-9d9d4ef62a37.jpg)

## B. Electron localization function

To visualize and quantify electron localization in Ti-N crystals discussed in this paper, one can make use of the dimensionless ELF given by [39, 54]

$$\text{ELF} ( \tilde { r } ) = \left [ 1 + \left \{ \frac { \text{D} _ { \sigma } ( \tilde { r } ) } { \text{D} _ { \sigma } ^ { 0 } ( \tilde { r } ) } \right \} ^ { 2 } \right ] ^ { - 1 }, \text{ \quad \ \ } ( 5 )$$

where D σ = τ σ - ∇ ( ρ σ ( ⃗ r )) 2 / (4 ρ σ ( ⃗ r ) is the electron localization measure determining the electron pair density curvature of an electron having identical spins at position ⃗ r and D 0 σ ( ⃗ r ) is the same for electron gas, which reads

$$D _ { \sigma } ^ { 0 } ( \vec { r } ) = \frac { 3 } { 5 } \left ( 6 \pi ^ { 2 } \right ) ^ { 2 / 3 } \rho _ { \sigma } ^ { 5 / 3 } ( \vec { r } ) \,,$$

where ρ σ ( ⃗ r ) is the spin electron density at ⃗ r . ELF ( ) ⃗ r at any point ⃗ r satisfies 0 ≤ ELF( ) ⃗ r ≤ 1 with ELF=1 indicating perfect localization. A system with a uniform charge distribution has an ELF value of 0.5. ELF is less than 0.5 if the electrons are delocalized, whereas the system with localized electrons has ELF ¿ 0.5.

ELF maps for Ti-N crystalline systems and 2 D -contur maps along a few crystallographic directions are shown in Figs. 4, 5, and 6. To gain an accurate and deep understanding of the spatial distribution of the IAEs, a color

FIG. 8. Work function (Φ S ) along (001) plane in TiN (a), Ti 2 N (b) and AFM Ti N 3 2 (c) and the same along (110) plane in TiN (d), Ti N (e) and AFM Ti N 2 3 2 (f) as a function of distance. The oscillating solid blue line for each titanium nitride is the planner average energy over a plane perpendicular to the crystal axis, whereas the solid red line is the macroscopic average energy. The dashed gray line represents the average electrostatic potential, and the dashed green line is the Fermi energy of the slab E F slab , . The same energy range is used in the vertical axis of all work function plots for better comparison.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00897v1-figure-72d3b980-c3b9-476b-b4f1-12a00f674861.jpg)

scale bar in the range of 0.0 to 0.7 was used with a continuous color gradient from blue (minimum ELF) to red (Maximum ELF). In TiN, electrons are strongly localized around the Ti atom, whereas they form a uniformly dense electron gas around N. The 3D iso-surface plot of ELF and 2 D contour plots in different directions confirm the absence of any anionic electrons in the TiN crystal. We witnessed four IAEs in a unit cell of Ti 2 N, filled with maximally localized anionic electrons with an identical shape of iso-surface with the same color depth (see Fig. 5). Further, the (100) side view of ELF in Ti 2 N indicates the one-dimensional distribution of IAEs within

the lattice crystal where the localized IAEs are interconnected with delocalized electron channel as observed in transition metal-rich oxide and chalcogenide electride materials [39]. Ti N 3 2 hosts a layer of IAEs between every two adjacent pairs of MXene layers. The IAE layers in different spatial locations receive a non-identical environment, making them distinct from each other. ELF of a Ti N 3 2 crystal depends not only on its magnetic configuration but also on its spin alignment. Anionic electrons in the Ti 3 N AFM configuration are comparatively more 2 localized at E A empty spheres but less localized at E B s as shown in Fig. 6. The 3 D iso-surfaces of the localized anionic electrons make almost a dumbbell shape with slightly different weights when viewed from the (110) plane; 2 D mapped ELF of the anionic electrons form a kite-like shape. The (001) plane (top view) ensures the two-dimensional distribution of IAEs in Ti N . 3 2 E A and E B exchange the localization/delocalization feature of anionic electrons as the spin is flipped. Interestingly, in the FM configurations of Ti N 3 2 the anionic electrons localize in both E A and E B spheres if electrons in the host atom species are at down spin, but these electrons delocalize if electrons in the host atom species are at up spin.

## C. Work function

The work function of a solid is the energy needed to remove an electron from the solid's Fermi level to a point where the electron becomes free from its influence, such that the removed electron has no kinetic energy. The work function provides deep insights into the solid's electronic, structural, physical, and optical properties, thereby guiding the solid's manipulation, design, and optimization for various electronic, optoelectronic, nanotechnological, and photovoltaic applications. The DFT calculations for determining work function involve evaluating vacuum energy, a constant electrostatic potential value along the surface's perpendicular direction, and the Fermi level of the material and finding the difference between them. Based on the calculation route, two work function calculation methods are common. The first method evaluates a slab's Fermi energy E F slab , and electrostatic potential in the vacuum region V vac using self-consistent calculations, and the work function is calculated as

$$\Phi _ { S } = V _ { \text{vac} } - E _ { F, \text{slab} } \,,$$

where subscript S in Φ indicates that it is the first method accounting for the slab's Fermi energy. The second method uses a two-step calculation: (i) computing the electrostatic potential in the vacuum region V vac and in the slab's interior region V Int slab , using a supercell model and (ii) evaluating the Fermi level of the bulk E F bulk , from a separate calculation. Finally, Φ B from the second method reads

$$\Phi _ { B } = V _ { \text{vac} } - E _ { F, \text{bulk} } ^ { * } \,,$$

where E ∗ F bulk , = E F bulk , +V Int slab , is the effective value of bulk's Fermi level. For further detail, we refer to Ref. [39].

TABLE II. Calculated work functions for Ti-N crystals along (001) and (110) planes. TiN-Ti and TiN-N represent the specified TiN Crystal face with the Ti and N atoms termination, respectively.

| system   | plane   |   E ∗ F , bulk (eV) |   E F , slab (eV) |   Φ S (eV) |   Φ B (eV) |
|----------|---------|---------------------|-------------------|------------|------------|
| TiN-Ti   | (001)   |              4.0625 |            4.0467 |     4.5087 |     4.4929 |
| TiN-N    | (001)   |              2.8787 |            2.8641 |     6.7717 |     6.7571 |
| TiN      | (110)   |              4.015  |            3.9755 |     2.904  |     2.8645 |
| Ti 2 N   | (001)   |              2.5569 |            2.5403 |     3.5622 |     3.5456 |
| Ti 2 N   | (110)   |              3.0394 |            3.0112 |     3.3055 |     3.2772 |
| Ti 3 N 2 | (001)   |              5.6409 |            5.6373 |     4.5733 |     4.5697 |
| Ti 3 N 2 | (110)   |              2.14   |            2.163  |     3.4979 |     3.5208 |

The calculated work functions Φ S and Φ B for Ti-N systems along with the E F slab , and E ∗ F bulk , for two different planes namely (001) and (110) are presented in table II. The energy as a function of distance plots to calculate the work function along (001) plane for Ti-N compounds are shown in Fig 8. The (001) surface of TiN is terminated with Ti or N atoms. The calculated value of workfunctions Φ S = 4 5087 and Φ . B = 4 4929 eV along (001) . surface with Ti atom termination are in good agreement with 4.40 -4.53 eV measured with Kelvin Probe Force Microscopy studying surface effect and controllable aging process on effective work function of TiN films [55] and 4.30 -4.65 eV measured on studying TiN's work function dependency on the film thickness and nitrogen flow during reactive sputter deposition [56]. The calculated Φ S and Φ B work functions for TiN crystal along the (001) plane are larger with N atom termination than the Ti atom termination. The calculated work function for (110) face of TiN Φ S = 2 9040 and Φ . B = 2 8645 eV are . also close to 3.17 eV calculated using a different computational package in Ref. [57]. The calculated values of the work function are as low as 3.2772 eV along (110) face for Ti 2 N and 3.3479 eV along (001) face for Ti N . 3 2 Being low work function MXene electrides, Ti 2 N and Ti N 3 2 can be used as efficient electron emitters in field emission displays that have better brightness, quicker response, and wider angular view than liquid crystal displays and organic light-emitting diodes.

## D. Bader charge analysis

The Bader charge analysis method [58, 59] in conjunction with DFT provides a clear picture of electronic charge distribution in a solid or a molecular system, the extent of charge transfer from one atom to the other in a compound, the oxidation state of elements in the compound, and the bonding characteristics. In nonelectrides, the absence of IAEs makes computation easy

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00897v1-figure-33ccd5dc-9417-42c5-9ac9-fd04f33160c1.jpg)

FIG. 9. 3 D Isosurfaces showing the magnetization density for AFM Ti N 3 2 crystal in the units of µ B / ˚ A 3 with isosurface value of 0.6 (a), 2 D map of MD along (110) side view (b), and 2 D maps of MD along (001) plane slicing through plane containing B empty spheres (c) and the plane containing A empty spheres. The color scale bar shows magnetization density in the range of -0 006 to 0 006 . . µ B / ˚ measured in color. A 3

FIG. 10. The same as in Fig. 9, but for FM Ti N 3 2 crystal.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00897v1-figure-35c6ad74-5eed-42a1-b1bc-bb11a103ed16.jpg)

FIG. 11. Spin structure models of optimized Ti N 3 2 crystal.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00897v1-figure-ea7dd79e-7463-4701-bcf9-394dafa55e27.jpg)

as electrons are concentrated only around the corresponding host atom species. However, the same is not true for electrides such as Ti 2 N and Ti N 3 2 as electrons are neither fully transferred nor shared, resulting in a positively charged lattice framework. To accurately account for the contribution of IAES in the electrides of our interest, the Bader charge analysis is performed by using total charge density with reference to the ELF, where the volume of the Bader basin around IAEs is approximated by estimating the region of the non-nuclear attractor (NNA), which can be done by locating the positions for electron density Laplacian's local maxima [60].

FIG. 12. Bandstructure of Ti N 3 2 monolayer along its high symmetry points and its atom resolved projected density of states. Fermi level (E F ) is shifted to zero energy level.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00897v1-figure-237ac7bf-2aa7-4399-b02f-9b00882ef127.jpg)

The NNAs have localized electronic charge density that is not associated with any atomic nucleus but is significant enough to affect electride materials' physical and chemical properties.

TABLE III. The average charge transfer or oxidation states on atomic species of Ti-N with reference to ELF. The plus sign denotes the charge lost, and the minus sign denotes the charge gained during the charge transfer process. The total charge is conserved.

| compounds           |    Ti |      N | E      | E A    | E B    |
|---------------------|-------|--------|--------|--------|--------|
| TiN                 | 2.566 | -2.566 |        |        |        |
| Ti 2 N              | 2.431 | -2.318 | -2.543 |        |        |
| Ti 3 N 2 AFM UP     | 2.405 | -2.525 | -1.083 | -1.880 | -0.285 |
| Ti 3 N 2 AFM DN     | 2.403 | -2.524 | -1.080 | -1.880 | -0.280 |
| Ti 3 N 2 FM UP      | 2.454 | -2.501 | -1.181 | -1.730 | -0.633 |
| Ti 3 N 2 FM DN      | 2.365 | -2.54  | -2.016 | -2.016 | 0.000  |
| Mono Ti 3 N 2 FM UP | 2.422 | -2.324 | -0.654 | -0.810 | -0.498 |
| Mono Ti 3 N 2 FM DN | 2.356 | -2.658 | -0.438 | -0.708 | -0.169 |

The average value of charge transfer for an i th atomic species is calculated using q i = N V E,i -Q B,i , where N V E,i is the number of pseudo-potential valence electrons and Q B,i is the Bader charge enclosed within the Bader basin of the i th species and summarized in Table III. The charge transfer value gives an insight into the degree of oxidation or reduction an atomic species undergoes. A positive value of charge transfer of Ti and a negative value of N and Es indicates that titanium gets oxidized while nitrogen and empty spheres get reduced. Empty spheres in the different A and B planes in Ti 3 N 2 do not equally reduce. Interesting things happen in the FM configuration of Ti N 3 2 up spin state, both A and B have negative oxidation values. For the down spin state, however, only the empty spheres in the A plane have negative oxidation values, while the empty spheres in the B plane have an oxidation value of zero, which means they neither oxidize nor reduce, as shown in Fig. 7B.

FIG. 13. The optimized structure of Ti N 3 2 monolayer viewing along the rotated (100) plane (A) and the ELF of Ti N 3 2 monolayer for up spin (B) and down spin(C) plotted with isosurface value of 0.6. Black color circles denote empty spheres. In each (A) and (B), 3 D isosurfaces of ELF (a) and 2 D contour plot mapping of the ELF viewing through (110) side plane (b) and the same taking a slice through (001) plane are presented. The color scale bar shows the ELF value measured in color with blue for minimum (0.0) and red for maximum (0.7).

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00897v1-figure-fec605cf-0202-4e18-a105-ae39473ab640.jpg)

FIG. 14. 3 D plot of MD for Ti N 3 2 monolayer viewing along rotated (100) plane (a) and 2 D contour maps of MD along the (110) plane. The color scale bar shows magnetization density in the range of -0 006 to 0 006 . . µ B / ˚ measured in color with A 3 minimum for blue and maximum for red colors.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00897v1-figure-46e6af6b-c37c-474b-ad6e-19aafec9b88b.jpg)

## V. MAGNETISM

Nonmagnetic electrides with low work function can be excellent electron emitters, whereas magnetic electrides with low work function are excellent for magnetic tunnel junctions, and magnetic electrides with IAEs having the same spin can be used in spintronic devices for spin injection [61]. The spin-polarized nature of IAEs provides a step towards the correlated electronic phases in the quantum regime [62, 63].

Non-electride TiN and electride Ti N 2 are nonmagnetic; however, Ti N 3 2 has a complex magnetic structure. The structural analysis shows that FM and AFM states have competing stability, with the AFM configuration being slightly more stable. To understand the magnetic dipole alignment and distribution within the Ti N 3 2 crystal, we calculated the magnetization density (MD) extracted from the difference in spin up and down components of pre-converged charge density in the unit of µ B / ˚ A 3 for both AFM and FM configurations with an isosurface value of 0.6 and presented along with 2D contour plots of MD map along the specified plan view in Figs. 9 and 10, respectively [64]. In both AFM and FM configurations, Ti N 3 2 has the magnetization density concentrated about NNA and has a vanishing value at the edge. A unit cell of Ti N 3 2 has two formula units, with 6 Ti atoms whose calculated values of magnetic moments are -0.229, 0.229, 0.229, -0.229, 0.000, and

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00897v1-figure-5ff12504-ddd4-49cd-91b2-45bb3e37a054.jpg)

FIG. 15. Phonon dispersion curves and phonon density of states (PHDOS) for (a) TiN, (b) Ti N, and (c) a monolayer of 2 Ti N . 3 2

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00897v1-figure-e6b6ea72-add1-4897-96c2-9492982387b2.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00897v1-figure-3ef40736-1f50-4b18-ab27-f6aee0333097.jpg)

0.000 µ B . Nitrogens have no magnetic moment, while the four empty spheres have magnetic moments of -0.029, 0.029, 0.217, -0.217 µ B .

Any change in the atomic coordination and local environment, symmetry, change in electronic structure, and surface and edge effect significantly alter the magnetism of a magnetic material [65-68]. The MD in each configuration is concentrated near and at the IAEs, so the Ti atoms at the corners of the crystal receive a completely different electronic environment than the Ti residing inside the crystal. This leads to a complete loss of the magnetic moment for atoms at the surface and corners, while non-zero magnetism for Ti resides in the inner part of the crystal.

the authors studied the impact of surface oxidation on 2D Ti N 3 2 MXene and analyzed its electron transport and optical properties and found that surface oxidation decreases metallic conductivity, absorption, and reflectivity in the visible range.

The vanishing magnetic moment for Ti atoms at the edge but the unidirectional alignment of spins in the remaining atoms in each layer results in each layer in Ti 3 N 2 being ferromagnetic. However, the adjacent Ti N 3 2 layers are ordered antiferromagnetically, resulting in zero magnetic moments in the bulk Ti 3 N 's ground state. Hexago2 nal symmetry and interlayer antiferromagnetism indicate the possibility of altermagnetism in Ti 3 N . However, ro2 tation and translation operations in Ti 3 N cannot remove 2 the spin degeneracy, confirming that the bulk Ti N 3 2 is only antiferromagnetic but not altermagnetic.

## VI. TI N 3 2 MONOLAYER

Competing AFM and FM configurations in Ti N 3 2 motivates a thorough study of a Ti N 3 2 single layer. The absence of inter-layer interactions, larger surface area, and electron confinement effect cause a single-layered 2 D structure to have different physical properties than a 3D bulk counterpart. Transition metal doped Ti N 3 2 monolayers are good HER catalysts, although the nitrogenvacancy and nitrogen-vacancy complex with most of the transition metal dopants are endothermic [69]. Pristine and passivated Ti N 3 2 monolayers remain metallic even after Li + , Na + , and Mg ++ ions adsorption with relatively low diffusion barriers providing faster charge and discharge rates if used as anode materials [70]. In Ref. [71],

Ti N 3 2 monolayer stabilizes in its ferromagnetic configuration and has a repeating pattern of equilateral triangles with reflections across certain axes when viewed from the c -axis as shown in Fig. 13(A) and falls in hexagonal p m 3 1 space group with the international number of 164. The optimized Ti N 3 2 monolayer has a lattice parameter of 2.9902 ˚ A in which two lattice vectors in a plane make 120 0 angle. To have an interlayer interaction-free 2 D structure, we introduced a vacuum space, which in turn optimized with a single layer with a c value of 31.496 ˚ . A Our first-principle DFT study on the electronic structure and magnetism of Ti N 3 2 monolayer indicates it to be a ferromagnetic electride, which has not yet been reported in the literature. The nominating spin-up IAEs form a two-dimensional hexatic pattern on the surfaces, whereas the spin-down IAEs are distributed in a regular triangular fashion, forming a novel Wigner crystal.

The Ti N 3 2 monolayer is metallic, and it has a nonoverlapping energy level for different spins, as shown in Fig. 12. The atomic species resolved density of states shows that the empty spheres contribute significantly to the density of states, revealing that Ti 3 N 2 monolayer is also an electride. The density of states has a significantly lower value for the up spin than the down spin. Anionic electrons localization in the plane above and below the 2 D FM Ti N 3 2 layer is larger for electrons' up spin than the down spin alignment (see Figs. 13(B) and (C)). The MD's 3 D isosurface plot and its 2 D contour map, as shown in Fig. 14, shows that the Ti residing in the inner part of the monolayer has non-zero magnetic moment while Tis at the corner have vanishing magnetic moment, which is consistent with what we have seen in its bulk configuration. Numerically, two Ti atoms in the inner part of the monolayer have a magnetic moment of 0.594 µ B each, vanishing magnetic moment for N atoms, and 0.170, 0.162, 0.141, and 0.135 µ B for four empty spheres.

## VII. PHONON'S DISPERSION AND DENSITY OF STATES

Using DFPT in conjunction with VASP and the Phonopy program, we calculated the phonon band structure, measuring the phonon frequency variation with their wave vector along the high symmetry point of irreducible Brillion Zone (BZ) boundaries and phonon density of states (PHDOS), showing the number of phonon states available at each frequency for TiN, Ti N, and 2 a monolayer of Ti N 3 2 and presented in Fig. 15. The phonon dispersion curves show nonzero values only for the real part of vibrational frequency for TiN, Ti 2 N, and a monolayer of Ti N , validating their dynamical stabil3 2 ity. The DFPT calculations reveal vibrational mode instability in Ti 3 N 2 crystal for as big as 4 × 4 × 2 supercell of 320 atoms as it has a nonzero imaginary component within the BZ. TiN, Ti 2 N, and a monolayer of Ti 3 N have 2 frequency gaps of 5.7219, 4.2261, and 4.1545 THz between their acoustic and optical phonon modes, which in eV reads 2 3664 . × 10 -2 , 1 7478 . × 10 -2 , and 1 7181 . × 10 -2 respectively.

## VIII. CONCLUSIONS

The electronic structure and magnetism of crystalline Ti-N systems, namely, TiN, Ti N, and Ti N , are studied 2 3 2

- [1] M. Naguib, M. W. Barsoum and Y. Gogotsi, Adv. Mater. , 2021, 33 , 2103393.
- [2] M. Naguib, M. Kurtoglu, V. Presser, J. Lu, J. Niu, M. Heon, L. Hultman, Y. Gogotsi and M. W. Barsoum, Adv. Mater. , 2011, 23 , 4248-4253.
- [3] M. P. Bilibana, Adv. Sens. Energy Mater , 2023, 2 , 100080.
- [4] Y. Gogotsi and B. Anasori, ACS Nano , 2019, 13 , 84918494.
- [5] K. Chaudhuri, M. Alhabeb, Z. Wang, V. M. Shalaev, Y. Gogotsi and A. Boltasseva, ACS Photonics , 2018, 5 , 1115-1122.
- [6] B. Soundiraraju and B. K. George, ACS Nano , 2017, 11 , 8892-8900.
- [7] H. Li, R. Fan, B. Zou, J. Yan, Q. Shi and G. Guo, J. Nanobiotechnology , 2023, 21 , 73.
- [8] F. Shahzad, M. Alhabeb, C. B. Hatter, B. Anasori, S. M. Hong, C. M. Koo and Y. Gogotsi, Science , 2016, 353 , 1137-1140.
- [9] D. Bury, M. Jakubczak, R. Kumar, D. ´ cie˙ zy´ska, J. BoS n gacki, P. Marcinowski and A. M. Jastrzebska, MRS Bulletin , 2023, 48 , 271-282.
- [10] Y. Pei, X. Zhang, Z. Hui, J. Zhou, X. Huang, G. Sun and W. Huang, ACS Nano , 2021, 15 , 3996-4017.
- [11] M. Zubair, M. M. Ul Hassan, M. T. Mehran, M. M. Baig, S. Hussain and F. Shahzad, Int. J. Hydrogen Energy. , 2022, 47 , 2794-2818.
- [12] A. Djire, H. Zhang, B. J. Reinhart, O. C. Nwamba and N. R. Neale, ACS Catalysis , 2021, 11 , 3128-3136.

using the first principle DFT calculating various physical quantities such as formation enthalpy, cohesive energy, exfoliation energy to extract layered MXene structure from Al-based MAX phase, band structure, DOS, charge transfer, work function, ELF, magnetic moment, and MD. Producing the results consistent with the available results for well-studied crystalline TiN systems and extending a detailed analysis to Ti 2 N and Ti N 3 2 MXene, we discover a new family of electride materials having low work functions at least along some crystallographic direction and variety in magnetic properties. TiN, which does not follow the general formula of MXene, is neither a magnetic nor an electride. Ti N is a non-magnetic MX2 ene electride. Bulk Ti N 3 2 Mxene has competing AFM and FM configurations with a slight preference for AFM configuration, which is a magnetic electride for both configurations. Calculations show a monolayer of Ti N 3 2 MXene is a ferromagnetic electride.

## ACKNOWLEDGEMENTS

The Department of Energy BES-RENEW award number DE-SC0024611 has supported this research work. D. T. and S. K. acknowledge the Center for Computational Assisted Science and Technology (CCAST) at North Dakota State University for providing computational resources for this study. .

- [13] D. Johnson, H.-E. Lai, K. Hansen, P. B. Balbuena and A. Djire, Nanoscale , 2022, 14 , 5068-5078.
- [14] K. Allen-Perry, W. Straka, D. Keith, S. Han, L. Reynolds, B. Gautam and D. E. Autrey, Materials , 2021, 14 , 694.
- [15] N. C. Frey, C. C. Price, A. Bandyopadhyay, H. Kumar and V. B. Shenoy, in Predicted Magnetic Properties of MXenes , ed. B. Anasori and Y. Gogotsi, Springer International Publishing, Cham, 2019, pp. 291-300.
- [16] Y. Limbu, G. C. Kaphle, A. L. Karn, N. K. Shah, H. Paudyal and D. Paudyal, J. Magn. Magn. Mater. , 2022, 563 , 169895.
- [17] T. Sakhraoui and F. Karlick´, y ACS Omega , 2022, 7 , 42221-42232.
- [18] N. Zhang, Y. Hong, S. Yazdanparast and M. A. Zaeem, 2D Materials , 2018, 5 , 045004.
- [19] L. Pauling, The Nature of the Chemical Bond and the Structure of Molecules and Crystals: An Introduction to Modern Structural Chemistry , Cornell University Press, 1960.
- [20] W. Research, ElementData , https://reference. wolfram.com/language/ref/ElementData.html , 2014, Accessed on 12-March-2024.
- [21] L. Nickelson, Electromagnetic Theory and Plasmonics for Engineers , Springer Singapore, 2019, p. 112.
- [22] C. Liu, S. A. Nikolaev, W. Ren and L. A. Burton, J. Mater. Chem. C , 2020, 8 , 10551-10567.
- [23] R. H. Huang and J. L. Dye, Chem. Phys. Lett. , 1990, 166 , 133-136.

- [24] X. Tang, A. E. Kuehster, B. A. DeBoer, A. D. Preston and K. Ma, Ceram. Int. , 2021, 47 , 16614-16631.
- [25] H. Hosono, S.-W. Kim, S. Matsuishi, S. Tanaka, A. Miyake, T. Kagayama and K. Shimizu, Philos. Transact. A Math. Phys. Eng. Sci. , 2015, 373 , 20140450.
- [26] J. Hu, B. Xu, S. A. Yang, S. Guan, C. Ouyang and Y. Yao, ACS Appl. Mater. Interfaces , 2015, 7 , 2401624022.
- [27] T.-N. Ye, J. Li, M. Kitano and H. Hosono, Green Chem. , 2017, 19 , 749-756.
- [28] X. Zhang, W. Meng, Y. Liu, X. Dai, G. Liu and L. Kou, J. Am. Chem. Soc. , 2023, 145 , 5523-5535.
- [29] A. Jain, S. P. Ong, G. Hautier, W. Chen, W. D. Richards, S. Dacek, S. Cholia, D. Gunter, D. Skinner, G. Ceder and K. A. Persson, APL Materials , 2013, 1 , 011002.
- [30] B. Wan, Y. Lu, Z. Xiao, J. Du, J. Wang and W. Hu, npj Computational Materials , 2018, 4 , 77.
- [31] Y. Zhang, H. Wang, Y. Wang, L. Zhang and Y. Ma, Phys. Rev. X , 2017, 7 , 011017.
- [32] R. L. Daake and J. D. Corbett, Inorganic Chemistry , 1977, 16 , 2029-2033.
- [33] K. Lee, S. W. Kim, Y. Toda, S. Matsuishi and H. Hosono, Nature , 2013, 494 , 336-340.
- [34] A. Walsh and D. O. Scanlon, J. Mater. Chem. C , 2013, 1 , 3525-3528.
- [35] L. M. McRae, R. C. Radomsky, J. T. Pawlik, D. L. Druffel, J. D. Sundberg, M. G. Lanetti, C. L. Donley, K. L. White and S. C. Warren, J. Am. Chem. Soc. , 2022, 144 , 10862-10869.
- [36] X. Zhang, Z. Xiao, H. Lei, Y. Toda, S. Matsuishi, T. Kamiya, S. Ueda and H. Hosono, Chem. Mater. , 2014, 26 , 6638-6643.
- [37] D. Y. Novoselov, D. M. Korotin, A. O. Shorikov, V. I. Anisimov and A. R. Oganov, J. Phys. Chem. C , 2021, 125 , 15724-15729.
- [38] J. Zhou, L. Shen, M. Yang, H. Cheng, W. Kong and Y. P. Feng, Chem. Mater. , 2019, 31 , 1860-1868.
- [39] D. Thapa, S. B. Sharma, D. Alf`, D. Kilin and S. Kilina, e Mater. Horiz. , 2024, 11 , 3613-3622.
- [40] J. Chae, J. Lee, Y. Oh and G. Kim, Phys. Rev. B , 2021, 104 , 125403.
- [41] A. Brager, Acta Physicochimica (USSR) , 1939, 9 , 617632.
- [42] G. Lobier and J. Marcon, Comptes Rendus Hebdomadaires des Seances de l'Academie des Sciences, Serie C, Sciences Chimiques (1966-) , 1969, 268 , 1132-1135.
- [43] J. H. Jung, C.-H. Park and J. Ihm, Nano Lett. , 2018, 18 , 2759-2765.
- [44] U. Yorulmaz, ˙ Ilker Demiro˘lu, D. Cakir, O. G¨lseren and g ¸ u C. Sevik, Journal of Physics: Energy , 2020, 2 , 032006.
- [45] G. Kresse and J. Furthm¨ller, u Comput. Mater. Sci. , 1996, 6 , 15-50.
- [46] G. Kresse and J. Furthm¨ller, u Phys. Rev. B , 1996, 54 , 11169-11186.
- [47] P. E. Bl¨chl, o Phys. Rev. B , 1994, 50 , 17953-17979.
- [48] J. P. Perdew, K. Burke and M. Ernzerhof, Phys. Rev. Lett. , 1996, 77 , 3865-3868.
- [49] K. Momma and F. Izumi, J. Appl. Crystallogr. , 2011, 44 , 1272-1276.
- [50] S. Baroni, P. Giannozzi and A. Testa, Phys. Rev. Lett. , 1987, 58 , 1861-1864.
- [51] S. Baroni, S. de Gironcoli, A. Dal Corso and P. Giannozzi, Rev. Mod. Phys. , 2001, 73 , 515-562.
- [52] A. Togo, L. Chaput, T. Tadano and I. Tanaka, J. Phys. Condens. Matter , 2023, 35 , 353001.
- [53] A. Togo, J. Phys. Soc. Jpn. , 2023, 92 , 012001.
- [54] A. D. Becke and K. E. Edgecombe, J. Chem. Phys. , 1990, 92 , 5397-5403.
- [55] Y. Zhuang, Y. Liu, H. Xia, Y. Li, X. Li and T. Li, AIP Advances , 2022, 12 , 125222.
- [56] S. A. Vitale, J. Kedzierski, P. Healey, P. W. Wyatt and C. L. Keast, IEEE Transactions on Electron Devices , 2011, 58 , 419-426.
- [57] A. Calzolari and A. Catellani, IEEE Access , 2020, 8 , 156308-156313.
- [58] W. Tang, E. Sanville and G. Henkelman, J. Phys. Condens. Matter , 2009, 21 , 084204.
- [59] E. Sanville, S. D. Kenny, R. Smith and G. Henkelman, J. Comput. Chem. , 2007, 28 , 899-908.
- [60] Q. Zhu, T. Frolov and K. Choudhary, Matter , 2019, 1 , 1293-1303.
- [61] X. Zhang, W. Meng, Y. Liu, X. Dai, G. Liu and L. Kou, J. Am. Chem. Soc. , 2023, 145 , 5523-5535.
- [62] G. Giuliani and G. Vignale, Quantum Theory of the Electron Liquid , Cambridge University Press, 2005.
- [63] S. Kim, J. Bang, C.-y. Lim, S. Y. Lee, J. Hyun, G. Lee, Y. Lee, J. D. Denlinger, S. Huh, C. Kim, S. Y. Song, J. Seo, D. Thapa, S.-G. Kim, Y. H. Lee, Y. Kim and S. W. Kim, Nat. Mater. , 2022, 21 , 1269-1274.
- [64] S. Kim, J. Bang, C. young Lim, S. Y. Lee, J. Hyun, G. Lee, Y. Lee, J. D. Denlinger, S. Huh, C. Kim, S. Y. Song, J. Seo, D. Thapa, S.-G. Kim, Y. H. Lee, Y. Kim and S. W. Kim, npj Quantum Materials , 2021, 6 , year.
- [65] J. M. D. Coey, Permanent Magnetism , Routledge, New York, 1st edn, 1999.
- [66] A. J. Freeman and R. E. Watson, Phys. Rev. , 1962, 127 , 2058-2075.
- [67] C. Kittel, Introduction to Solid State Physics , Wiley, Hoboken, NJ, 8th edn, 2004.
- [68] P. Gambardella, S. Rusponi, M. Veronese, S. S. Dhesi, C. Grazioli, A. Dallmeyer, I. Cabria, R. Zeller, P. H. Dederichs, K. Kern, C. Carbone and H. Brune, Science , 2003, 300 , 1130-1133.
- [69] I. C. Onyia, S. O. Ezeonu, D. Bessarabov and K. O. Obodo, Comput. Mater. Sci. , 2021, 197 , 110613.
- [70] H. Yu, W. Lin, Y. Zhang, Y. Li, K. Ding, S. Huang and W. Chen, RSC Adv. , 2019, 9 , 40340-40347.
- [71] S. Aghaei, I. Torres, A. Baboukani, I. Khakpour and C. Wang, Comput. Condens. Matter , 2019, 20 , e00382.
Chandra M. Adhikari Dinesh Thapa Talon D. Alexander Christopher K. Addaman Shubo Han Bishnu P. Bastakoti Daniel E. Autrey Svetlana Kilina Binod K. Rai Bhoj R. Gautam
Confinement of quasi-atomic structures in Ti$_2$N and Ti$_3$N$_2$ MXene Electrides
Metal carbides, nitrides, or carbonitrides of early transition metals, better
known as MXenes, possess notable structural, electrical, and magnetic
properties. Analyzing electronic structures by calculating structural
stability, band structure, density of states, Bader charge transfer, and work
functions utilizing first principle calculations, we revealed that titanium
nitride Mxenes, namely Ti$_2$N and Ti$_3$N$_2$, have excess anionic electrons
in their pseudo-atomic structure inside the crystal lattice, making them MXene
electrides. Bulk Ti$_3$N$_2$ has competing antiferromagnetic (AFM) and
ferromagnetic(FM) configurations with slightly more stable AFM configurations,
while the Ti$_2$N MXene is nonmagnetic. Although Ti$_3$N$_2$ favors AFM
configurations with hexagonal crystal systems having $6/mmm$ point group
symmetry, Ti$_3$N$_2$ does not support altermagnetism. The monolayer of the
Ti$_3$N$_2$ MXene is a ferromagnetic electride. These unique properties of
having non-nuclear interstitial anionic electrons in the electronic structure
of titanium nitride MXene have not yet been reported in the literature. Density
functional theory calculations show TiN is neither an electride, MXene, or
magnetic.