Persistent spin helix on a diamond surface

Using first-principles calculations, we determine the impact of hydroxyl (OH) termination on the spin–orbit coupling (SOC) coefficient and spin texture of diamond surfaces. The persistent spin helix (PSH) state is revealed in the OH-terminated diamond (111) surface. The value of coefficient PSH α PSH was 14.2 meV·Å, which is larger than that of the zinc-blende n-type quantum well of GaAs/AlGaAs. The large α PSH is caused by the asymmetric charge density distribution near the surface rather than the size of the atomic SOC. The large α PSH will ensure the wavelength of PSH (λ PSH ), which will aid in the miniaturization of spintronic devices.

S emiconductor surfaces are unstable owing to their dangling bond, resulting in chemical reactivity. To overcome this, surface termination is required. Previously, hydrogen-terminated silicon surfaces were considered for termination semiconductor surfaces by hydrogen, where a silicon substrate with a surface that has been chemically passivated has hydrogen atoms bound to it. 1) Terminating the silicon surface using hydrogen is the first step in upcoming functionalization devices considering that it is submissive to chemical modification. 2) Additionally, the surface termination of semiconductors such as Ge, 3) GaAs and InP, 4) as well as diamond surfaces, 5) have been studied.
Termination of the diamond surface strongly affects the surface properties of the diamond. Hydrogen-or oxygenterminated diamond surfaces have been used to build diamond devices. [6][7][8] In the case of a hydrogen-terminated diamond, a hydrogenated diamond increases its electrical conductivity. 9,10) In contrast, an oxygen-terminated diamond produces non-conductive electrical properties, shows positive electron affinity, and has a larger Schottky barrier height in metal/diamond interactions compared to a hydrogenated diamond. 11) Furthermore, the oxygen-terminated diamond (111) surface tends to be more reactive, and the roughness 12) remains on its surface. Therefore, the combination of hydrogen and oxygen (OH) termination may solve the oxygenated diamond roughness issues on the surface. 13) Surface termination affects the electronic properties of materials. A hydrogenated diamond surface possesses strong Rashba-type spin-orbit coupling (SOC), which is approximately 9.74 ± 0.1 meV 14) and 4.6-24.5 meV, 15) which is yielded from the highly asymmetric confinement potential. Therefore, a diamond surface with hydrogen termination offers enormous promise for the investigation and use of spin transport phenomena. 15) These findings guide us to study the effect of OH termination of a diamond surface on its electronic properties, especially in the SOC property.
Unlike the hydrogen-terminated (H-terminated) diamond case, which possesses Rashba-type SOC, this study demonstrated that the OH-terminated diamond yields persistent spin helix (PSH) SOC using first-principles density functional calculations. The PSH has a one-dimensional orientation of spin textures, which is predicted to induce a strongly enhanced spin relaxation time. 16) We found that the coupling constant that defines the SOC strength, α PSH , for an OHterminated diamond, is larger than the other α PSH values observed in the zinc-blende n-type quantum well (QW) of GaAs/AlGaAs, 17,18) which is promising as a material for miniaturized spintronic devices.
First-principles electronic structure calculations were performed based on the density functional theory using the OpenMX code. 19) We calculated the OH-terminated diamond (111) surface and focused on the top diamond surface of termination, as shown in Fig. 1(a). The OH-terminated diamond surface (111) was conducted using a slab model comprising 9 C-C (carbon and carbon) bonding layers, as shown in Fig. 1(b). The lattice constant of the diamond (111) surface hexagonal supercell a 111 equal to b 111 was set to 2.527 Å, and c 111 was 39.285 Å, which included the vacuum region to desist interactions between neighboring slabs. a 111 was calculated from a a 2 111 = , where a is 3.574 Å, the theoretical lattice constant of the diamond. 20) We applied the generalized gradient approximation 21) as the exchange-correlation functional, with a k-point grid of 20 × 20 × 1, to discretize the first Brillouin zone in the kspace, as shown in Fig. 1(c), and a cutoff energy of 500 Rydberg for the numerical integration as well as the solution of Poisson's equation. 22) Furthermore, we applied the effective screening medium approach 23) to eliminate the dipoledipole interaction between slab models. The bottom side, hydrogen and four carbon atoms were fixed during the geometrical optimization to maintain the bulk structure, until the total force working on each atom was less than 0.005 eV/ Å. The linear combination of multiple pseudo-atomic orbitals 24)-26) was used to expand the wave functions for each atom, defined as C6.0-s3p2d2, H7.0-s2p2d1 and O7.0-s3p2d2. For example, C6.0-s3p2d2 means the cutoff radius of the carbon atom is 6.0 Bohr, and three primitive orbitals for s and two primitive orbitals for the p and d components are used. We used the norm-conserving pseudopotentials including the 2s and 2p electrons for C and O and the 1s electron for H. Through the fully relativistic total-angularmomentum-dependent pseudopotentials 27) the SOC was incorporated in this computation, and the spin textures in kspace were derived using the k-space spin density matrix of the spinor wave function. 28,29) In our OH-terminated diamond (111) surface system, as shown in Fig. 1(b), the bond between the carbon atom and oxygen was 1.425 Å, which was slightly larger than that in the previous work (1.39 Å). 30) The bond between the oxygen and hydrogen atoms was 0.99 Å, which was slightly smaller than that in the previous study (1.01 Å). 30) The bonds between the carbon atoms were between 1.526 and 1.550 Å, which was consistent with the previous theoretical work (between 1.51 and 1.55 Å). 30) On the bottom side, the hydrogenated one, the calculated distance between the hydrogen and carbon atoms was 1.11 Å, which was consistent with that in the previous theoretical study work, 1.11 Å. 30) We calculated the energy band dispersion, as shown in Fig. 2(a). The purple and black lines indicate the energy bands of the diamond surface (111) with OH termination and using a supercell, respectively. The energy gap for these systems was 2.94 eV at the Γ point, direct band gap, for the OH-terminated diamond (111) surface system, and 4.4 eV for the bulk diamond (111) system. The OH termination yielded a narrower band gap owing to several bands within the bulk band gap. The occupied surface state was identified by the 2p orbital of an oxygen atom, based on the partial density of states, as shown in Fig. 2(b).
where α R satisfies α R = 2E R /k R , k x and k y are wave vector components in the x and y directions, and σ x and σ y are the Pauli matrix vectors. The Rashba parameter α R suggested by a theoretical work 32) can be written as where c is the velocity of light, f(z) is the surface wave function, and V(z) the crystal potential averaged for the xy plane. The light elements can produce strong SOC with a strong electric field and charge asymmetry as presented in the previous theoretical works on graphene 33) and SrTiO 3 (001). 34) In the OH-terminated diamond (111) surface system, the OH termination induced the in-plane electric field, resulting in electric polarization along the y-direction. There is a mirror plane at the yz plane. 35) To comprehend the origin of the spin texture, we considered the SOC of surface states based on the system's symmetry. Applying mirror reflection symmetry operation to this material resulting in transformation from (x, y, z) to (−x, y, z), the Pauli matrices σ x , σ y , σ z transforms to σ x , −σ y , −σ z . We should always remember that the symmetry operation should not alter the Hamiltonian (see supplementary data, Sect. 2). By keeping R −1 HR = H, the SOC Hamiltonian can be given as H SOC = α 1 k x σ z + α 2 k x σ y + α 3 k y σ x , where k x and k y are wave vectors in the x and y (b) (c) (a)  Fig 2(a). The energy difference (ΔE) and k difference (Δk) are shown to calculate the large spin splitting (α). The origin of the energy axis is set to the VBM. directions, respectively, σ x , σ y and σ z are Pauli matrices, and α 1 , α 2 , α 3 are coupling constants that define the SOC strength. Here, α 1 is characterized by the in-plane electric field E y , whereas α 2 and α 3 relate to the out-of-plane electric field E z that is owing to the surface effect. The expression of H SOC is the same as in the case of a ZnO(10 1 0) surface. 35) The SOC effect on the surface state of the OH-terminated diamond (111) surface is shown in Fig. 2(c). For the Γ-M direction, k x ≈ 0, and H SOC = α 3 k y σ x . Because the calculated bands are degenerated on the Γ-M line, α 3 ≈ 0. For the Γ-K′ direction, we found that the bands were splitting, as shown in Fig. 2(c). Along this direction, k y = 0, and hence, only the first and second terms remain in the spin-orbit Hamiltonian, α 1 k x σ z + α 2 k x σ y . Consequently, the expected spin value of the spin component mostly comprised 〈S y 〉 and 〈S z 〉 as shown in Fig. 3(b). The SOC strength parameter K a G-¢ can be evaluated in the same way as for the Rashba system: Through the band dispersion method, as shown in Fig. 2(c), we get K a G-¢ is 14.2 meV·Å. This is larger than other PSH coefficients observed in GaAs/AlGaAs, which were 3.5 to 4.9 meV·Å 17) and 2.77 meV·Å. 18) In addition, the α PSH of the OH-terminated diamond is significantly larger than the Rashba parameter α R of the H-terminated diamond. The calculated results of the H-terminated diamond are available in Sect. 3 of the supplementary material.
We calculated the spin texture projected to the k x -k y plane and analyzed the spin texture as shown in Fig. 3(a) in the kspace at two meV below the degenerate point. We found the unidirectional spin texture in the k-space at the band around the VBM at the Γ point. Contrary directions can be seen in the spin textures from the inner and outer Fermi arcs, as shown in Fig. 3(a). Except for the Fermi arc surrounding the degenerate points at θ = π/2 and θ = -π/2, the spin components 〈S x 〉 and 〈S z 〉 were nearly constant. Performing the global spin rotation U, 16) when α 1 is equal to α 2 , we could perform the global spin rotation π/4 around the x-axis 36) to the present spin Hamiltonian using  Fig. 3(a). The spin textures exhibited a quasi-one-dimensional orientation in the in-plane y-direction, as shown in Fig. 3(a), despite possessing an outof-plane component, as shown in Fig. 3(b), thereby indicating the formation of PSH. 16) The PSH state, having unidirectionally oriented spin components, yielded a long spin lifetime. 16) Yoshida et al. experimentally observed the OH-terminated diamond (111) surfaces by Fourier transform infrared spectroscopy. 13) Because H SOC is mostly affected by the electric field, the origin of the spin textures was further clarified by evaluating the electric polarization. As shown in Fig. 3(c), given that the spin-split surface state is strongly localized in the first two bilayers, a strong electric polarization is expected to occur in the C-C-OH bilayer, as shown in Fig. 3(d). Here, the first bilayer comprising H and O atoms has strong polarity. Moreover, the expected values of spin are dominated by the oxygen atom compared to the hydrogen atom. Furthermore, the oxygen atoms dominated the occupied surface state, as shown in the partial density of the state in Fig. 2(b).
In summary, we performed first-principles density functional calculations for the OH-terminated diamond (111) surface and demonstrated that the spin textures have a quasione-dimensional orientation. Furthermore, we found that the persistent spin helix occurred owing to the OH termination on the surface of the diamond (111). The coefficient PSH α PSH of the OH-terminated diamond in the Γ-K′ direction of 14.2 meV·Å was larger than that observed in the zinc-blende n-type QW structure of GaAs/AlGaAs, of 3.5 to 4.9 meV·Å 17) and 2.77 meV·Å. 18) These findings suggest that OH-terminated diamonds with in-plane electric polarization and mirror symmetry can be used to create a persistent spin helix state, which is a promising technique for spintronic devices.