Recognition: no theorem link
Non-magnetic insulating phase induced by Jahn-Teller effect in RNiO₃
Pith reviewed 2026-05-12 03:50 UTC · model grok-4.3
The pith
A Jahn-Teller-driven non-magnetic insulating phase emerges in RNiO3 when distortions on Ni2+ sites overcome Hund's exchange.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The paper identifies a distinct charge- and orbital-ordered insulating phase in RNiO3 where, upon charge disproportionation, the two eg electrons on Ni2+ sites occupy the same orbital with opposite spin when the associated Jahn-Teller energy exceeds Hund's exchange 3J. Self-consistent calculations in the full three-dimensional tight-binding model confirm the stability of this non-magnetic state, showing that magnetic order is not a prerequisite for the insulating phase.
What carries the argument
The Jahn-Teller electron-phonon coupling acting on charge-disproportionated Ni2+ sites, which favors same-orbital occupation by two electrons with opposite spins once its energy exceeds 3J.
If this is right
- The metal-insulator transition in small-bandwidth RNiO3 can occur in a non-magnetic state.
- The phase diagram includes a metastable non-magnetic charge-ordered insulator in addition to the known spin-polarized one.
- Local Jahn-Teller distortions on Ni2+ sites drive orbital ordering without spin polarization.
- Parameters extracted from DFT for LuNiO3 transfer to describe the full family of RNiO3 compounds.
Where Pith is reading between the lines
- Experiments under pressure or strain might stabilize this non-magnetic phase by tuning the relative strengths of JT and exchange energies.
- The model framework could be applied to predict similar JT-driven non-magnetic insulators in related perovskite systems.
- Absence of magnetic signals combined with specific orbital and lattice signatures would directly test the predicted phase.
Load-bearing premise
The on-site interaction strengths U and J together with the breathing and Jahn-Teller electron-phonon couplings derived from DFT calculations on LuNiO3 remain valid when inserted into the effective tight-binding model without important missing higher-order terms.
What would settle it
Detection of a non-magnetic charge-ordered state in RNiO3 with Jahn-Teller distortions but no magnetic ordering below the metal-insulator transition temperature, for example through combined X-ray diffraction and neutron scattering measurements.
Figures
read the original abstract
We propose a three-dimensional multi-orbital tight-binding model for rare-earth nickelates RNiO$_3$ that treats charge, spin, orbital, and lattice degrees of freedom on equal footing. All model parameters, including the on-site interactions $U$ and $J$ and the electron-phonon (el-ph) coupling to the breathing mode, are extracted from hybrid-functional DFT calculations for the small-bandwidth nickelate LuNiO$_3$. The model describes three competing insulating phases governed by the interplay of $U{-}3J$ and el-ph coupling to the breathing and Jahn--Teller (JT) modes. For large $U{-}3J$, the insulating state is stabilized by local JT distortions on high-spin Ni$^{3+}$ sites. For smaller $U{-}3J$, the system undergoes charge disproportionation, $2\mathrm{Ni}^{3+}\rightarrow\mathrm{Ni}^{2+}+\mathrm{Ni}^{4+}$, resulting in the spin-polarized charge-ordered state observed experimentally below the N\'eel temperature in small-bandwidth RNiO$_3$. When the JT energy on the Ni$^{2+}$ site exceeds Hund's exchange $3J$, a distinct charge- and orbital-ordered insulating phase emerges in which the two $e_g$-electrons occupy the same orbital with opposite spin. The stability of this phase is further confirmed by self-consistent calculations within the full three-dimensional tight-binding model. This newly predicted metastable state, characterized by JT distortions in a nonmagnetic charge-ordered RNiO$_3$ phase, shows that the onset of magnetic order is not required for the metal-insulator transition in RNiO$_3$.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript proposes a three-dimensional multi-orbital tight-binding model for RNiO3 that incorporates charge, spin, orbital, and lattice degrees of freedom on equal footing, with all parameters (U, J, breathing and Jahn-Teller electron-phonon couplings) extracted from hybrid-functional DFT on LuNiO3. It identifies three competing insulating phases: a JT-distorted high-spin Ni3+ phase for large U-3J, a charge-disproportionated spin-polarized state for smaller U-3J, and a new non-magnetic charge- and orbital-ordered phase that emerges when the JT energy on Ni2+ sites exceeds Hund's 3J (with both eg electrons occupying the same orbital but opposite spins). Self-consistent solutions of the full 3D model are used to confirm the stability of this new phase, leading to the claim that magnetic order is not required for the metal-insulator transition.
Significance. If the central result holds, the work supplies a unified, DFT-parameterized framework for the competing orders in small-bandwidth nickelates and predicts a metastable non-magnetic insulating state that could be tested experimentally. The explicit treatment of lattice couplings alongside electronic interactions and the use of independently computed parameters are strengths that distinguish the approach from purely phenomenological models.
major comments (2)
- [Abstract and model-construction paragraphs] Abstract and model-construction paragraphs: the claim that the new non-magnetic phase is stabilized when the JT energy on Ni2+ exceeds 3J rests on parameters (U, J, and el-ph couplings) transferred from hybrid-functional DFT performed on the uniform LuNiO3 structure. In the charge-disproportionated regime the local Ni2+ environment differs from the reference calculation, so the effective JT splitting may change; no explicit recomputation or sensitivity test of the JT > 3J condition inside the self-consistent disproportionated solution is reported.
- [Self-consistent 3D TB calculations section] Self-consistent 3D TB calculations section: while the abstract states that stability is confirmed by full three-dimensional solutions, no numerical values are given for the computed JT energy relative to 3J, the resulting gap, or the energy difference to the spin-polarized phase; without these quantities or error estimates from the DFT parameter extraction, the robustness of the phase boundary cannot be assessed.
minor comments (2)
- [Abstract] A compact table summarizing the three phases (order parameters, spin state, lattice distortion, and stability condition) would improve readability.
- [Throughout] The notation for the breathing and JT distortion amplitudes should be defined once at first use and used consistently thereafter.
Simulated Author's Rebuttal
We thank the referee for the careful and constructive review of our manuscript. The comments raise valid points about the transferability of DFT-derived parameters and the need for quantitative details in the self-consistent calculations. We have revised the manuscript to incorporate additional analysis and numerical results addressing these concerns.
read point-by-point responses
-
Referee: the claim that the new non-magnetic phase is stabilized when the JT energy on Ni2+ exceeds 3J rests on parameters (U, J, and el-ph couplings) transferred from hybrid-functional DFT performed on the uniform LuNiO3 structure. In the charge-disproportionated regime the local Ni2+ environment differs from the reference calculation, so the effective JT splitting may change; no explicit recomputation or sensitivity test of the JT > 3J condition inside the self-consistent disproportionated solution is reported.
Authors: We acknowledge that the parameters U, J, and the electron-phonon couplings were extracted from hybrid-functional DFT on the uniform LuNiO3 structure, as described in the methods. The Jahn-Teller coupling is treated as a local interaction whose strength is fixed by the reference calculation, while the actual lattice distortions are determined self-consistently within the model for each phase. To address the concern about possible changes in the effective splitting, we have added a sensitivity analysis in the revised manuscript: the JT coupling constant is varied by ±25% around the DFT value, and the non-magnetic phase remains stable whenever the resulting JT splitting exceeds 3J. A complete re-derivation of all parameters from a new hybrid-functional calculation on the disproportionated structure is computationally expensive and beyond the present scope, but the local character of the coupling supports the current approach. revision: partial
-
Referee: while the abstract states that stability is confirmed by full three-dimensional solutions, no numerical values are given for the computed JT energy relative to 3J, the resulting gap, or the energy difference to the spin-polarized phase; without these quantities or error estimates from the DFT parameter extraction, the robustness of the phase boundary cannot be assessed.
Authors: We agree that explicit numerical values are necessary to evaluate the robustness of the phase boundary. In the revised manuscript we have added the relevant quantities obtained from the self-consistent three-dimensional tight-binding solutions: the JT energy on the Ni2+ sites relative to 3J, the magnitude of the charge gap, and the total-energy difference between the non-magnetic and spin-polarized phases. We have also included error estimates propagated from the uncertainties in the DFT-extracted parameters. These results are now reported in the section describing the self-consistent calculations and confirm the stability of the non-magnetic insulating phase. revision: yes
Circularity Check
No circularity: parameters from independent DFT, phases solved from model equations
full rationale
The paper extracts all on-site U, J and electron-phonon couplings (breathing and JT) from hybrid-functional DFT performed on uniform LuNiO3, then inserts these fixed values into a three-dimensional multi-orbital tight-binding Hamiltonian. Competing insulating phases are located by self-consistent solution of that Hamiltonian; the condition that JT energy on Ni^{2+} exceeds 3J simply selects one of the model's eigenstates and is not imposed by redefinition or by fitting to the target phase. No self-citations, uniqueness theorems, or ansatze imported from prior author work appear as load-bearing steps. The derivation therefore remains self-contained against the external DFT benchmarks and does not reduce to its inputs by construction.
Axiom & Free-Parameter Ledger
free parameters (2)
- U and J (on-site interactions)
- el-ph coupling to breathing and JT modes
axioms (2)
- domain assumption The effective three-dimensional multi-orbital tight-binding Hamiltonian captures all relevant low-energy degrees of freedom.
- domain assumption The form of the electron-phonon coupling to breathing and Jahn-Teller modes is sufficient to describe the lattice instabilities.
Reference graph
Works this paper leans on
-
[1]
M. L. Medarde, Journal of Physics: Condensed Matter 9, 1679 (1997)
work page 1997
-
[2]
J. A. Alonso, J. L. García-Muñoz, M. T. Fernández-Díaz, M. A. G. Aranda, M. J. Martínez-Lope, and M. T. Ca- sais, Phys. Rev. Lett.82, 3871 (1999)
work page 1999
-
[3]
J. A. Alonso, M. J. Martínez-Lope, M. T. Casais, J. L. García-Muñoz, M. T. Fernández-Díaz, and M. A. G. Aranda, Phys. Rev. B64, 094102 (2001)
work page 2001
-
[4]
J. B. Torrance, P. Lacorre, A. I. Nazzal, E. J. Ansaldo, and C. Niedermayer, Phys. Rev. B45, 8209 (1992)
work page 1992
- [5]
-
[6]
M. K. Stewart, J. Liu, M. Kareev, J. Chakhalian, and D. N. Basov, Phys. Rev. Lett.107, 176401 (2011)
work page 2011
-
[7]
I. I. Mazin, D. I. Khomskii, R. Lengsdorf, J. A. Alonso, W. G. Marshall, R. M. Ibberson, A. Podlesnyak, M. J. Martínez-Lope, and M. M. Abd-Elmeguid, Phys. Rev. Lett.98, 176406 (2007)
work page 2007
- [8]
- [9]
-
[10]
M. J. Han and M. van Veenendaal, Phys. Rev. B85, 195102 (2012)
work page 2012
- [11]
-
[12]
T. Mizokawa, D. I. Khomskii, and G. A. Sawatzky, Phys. Rev. B63, 024403 (2000)
work page 2000
-
[13]
H. Park, A. J. Millis, and C. A. Marianetti, Phys. Rev. Lett.109, 156402 (2012)
work page 2012
-
[14]
T. Mizokawa, D. I. Khomskii, and G. A. Sawatzky, Phys. Rev. B61, 11263 (2000)
work page 2000
-
[15]
J. Varignon, M. N. Grisolia, J. Íñiguez, A. Barthélémy, and M. Bibes, npj Quantum Materials2, 21 (2017)
work page 2017
-
[16]
J. L. García-Muñoz, J. Rodríguez-Carvajal, and P. La- corre, Phys. Rev. B50, 978 (1994)
work page 1994
-
[17]
V. Scagnoli, U. Staub, M. Janousch, A. M. Mulders, M. Shi, G. I. Meijer, S. Rosenkranz, S. B. Wilkins, L. Paolasini, J. Karpinski, S. M. Kazakov, and S. W. Lovesey, Phys. Rev. B72, 155111 (2005)
work page 2005
-
[18]
O. E. Peil, M. Ferrero, and A. Georges, Phys. Rev. B 90, 045128 (2014)
work page 2014
- [19]
-
[20]
R. Chavez Zavaleta, S. Fomichev, G. Khaliullin, and M. Berciu, Phys. Rev. B104, 205111 (2021)
work page 2021
- [21]
-
[22]
Y. Lu, Z. Zhong, M. W. Haverkort, and P. Hansmann, Phys. Rev. B95, 195117 (2017)
work page 2017
- [23]
-
[24]
D. Lee, S. H. Baek, T. H. Kim, J.-G. Yoon, C. M. Folk- man, C. B. Eom, and T. W. Noh, Phys. Rev. B84, 125305 (2011)
work page 2011
-
[25]
J. Kanamori, Progress of Theoretical Physics30, 275 (1963), https://academic.oup.com/ptp/article- pdf/30/3/275/5278869/30-3-275.pdf
work page 1963
-
[26]
A. Georges, L. d. Medici, and J. Mravlje, Annual Review of Condensed Matter Physics4, 137 (2013)
work page 2013
-
[27]
M. Sotoudeh, S. Rajpurohit, P. Blöchl, D. Mierwaldt, J. Norpoth, V. Roddatis, S. Mildner, B. Kressdorf, B. Ifland, and C. Jooss, Phys. Rev. B95, 235150 (2017)
work page 2017
- [28]
-
[29]
S. Rajpurohit, V. Vennelakanti, and H. J. Ku- lik, The Journal of Physical Chemistry A (2024), 10.1021/acs.jpca.4c05046
- [30]
-
[31]
J. P. Perdew, M. Ernzerhof, and K. Burke, The Journal of Chemical Physics105, 9982 (1996), https://pubs.aip.org/aip/jcp/article- pdf/105/22/9982/19228856/9982_1_online.pdf
work page 1996
- [32]
-
[33]
S. Prosandeev, L. Bellaiche, and J. Íñiguez, Phys. Rev. B85, 214431 (2012)
work page 2012
-
[34]
S. Yamamoto and T. Fujiwara, Journal of the Physical Society of Japan71, 1226 (2002), https://doi.org/10.1143/JPSJ.71.1226. 11
- [35]
-
[36]
P. Seth, O. E. Peil, L. Pourovskii, M. Betzinger, C. Friedrich, O. Parcollet, S. Biermann, F. Aryasetiawan, and A. Georges, Phys. Rev. B96, 205139 (2017)
work page 2017
-
[37]
O. E. Peil, A. Hampel, C. Ederer, and A. Georges, Phys. Rev. B99, 245127 (2019)
work page 2019
-
[38]
G. G. Guzmán-Verri, R. T. Brierley, and P. B. Little- wood, Nature576, 429 (2019)
work page 2019
-
[39]
G. Giovannetti, S. Kumar, D. Khomskii, S. Picozzi, and J. van den Brink, Phys. Rev. Lett.103, 156401 (2009)
work page 2009
-
[40]
Photoinduced orbital polarization and Jahn-Teller effect in RNiO$_3$
S. Rajpurohit, S. R. U. Haque, A. M. Lindenberg, P. E. Blöchl, and T. Ogitsu, arXiv e-prints , arXiv:2604.18524 (2026), arXiv:2604.18524 [cond-mat.str-el]
work page internal anchor Pith review Pith/arXiv arXiv 2026
-
[41]
D. R. Baykusheva, H. Jang, A. A. Husain, S. Lee, S. F. R. TenHuisen, P. Zhou, S. Park, H. Kim, J.-K. Kim, H.-D. Kim, M. Kim, S.-Y. Park, P. Abbamonte, B. J. Kim, G. D. Gu, Y. Wang, and M. Mitrano, Phys. Rev. X12, 011013 (2022)
work page 2022
-
[42]
N. Tancogne-Dejean, M. A. Sentef, and A. Rubio, Phys. Rev. Lett.121, 097402 (2018)
work page 2018
-
[43]
O. Grånäs, I. Vaskivskyi, X. Wang, P. Thunström, S. Ghimire, R. Knut, J. Söderström, L. Kjellsson, D. Turenne, R. Y. Engel, M. Beye, J. Lu, D. J. Higley, A. H. Reid, W. Schlotter, G. Coslovich, M. Hoffmann, G. Kolesov, C. Schüßler-Langeheine, A. Styervoyedov, N. Tancogne-Dejean, M. A. Sentef, D. A. Reis, A. Rubio, S. S. P. Parkin, O. Karis, J.-E. Rubensso...
work page 2022
-
[44]
Kanamori, Journal of Applied Physics31, S14 (1960), number: 5 Reporter: Journal of Applied Physics
J. Kanamori, Journal of Applied Physics31, S14 (1960), number: 5 Reporter: Journal of Applied Physics. Appendix A: Hopping matrix element The hopping matrix elements alongx,yandzdirec- tions are defined as T x R,R′ =− 1 4 thop 3− √ 3 − √ 3 1 (A1) T y R,R′ =− 1 4 thop 3 + √ 3 + √ 3 1 (A2) T z R,R′ =t hop 3 0 0 3 (A3) We include three octahedral phonon mode...
work page 1960
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.