Next Article in Journal
The Role of Orbital Nesting in the Superconductivity of Iron-Based Superconductors
Next Article in Special Issue
Connecting Complex Electronic Pattern Formation to Critical Exponents
Previous Article in Journal
Stress and Strain Prediction of Zirconium Nitride under Oxygen Doping and Vacancy Introduction
Previous Article in Special Issue
Charge Order and Suppression of Superconductivity in HgBa2CuO4+d at High Pressures
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Character of Doped Holes in Nd1−xSrxNiO2

by
Tharathep Plienbumrung
1,2,
Michael Thobias Schmid
3,
Maria Daghofer
1,2 and
Andrzej M. Oleś
4,5,*
1
Institute for Functional Matter and Quantum Technologies, University of Stuttgart, Pfaffenwaldring 57, 70550 Stuttgart, Germany
2
Center for Integrated Quantum Science and Technology, University of Stuttgart, Pfaffenwaldring 57, 70550 Stuttgart, Germany
3
Waseda Research Institute for Science and Engineering, Waseda University, Shinjuku, Tokyo 169-8555, Japan
4
Institute of Theoretical Physics, Jagiellonian University, ul. Prof. S. Łojasiewicza 11, 30-348 Kraków, Poland
5
Max Planck Institute for Solid State Research, Heisenbergstrasse 1, 70569 Stuttgart, Germany
*
Author to whom correspondence should be addressed.
Condens. Matter 2021, 6(3), 33; https://doi.org/10.3390/condmat6030033
Submission received: 29 July 2021 / Revised: 16 August 2021 / Accepted: 17 August 2021 / Published: 19 August 2021
(This article belongs to the Special Issue Quantum Complex Matter from Charge Density Waves to Superstripes)

Abstract

:
We investigate charge distribution in the recently discovered high- T c superconductors, layered nickelates. With increasing value of charge-transfer energy, we observe the expected crossover from the cuprate to the local triplet regime upon hole doping. We find that the d p Coulomb interaction U d p makes Zhang-Rice singlets less favorable, while the amplitude of local triplets at Ni ions is enhanced. By investigating the effective two-band model with orbitals of x 2 y 2 and s symmetries we show that antiferromagnetic interactions dominate for electron doping. The screened interactions for the s band suggest the importance of rare-earth atoms in superconducting nickelates.

1. Introduction: Superconducting Infinite-Layered Nickelates Nd1−xSrxNiO2

The discovery of Bednorz and Muller [1] started an intense search for novel superconductors with high values of the critical temperature T c and gave very interesting open questions [2]. However, in spite of a tremendous effort in the theory, the mechanism responsible for the pairing in cuprates is still unknown [3]. Yet, this is one of the fundamental open problems in modern physics.
Perhaps less spectacular was the recent discovery of superconductivity in infinite-layered NdNiO 2 doped by Sr [4] as the values of T c are “only” close to 15 K [5]. Nevertheless, it gave a new impulse to the theory of high- T c superconductivity at large. To some extent, the nickelate superconductor family is rather similar to cuprate superconductors [6], as once again two-dimensional (2D) planes of transition metal and oxygen ions play a central role here [7]. With Ni + ions, one again has d 9 electronic configuration and similar lattice structure, but the apical oxygens are absent [8]. By following the same analysis for a NiO 2 layer as the one performed for a CuO 2 layer, we expect Ni 1 + to exhibit antiferromagnetic (AFM) order. One finds that the charge-transfer energy and the gap in nickelate are larger than those of cuprates [9,10]. The doped holes reside on oxygen sites in cuprates forming the Zhang-Rice singlets [11]. On the contrary, the doped holes will likely reside on Ni sites in doped Nd 1 x Sr x NiO 2 .
Usually, nickelates are investigated by evaluating their electronic structures by density functional theory (DFT) and comparing them with those of cuprates [6,7,8,12,13,14]. This approach has its well known weak points—here we consider instead the tight-binding charge-transfer model [10], similar to that for cuprates [15], and shall concentrate on the electronic properties of nickelates. We shall also address the question of magnetic exchange and the absence of AFM order in nickelates [16].

2. Charge-Transfer Model Revisited: Charge-Transfer Model for an NiO2 Plane

We introduced the multiband d p Hamiltonian for a NiO 2 plane [10] starting from a 2D Ni 4 O 8 2 × 2 cluster with periodic boundary conditions (PBCs). The basis set includes four orbitals per NiO 2 unit cell: two e g orbitals { 3 z 2 r 2 , x 2 y 2 } { z , z ¯ } at each Ni + ion and one bonding 2 p σ orbital (either 2 p x or 2 p y ) at each oxygen ion in the 2D plane,
H = H d p + H p p + H diag + H int d + H int p .
Here, the first two terms in the Hamiltonian (1) stand for the kinetic energy: H d p includes the d p hybridization t p d and H p p includes the inter-oxygen p p hopping t p p ,
H d p = { m α ; j ν } , σ t p d d ^ m α σ p ^ j ν σ + H . c . ,
H p p = { i μ ; j ν } , σ t p p p ^ i μ σ p ^ j ν σ + H . c . ,
where d ^ m α σ ( p ^ j ν σ ) is the creation operator of an electron at nickel site m (oxygen site i) in an orbital α { z , z ¯ } ( μ { x , y } ). Here, z and z ¯ stand for 3 z 2 r 2 and x 2 y 2 orbitals, while { x , y } stand for p x and p y orbitals. The elements { t p d , t p p } are accompanied by the phase factors, which follow from orbital phases [15].
The one-particle (level) energies are included in H diag , where we introduce the charge-transfer energy between d and p orbitals,
Δ = ε d ε p .
The model is completed by the Coulomb interactions in the d and p orbitals. For the d electrons,
H int d = m α U α n m α n m α + U 1 2 J H i n m 1 n m 2 2 J H i S m 1 · S m 2 + J H i d m 1 d m 1 d m 2 d m 2 ,
where α = z , z ¯ stands for z 2 and x 2 y 2 orbitals. d m α σ ( d m α σ ) is electron creation (annihilation) operator at site m in orbital α with spin σ . n m α σ is a number operator. S m α stands for spin operator in α orbital at site m. The two Kanamori parameters to describe the interactions between 3 d electrons are { U , J H } . U α is Coulomb repulsion element for α orbital. J H and U = U 2 J H are Hund’s exchange and interorbital Coulomb interaction [17]. Similar interactions with { U p , J H p } are written for p electrons.
Figure 1 shows different weight distributions for hole occupations at Ni and O sites for the parameters of Table 1. The Ni-O charge-transfer model has been studied via impurity [9] as well as lattice approach [10]. It has been established that the holes in undoped compounds remain within d x 2 y 2 orbitals. Thereby, we have assumed that Nd does not contribute to the electronic structure, and the system without Sr is a Mott insulator. Doping by Sr gives a doped hole, which tends to reside at nickel sites rather than at oxygen sites. The nickelate (cuprate) regime [9] is highlighted in blue (orange) in Figure 1. In the nickelate regime, the holes reside predominantly at Ni sites. This is the essential difference with cuprates where a doped hole (for hole doping) resides predominantly at oxygen and forms a Zhang-Rice singlet [11].
Including intersite Coulomb repulsion U d p enhances the hole occupancy at Ni sites [10] and shifts the doping crossover to lower values of the charge transfer energy Δ . All the on-site energies of the Ni( 3 d ) orbitals have been included in the Ni-O hybridization terms t p d [9]. Similar results were obtained for finite e g orbital splitting, where Δ z = 1 eV should be considered the upper limit.
The NiO 2 compound is a Mott–Hubbard insulator. It is then possible to replace the charge-transfer model by the d-only Hubbard model. The effective Ni–Ni hoppings can be derived from second-order perturbation theory [18]. The next question is on which d-orbitals the doped holes preferably reside?
The asymmetric distribution of holes suggests that one could replace the Ni-O model (1) with a Ni d-only model as oxygen p-orbitals become unimportant. The DFT calculation shows the band structure of NdNiO 2 that two bands are crossing Fermi level. In the orbital-resolved band structure, the lower band has d x 2 y 2 character and the upper band contains both Nd and Ni contributions. The large charge transfer energy as well as the presence of electron pocket at Γ are the two striking features of the nickelate compound. The empty 5d states of Nd are responsible for providing electron pocket. The empty 5 d states are below the Fermi level, in other words, these states provide the hole states into Ni band by the so-called ’self-doped’ effect [19]. Furthermore, the 5 d states was shown to be hybridized with Ni apical orbitals, i.e., 3 d z 2 and 4s orbitals [20] lead us to construct the effective two-band model consisting of Ni in-plane orbital, d x 2 y 2 , and Ni off-plane orbital, the modified s orbital. In this work, we present the character of doped holes in the realistic two-band model of NdNiO 2 compound.

3. Electronic Structure Calculations

The ab initio electronic structure calculations using DFT were performed with the Quantum Espresso code [21,22,23] using a plane-wave pseudopotential method, combining a projector augmented wave method [24] and a specific choice of pseudopotentials [25]. Within this work, we chose an energy cutoff of 600 eV and a Γ point centered Brillouin zone mesh of size 16 × 16 × 16 . For all calculations we used the same crystal structure published in [12].
The two-band model is derived by performing a Wannier projection onto the DFT band structure as implemented within the Wannier90 interface [26], including onsite energies and hopping parameters of each Wannier orbital. The projection was performed onto a Ni( x 2 y 2 ) and an apical Ni(s) Wannier orbital within an energy window ranging from 2.0 to 4.0 eV. The DFT band structure, represented here by the projected Wannier bands and its density of states (DOS), are shown in Figure 2. Panel Figure 2a shows the DFT band structure (black lines) and Wannier bands (red dashed lines). The right panel Figure 2b contains the resulting Wannier DOS around the Fermi energy for the bands of two symmetries: x 2 y 2 and s. It is clear that the fit cannot be perfect and the electron/hole estimation is wrong, but this plays a minor role fr the questions asked here and our result qualitatively matches previous studies [12,20].
The Wannier Hamiltonian is given within the real space notation and is of the form,
H ( R ) = i j α σ t i j α ( R ) c i α σ c j α σ .
The basis is ordered following the convention { d x 2 y 2 , s } . Explicit numerical values for the hopping parameters are given in Table 2 for the terms with leading contributions, i.e., terms larger than 0.001 (the other terms were neglected). Note that the relation H ( R ) = H T ( R ) holds for Wannier models. As a consequence, both terms (for distances R and R ) need to be included in further calculations. From H ( R ) , a tight-binding Hamiltonian can be constructed by applying a Fourier transformation of the form,
H i j ( k ) = R e i k R H i j ( R ) ,
where R describes the distances of the Wannier orbitals | i j | and is typically represented in terms of the lattice vectors.

4. Effective Two-Band Hamiltonian

Following the idea of neglecting oxygen orbitals, we construct an effective model containing only Ni( 3 d ) orbitals. The band structure calculation shows that only two bands contribute at Fermi level. Therefore, the two-band model of nickelates is capable of reproducing the physics of nickelates. It consists of x 2 y 2 orbital and the s orbital, which includes rare-earth 5 d states and Ni apical orbitals.
We consider the two-band Hamiltonian H = H kin + H int . The kinetic part is
H kin = i α σ ϵ α a i α σ a i α σ + i j α β σ t i j α β a i α σ a j β σ ,
while the interactions are given in a similar way to Equation (5). The orbitals x 2 y 2 and s have the diagonal energies ϵ α , being 0 and ϵ . The bands are constructed following Ref. [20]. The oxygen 2 p orbitals are included implicitly in x 2 y 2 ones. Note that here, one should not confuse x 2 y 2 symmetry with Ni( x 2 y 2 ) orbital. The former is an effective orbital including in-plane oxygen, while the latter is solely a Ni( 3 d ) orbital. The s orbital contains contributions from Nd( 5 d ), Ni( 3 d z 2 ) and Ni( 4 s ) orbitals. The symmetry of this hybrid orbital is the same as that of atomic s orbital.
Instead of Nd( 4 f ) electronic states, the empty Nd( 5 d ) orbitals are responsible for the striking electron pocket at the Γ point in the band structure [6]. The Nd atoms are originally located in the off-plane direction of Ni-O plane. In this effective model, the rare-earth atoms are included into Ni atoms via s orbital, indicating the model is three-dimensional and the Coulomb interaction as well as Hund’s coupling need to be screened. For instance, the two electrons sitting within s orbital could be located either at Nd or at Ni atom. Therefore, the Coulomb repulsion between these two electrons within the s orbital is reduced. We introduce a parameter α [ 0 , 1 ] to represent the reduced Coulomb interaction and Hund’s coupling, i.e., U 2 = α U 1 and J = α J H . α = 1 ( 0 ) stands for weak (strong) screening effect from rare-earth atoms, and we consider two parameter sets A and B, given in Table 3. We study the effective two-band model via Lanczos algorithm [27] on a 2 × 2 × 2 unit cell.
The undoped nickelate corresponds to quarter-filling, i.e., 8 electrons. The stoichimetric compound has d 9 configuration, where Ni( x 2 y 2 ) orbital is half-filled, and the weak hybridization of Ni and Nd causes Ni( x 2 y 2 ) orbital to be away from half-filling and creates the self-doping effect [19,28]. The Ni–Ni hopping integrals are obtained via fitting Wannier functions on DFT band calculation. We aim to address the question: where are the doped holes located in the two-band model?
In the absence of an electron hoppings Hamiltonian, the ground state is a trivial antiferromagnet where x 2 y 2 orbital is half-filled and the two canonical AFM phases, C-AFM and G-AFM, are degenerate. Including hopping elements to further neighbors leads to metallic behavior. In one-band Hubbard model, the metal–insulator transition occurs when U = 2 z t , where z = 4 is the number of neighbors in the 2D plane. Similarly, the two-band description can lead to partial orbital-selective Mott transition, where one band is insulating and the other one remains metallic. In the one-band version, the two relevant configurations are singly occupied or form double occupation within a single site. The quarter filling two-band version, however, contains several possible configurations where the lowest energy is still a single occupancy followed by a local triplet state with energy ( U 3 J H ) [17].

5. Results and Discussion

We begin with electronic density distribution on 2 × 2 × 2 clusters obtained by exact diagonalization with twisted boundary condition (TBC). The hopping parameters used here are given in Table 2. With PBC, one requires hopping integrals t i j α β to be scaled by the factor of 1/2 due to double-counting at the boundary. To avoid this additional factor in the hopping terms, we replace the PBC with TBC. Instead of having a constant phase t N + i = t i as in PBCs, the hopping terms at the boundary are modified by twisted angles ( ϕ x , ϕ y , ϕ z ) , giving t N + i = e i ϕ · r i t i ; for more details, see Refs. [29,30]. The PBC corresponds to ( 0 , 0 , 0 ) , while ( π , π , π ) is obtained for the anti-periodic boundary condition. In what follows, the observables are obtained by averaging over several twisted angles.
In Figure 3a,b, we show the undoped orbital-resolved densities, i.e., the electron occupancy within each orbital: x 2 y 2 (predominantly occupied) and s (it usually has decent amount of electrons). Finite occupancy of s orbital is expected due to self-doping effect. We further see that increasing of the parameter α causes the occupation on x 2 y 2 orbital to approach half-filling. Once the x 2 y 2 orbitals are almost half-filled, we show below that the G-AFM ground state is the only dominant magnetic ground state. Including screening effect into x 2 y 2 orbital by lowering its Coulomb interaction gives this orbital as slightly more favorable, see also Figure 3a,b.
By adding one hole, we see that both orbitals can be occupied by holes depending on how strong the screening effects of s orbital are, see Figure 3c,d. The screening effect on x 2 y 2 orbitals makes little changes in the one-hole case. An interesting scenario arises when two holes are added into the undoped Ni-O plane. The unscreened x 2 y 2 remains roughly the same, as s orbital is effectively favorable for holes at small α . When taking the screening effect into account, holes are effectively occupying the x 2 y 2 orbital regardless of the parameter α .
Theoretical studies based on Ni( e g ) bands [10,19,28] have suggested that high-spin S = 1 state Ni 2 + is favorable for hole doping. In contrast, RIXS measurements [31] show that two holes are residing mainly within Ni( x 2 y 2 ). Due to the limitation of the measurement, we cannot determine the occupation of the rare-earth 5 d states. The difference between the predictions of the theoretical model and the experiment arises from the rare-earth atoms. In the Ni( e g ) with d 8 configuration, the x 2 y 2 is half-filled and forming a high-spin state together with another 3 d orbital with energy ( U 3 J H ) [10]. This energy is significantly smaller than the energy U for a low-spin d 8 state.
On the contrary, the s orbital, in the hole picture, has lower energy than x 2 y 2 orbital and the quarter-filling of electrons corresponds to 3 4 -filling by holes via particle-hole transformation. In this case, the s orbital is filled by holes and leaves half-filled x 2 y 2 orbital—then holes reside mainly on x 2 y 2 orbital. This finding suggests that for hole doping, the one-band Hubbard model may be sufficient [32]. This intuitive picture of hole configuration requires the x 2 y 2 to be nearly half-filled, in other words, it suggests that the ground state is a strong antiferromagnet as in cuprates [11,33]. However, NMR experiments [34,35] report no observation of long-range AFM order in the RNiO 2 (R = Nd,La) down to 2 K.
According to the DFT+sicDMFT approach [28], the paramagnetic ground state has the lowest energy followed by C–AFM with energy difference about 20 meV/atom and G-AFM with energy 105 meV/atom. The C–AFM state has parallel spins along the c axis, while G–AFM phase has antiparallel spin alignment to its nearest neighbor in all directions. To address the magnetic order in this two-band model, the spin structure factor of undoped profile,
S ( k ) S ( k ) = i j e i k ( r i r j ) S i · S j ,
is shown in Figure 4. With this small cluster size we can observe only six relevant k -points in the Brillouin zone and only simple types of AFM order can be investigated. As a universal feature, the spin structure factor at ( π , π , π ) is enhanced as the x 2 y 2 orbital is close to forming a half-filled band; it indicates that G–AFM order is favored for electron doping [13], see Figure 4.
At small α , the spin structure factor (9) shows a competition between several magnetic ordered states, dominated by ( π , π , 0 ) and ( π , π , π ) spin correlations, which stand for C–AFM and G–AFM order. Other types of order, in particular A–AFM order, are not favored. In the recent x-ray scattering experiment [36], the existence of AFM correlations was confirmed, but the small cluster size prevents us from concluding whether AFM long-range order could be stable or not in this parameter regime.
The screening effect from rare-earth atoms [19,37] could be responsible for the competition of magnetic orders. At small α where the screening effect from rare-earth is strong, the electrons tend to form on-site triplets with energy ϵ + α ( U 1 3 J H ) [10], coexisting with singly occupied x 2 y 2 orbitals. Then, increasing α enhances the interactions on s orbitals and the on-site triplet energy surpasses the bandwidth of s orbital, i.e., ϵ + α ( U 1 3 J H ) > 2 z t ( z = 4 ). The largest hopping elements are the in-plane hoppings along x and y direction of x 2 y 2 orbital. Another transition then occurs when the on-site triplet energy overcome the bandwidth of x 2 y 2 orbital, becoming a Mott insulator.
Figure 5 shows the competition between low-spin and high-spin states in NiO 2 planes of the novel Ni-layered superconductors for decreasing number of holes, i.e., for electron doping. Here we use D α = 1 N i n i α n i α and T = 1 N i 3 4 + S i 1 · S i 2 , where N is the number of Ni sites. We show here that the double occupancy D α is enhanced with increasing number of holes in the NiO 2 plane. Increasing α increases the double occupancy D x 2 y 2 , while the double occupancies in s orbitals are almost absent. Thus, the model reduces in this regime to the one-band Hubbard model [32].
Simultaneously, the amplitude of high-spin states T at Ni sites is reduced in the regime of large α . Therefore, we conclude that electron doped materials have preferably low-spin configuration. On the contrary, hole doping may favor locally high-spin ( S = 1 ) states [10] instead of singlets ( S = 0 ) for the double occupancies of x 2 y 2 orbitals. Whether such local triplets could play a role in nickelate superconductivity is still an open question.

6. Summary and Conclusions

In summary, we have replaced the half-filled d p charge-transfer model by the effective two-band model at quarter-filling. The symmetries of each band are given by x 2 y 2 and s symmetry. The contributions of Nd and Ni atoms to s orbital lead to the reduction of Hund’s exchange and Coulomb repulsion, with its strength being scaled by the parameter α . The model shows the three competing phases: metal, orbital-selective, and Mott insulator. The Mott insulator is realized when the splitting between onsite triplet and singly occupied state is larger than the size of the bandwidth, similar to one-band Hubbard model. It is followed by the orbital-selective Mott insulator where the orbital splitting separates the x 2 y 2 and s bands from each other.
The nonmagnetic state is unlikely within this small unit cell and only AFM configuration can be realized. While the G–AFM phase is clearly dominating at large α , the competition between C–AFM and G–AFM is found at low α where the on-site triplet competes with AFM ground state, indicating the tendency toward AFM ordering. The holes are doped differently among the three phases. In metallic phase, the itinerant s orbital is favorable for holes. On the other hand, in the orbital-selective phase, the situation is slightly complicated while adding one hole favors s orbital, but adding two holes can effectively lead again to x 2 y 2 occupation, as in the Mott insulating phase.
When the screening effect on x 2 y 2 orbital is included, the stability of the orbital-selective phase is enhanced and the metallic phase vanishes. The two-band model at quarter-filling, therefore, connects the two controversial scenarios of which orbitals preferred by doped-holes via the parameter α . The parameter somehow represents the screening effect from rare-earth orbital, showing the importance of rare-earth atoms for the electronic structure of superconducting nickelate.
Finally, we remark that the present approach has several limitations, which could be seen as open questions at present, to be resolved in the future. First of all, electronic phase separation could occur in nickelates as in other superconductors, here we give the data for iron superconductors [38]. A very interesting open problem is how many electronic components contribute to the superconducting state and whether orbital selective superconductivity may occur in nickelates [20]. One could address these questions in a more satisfactory way only by extending the two-band model by spin-orbit coupling.

Author Contributions

T.P. performed numerical analysis of the two-band model. The two-band model was derived from the d p charge-transfer model by M.T.S. Conceptualization, M.D. and A.M.O.; writing—original draft preparation, A.M.O.; supervision, M.D.; writing—review and editing, A.M.O.; T.P., M.T.S., M.D. and A.M.O. selected the relevant information, analyzed the predictions of the theory versus experimental data, developed the interpretation of the numerical results, and wrote the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

T.P. acknowledges Development and Promotion of Science and Technology Talents Project (DPST). A. M. Oleś kindly acknowledges Narodowe Centrum Nauki (NCN, National Science Centre, Kraków, Poland) Project No. 2016/23/B/ST3/00839.

Acknowledgments

We would like to thank Andres Greco, Peter Horsch, Krzysztof Rościszewski, and George A. Sawatzky for many insightful discussions. A. M. Oleś is grateful for the Alexander von Humboldt Foundation Fellowship (Humboldt-Forschungspreis).

Conflicts of Interest

We declare no conflict of interest.

References

  1. Bednorz, J.G.; Muller, K.A. Possible high-Tc superconductivity in the Ba-La-Cu-O system. Z. Phys. B 1986, 64, 189–193. [Google Scholar] [CrossRef]
  2. Campi, G.; Bianconi, A.; Poccia, N.; Bianconi, G.; Barba, L.; Arrighetti, G.; Innocenti, D.; Karpinski, J.; Zhigadlo, N.D.; Kazakov, S.M.; et al. Inhomogeneity of charge-density-wave order and quenched disorder in a high-Tc superconductor. Nature 2015, 525, 359–362. [Google Scholar] [CrossRef]
  3. Keimer, B.; Kivelson, S.A.; Norman, M.R.; Uchida, S.; Zaanen, J. From quantum matter to high-temperature superconductivity in copper oxides. Nature 2015, 518, 179–186. [Google Scholar] [CrossRef] [PubMed]
  4. Li, D.; Lee, K.; Wang, B.Y.; Osada, M.; Crossley, S.; Lee, H.R.; Cui, Y.; Hikita, Y.; Hwang, H.Y. Superconductivity in an infinite-layer nickelate. Nature 2019, 572, 624–627. [Google Scholar] [CrossRef]
  5. Li, D.; Wang, B.Y.; Lee, K.; Harvey, S.P.; Osada, M.; Goodge, B.H.; Kourkoutis, L.F.; Hwang, H.Y. Superconductiviting Dome in Nd1-xSrxNiO2 Infinite Layer Films. Phys. Rev. Lett. 2020, 125, 027001. [Google Scholar] [CrossRef] [PubMed]
  6. Botana, A.S.; Norman, M.R. Simitarities and differences between LaNiO2 and CaCuO2 and implications for superconductivity. Phys. Rev. X 2020, 10, 011024. [Google Scholar]
  7. Jarlborg, T.; Bianconi, A. Multiple Electronic Components and Lifshitz Transitions by Oxygen Wires Formation in Layered Cuprates and Nickelates. Condens. Matter 2019, 4, 15. [Google Scholar] [CrossRef] [Green Version]
  8. Nomura, Y.; Arita, R. Superconductivity in infinite-layer nickelates. arXiv 2021, arXiv:2107.12923. [Google Scholar]
  9. Jiang, M.; Berciu, M.; Sawatzky, G.A. Critical Nature of the Ni Spin State in Doped NdNiO2. Phys. Rev. Lett. 2020, 124, 207004. [Google Scholar] [CrossRef] [PubMed]
  10. Plienbumrung, T.; Daghofer, M.; Oleś, A.M. Interplay between Zhang-Rice singlet and high-spin states in a model for doped NiO2 planes. Phys. Rev. B 2021, 103, 104513. [Google Scholar] [CrossRef]
  11. Zhang, F.C.; Rice, T.M. Effective Hamiltonian for the superconducting Cu Oxides. Phys. Rev. B 1988, 37, 3759–3761. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Gu, Y.; Zhu, S.; Wang, X.; Hu, J.; Chen, H. A substantial hybridization between correlated Ni–d orbital and itinerant electrons in infinite-layer nickelates. Commun. Phys. 2020, 3, 84. [Google Scholar] [CrossRef]
  13. Been, E.; Lee, W.-S.; Hwang, H.Y.; Cui, Y.; Zaanen, J.; Devereaux, T.; Moritz, B.; Jia, C. Electronic Structure Trends Across the Rare-Earth Series in Superconducting Infinite-Layer Nickelates. Phys. Rev. X 2021, 11, 011050. [Google Scholar]
  14. Zhang, R.; Lane, C.; Singh, B.; Nockeleinen, J.; Barbiellini, B.; Markiewicz, R.S.; Bansil, A.; Sun, J. Magnetic and f-electron effects in LaNiO2 and NdNiO2 nickelates with cuprate–like 3dx2-y2 band. Commun. Phys. 2021, 4, 118. [Google Scholar] [CrossRef]
  15. Arrigoni, E.; Aichhorn, M.; Daghofer, M.; Hanke, W. Phase diagram and single-particle spectrum of CuO2 high-Tc layers: Variational cluster approach to the three-band Hubbard model. New J. Phys. 2009, 11, 055066. [Google Scholar] [CrossRef] [Green Version]
  16. Nomura, Y.; Nomoto, T.; Hirayama, M.; Arita, R. Magnetic exchange coupling in cuprate–analog d9 nickelates. Phys. Rev. Res. 2020, 2, 043144. [Google Scholar] [CrossRef]
  17. Oleś, A.M. Antiferromagnetism and correlation of electrons in transition metals. Phys. Rev. B 1983, 28, 327–339. [Google Scholar] [CrossRef]
  18. Hu, L.-H.; Wu, C. Two-band model for magnetism and superconductivity in nickelates. Phys. Rev. Res. 2019, 1, 032046. [Google Scholar] [CrossRef] [Green Version]
  19. Zhang, G.-M.; Yang, Y.F.; Zhang, F.C. Self-doped Mott insulator for parent compounds of nickelate superconductors. Phys. Rev. B 2020, 101, 020501. [Google Scholar] [CrossRef] [Green Version]
  20. Adhikary, P.; Bandyopadhyay, S.; Das, T.; Dasgupta, I.; Saha-Dasgupta, T. Orbital selective superconductivity in a two-band model of infinite-layer nickelates. Phys. Rev. B 2020, 102, 100501. [Google Scholar] [CrossRef]
  21. Giannozzi, O.; Baroni, S.; Bonini, N.; Calandra, M.; Car, R.; Cavazzoni, C.; Ceresoli, D.; Chiarotti, G.L.; Cococcioni, M.; Daboet, I.; et al. QUANTUM ESPRESSO: A modular and open-source software project for quantum simulations of materials. J. Phys. Condens. Matter 2009, 21, 395502. [Google Scholar] [CrossRef] [PubMed]
  22. Giannozzi, O.; Andreussi, O.; Brumme, T.; Bunau, O.; Buongiorno Nardelli, M.; Calandra, M.; Car, R.; Cavazzoni, C.; Ceresoli, D.; Cococcioni, M.; et al. Advanced capabilities for materials modelling with Quantum ESPRESSO. J. Phys. Condens. Matter 2017, 29, 465901. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Giannozzi, O.; Baseggio, O.; Bonfà, P.; Brunato, D.; Car, R.; Carnimeo, I.; Cavazzoni, C.; de Cironcoli, S.; Delugas, P.; Ruffino, F.F.; et al. Quantum ESPRESSO toward the exascale. J. Chem. Phys. 2020, 152, 154105. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Lejaeghere, K.; Bihlmayer, G.; Björkman, T.; Blaha, P.; Blügel, S.; Blum, V.; Caliste, D.; Castelli, I.E.; Clark, S.J.; Dal Corso, A.; et al. Reproducibility in density functional theory calculations of solids. Science 2016, 351, aad3000. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Dal Corso, A. Pseudopotentials periodic table: From H to Pu. Comput. Mater. Sci. 2014, 95, 337–350. [Google Scholar] [CrossRef]
  26. Pizzi, G.; Vitale, V.; Arita, R.; Blügel, S.; Freimuth, F.; Géranton, G.; Gibertini, M.; Gresch, D.; Johnson, C.; Koretsune, T.; et al. Wannier90 as a community code: New features and applications. J. Phys. Condens. Matter 2020, 32, 165902. [Google Scholar] [CrossRef]
  27. Koch, E. The Lanczos Method. In The LDA+DMFT Approach to Strongly Correlated Materials; Pavarini, E., Koch, E., Vollhardt, D., Lichtenstein, A., Eds.; Forschungzentrum Jülich: Jülich, Germany, 2011. [Google Scholar]
  28. Lechermann, F. Doping-dependent character and possible magnetic ordering of NdNiO2. Phys. Rev. Mater. 2021, 5, 044803. [Google Scholar] [CrossRef]
  29. Shiroshi, M.; Wadati, M. Integrable boundary conditions for the one-dimensional Hubbard model. J. Phys. Soc. Jpn. 2012, 66, 2288–2301. [Google Scholar] [CrossRef] [Green Version]
  30. Poilblanc, D. Twisted boundary conditions in cluster calculations of the optical conductivity in the two-dimensional lattice models. Phys. Rev. B 1991, 44, 9562–9581. [Google Scholar] [CrossRef] [PubMed]
  31. Rossi, M.; Lu, H.; Nag, A.; Li, D.; Osada, M.; Lee, K.; Wang, B.Y.; Agrestini, S.; Garcia-Fernandez, M.; Chuang, Y.-D.; et al. Orbital and Spin Character of Doped Carriers in Infinite-Layer Nickelates. arXiv 2020, arXiv:2011.00595. [Google Scholar]
  32. Kitatani, M.; Si, L.; Janson, O.; Arita, R.; Zhong, Z.; Held, K. Nickelate superconductors—A renaissance of the one–band Hubbard model. NPJ Quantum Mater. 2020, 5, 59. [Google Scholar] [CrossRef]
  33. Anisimov, V.I.; Bukhvalov, D.; Rice, T.M. Electronic structure of possible nickelate analogs to the cuprates. Phys. Rev. B 1999, 59, 7901–7906. [Google Scholar] [CrossRef]
  34. Hayward, M.A.; Green, M.A.; Rosseinsky, M.J.; Sloan, J. Sodium Hydride as a Powerful Reducing Agent for Topotactic Oxide Deinteraction: Synthesis and Characterization of the Nickel(I) Oxide LaNiO2. J. Am. Chem. Soc. 1999, 121, 8843–8854. [Google Scholar] [CrossRef]
  35. Hayward, M.A.; Rosseinsky, M.J. Synthesis of the infinite layer Ni(I) phase NdNiO2+x by low temperature reduction of NdNiO3 with Sodium Hydride. Solid State Sci. 2003, 5, 839–850. [Google Scholar] [CrossRef]
  36. Lu, H.; Rossi, M.; Nag, A.; Osada, M.; Li, D.F.; Lee, K.; Wang, B.Y.; Garcia-Fernandez, M.; Agrestini, S.; Shen, Z.X.; et al. Magnetic excitations in infinite-layer nickelates. Science 2021, 373, 213–216. [Google Scholar] [CrossRef] [PubMed]
  37. Sakakibara, H.; Usui, H.; Suzuki, K.; Kotani, T.; Aoki, H.; Kuroki, K. Model Construction and a Possibility of Cupratelike Pairing in a New d9 Nickelate Superconductor (Nd,Sr)NiO2. Phys. Rev. Lett. 2020, 125, 077003. [Google Scholar] [CrossRef]
  38. Ricci, A.; Poccia, N.; Joseph, B.; Arrighetti, G.; Barba, L.; Plaisier, J.; Campi, G.; Mizuguchi, Y.; Takeya, H.; Takano, Y.; et al. Intrinsic phase separation in superconducting K0.8Fe1.6Se2 (Tc=31.8 K) single crystals. Superconduct. Sci. Technol. 2011, 24, 082002. [Google Scholar] [CrossRef]
Figure 1. Doping tendency of NiO 2 model as a function of crystal-field splitting Δ : (a) U d p = 0 ; (b) U d p = 1 eV. (Blue) undoped, N = 2 and N = 2 ; (Orange) one hole doped, N = 3 , N = 2 ; (Green) two holes doped, N = 3 , N = 3 in Ni 4 O 8 (Cu 4 O 8 ) cluster. The nickelate (cuprate) regime is highlighted in blue (orange). The model parameters are given in Table 1.
Figure 1. Doping tendency of NiO 2 model as a function of crystal-field splitting Δ : (a) U d p = 0 ; (b) U d p = 1 eV. (Blue) undoped, N = 2 and N = 2 ; (Orange) one hole doped, N = 3 , N = 2 ; (Green) two holes doped, N = 3 , N = 3 in Ni 4 O 8 (Cu 4 O 8 ) cluster. The nickelate (cuprate) regime is highlighted in blue (orange). The model parameters are given in Table 1.
Condensedmatter 06 00033 g001
Figure 2. DFT band structure and DOS of the Wannier Hamiltonian: (a) DFT band structure (black solid lines) compared with Wannier band structure (red dashed lines), and (b) Wannier DOS for the Ni( x 2 y 2 ) (red) and apical Ni(s) orbital (blue).
Figure 2. DFT band structure and DOS of the Wannier Hamiltonian: (a) DFT band structure (black solid lines) compared with Wannier band structure (red dashed lines), and (b) Wannier DOS for the Ni( x 2 y 2 ) (red) and apical Ni(s) orbital (blue).
Condensedmatter 06 00033 g002
Figure 3. Orbital-resolved electron densities as obtained for: (a) set A, and (b) set B, see Table 3. Strongly anisotropic electron distribution over the d x 2 y 2 and s orbitals is favored when the interactions are screened by α > 0.5 , as shown for: (c) set A, and (d) set B.
Figure 3. Orbital-resolved electron densities as obtained for: (a) set A, and (b) set B, see Table 3. Strongly anisotropic electron distribution over the d x 2 y 2 and s orbitals is favored when the interactions are screened by α > 0.5 , as shown for: (c) set A, and (d) set B.
Condensedmatter 06 00033 g003
Figure 4. Spin structure factor (9) of the three-dimensional two-band model, as obtained for: (a) set A, and (b) set B, see Table 3.
Figure 4. Spin structure factor (9) of the three-dimensional two-band model, as obtained for: (a) set A, and (b) set B, see Table 3.
Condensedmatter 06 00033 g004
Figure 5. Competition between low-spin and high-spin states as obtained for the parameters of: (a,c) set A, and (b,d) set B, see Table 3. Panels (a,b) show local double occupancy D α in each α = x 2 y 2 , s orbital, while panels (c,d) show local triplet states T.
Figure 5. Competition between low-spin and high-spin states as obtained for the parameters of: (a,c) set A, and (b,d) set B, see Table 3. Panels (a,b) show local double occupancy D α in each α = x 2 y 2 , s orbital, while panels (c,d) show local triplet states T.
Condensedmatter 06 00033 g005
Table 1. Parameters of the NiO 2 charge-transfer model (all in eV) used in [10].
Table 1. Parameters of the NiO 2 charge-transfer model (all in eV) used in [10].
t p d t p p Δ U z = U z ¯ J H U p J H p
1.300.557.08.41.24.40.8
Table 2. Hoppings parameters for the two-band model (all in eV) used in ED calculations.
Table 2. Hoppings parameters for the two-band model (all in eV) used in ED calculations.
(x,y,z) t ijk α β
(0,0,0) 0.2 0 0 1.2
(1,0,0) 0.380 0.050 0.050 0.031
(0,1,0) 0.380 0.050 0.050 0.031
(0,0,1) 0.039 0 0 0.076
(1,1,0) 0.088 0 0 0.111
(1,0,1) 0.001 0.009 0.009 0.252
(0,1,1) 0.001 0.009 0.009 0.252
(1,1,1) 0.015 0 0 0.056
Table 3. Parameters of the two-band model (all in eV) used in exact diagonalization calculations. The reference energies for the two bands of x 2 y 2 and s symmetry are 0 and ϵ , respectively.
Table 3. Parameters of the two-band model (all in eV) used in exact diagonalization calculations. The reference energies for the two bands of x 2 y 2 and s symmetry are 0 and ϵ , respectively.
Set ϵ t U 1 J H
A1.210.388.01.2
B1.210.384.00.6
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Plienbumrung, T.; Schmid, M.T.; Daghofer, M.; Oleś, A.M. Character of Doped Holes in Nd1−xSrxNiO2. Condens. Matter 2021, 6, 33. https://doi.org/10.3390/condmat6030033

AMA Style

Plienbumrung T, Schmid MT, Daghofer M, Oleś AM. Character of Doped Holes in Nd1−xSrxNiO2. Condensed Matter. 2021; 6(3):33. https://doi.org/10.3390/condmat6030033

Chicago/Turabian Style

Plienbumrung, Tharathep, Michael Thobias Schmid, Maria Daghofer, and Andrzej M. Oleś. 2021. "Character of Doped Holes in Nd1−xSrxNiO2" Condensed Matter 6, no. 3: 33. https://doi.org/10.3390/condmat6030033

APA Style

Plienbumrung, T., Schmid, M. T., Daghofer, M., & Oleś, A. M. (2021). Character of Doped Holes in Nd1−xSrxNiO2. Condensed Matter, 6(3), 33. https://doi.org/10.3390/condmat6030033

Article Metrics

Back to TopTop