Giant Zeeman-type spin splitting of free electron/hole states on quasi-2D perovskite niobates: a theoretical prediction

We study the spin-orbit interaction of two-dimensional electron/hole gas (2DEGs/2DHGs) on quasi-2D potassium niobates (KNs) via first-principles calculations. The strong surface polarity changes the free surface states from 2DEGs to 2DHGs. The in-plane dipole maintained on 2D models leads to giant Zeeman-type spin splitting, as high as 566 meV for the (001)c facet KN and 1.21 eV for the (111)c facet KN. The thickness-dependent Zeeman-type spin splitting shows a linear relation with respect to 1/r, while the corresponding in-plane polarization quantum has a linear relation of 1/(2^0.5)with respect to a decrease in thickness. Interestingly, the 2DHGs with molecular-like orbital character is solely constituted by O 2p states, showing logic switchable behavior at extremely thin samples with enormous Zeeman-type splitting that can switch between insulator and conductor by opposite spin polarization.

Here, we first propose an intrinsic Zeeman-type spin splitting based on the 2DEGs/2DHGs system of a quasi-2D perovskite niobate in which its inversion symmetry is broken by maintaining an additional dipole in its 2D plane, as depicted in Fig. 1a.
The formulas of K n+1 Nb n O 3n+1 and K n+1 Nb n O 3n+3 are theoretical 2D models used in the ab initio calculation, with facets of (001) c and (111) c (oriented by 0° and 60° tilted oxygen octahedrons, respectively). 2D crystal models are relaxed through geometry optimization according to system energy, force and stress by density functional theory (DFT, implemented in the Cambridge Serial Total Energy Package 12 ) with the Perdew, Burke, and Ernzerhof (PBE) functional 13 . Both of the models are cut off by a K-O plane as the terminating surface. Hydrogen atoms are used to cap the dangling bonds of both surfaces to avoid the polar catastrophe 14 .
Dissimilar to traditional 2DEGs materials based on SrTiO 3 15 , in which the outer-shell valence is composed of 5s 2 for Sr and 3d 2 4s 2 for Ti, the niobate perovskite shows an odd number of valence electrons; for instance, KNbO 3 is composed of 4s 1 for K and 4d 4 5s 1 for Nb, giving rise to an out-of-plane effective magnetic momentum (B eff ) on both quasi-2D models due to the breaking of time-reversal symmetry, as shown in Fig. 1a. The quantum-confined spin-up electrons construct the out-of-plane magnetic order, transforming such a band insulator into a conductor and magnets, as sketched in Fig. 1a. After hydrogen capping of their surfaces, two different crystal facet configurations were evaluated with distinct origins of conductivity. The (001) c facet shows the traditional n-type free electron conduction (2DEGs), while the (111) c facet maintains abnormal p-type free hole (2DHGs) conduction. The reversed behavior is ascribed to the massive polar anisotropy of the varied facets, and more than 200 meV/Å 2 differences in surface energy (calculated by the PBE functional) can be obtained between (111) c and (001) c . Figure 1b,c shows the calculated electronic structures of the (001) c facet (0.80 nm sample) and (111) c facet (0.46 nm sample) oriented quasi-2D KN, respectively. Strong spin-polarized free electron or hole states (approximately 1 eV in energy) can be observed in both cases (marked by light gray areas). Compared with the confined free states around the Fermi level (E F ), subbands distant from E F show small energy splitting [details can be seen in the right panels of Fig. 1b,c, which quantitively give the distribution of the spin density of states (SDOS) of both models], indicating that the spin-polarized free states are mainly contributed by 2D confined Fermi gas.
Under Zeeman-type SOI, the Hamiltonian of an electron/hole in a 2D quantum well is described by , where H ĉ and H ẑ describe the contributions of the quantized well potential and the Zeeman-type SOI, respectively. Hence, their intrinsic energy dispersion can be written as , where M is the effective mass-related parameter and g factors define the power of the Zeeman-type SOI. Treating the H ẑ perturbatively 17 , = + ′ g g g k z 0 2 , g 2 0 ≈ is the Landè factor, and ′ g k 2 corresponds to the second-order perturbative term 18 . Considering that the origin of spin splitting is symmetry breaking in 2D KN, the in-plane polarization quantum (P x,y ) along orthogonal directions is calculated with various thicknesses.
= Ω P e R x,y i , where e is the electronic charge, R i is the i′th lattice vector, and Ω is the unit cell volume, this equation is given by the self-returned Hamiltonian of the change in polarization per unit volume for momentum k 19 . As shown in Fig. 2a,c, the values of P x,y for both the (100) c and (111) c facet 2D KN show an evolutionary www.nature.com/scientificreports www.nature.com/scientificreports/ tendency of r −1/2 with respect to a decrease in thickness and show almost identical values for orthogonal x and y directions of each thickness, which means that P x,y is isotropic in orthogonal directions of the 2D plane. The Zeeman splitting energy for an electron or hole can be described by eff . With a decrease in well thickness and a small B eff , g k g Taking into account the constant effective magnetic order of /2  and 3 /2  for constructed (001) c and (111) c 2D models, respectively, they would give rise to a 1/r relationship of g z with respect to the well thickness, which is verified by the result of the thickness dependence of E z in Fig. 2b,d. In the (001)c facet case, the 2DEGs is contributed by hybridization of the Nb 4d and K 3p states. Figure  3a shows the calculated conduction bands of samples with thicknesses varying from 0.40 nm to 4.00 nm.
Comparing the spin-polarized (dark and red dotted lines) and off-spin-polarized (cyan dotted lines) bands, the Zeeman-type SOI dramatically increases the absolute value of occupied energy (almost twice in energy) of free states, although this energy also increases with respect to an increase in dimension confinement (see off-spin-polarized subbands of samples with a reduction in thickness). To determine the contribution of the Zeeman-type spin splitting to such an enhanced free state, the splitting energy is obtained by calculating the difference of the CBM with the opposite spin. Conventionally, the CBM is constituted by the t 2g subband of Nb 4d 20 ; however, our results show an abnormal contribution by K 3p electrons to the first spin-down subband (marked by green balls in the subbands of the 0.40 nm sample in Fig. 3a), which is ascribed to the subband distortion due to the strong quantum confinement at the extremely thin sample and explains the abnormal E z of the 0.40 nm sample (Figs. 2b and 3a). The distribution detail of these abnormal K 3p states can be seen at the calculated partial density of states (PDOS) of the 0.40 nm sample (Fig. 3b). Furthermore, the PDOS of samples shows that the quantized states below E F (nondegenerate K 3p or Nb 4d) would become degenerate by increasing the 2D well width and would hybridize into the condensed 2DEGs. In addition, the E z of other samples (0.80-4.00 nm, shown as the red diamonds in the inset of Fig. 2b) present a 1/r relation with respect to a decrease in thickness. This relation is different from the evolutionary tendency of energy gap E g (red dotted line in the inset of Fig. 3a) without the perturbation of the Zeeman-type SOI, which has a quadratic relation with respect to a decrease in thickness. In a 2D square well, E g is given by Ε Ε = + is the reduced effective mass. The black dotted line is obtained by fitting a best-fit line to calculate E g values of samples with varied thicknesses, where the fitting result of m * ~ 0.9 is obtained for such an indirect band gap. Taking into account that the initialized spin momentum is /2 in (001) c facet system, the calculated maximum intrinsic E Z-e of 566 meV can be obtained for a 0.80 nm sample, twice the largest values reported for either the Rashba or Zeeman SOI systems (approximately 150 ~ 300 meV) 7,22 .
By tilting the angle of the oxygen octahedron to 60°, the free states switch from CBM to the VBM. The hole conductance is merely contributed by the O 2p states and shows 100% spin polarization around E F when the thickness is smaller than 1.15 nm (see Fig. 4a,b).
This phenomenon is of critical importance for realizing the application of spin-logic-circuit by using perovskite oxides, where a spin-up state leads to an insulator (presenting as |0), while the spin-down state gives rise www.nature.com/scientificreports www.nature.com/scientificreports/ to a conductor (presenting as |1). Compared with 2DEGs on the (001) c facet KN, the 2DHGs on the (111) c facet shows a stronger thickness dependence, regardless of the purely quantum confined 2DHGs that are insensitive to the spin polarization [see the cyan dotted line marking subbands of 0.46 nm to 4.60 nm samples in Fig. 4a, where the strong confined potential induced splitting of one heavy hole and two light hole subbands (0.46 nm sample), gradually becoming degenerate and decreasing to the E F (4.60 nm sample) at a high symmetric Γ point]. The size dependence of the Zeeman-type SOI leads to a much stronger enhancement of the occupied energy of free holes, almost 5 times that of the off-spin polarized counterpart, compared to the VBM of the 0.46 nm sample in Fig. 4a. The calculated E Z-h values of samples with thicknesses varying from 0.46 nm to 4.60 nm show a 1/r dependence with respect to a decrease in well width (Figs. 2d and 4a). The confined energy gap shows a quadratic relation with respect to a decrease in thickness; it is comparable to the fitted dotted line from each energy gap (m* ~ 0.6) without regard to spin polarization. The slight difference in the fitting m * between the (001) c and (111) c facet KN might be ascribed to the transition of the band gap from indirect to direct, which is induced by the 2D symmetry change (see Fig. 1b,c). The maximum E Z-h of ~1.21 eV is obtained for the 0.46 nm sample, which is 36 times that of the 4.60 nm sample (E Z-h ~ 33 meV). Furthermore, the B eff on energy dispersion is remarkably different for 2D holes than for electrons, leading to a different E Z-h at various k points (samples of thickness smaller than 1 nm in Fig. 4a). The different orbital angular momenta of L = 1 for the free holes (O 2p) and L = 2 for the free electrons (Nb 4d), in conjunction with the strong lattice distortion induced by strong quantum confinement, would contribute to the momentum dependent behavior of E Z-h It is similar to the case of k dependent g factor of holes in GaAs 2D quantum well, which the effective mass m ± for spin splitting subbands of heavy holes diverge further with an increase in B 23 .
In conclusion, our theoretical calculations show that the in-plane dipole can be maintained in both the (001) c and (111) c facet-oriented quasi-2D KN models. The constructed 2D model with time-reversal symmetry breaking along the z direction gives rise to an effective magnetic moment out of the 2D plane. In addition to the thickness dependent in-plane dipole, giant Zeeman-type spin splitting are observed. These free electrons and holes are susceptible to the Zeeman-type SOI, and the maximum E z of 566 meV for 2DEGs and 1.21 eV for 2DHGs can be obtained. Furthermore, these values all follow a 1/r rule with respect to the dimension reduction. Interestingly, taking advantage of giant Zeeman-type spin splitting, 100% spin-polarized O 2p states can be maintained near E F , which provides a new route toward the application of all-oxide spintronic devices in on-chip logic circuits within an extreme scale. Experimentally, we have developed a novel synthesis method that allows the topotactic growth of freestanding perovskite niobates into specific 2D facets, taking the first step toward the realization of spin computing in the near future.

Methods
The first-principles calculations are performed with the Cambridge Serial Total Energy Package (known as CASTEP). A bulk KNbO 3 crystal is established first and then cut off by a K-O plane as the terminating surface to build 2D crystal slabs with different crystal surfaces exposed. There should be a vacuum spacing layer with a value of 1 nm between neighboring repeat units of 2D crystal slabs. Hydrogen atoms are set to balance the surface www.nature.com/scientificreports www.nature.com/scientificreports/ charge of oxygen bonds 14 . Geometry optimization based on DFT with the Perdew, Burke, and Ernzerhof functional is performed to relax the crystal models 13 . The plane-wave cut off energy is set to 340 eV, and pseudopotentials that describe the ionic core are set to ultrasoft pseudopotentials. The k-point sampling grid for bulk KNbO 3 is 4 × 4 × 2, while it is 3 × 3 × 1 for 2D crystal slabs.
The in-plane polarization properties are calculated by the Perdew, Burke, and Ernzerhof functional available in commercial software Quantum ATK with PseudoDojo pseudopotential. The k-point sampling grid for 2D crystal slabs is 6 × 6 × 1 for 2D crystal slabs. All of the 2D slabs are fully relaxed.

Data availability
The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.