Nonsymmorphic symmetry protected node-line semimetal in the trigonal YH3.

Using ab initio calculations based on density-functional theory and effective model analysis, we propose that the trigonal YH3 (Space Group: P[Formula: see text]c1) at ambient pressure is a node-line semimetal when spin-orbit coupling (SOC) is ignored. This trigonal YH3 has very clean electronic structure near Fermi level and its nodal lines locate very closely to the Fermi energy, which makes it a perfect system for model analysis. Symmetry analysis shows that the nodal ring in this compound is protected by the glide-plane symmetry, where the band inversion of |Y+, d xz 〉 and |H1-, s〉 orbits at Γ point is responsible for the formation of the nodal lines. When SOC is included, the line nodes are prohibited by the glide-plane symmetry, and a small gap (≈5 meV) appears, which leads YH3 to be a strong topological insulator with Z2 indices (1,000). Thus the glide-plane symmetry plays an opposite role in the formation of the nodal lines in cases without and with SOC. As the SOC-induced gap is so small that can be neglected, this P[Formula: see text]c1 YH3 may be a good candidate for experimental explorations on the fundamental physics of topological node-line semimetals. We find the surface states of this P[Formula: see text]c1 phase are somehow unique and may be helpful to identify the real ground state of YH3 in the experiment.

Introduction

Topological semimetals (TSMs) have attracted great attention for both theoretical interests and experimental applications in recent years. Different from time-reversal symmetry (TRS) protected Z2 topological insulators (TIs)[1,2] which are insulating in the bulk, TSMs are materials where the conduction and the valence bands cross with each other at certain locations in the Brillouin Zone (BZ). Usually, the band crossing is protected by certain symmetries, i.e., perturbations on the Hamiltonian which respect the symmetries can not break the crossing. Recently, several types of TSMs have been proposed to investigate the fermion-like excitations, including Dirac fermions[3-6], Weyl fermions[7-10] and nodal lines[11-16]. These compounds are named as: Dirac semimetals (DSMs), Weyl semimetals(WSMs) and node-line semimetals (NLSMs), respectively.

Up to now, there have been a lot of reports about the progress in DSMs and WSMs, for example, 3D Dirac semimetals have recently been identified experimentally in Cd3As2[17-21] and Na3Bi systems[22,23]. Similarly, TaAs[9,24-26], NbAs[27], TaP[28], WTe2[29,30], and MoTe2[31] etc are verified to be WSMs experimentally in recent years. Different from DSMs and WSMs, in which the conduction and valence bands touch at discrete points, the crossings of NLSMs form a closed loop in the BZ. Although many candidates of NLSMs have been proposed and much efforts has been made to investigate them, the corresponding progress in the experiment is slow, because an open surface usually breaks the inversion or some mirror symmetries which are important to the formation of nodal lines[32].

Materials experimentally confirmed (or partially confirmed) to host nodal line include Be metal[16], ZrSiS[33], PbTaSe2[14] and ZrSiSe/ZrSiTe[34]. Therefore, theoretical predictions on more candidates of node-line semimetal are still in demand. It is well known that symmetries play important roles in identifying various of TIs and topological superconductors (TSCs)[35-39]. In fact, symmetries are also important in classifying NLSMs, for instance, three types of NLSMs protected by different symmetries have been proposed: (a) mirror symmetry protected NLSMs[8,14,16,40], (b) coexistence of TRS and inversion symmetry (IS) protected NLSMs[11-13,32] and (c) nonsymmorphic symmetry protected NLSMs[15,32].

Hydrides is a large class of materials and has been extensively investigated in many aspects, including energy storage[41] and superconductivity[42-46], etc. Since Ashcroft proposed that high T superconductivity can be obtained in hydrogen and hydrides under high pressure[47,48], many hydrides have been investigated. Yttrium-hydrogen system becomes interesting due to the same reason. For instance, a fcc YH3 has been predicted to be a superconductor with T ~ 40 K at 17.7 GPa[49]. Later work predicts that two new yttrium hydrides, i.e., YH4 and YH6, are also superconductors with T ~ 84–95 K and T ~ 251–264 K at 120 GPa, respectively[50]. Very recently, YH10 in the space group of both and has been predicted to be a room-temperature superconductor under very high pressure[51,52].

Though many works about superconductivity of yttrium-hydrogen system under pressure have been implicated, very few explorations on the topological electronic properties of hydrides have been reported so far[53-55]. In this work, we predicted that YH3 in the space group of Pc1 at ambient pressure is a node-line semimetal when spin-orbit coupling (SOC) is ignored. Especially, the YH3 system we studied has extremely clean electronic structures near the Fermi level; i.e., there are no other pockets. The energy of the crossing points along the nodal loop varies within a very small energy range, from around −5 to 35 meV. Therefore, this nodal loop is very “flat” in the energy and momentum space, which makes YH3 a perfect model system for NLSMs. In general, NLSMs without SOC will transform into either insulators, DSMs, WSMs or even double NLSMs when SOC is considered, which is much related to the symmetries in the corresponding systems[56]. While in this work, when SOC is included, the three nodal lines around Γ point will be gapped out with a small gap (≈5 meV), making YH3 a topological insulator with Z2 = (1,000). Nevertheless, further calculation shows that the gap induced by SOC along the nodal ring is very small (about 5 meV), which indicates that the effect of SOC is negligible and the characteristic of the nodal ring can be preserved.

Methods

Calculations of the band structures are performed using the full-potential linearized augmented plane-wave (FP-LAPW) method[57,58] implemented in the WIEN2k[59] package. We use 13 × 13 × 11 k-mesh for the BZ sampling and −7 for the plane wave cut-off parameter RK for the electronic structure calculation, where the R is the minimum muffin-tin radius and K is the plane-wave vector cut-off parameter. SOC is taken into consideration by a second-variation method[60]. The tight-binding models are constructed with the maximally localized Wannier functions (MLWFs) method[61-63], the corresponding hopping parameters are determined from the projections of the bulk Bloch wave functions. The projected surface states are calculated using surface Green’s function in the semi-infinite system[64,65].

Results and Discussions

The crystal structure of YH3

Historically, three different phases of YH3 have been reported to be the ground state potentially. Two of them are experimentally favoured with trigonal Pc1 and hexagonal P63cm symmetry[66-71], while the third candidate is in the space group of P63 which was predicted theoretically[72]. It seems that the ground state of YH3 at ambient pressure is still under debate because the three candidates have very tiny total energy difference (0.001 eV/f.u.) from theoretical point of view. First, the hexagonal P63 structure is only theoretically proposed and seems to disagree with the neutron-diffraction results[73,74]. Second, later neutron-diffraction experiments[71,75] identify that the Pc1 structure is stable from the ambient pressure up to 12 GPa. Thus, in the following, we only focus on the phase of YH3 with the Pc1 symmetry.

The crystal structure and corresponding BZ of Pc1 (Space Group No. 165) YH3 is shown in Fig. 1(a,b) respectively. We use the experimental lattice parameters from literature[71] in our calculations, which are listed in Table 1.

Figure 1

(a) Crystal structure of YH3 at the ambient pressure with Pc1 symmetry. H1, H2 and H3 atoms occupy the 2a , 4d and 12 g (0.348, 0.025, 0.093) sites, respectively, while Y atoms lie at the 6f sites. (b) The corresponding BZ and its projection onto the (010) direction. The red ring on the shadow plane surrounding the Γ point represents the node-line structure in the BZ.

Table 1

The experimental lattice parameters of YH3 with Pc1 symmetry[71].

Phase	a = b	c	α = β	γ
P3¯c1	6.359 (Å)	6.607 (Å)	90°	120°

Band structures without SOC and the corresponding model analysis

From the band structures of Pc1 YH3 without SOC shown in Fig. 2(a), we can find three Dirac crossings composed of the conduction band minimum (CBM) and the valence band maximum (VBM) near Γ along the high-symmetry path in the BZ at the first sight. Detailed first-principle calculations indicates that Dirac crossings lying in the plane mΓ– is protected by the G symmetry, while Dirac crossing along K → Γ is not symmetry-protected (both the irreducible representations of CBM and VBM are Γ2) i.e., hybridization between the CBM and the VBM will open a gap in-between. Further calculations indicates that the CBM and the VBM are contributed mainly by Y-d and H1-s orbits, respectively, as the fat-band showed in Fig. 2(a). The band inversion (not caused by SOC) of Y-d and H1-s at Γ point leads to a gap of 0.302 eV.

Figure 2

(a) Corresponding fat-band structure of YH3 in the space group of Pc1 without SOC. The bands between CBM and VBM is gapped along K → Γ as the irreducible representations have showed. (b) Corresponding band structure of YH3 in the space group of Pc1 with SOC.

To give more insights of the nodal lines surrounding the Γ point, we have established an effective Hamiltonian by k · p method. Taking the crystal symmetry and TRS into consideration, the effective Hamiltonian can be written as follows:Here, the τ and τ are Pauli matrices, τ0 is a 2 × 2 unit matrix. This system has both TRS and IS, thus, the component of τ must be zero[11]. We can obtain the eigenvalues of the two-level system by diagonalizing the 2 × 2 effective Hamiltonian and the results are . Band crossings of the nodal lines will occur when g = g = 0. It’s clear that gives us MB > 0 and MC > 0. We find that MB > 0 and MC > 0 are exactly the condition of band inversion. Furthermore, confines the band crossings of the node-line in the k = 0, ± planes, i.e., there are three nodal rings surrounding the Γ point and lying in mirror-invariant planes mΓ–, as shown with the shadow sector in Fig. 1(b). It’s obvious that these three nodal rings are related to each other by R3 symmetry.

Band structures with SOC and the corresponding model analysis

When SOC is considered, band crossings of the three nodal lines will disappear, as the corresponding band structure with SOC shown in Fig. 2(b). We will further explain the above-mentioned phenomenon in the following. Taking SOC into account, spin and orbital angular momentum are coupled together, which generates a group of new eigenstates with certain total angular momentum quantum numbers. Then we can mark these new eigenstates of the CBM and VBM as , and . Here subscripts s and d denote corresponding orbits consisting of the new eigenstates and the superscripts ± represent the parities of corresponding eigenstates, respectively.

According to the analysis of irreducible representations and projected orbits, the CBM and VBM at the Γ point (denoted as and ) are mainly composed of and basis, respectively. If we arrange the 4 basis in the order of , , , , and then take the time-reversal and D3 point-group symmetries at the Γ point into consideration, we can give the character table of Γ matrices and the polynomials of momentum for this system as shown in Table 2.

Table 2

The character table for the Pc1 YH3.

Γ	Representation	T/P	k→
Γ1	Γ˜1+/A_1g	+	1,kx2+ky2,kz2
Γ2	Γ˜2−/A_2u	−	k_z,kz3,3kx2k_y−ky3
Γ3	Γ˜1−/A_1u	−	kx3−3k_xky2
{Γ4, Γ5}	Γ˜3−/E_u	−	{k_x,k_y},{kx3+k_xky2,kx2k_y+ky3}

From Table 2, our model Hamiltonian yields aswhich describes the dispersion of the CBM and VBM around the Γ point. Here we use the following Γ matrices:which satisfy the Clifford algebra {Γ, Γ} = 2δ. While the other ten Γ matrices are given by . Presence of both TRS and IS will forbid the existence of these ten Γ terms in our model Hamiltonian. In Equation (2),  − ,  − ,  +  − ,  −  −  and  +  with k− = k − ik.

From the model Hamiltonian given in Equation (2) together with the band structure shown in Fig. 1(b), we can draw some conclusions as the following. First of all, will break the particle-hole symmetry for the CBM and VBM around the Γ point. Secondly, D1 in will lead to a gap at the Γ point. Thirdly, to reproduce band inversion, we need that . More importantly, the dispersions of the model Hamiltonian given by Equation (2) are and both dispersions are doubly degenerate. As a result, a band crossing of this model requires f1 = f2 = f3 = f4 = f5 = 0. Several discrete points near the Γ point may satisfy the above-mentioned conditions. For example, when stands. It means that we may find some Dirac crossings at the first sight. However, on one hand, is a very rigorous condition and can not be obtained without other symmetries. More importantly, on the other hand, we can explain that points lying in the plane mΓ– must induce a gap in the following. There are three generator operators existing in the nonsymmorphic space group Pc1, we sign the three-fold rotation axis around z axis, the glide plane located at x = 0, and the inversion symmetry as R3, G and P, respectively. The operation of G acts in both the real space (x, y, z) and the spin space as

Similar with the analysis in the work by Fang et al.[32], we can easily find

On the mirror invariant plane k = 0, the bands can be labeled by its G eigenvalues. When SOC is considered, we havethe minus sigh is because includes a 2Π rotation in the spin space, which gives a −1 for a spin − system. So the eigenvalue of G is ± . The existence of both P and T will ensure all bands locally degenerated at every point in the BZ when SOC is considered, and the degenerated bands are related to each other by P  T. Suppose at (0, k, k), a Bloch function is an eigenstate of G with eigenvalue , then its degenerate partner under G,It means that the degenerated bands on the k = 0 plane have opposite G eigenvalues, and two sets of such doublet bands generally anticross, i.e., the bands with the same G eigenvalues hybridize and avoid crossing. In other words, nodal lines near the Γ point in the k = 0 plane (without SOC) will disappear in the whole BZ when SOC is considered, and G symmetry plays the key role of prohibiting the band crossing between CBM and VBM, even though the gap is very small (≈5 meV) as the k · p Hamiltonian in Equation (2) has showed. As a result, this node-line semimetal will become an insulator when SOC is considered. With both TRS and IS in this system, we can easily calculate the Z2 index by multiplying all the parities of all the occupied bands at all time-reversal-invariant momenta (TRIMs) using the method by Fu and Kane[76]. The results are shown in Table 3, which indicates the Pc1 YH3 is a strong TI with Z2 = (1,000) when SOC is taken into consideration. Nevertheless, the effect of the SOC is negligible because the H atom is small and the Y atom is also not very heavy; there is only one d electron in the Y atom. We have performed calculations with SOC and found that the SOC induced gap along the nodal ring is very small (about 5 meV), which indicates that the effect of SOC could be ignored and the characteristic of the nodal ring can be preserved.

Table 3

The parities product of all the occupied bands at the eight TRIMs for the Pc1 phase of YH3.

TRIM	Γ	3 M	A	3 L;	Total
Parity	+	−	−	−;	−

Surface states with and without SOC

Exotic topological surface states are an important property to identify various topological phases. Based on the tight-binding model constructed with WLWFs and surface Green function methods[64,65,77], we have calculated the 〈010〉 surface states of the Pc1 YH3 without SOC and the 〈001〉 surface states with SOC, as shown in Fig. 3(a,b), respectively. In particular, as the nice picture shows in Fig. 3(a), topological protected surface states without SOC signed with a bright curve connects the nodal points across the boundary of BZ. When SOC is included, the gap along Γ → A is so small (≈5 meV) that we may consider the CBM and VBM are nearly touched in the bulk band structure, this phenomenon can be proved that we can find the corresponding 〈001〉 surface states connecting the psudo-touch points from Fig. 3(b). We think these topological protected surface states may be helpful to identify the real ground state of YH3 from those three candidates. For example, we propose that angle-resolved photoemission spectroscopy (ARPES) technique can be used to investigate these surface states of this node-line semimetal candidate. If a bright surface state could be found in the 〈010〉 direction, and two parabolic bright curves with negative mass touching at the point could be found in the 〈001〉 direction, then the YH3 sample should be in the Pc1 symmetry.

Figure 3

(a) The surface states without SOC of Pc1 YH3 terminated with H atoms in the 〈010〉 direction, (b) The surface states with SOC of Pc1 YH3 terminated with H atoms in the 〈001〉 direction.

Conculsion

In conclusion, based on first-principles calculations and effective model analysis, we propose that the Pc1 YH3 is a nonsymmorphic symmetry protected node-line semimetal when SOC is ignored. This system has very clean electronic structures, there are no other pockets except the ones composing the node line near the Fermi level. The energy of the crossing points along the nodal loop has a very small energy range, from around −5 to 35 meV. Therefore, this nodal loop is very “flat” in the energy and momentum space, which makes YH3 a perfect system for model analysis. There are three node-lines related to each other by R3 symmetry locating on three planes (signed as mΓ–) surrounding the Γ point. While SOC is taken into consideration, the band crossings consisting of these node-lines will be gapped, and the Pc1 YH3 becomes a strong topological insulator with Z2 indices (1,000). At last, we have calculated the surface states of this system to verify its topological properties. We think our predictions should be helpful to identify the real ground state of YH3 in experiments in the future.