ย้ถนำฐา๔

Skip to main content
SHARE
Publication

Magnetic phase diagram of a two-orbital model for bilayer nickelates with varying doping

Publication Type
Journal
Journal Name
Physical Review B
Publication Date
Page Number
195135
Volume
110

Motivated by the recently discovered high-๐‘‡๐‘ bilayer nickelate superconductor La3โขNi2โขO7, we comprehensively research a bilayer 2ร—2ร—2 cluster for different electronic densities ๐‘› by using the Lanczos method. We also employ the random-phase approximation to quantify the first magnetic instability with increasing Hubbard coupling strength, also varying ๐‘›. Based on the spin structure factor ๐‘†โก(๐‘ž), we have obtained a rich magnetic phase diagram in the plane defined by ๐‘› and ๐‘ˆ/๐‘Š, at fixed Hund coupling, where ๐‘ˆ is the Hubbard strength and ๐‘Š the bandwidth. We have observed numerous states, such as A-AFM, Stripes, G-AFM, and C-AFM. At half-filling, ๐‘›=2 (two electrons per Ni site, corresponding to ๐‘=16 electrons), the canonical superexchange interaction leads to a robust G-AFM state (๐œ‹,๐œ‹,๐œ‹) with antiferromagnetic couplings both in-plane and between layers. By increasing or decreasing electronic densities, ferromagnetic tendencies emerge from the โ€œhalf-emptyโ€ and โ€œhalf-fullโ€ mechanisms, leading to many other interesting magnetic tendencies. In addition, the spin-spin correlations become weaker both in the hole or electron doping regions compared with half-filling. At ๐‘›=1.5 (or ๐‘=12), density corresponding to La3โขNi2โขO7, we obtained the โ€œStripe 2โ€ ground state (antiferromagnetic coupling in one in-plane direction, ferromagnetic coupling in the other, and antiferromagnetic coupling along the ๐‘ง axis) in the 2ร—2ร—2 cluster. In addition, we obtained a much stronger AFM coupling along the ๐‘ง axis than the magnetic coupling in the ๐‘ฅโข๐‘ฆ plane. The random-phase approximation calculations with varying ๐‘› give very similar results as Lanczos, even though both techniques are based on quite different procedures. Additionally, a state with ๐‘ž/๐œ‹=(0.6,0.6,1) close to the E-phase wavevector is found in our RPA calculations by slightly reducing the filling to ๐‘›=1.25, possibly responsible for the E-phase SDW recently observed in experiments. Our predictions can be tested by chemically doping La3โขNi2โขO7.