Search tips
Search criteria 


Logo of scirepAboutEditorial BoardFor AuthorsScientific Reports
Sci Rep. 2015; 5: 12285.
Published online 2015 July 20. doi:  10.1038/srep12285
PMCID: PMC5387395

Electronic, Dielectric, and Plasmonic Properties of Two-Dimensional Electride Materials X2N (X=Ca, Sr): A First-Principles Study


Based on first-principles calculations, we systematically study the electronic, dielectric, and plasmonic properties of two-dimensional (2D) electride materials X2N (X = Ca, Sr). We show that both Ca2N and Sr2N are stable down to monolayer thickness. For thicknesses larger than 1-monolayer (1-ML), there are 2D anionic electron layers confined in the regions between the [X2N]+ layers. These electron layers are strongly trapped and have weak coupling between each other. As a result, for the thickness dependence of many properties such as the surface energy, work function, and dielectric function, the most dramatic change occurs when going from 1-ML to 2-ML. For both bulk and few-layer Ca2N and Sr2N, the in-plane and out-of-plane real components of their dielectric functions have different signs in an extended frequency range covering the near infrared, indicating their potential applications as indefinite media. We find that bulk Ca2N and Sr2N could support surface plasmon modes in the near infrared range. Moreover, tightly-bounded plasmon modes could exist in their few-layer structures. These modes have significantly shorter wavelengths (few tens of nanometers) compared with that of conventional noble metal materials, suggesting their great potential for plasmonic devices with much smaller dimensions.

Electrides are a special kind of ionic solids with cavity-trapped electrons serving as the anions1,2. These electrons are spatially separated from the cations in a regular crystalline array, and are not bound to any particular atom, molecule, or bond. The first crystalline electride, Cs+·(18-crown-6)2e, was synthesized in 19833, and several other electrides have been successfully discovered or predicted later on4,5,6,7,8,9,10,11,12,13,14. The properties of electrides are closely connected to the topology and geometry of the cavities which confine the anionic electrons2. The early examples of electride demonstrate confinement of zero-dimensional cavities or one-dimensional weakly linked channels2. Recently, a new type of electride with two-dimensional (2D) confinement of anionic electrons was discovered in dicalcium nitride (Ca2N) which has a layered structure15. From formal valence consideration, each unit cell of Ca2N should have one excess electron. Experimental measurements on properties such as the transport parameters and work functions, combined with first-principles calculations, indeed proved that there is a built-in anionic electron layer confined between the calcium layers, which agrees well with the chemical formula [Ca2N]+·e15. This discovery of 2D electride materials (here 2D refers to the anionic electron confinement topology) has generated great research interest. First-principles calculations have predicted several other 2D electrides16,17, including the other two alkaline earth sub-nitrides Sr2N and Ba2N18. Recent theoretical work also indicated that the monolayer Ca2N might be mechanically exfoliated from the bulk, meanwhile, its electron confined layer could be maintained, and suitable encapsulation layers were proposed to protect it in ambient environment19.

Motivated by these recent experimental and theoretical progress on 2D electrides, and also by the surge of research activities on 2D materials in recent years initiated by the discovery of graphene, in this work, we conduct a systematically investigation of the electronic, dielectric, and plasmonic properties of monolayer and few-layer alkaline earth sub-nitrides Ca2N and Sr2N. We find that besides Ca2N, Sr2N is also dynamically stable down to monolayer thickness. Their phonon spectra exhibit characteristic features of 2D materials. For thickness greater than 1-ML, besides the two surface electron bands, there are additional bands crossing the Fermi energy which are due to the anionic electron layers confined in the 2D interlayer regions between the [X2N]+ layers. These interlayer bands are lower in energy than the surface bands and are nearly degenerate, indicating that they are only weakly coupled. The thickness dependence of properties such as surface energy, work function, and dielectric function are analyzed in detail. The change in property with thickness is most dramatic when going from 1-ML to 2-ML, which is associated with the appearance of the first interlayer electronic band. We find that Sr2N has a lower surface energy and lower work function compared with Ca2N. Due to their intrinsic structural anisotropy, these layered materials have highly anisotropic dielectric functions. In particular, the in-plane and out-of-plane components of dielectric function can have different signs, which occurs already for the bulk form in the near infrared frequency range. Moreover the dissipation is low in the range for the bulk materials. This shows that bulk Ca2N and Sr2N could be ideal low-loss indefinite media20. We further investigate the possibility of guided surface plasmon modes for these metallic materials and show that bulk Ca2N and Sr2N have good performance comparable to the noble metals but operating at lower frequencies, while their thin film structure can have strongly bounded plasmon modes which have advantage of much shorter wavelength compared to the conventional noble metals. Our findings thus identify a class of promising 2D electride materials and reveal their great potentials for future electronics and plasmonics applications.


Crystal structure, dynamical stability, and surface energy

The alkaline earth subnitrides X2N (X = Ca, Sr, and Ba) in their bulk form can be synthesized by a direct solid-state reaction of the elements or through chemical reduction of the corresponding sesquinitrides21. Their structures consist of alternating ABC-stacked (X-N-X) hexagonal layers, which are in the An external file that holds a picture, illustration, etc.
Object name is srep12285-m1.jpg space group (anti-CdCl2-type) with a high c/a ratio, where a and c are the in-plane and out-of-plane unit cell dimensions22,23. Each (X-N-X) layer is closely packed. For example, the [Ca2N] layer has a thickness of 2.51 Å which is much smaller than that of the ordinary fcc Ca (111) layers (3.18 Å). This is typically attributed to the ionic bounding in the layer. Meanwhile, the separation between two neighboring layers is relatively large, leading to 2D confined spaces for the anionic electrons15.

In this work, we investigate the properties of monolayer and few-layer thin films of these subnitrides (here one monolayer (ML) refers to one (X-N-X) unit layer). Since large interlayer separations (>3.5 Å) usually signal a possible role of long-range van der Waals interaction, the van der Waals corrections was included in the first-principles calculation24. The calculation details are described in the Method section. As an example, the crystal structure of 3-ML X2N (X = Ca, Sr) is schematically shown in Fig. 1. The optimized lattice parameters for isolated Ca2N and Sr2N thin films with thickness from 1-ML to 5-ML are listed in Table 1. (A comparison between DFT and experimental results for their bulk materials are shown in Supplementary Table S1, which exhibit good agreement between the two.) One observes that despite slight variations with respect to the number of layers, the calculated in-plane lattice constants, layer thicknesses, and interlayer distances are close to the experimental values of the corresponding bulk structures (for bulk Ca2N, a = 3.62 Å, layer thickness = 2.51 Å, and layer separation = 3.86 Å; for bulk Sr2N, a = 3.85 Å, layer thickness = 2.71 Å, and layer separation = 4.19 Å)23. Generally, the interlayer spacing is much larger than the layer thickness by about 50%. Sr2N has a larger interlayer spacing (around 4 Å) than Ca2N due to the larger atomic number of Sr.

Figure 1
a Top and b side view of crystal structure of 3-ML X2N (X = Ca, Sr). It consists of three (X-N-X) unit layers with ABC-stacking. The symbols L1, L2, and G12 refer to the thicknesses of outermost layer, the next outermost layer, and the ...
Table 1
Lattice parameters of Ca2N and Sr2N thin films with thicknesses from 1-ML to 5-ML.

The dynamical stability of 1-ML Ca2N has been studied in Ref. 19. Here we find that few-layer Ca2N and Sr2N (including Sr2N monolayer) are also dynamically stable through analysis of their phonon spectra. In contrast, the similar structures of 1-ML and 2-ML Ba2N are not stable due to the presence of imaginary frequencies in their phonon spectra (see Supplementary Fig. S1 for the case of monolayer Ba2N). Therefore, we only focus on Ca2N and Sr2N in the present work. As representative examples, the phonon dispersions of monolayer Ca2N and Sr2N are shown in Fig. 2. The absence of imaginary frequencies in the whole Brillouin zone demonstrates the dynamical stability of the corresponding structures. The in-plane and out-of-plane transverse acoustic modes are not degenerate due to the structural anisotropy which is typical for layered materials. More importantly, while the in-plane transverse acoustic modes have a linearly dependence on their wavevectors in the vicinity of Γ-point, the out-of-plane acoustic (ZA) phonons exhibit a quadratic dispersion around Γ-point. The parabolic dispersion of ZA modes is a characteristic feature of layered materials, which is consistent with the macroscopic elastic theory of thin plates25. This feature has been frequently observed in other layered 2D materials, such as graphene26, layered transition metal dichalcogenides27, and phosphorene28.

Figure 2
Phonon dispersions.

We then calculate the surface energies of 1-ML to 5-ML Ca2N and Sr2N and the results are plotted in Fig. 3. Because of their layered structure and the large interlayer separations, the surface energies are relatively low. The calculated values are comparable to the surface energy of graphene (~12 meV/Å2)29. Previous calculations on 1-ML Ca2N have suggested its possibility to be mechanically exfoliated from the bulk. From our result in Fig. 3, the surface energy decreases with the increasing number of layers. Hence few-layer Ca2N would be more easily exfoliated compared with monolayer Ca2N, as naturally expected. Sr2N has an even larger interlayer spacing. Surface energy for 1-ML Sr2N is about 29 meV/Å2, and it decreases to about 20 meV/Å2 for 5-ML. With the same number of layers, Sr2N’s surface energy is smaller than that of Ca2N by about 6 ~ 8 meV/Å2. This implies that Sr2N monolayer and few-layers could be more easily obtained by mechanical exfoliation using Scotch tapes or AFM tips.

Figure 3
Surface energies for few-layer Ca2N and Sr2N as a function of thickness d from 1-ML to 5-ML.

Electronic structure and 2D confined electron layers

Bulk Ca2N and Sr2N have conduction electrons confined in the 2D interlayer regions playing the role of anions15,18. One key question is whether these anionic electron layers would be maintained in their thin film structures. We find that the answer is positive. In Fig. 4, we show the electronic band structures of X2N (X = Ca, Sr) from 1-ML to 3-ML. First, one notices that all these systems are metallic, with Fermi level lying in partially filled dispersive energy bands. In 1-ML Ca2N (and Sr2N), there are two bands crossing Fermi level (see Fig. 4a,d). By inspecting the charge density distribution, one can verify that these two bands are from the 2D confined electron layers residing on the two sides of the [Ca2N]+ layer, and were referred to as the 2D electron gas in free space states19. The two bands are energy splitted due to the coupling between the two sides. When a second [Ca2N] layer is added, i.e. for a 2-ML structure, there appears an additional band crossing Fermi level, which has similar dispersion as the two surface bands but is lower in energy (see Fig. 4b). This band is from the confined anionic electron layer between the two [Ca2N]+ layers, as we will show in the following. The energy for this 2D band is lower because the confinement in the interlayer region is stronger than that for the surface, leading to a deeper potential well. One also observes that the splitting between the two surface bands is reduced, as a result of the reduced coupling between the two with increasing structural thickness. When going to 3-ML, there are two interlayer gap regions and indeed there is one more band crossing Fermi level (see Fig. 4c). One observes that the two bands from the interlayer-confined 2D electrons almost coincide, indicating that these states are strongly confined and the coupling between neighboring interlayer regions is very small. This clearly demonstrates the 2D character of the confinement topology. The above features persist when more [Ca2N] layers are added. Then more interlayer-confined anionic electron bands will appear and are nearly degenerate in energy. Their dispersion and bandwidth are almost independent on the thickness. For Sr2N, the band structure is quite similar to that of Ca2N with the same number of layers (with a slight decrease of the bandwidth, c.f. Fig. 4d–f), indicating that the 2D electride character is also maintained in its thin film form.

Figure 4
Electronic band structures of few-layer Ca2N and Sr2N.

To visualize the real space distribution of the electronic states, we take 2-ML Ca2N as an example, and plot their partial electron densities for states within 0.05 eV around Fermi energy EF in Fig. 5a. It is clear that the states around Fermi level are mainly located in three confined 2D regions: outside the two surfaces and in the interlayer space between the [Ca2N]+ layers. We further plot the total charge densities of each of the three conduction bands for 2-ML Ca2N. As shown in Fig. 5b–d, the lower band is from the 2D states in between the [Ca2N]+ layers, while the two higher bands are from the states confined on the surfaces. Similar result holds for Sr2N. This analysis validates our claims in the previous discussion.

Figure 5
Electron density distribution in 2-ML Ca2N.

Furthermore, electron-localization function (ELF) is useful for the analysis of the degree of electron localization and the bonding character30,31. From the ELF maps of 2-ML Ca2N and 2-ML Sr2N (see Supplementary Fig. S3), one finds that the bonding between the confined electron layers and the [Ca2N]+ ([Sr2N]+) layers is of ionic type. When one valence electron is removed, the surface confined layers will be vacated while the interlayer anionic electron layer remain largely intact (c.f. Supplementary Fig. S3).

Previous studies found that bulk Ca2N has a highly anisotropic work function ranging from 2.6 eV (for (100) surface) to 3.4 eV (for (001) surface)15,32. In Fig. 6, we plot our results of the out-of-plane work function as a function of the number of layers for both Ca2N and Sr2N. One notes that there is a sharp decrease of the work function from 1-ML to 2-ML for both materials, due to the appearance of the strongly confined interlayer states that decouples the two surfaces. For Ca2N, with increasing layer number, the work function approaches a value around 3.4 eV, which is consistent with previous studies. The work function of Sr2N is less than that of Ca2N by about 0.2 eV, showing that the electrons in Sr2N are more loosely bound. We average the partial electron density around Fermi level in the layer plane. The resulting 1D density profiles for 2-ML with and without one electron removed are plotted in Fig. 7. One observes that there are three peaks corresponding to the three anionic electron layers. After one electron is taken away, the two side peaks disappear while the central peak remains. The two side peaks of Sr2N have larger width than those of Ca2N, implying a higher degree of delocalization of these surface states hence a lower work function.

Figure 6
Work functions for few-layer Ca2N and Sr2N as a function of thickness d from 1-ML to 5-ML.
Figure 7
Partial density of states for the energy range |E  Ef| < 0.05 eV averaged over the ab plane for a 2-ML Ca2N and b 2-ML Sr2N. In each figure, the blue curve and the red curve are for the distributions ...

Dielectric function

The first-principles calculation based on DFT has proved to be a powerful tool for the study of dielectric function for metals, including ultrathin metallic films down to few-layer thickness33,34,35,36,37,38,39. In crystalline solids, the dielectric function ε(ω) consists of two contributions: a Drude-like intraband contribution and an interband contribution. The imaginary part of the interband contribution involves the interband matrix elements of the momentum operators, and can be evaluated directly in DFT34. Its real part can then be calculated via the Kramers-Kronig relation. The intraband contribution is typically treated by the Drude model, in which the plasma frequency (tensor) ωp can be evaluated by DFT34. The calculation details are presented in the Method section.

We first consider the dielectric properties of bulk Ca2N and Sr2N. Simple metals such as Au, Ag, and Al generally have a large isotropic bulk plasma frequency around 10 eV. In contrast, Ca2N and Sr2N are of layered structure, hence their dielectric properties are expected to exhibit intrinsic anisotropy. We find that bulk Ca2N has an in-plane (i.e. in the layer plane) plasma frequency ωp,xx of 3.14 eV and an out-of-plane plasma frequency ωp,zz of 0.95 eV; bulk Sr2N has an in-plane plasma frequency of 2.94 eV and an out-of-plane plasma frequency of 1.05 eV. These frequencies are significantly lower than those for simple metals, partly because the conducting electron density is lower. For each material, the out-of-plane plasma frequency is much lower compared with the in-plane plasma frequency, showing that the neighboring anionic electron layers are only weakly coupled. In Fig. 8, we plot the real and imaginary parts of in-plane and out-of-plane dielectric functions for the two materials. They exhibit the character of metallic behavior, i.e. Drude peaks at low energy due to intraband contribution and the real part of ε(ω) crossing from negative value to positive value with increasing frequency. The strong anisotropy in the dielectric function is also quite obvious. Remarkably, because of this anisotropy, the real part of in-plane dielectric function, Reεxx(ω), changes sign at a frequency which is different from that for Reεzz(ω), leading to an extended frequency window in which the two components have different signs (marked as the blue shaded regions in Fig. 8). This sign difference is the characteristic feature of so-called indefinite media20, which was proposed in the study of metamaterials and has important potential applications such as near-field focusing and building hyperlenses that can transform evanescent fields into propagating modes40,41. Previous realizations of indefinite media are mostly in artificially assembled structures which require complicated fabrication process and usually have high dissipation. Our results suggest that crystalline solids Ca2N and Sr2N in their bulk form would just be indefinite materials for a frequency range spanning the near infrared. Moreover, one notes that the imaginary parts of εxx and εzz in this frequency range are very small (~0.1), implying that they could be ideal low-loss indefinite materials.

Figure 8
Anisotropic dielectric functions for bulk X2N: a for Ca2N and b for Sr2N. The real and imaginary parts of in-plane component εxx and out-of-plane component εzz are plotted using different colors. The blue shaded region in each figure indicates ...

Next we focus on the dielectric functions of Ca2N and Sr2N few-layer structures. With decreasing film thickness, one expects that the anisotropy effects would be even stronger. The intraband contribution to the out-of-plane dielectric response εzz becomes negligible, and it has been shown that for simple metals which are isotropic in bulk can exhibit strongly anisotropic ε(ω) down to few-layer thickness38. In Fig. 9a, we plot the results of in-plane plasma frequency ωp,xx as a function of the number of layers for both materials. One observes that the frequency slowly decreases with increasing number of layers and approaches its bulk value. The plasma frequency of Sr2N is lower than that of Ca2N by about 0.2 eV for each thickness. The square of plasma frequency is roughly proportional to the product of the carrier density and the inverse of effective mass. In Fig. 9b, one can see that from 1-ML to 2-ML, there is in fact a sharp decrease in the carrier density, as characterized by the DOS per unit volume around Fermi level, primarily due to the appearance of the large interlayer space of 2-ML. However, from the band structure in Fig. 4b, one notices that the additional 2D electron band (marked in red color) has a smaller effective mass around Fermi level than the two surface bands, hence decreasing the average effective mass. This compensates the decrease in carrier density, resulting in an overall small change in the plasma frequency between 1-ML and 2-ML.

Figure 9
a, In-plane plasma frequency ωp,xx for Ca2N and Sr2N thin films as a function of thickness d from 1-ML to 5-ML. b, electron density of states (DOS) of Ca2N and Sr2N at Fermi energy as a function of film thickness. In each figure, the dashed lines ...

In Fig. 10, we show the interband contribution to the imaginary part of the in-plane dielectric function, Imεxx(ω), which is closely connected to band structure features. One noticeable feature is the sharp peak around 0.3 eV for both materials at 1-ML thickness. This peak drops from 1-ML to 2-ML and its position shifts to a lower energy. Its appearance can be attributed to the (almost) parallel sections of the two surface bands for 1-ML, as indicated in the red shaded regions in Fig. 4a,d. From 1-ML to 2-ML, the splitting between two surface bands decreases and the volume increases a lot due to the appearance of the large interlayer region, resulting in the observed change of the peak. For thicknesses larger than 1-ML, there is another peak around 1 eV, which is due to the transition between the interlayer 2D bands and the surface bands (c.f. Fig. 4). From the band structures, it is clear that the response below 1eV is mainly contributed by the states from the 2D anionic electron layers.

Figure 10
The interband contribution to Imεxx(ω) for a Ca2N and b Sr2N thin film structures from 1-ML to 5-ML. Results for different thicknesses are plotted using different colors.

The total in-plane and out-of-plane dielectric function with their real and imaginary parts for Ca2N few-layers are shown in Supplementary Fig. S4. For εxx, the low energy part is dominated by the intraband Drude-like contribution. Reεxx(ω) is negative below about 1eV. Meanwhile, the intraband contribution is negligible for εzz(ω). Hence Reεzz(ω) is largely positive for the low energy part. The sign difference between Reεzz and Reεxx in the low energy range again signals a possible indefinite material. Similar features are also demonstrated in the results for Sr2N few-layer structures (see Supplementary Fig. S5).

Surface plasmon modes

Surface plasmon modes are confined electromagnetic excitations propagating at an interface between a conductor and a dielectric42. It typically requires a sign change of Reε across the interface. Therefore the most commonly used materials for plasmonic applications are metals such as Au, Ag, and Al which have a range of frequencies with negative Reε(ω)43. Since Ca2N and Sr2N are also metals, one may naturally wonder whether they could also support surface plasmon modes. This question is particularly interesting when considering the frequency range in which the material shows indefinite medium property. In such case, one expects that the positive Reεzz component may compete with the negative Reεxx component and tend to destroy the bounded plasmon modes.

Let’s first consider the interface between bulk Ca2N or Sr2N and a dielectric medium characterized by a frequency-independent dielectric constant εd > 0. The interface is parallel to the layer plane. As we have shown, Ca2N and Sr2N in the bulk form already have strong anisotropy in their dielectric functions. Assume that the interface supports a transverse magnetic (TM) plasmon mode travelling along the interface with a wave vector β. Following standard derivation using Maxwell’s equations42, one obtains that

An external file that holds a picture, illustration, etc.
Object name is srep12285-m2.jpg

Here k0 = ω/c is the wave vector in vacuum, c is the speed of light, and εxx and εzz are the two components of the dielectric function for our conducting material. If εxx = εzz, the result in Eq.(1) reduces to the familiar result for isotropic metal42. The dispersion characteristics of the surface plasmon modes for both materials are shown in Fig. 11, with εd = 2.25 (appropriate for SiO2). In Fig. 11a, the modes lying to the right of the light line (in the dielectric medium) are bounded to the surface. One observes that the results show characteristic surface plasmon peaks similar to simple metals and confined plasmon modes still exist within the frequency range where the material shows indefinite medium property. A major difference is that for simple metals the peak position, known as surface plasmon frequency ωsp, occurs at higher energies, e.g. ωsp ~ 3.4 eV for Ag/SiO2 interface44; while ωsp here is much lower, around 1.1 eV–1.2 eV in the near infrared range. The Imβ shown in Fig. 11b is connected to the energy damping during the mode propagation. Since high dissipation occurs around ωsp, for practical applications, modes with frequencies less than ωsp are used. Here for Ca2N, if we take ω = 1.02 eV, the corresponding surface plasmon wavelength λsp [similar, equals] 656 nm, and we have a long propagation length L = 1/(2Imβ) [similar, equals] 5.15 μm and the decay length in the dielectric is An external file that holds a picture, illustration, etc.
Object name is srep12285-m3.jpg nm. Longer propagation length can be achieved at lower frequencies, e.g. at ω = 0.64 eV, we can have L [similar, equals] 173 μm but the mode confinement is reduced, with ld [similar, equals] 737 nm. These values are comparable to those for the noble metals (in the visible or UV frequency range) which are the usual building blocks for plasmonic devices. However, the operating frequency here is much lower. From above discussion, we see that despite the intrinsic anisotropy, bulk Ca2N and Sr2N can still support surface plasmons and could be suitable plasmonic materials in the near infrared frequency range. Experimentally, these surface plasmon modes may be excited, e.g. by prism coupling or by near-field excitation, and probed by standard optical techniques such as near-field optical microscopy, leakage radiation imaging, or scattered light imaging42.

Figure 11
Dispersion relation of surface plasmon modes for an interface between a dielectric medium with εd = 2.25 and bulk X2N (X = Ca, Sr). The light line in the dielectric medium is also shown in grey color.

We then turn to the thin films of Ca2N and Sr2N sandwiched between dielectric materials. Besides the change in the thickness-dependent dielectric function, an important effect in thin films is that the surface plasmon modes at two interfaces could couple and form two modes with opposite parity: a symmetric mode (L−) and an antisymmetric mode (L+)44. Their dispersions have been derived before in the study of ultrathin metallic films38, and are quoted in our Method section. In Fig. 12, we plot the surface plasmon dispersion characteristics for Ca2N with different film thicknesses. One observes that the antisymmetric modes are lying on the light line of the dielectric medium, indicating that they are squeezed out of the metal region forming unbounded modes propagating in the dielectric material, which is similar to the case of Au ultrathin films. Meanwhile, pronounced plasmonic peaks do show up for the symmetric modes, clearly indicating that these are bounded surface plasmon modes. Again, one observes that the variation with thickness is most dramatic between 1-ML and 2-ML. For 1-ML, the bounded modes occur around 1 eV. The high peak in Reβ vs. ω shows the modes are strongly bounded to the metallic layer. For thicknesses of 2-ML to 5-ML, the dispersions are quite close. The corresponding surface plasmon frequencies shift to around 1.2 eV. As for the imaginary part of β, 1-ML structure has relatively large values, whereas Imβ for 2-ML to 5-ML almost collapse on a single curve and deceases rapidly in the range below 0.8 eV. For Sr2N few-layers, the results for surface plasmon dispersion show similar features. The dispersions for symmetric (L−) modes are shown in Supplementary Fig. S6. The antisymmetric (L+) modes are again unbounded hence are not shown.

Figure 12
Dispersion characteristics of surface plasmon modes for Ca2N few-layers in a dielectric medium with εd = 2.25. a and b are for the antisymmetric mode (L+ mode). c and d are for the symmetric mode (L− mode). a and c show ...

Compared with their bulk results, one observes that the surface plasmon frequency ωsp is more or less the same, but both Reβ and Imβ are greatly increased by two orders of magnitude. This means that the plasmon wavelength and the confinement scale are much smaller, which are desired features. However, the propagation length is also decreased at the same time. Therefore, it is more meaningful to consider the dimensionless ratio Reβ/Imβ (which measures how many surface plasmon wavelengths can be covered before the wave loses most of its energy) as well as the wave localization (or wave shrinkage) quantified by λairsp, where λair = 2πc/ω is the wavelength in air. In Fig. 13, we show these two dimensionless characteristics for the two materials as functions of λair. The results for 3-ML and 4-ML are similar to those of 5-ML hence are not shown. One observes that for 1-ML Reβ/Imβ is small for most frequencies, due to the relatively high dissipation associated with Imβ (c.f. Fig. 12d). Generally, the wave localization reaches its peak near the surface plasmon resonance. However, there the ratio Reβ/Imβ is small due to the enhanced dissipation. While larger Reβ/Imβ can be achieved for longer wavelengths, the wave localization becomes poor in that range. This tradeoff is a typical feature for surface plasmons. For application purposes, a compromise may be reached somewhere in between when the two are comparable.

Figure 13
Characteristics Reβ/Imβ and λairsp for the symmetric surface plasmon modes in few-layer a Ca2N and b Sr2N, plotted as functions of the wavelength in air λair. λsp is the surface plasmon wavelength. The ...

Since 1-ML Ca2N and Sr2N have considerably larger dissipation, in the following we take 2-ML Ca2N and Sr2N as examples. For λair = 2 μm (0.64 eV) in the near infrared range, 2-ML Ca2N has a decay length ld = 5.4 nm in the dielectric, showing that the plasmon mode is strongly confined to the metallic layer. The ratios λairsp = 56 and Reβ/Imβ = 61. These two values are considerably larger than what noble metals could achieve at their optimal working frequencies. The corresponding values for Sr2N at λair = 2 μm (ω = 0.64 eV) are given by ld = 4.4 nm, λairsp = 68 and Reβ/Imβ = 17. In comparison, the surface plasmon modes for noble metals such as Au and Ag are unbounded at this wavelength. Noticeably, the surface plasmon wavelengths in these materials are very small, λsp = 34.9 nm for Ca2N and λsp = 28.9 nm for Sr2N (at ω = 0.64 eV), while the shortest λsp for Au and Ag (occur in the visible light range) would be at least larger than 100 nm. These suggest that Ca2N and Sr2N thin films have great potential for making plasmonic devices operating in the near infrared range with much smaller scales.


From our previous analysis of the electronic structures of these electride materials, one expects that the loosely bound surface anionic electrons are highly reactive in ambient conditions. They may act as good electron donors and as excellent catalysts for chemical reactions11. For physical applications, one needs to stabilize its property, e.g. by effective encapsulation. For 1-ML Ca2N, a possible encapsulation scheme using 2D insulating layers of graphane was proposed19. In the case of few-layer alkaline earth subnitrides, since as we discussed each interlayer anionic electron layer are strongly confined in 2D regions and the coupling between them is small, one expects that the most reactive electrons are from the surface layers, while the electron layers inside should be less reactive. A detailed study of this point and possible encapsulation schemes for Sr2N will be deferred to a future work.

For electride materials, the topology of the cavities confining anionic electrons is one key factor determining their properties. Although the cavity topologies for both Ca2N and Sr2N are similar, compared with Ca2N, Sr2N has a smaller electrostatic potential associated with the larger atomic number of Sr. This was reflected in its lower work function, and may also lead to a high electron mobility for Sr2N few-layers18. Previous experimental studies have shown that bulk Ca2N has high mobility of 520 cm2/(V·s)15. One expects that the mobility for Sr2N may be even higher. Hence Sr2N few-layers as new 2D conducting materials could have a good potential for electronics applications.

As having been demonstrated for other 2D layered materials, strain engineering has proved to be a powerful tool to modify and control the material properties. The layered materials discussed here have a large interlayer spacing (>3.5 Å) between the [X2N]+ layers. We also expect that applying strain could be a good method to tune the material properties such as the dielectric function, the plasmonic dispersion, and the carrier mobility. A systematic study of the strain effects is currently underway.


First-principles calculations

First-principles calculations are carried out using the Vienna ab-initio simulation package (VASP)45,46, based on the density functional theory (DFT). The exchange-correlation functional is treated using Perdew-Burke-Ernzerhof generalized gradient approximation47. The projector augmented wave (PAW) method48 is employed to model interactions between electrons and ions. The treated valence electrons are the 3p4s, 4s4p5s, and 2s2p for the Ca, Sr, and N atoms, respectively. The cutoff for plane-wave expansion is set to be 600 eV. The vertical distance between thin films (the thickness of the vacuum gap) is at least 18 Å, which is large enough to avoid artificial interactions between the film and its periodic images. Both the atomic positions and lattice constant were fully relaxed using conjugate gradient method. The convergence criteria for energy and force were set to be 10−5 eV and 0.01 eV/Å, respectively. DFT-D2 method was applied to describe the long-range van der Waals interaction. The Brillouin zone integrations have been carried out on a Γ-centered k-mesh. Monkhorst-Pack k-point meshes49 with sizes of 15 × 15 × 1 and 31 × 31 × 1 were used for geometry optimization and static electronic structure calculation, respectively. In the later study of optical properties, the sizes of k-mesh are significantly increased to 61 × 61 × 1 and 41 × 41 × 9 for thin films and bulk respectively, to achieve highly converged results. For the integration over the Brillouin zone in calculating dielectric functions, we used the first order Methfessel-Paxton method50 with a value of 0.1 eV. The phonon dispersions of the structures were calculated by using density functional perturbation theory as implemented in the PHONOPY code51,52.

Calculation of dielectric functions

The optical properties of solids are mainly due to the response of the electron system to a time-dependent electromagnetic perturbation. For metals, the optical complex dielectric function consists of interband and Drude-like intraband contributions:

An external file that holds a picture, illustration, etc.
Object name is srep12285-m4.jpg

The imaginary part of the interband part can be calculated using the results from DFT calculations as34

An external file that holds a picture, illustration, etc.
Object name is srep12285-m5.jpg

where α and β refer to Cartesian coordinates, eα(β) are unit vectors, V is the volume of the unit cell, An external file that holds a picture, illustration, etc.
Object name is srep12285-m6.jpg and En,k are the periodic part of the Bloch wave function and the corresponding eigenenergy for band n and wave vector k, and fnk is the Fermi-Dirac distribution function. The real part of interband contribution can be obtained through the Kramers-Kronig relation. The intraband contribution is usually modeled by the Drude model:

An external file that holds a picture, illustration, etc.
Object name is srep12285-m7.jpg
An external file that holds a picture, illustration, etc.
Object name is srep12285-m8.jpg

Here γ is a life-time broadening obtained either from a higher-order calculation or from experiments. In our calculation, we used the experimental determined electron life-time (of 0.6 ps)15 for bulk Ca2N to estimate γ (~1.1 meV). The same value was also used for the calculations of Sr2N due to their similar electronic structures. We have checked the sensitivity of our results’ dependence on γ by repeating the calculations of ε(ω) with γ value varying from γ/10 to 5γ. The obtained results are of little difference. This is because the value of γ is already quite small (reflecting the fact that these materials are good metals). The ωp,αβ is the plasma frequency tensor which can be calculated using

An external file that holds a picture, illustration, etc.
Object name is srep12285-m9.jpg

The dielectric functions and plasma frequencies were suitably renormalized to exclude the vacuum region from the unit cell in our calculations.

Surface plasmon modes calculations

For thin film structures, the plasmon modes at two surfaces would couple and form two modes with different parity44. The equations governing their dispersions have been derived before38:

An external file that holds a picture, illustration, etc.
Object name is srep12285-m10.jpg
An external file that holds a picture, illustration, etc.
Object name is srep12285-m11.jpg

The first equation above is for the antisymmetric (L+) mode and the second equation is for the symmetric (L−) mode. Here L in the equations is the film thickness. Other quantities in these two equations are defined in the main text. In our calculation, we solve the two equations numerically using a two-dimensional unconstrained Nelder-Mead minimization algorithm53 with a tolerance of 10−13 nm−1 in the complex wave vectors.

Additional Information

How to cite this article: Guan, S. et al. Electronic, Dielectric, and Plasmonic Properties of Two-Dimensional Electride Materials X2N (X=Ca, Sr): A First-Principles Study. Sci. Rep. 5, 12285; doi: 10.1038/srep12285 (2015).

Supplementary Material

Supplementary Information:


The authors thank D.L. Deng for helpful discussions. This work was supported by the MOST Project of China (Nos. 2014CB920903, 2011CBA00100), the NSF of China (Nos. 11174337, 11225418) , the SRFDPHE of China (No. 20121101110046) and the SUTD-SRG-EPD2013062.


The authors declare no competing financial interests.

Author Contributions Y.Y. and S.A.Y. conceived the idea and supervised the work. S.G. performed the calculation and the data analysis. S.G., L.Z., Y.Y. and S.A.Y. contributed to the interpretation of the results and wrote the manuscript. All authors contributed in the discussion and reviewed the manuscript.


  • Dye J. L. Electrons as Anions. Science 301, 607–608 (2003). [PubMed]
  • Dye J. L. Electrides: early examples of quantum confinement. Acc.Chem.Res. 42, 1564–1572 (2009). [PubMed]
  • Ellaboudy A., Dye J. L. & Smith P. B. Cesium 18-Crown-6 Compounds. A Crystalline Ceside and a Crystalline Electride. J. Am. Chem. Soc. 105, 6490–6491 (1983).
  • Huang R. H., Faber M. K., Moeggenborg K. J., Ward D. L. & Dye J. L. Structure of K+ (cryptand[2.2.2J) electride and evidence for trapped electron pairs. Nature 331, 599–601 (1988).
  • Dye J. L. Electrides: Ionic Salts with Electrons as the Anions. Science 247, 663–668 (1990). [PubMed]
  • Matsuishi S. et al. . High-Density Electron Anions in a Nonaporous Single Crystal: [Ca24Al28O64]4+(4e−). Science 301, 626 (2003). [PubMed]
  • Li Z. Y., Yang J. L., Hou J. G. & Zhu Q. S. Inorganic Electride: Theoretical Study on Structural and Electronic Properties. J. Am. Chem. Soc. 125, 6050–6051 (2003). [PubMed]
  • Toda Y., Yanagi H., Ikenaga E., Kim J. J., Kobata M., Ueda S., Kamiya T., Hirano M., Kobayashi K. & Hosono H. Work Function of a Room-Temperature, Stable Electride [Ca24Al28O64]4+(e)4. Adv. Mater 19, 3564–3569 (2007).
  • Kim S. W. & Hosono H. Synthesis and Properties of 12CaO·7Al2O3 Electride: Review of Single Crystal and Thin Film Growth. Philos. Mag. 92, 2596 (2012).
  • Fang C. M., de Wijs G. A., de Groot R. A., Hintzen H. T. & de With G. Bulk and Surface Electronic Structure of the Layered Sub-Nitrides Ca2N and Sr2N. Chem. Mater 12, 1847–1852 (2000).
  • Kitano M., Inoue Y., Yamazaki Y., Hayashi F., Kanbara S., Matsuishi S., Yokoyama T., Kim S. W., Hara M. & Hosono H. Ammonia synthesis using a stable electride as an electron donor and reversible hydrogen store. Nat Chem 4, 934–940 (2012). [PubMed]
  • Pickard C. J. & Needs R. J. Aluminium at Terapascal Pressures. Nat. Mater 9, 624 (2010). [PubMed]
  • Gatti M., Tokatly I. V. & Rubio A. Sodium: A Charge-Transfer Insulator at High Pressures. Phys. Rev. Lett. 104, 216404 (2010). [PubMed]
  • Miao M.-S. & Hoffmann R. High Pressure Electrides: A Predicative Chemical and Physical Theory. Acc. Chem. Res. 47, 1311 (2014). [PubMed]
  • Lee K., Kim S. W., Toda Y., Matsuishi S & Hosono H. Dicalcium nitride as a two-dimensional electride with an anionic electron layer. Nature 494, 336–340 (2013). [PubMed]
  • Inoshita T., Jeong S., Hamada N. & Hosono H. Exploration for Two-Dimensional Electrides via Database Screening and Ab Initio Calculation. Phys. Rev. X 4, 031023 (2014).
  • Tada T., Takemoto S., Matsuishi S. & Hosono H. High-Throughput ab Initio Screening for Two-Dimensional Electride Materials. Inorg. Chem. 53, 10347–10358 (2014). [PubMed]
  • Walsh A. & Scanlon D. Electron excess in alkaline earth sub-nitrides: 2D electron gas or 3D electride? J. Mater. Chem. C 1, 3525–3528 (2013).
  • Zhao S. T., Li Z. Y. & Yang J. L. Obtaining two-dimensional electron gas in free space without resorting to electron doping: an electride based design. J. Am. Chem. Soc. 136, 13313–13318 (2014). [PubMed]
  • Smith D. R. & Schurig D. Electromagnetic Wave Propagation in Media with Indefinite Permittivity and Permeability Tensors. Phys. Rev. Lett. 90, 077405 (2003). [PubMed]
  • Gregory D. H., Bowman A., Baker C. F. & Weston D. P. Dicalcium Nitride, Ca2N-a 2D “Excess Electron” Compound: Synthetic Routes and Crystal Chemistry. J. Mater. Chem. 10, 1635 (2000).
  • Baker C. F., Barker M. G. & Blake A. J. Calcium nitride (Ca2N), a redetermination. Acta Cryst. E57, i6–i7 (2001).
  • Reckeweg O. & DiSalvo F. J. Alkaline earth metal nitride compounds with the composition M2NX (M = Ca, Sr, Ba; X = [open square], H, Cl or Br). Solid State Sciences 4, 575–584 (2002).
  • Grimme S. Semiempirical GGA-type density functional constructed with a long-range dispersion correction. J Comput Chem. 27, 1787 (2006). [PubMed]
  • Zabel H. Phonons in layered compounds. J. Phys. Condens. Matter 13, 7679 (2001).
  • Liu F., Ming P. & Li J. Ab initio calculation of ideal strength and phonon instability of graphene under tension. Phys. Rev. B 76, 046120 (2007).
  • Molina-Sanchez A. & Wirtz L. Phonons in single-layer and few-layer MoS2 and WS2. Phys. Rev. B 84, 155413 (2011).
  • Zhu L., Zhang G. & Li B. Coexistence of size-dependent and size-independent thermal conductivities in phosphorene. Phys. Rev. B 90, 214302 (2014).
  • Zacharia R., Ulbricht H. & Hertel T. Interlayer cohesive energy of graphite from thermal desorption of polyaromatic hydrocarbons. Phys. Rev. B 69, 155406 (2004).
  • Becke A. D. & Edgecombe K. E. A simple measure of electron localization in atomic and molecular systems. The Journal of Chemical Physics 92, 5397–5403 (1990).
  • Silvi B. & Savin A. Classification of chemical bonds based on topological analysis of electron localization functions. Nature 371, 683–686 (1994).
  • Uijttewaal M. A., Wijs G. A. de. & Groot R. A. de. Low work function of the (1000) Ca2N surface. J. Appl. Phys. 96, 1751–1753 (2004).
  • Harl J., Kresse G., Sun L. D., Hohage M. & Zeppenfeld L. Ab initio reflectance difference spectra of the bare and adsorbate covered Cu(110) surfaces. Phys. Rev. B 76, 035436 (2007).
  • Harl J. The linear response function in density functional theory: Optical spectra and improved description of the electron correlation. Ph.D Dissertation submitted to Universitat Wien (2008).
  • He Y. & Zeng T. First-principle study and model of dielectric functions of silver nanoparticles. J. Phys. Chem. C 114 (42), 18023–18030 (2010).
  • Glantschnig K. & Ambrosch-Draxl C. Relativistic effects on the linear optical properties of Au, Pt, Pb and W. New J. Phys. 12, 103048 (2010).
  • Yan J., Jacobsen K. W. & Thygesen K. Y. First principles study of surface plasmons on Ag(111) and H/Ag(111). Phys. Rev. B 84, 235430 (2011).
  • Laref S., Cao J. R., Asaduzzaman A., Runge K., Deymier P., Ziolkowski R. W., Miyawaki M. & Muralidharan K. Size-dependent permittivity and intrinsic optical anisotropy of nanometric gold thin films: a density functional theory study. Opt. Express 21, 11827–11838 (2013). [PubMed]
  • Ming W. M., Blair S. & Liu F. Quantum size effect on dielectric function of ultrathin metal film: a first-principles study of Al(111). J. Phys. Condens. Matter 26, 505302 (2014). [PubMed]
  • Jacob Z., Alekseyev L. V. & Narimanov E. Optical hyperlens: Far-field imaging beyond the diffraction limit. Opt. Express 14, 8247–8256 (2006). [PubMed]
  • Liu Z., Lee H., Xiong Y., Sun C. & Zhang X. Far-field optical hyperlens magnifying sub-diffraction-limited objects. Science 315, 1686 (2007). [PubMed]
  • Maier S. A. Plasmonics: Fundamentals and Applications Springer, New York, 2007.
  • West P. R., Ishii S., Naik G. V., Emani N. K., Shalaev V. M. & Boltasseva A. Searching for better plasmonics materials. Laser & Photon. Rev. 4, 795–808 (2010).
  • Dionne J. A., Sweatlock L. A., Atwater H. A. & Polman A. Planar metal plasmon waveguides: frequency-dependent dispersion, propagation, localization, and loss beyond the free electron model. Phys. Rev. B 72, 075405 (2005).
  • Kresse G. & Hafner J. Ab initio molecular dynamics for open-shell transition metals. Phys. Rev. B 48, 13115–13118 (1993). [PubMed]
  • Kresse G. & Furthmüller J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. Rev. B 54, 11169–11186 (1996). [PubMed]
  • Perdew J., Burke K. & Ernzerhof M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 77, 3865–3868 (1996). [PubMed]
  • Blöchl P. Projector augmented-wave method. Phys. Rev. B 50, 17953–17979 (1994). [PubMed]
  • Monkhorst H. & Pack J. Special points for Brillouin-zone integrations. Phys. Rev. B 13, 5188–5192 (1976).
  • Methfessel M. & Paxton A. High-precision sampling for Brillouin-zone integration in metals. Phys. Rev. B 40, 3616–3621 (1989). [PubMed]
  • Togo A., Oba F. & Tanaka I. First-principles calculations of the ferroelastic transition between rutile-type and CaCl2type SiO2 at high pressures. Phys. Rev. B 78, 134106 (2008).
  • Gonze X. & Lee C. Y. Dynamical matrices, Born effective charges, dielectric permittivity tensors, and interatomic force constants from density-functional perturbation theory. Phys. Rev. B 55, 10355–10368 (1997).
  • Nelder J. A. & Mead R. A Simplex Method for Function Minimization. The Computer Journal 7, 308–313 (1965).

Articles from Scientific Reports are provided here courtesy of Nature Publishing Group