{"id":"783b2c57-09c3-4cd4-b98b-39184351f851","arxiv_id":"1909.02557","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"An impurity exact-diagonalization model places NdNiO2's NiO2 layers at a singlet-triplet crossover, with a cuprate-like singlet hole state and superexchange about ten times smaller than in cuprates.","lead":"This paper models the nickel-oxide layers of the newly discovered nickelate superconductor NdNiO2. It argues that added holes can form a cuprate-like spinless singlet, but that the magnetic coupling is about ten times too weak to explain the superconductivity.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The predicted 1A1 singlet ground state hinges on Δ and t_pd estimates from an unpublished fit; the paper's own phase diagram shows the expected parameter box straddles the 1A1/3B1 boundary, so a modest shift in either parameter flips the central claim.","rationale":"I agree with the reader's weakest assumption. The paper's own figures and text make the parameter sensitivity explicit: Fig. 4 gives Δ_c=8.1 eV for t_pd=1.5 eV, and the quoted Δ range up to 9 eV already crosses the boundary; for t_pd=1.3 eV the boundary moves to lower Δ, making the 3B1 state likely over much of the quoted range. The reliance on an unpublished fit (Ref. 21) with no error bars is the key gap. I do not see a more fundamental flaw: the impurity ED is internally consistent, the full d8 multiplet is included, the O sublattice is treated properly, and the authors explicitly caveat the Nd-layer issue. The superexchange estimate is a standard order-of-magnitude argument and is not the shaky part. Therefore the reader's CONDITIONAL verdict is appropriate, and my stress-test does not change it.","tokens_in":8463,"tokens_out":9225,"duration_ms":94381,"concrete_test":"Perform an independent first-principles extraction of Δ (Ni d10 – O 2p charge-transfer energy) and t_pd for NdNiO2 using e.g. cRPA or DFT+U with a published, verifiable method, and report uncertainties. Then run the same impurity exact diagonalization for the best-estimate (Δ,t_pd) point and the endpoints of the uncertainty box. If any plausible endpoint yields a 3B1 ground state—or if the best point lies at Δ>8.1 eV for t_pd=1.5 eV—the abstract's statement that the singlet 1A1 state is favored is not supported for those parameters.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim—that hole-doped NdNiO2 has a 1A1 Zhang-Rice-singlet-like ground state rather than a 3B1 triplet—depends entirely on where NdNiO2 sits in the Fig. 2 phase diagram. The paper does not compute Δ and t_pd for NdNiO2 here; it relies on 'fits to ab-initio results [21–23]' (Ref. 21 is an unpublished manuscript by the authors) and quotes t_pd≈1.3–1.5 eV, Δ≈7–9 eV without uncertainty estimates. In Fig. 4, for the central values t_pd=1.5 eV and A=6.0 eV, the 1A1→3B1 transition occurs at Δ_c=8.1 eV. The quoted range extends to Δ=9 eV, which is already on the 3B1 side. For t_pd=1.3 eV, the boundary shifts to smaller Δ, so a large part of the quoted parameter range yields a 3B1 triplet ground state. The paper explicitly concedes (Summary) that 'the triplet 3B1 state is close in energy so that small changes in the parameters ... could stabilize the triplet state, making SC unlikely.' Thus the positive claim that large pd hybridization favors the singlet is a fine-tuned statement, not a robust finding. If the true parameters lie at the lower t_pd / higher Δ end, the central claim collapses to the weaker statement that NdNiO2 is close to a crossover—which the abstract does not convey. The negative superexchange claim is less affected because it depends on the same Δ but would survive a switch to 3B1.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the nature of a hole doped into a NiO2 layer of NdNiO2, modeled as a Ni 3d9 impurity embedded in an infinite O 2p6 square lattice. Using exact diagonalization with the full Ni 3d8 multiplet, O 2p orbitals, and pd and pp hoppings, the authors compute the ground-state symmetry of the added hole as a function of the charge-transfer energy Δ and the Racah parameter A. They find a 1A1 Zhang-Rice-singlet-like ground state in the parameter range they associate with NdNiO2, with a 3B1 triplet state close in energy, and they argue that the large Δ makes the superexchange about one order of magnitude smaller than in cuprates, casting doubt on magnon-mediated pairing. The paper also notes that the parent compound's metallicity and the possible presence of Nd bands remain unresolved, and it proposes a testable sensitivity of the ground state to strain, pressure, and chemical substitution.","tokens_in":8731,"tokens_out":9496,"duration_ms":99372,"significance":"If the parameter assignment is correct, the result is significant: it would shift the description of doped NdNiO2 from a simple S=1 Mott system to a cuprate-like Zhang-Rice-singlet picture while simultaneously arguing against spin-fluctuation pairing. The exact-diagonalization calculation is exact for the stated model, the full 3d8 multiplet is included, and the parameters are not fitted to the target ground-state symmetry, so the calculation is not circular. The superexchange estimate is based on a standard formula and is robust to the singlet-triplet ambiguity, since it follows from the larger Δ. However, the central 'critical region' claim is only as reliable as the t_pd and Δ estimates, which carry no quantified uncertainty, and one key source is an unpublished self-cited fit.","major_comments":[{"comment":"The central claim that NdNiO2 falls inside the 1A1 region is not robust to the stated parameter uncertainty. The quoted ranges t_pd ≈ 1.3–1.5 eV and Δ ≈ 7–9 eV are taken from fits to ab initio results, one of which (Ref. 21) is an unpublished manuscript by the authors, and no error bars or sensitivity analysis are given. For the central values t_pd = 1.5 eV and A = 6.0 eV, Fig. 4 shows the 1A1–3B1 transition at Δ = 8.1 eV, which lies inside the quoted Δ range; for t_pd = 1.3 eV, the boundary moves to smaller Δ, so a large portion of the quoted parameter box yields a 3B1 triplet. The abstract states that the layers 'fall inside' the critical region without this caveat. Please either provide quantitative uncertainty estimates and a sensitivity analysis, or reformulate the central claim as conditional on a specific parameter box.","section":"Model and Results (Figs. 2 and 4)"},{"comment":"There is an internal inconsistency in the stated relation between A and Udd. The text states Udd = A + 8B + 3C ≈ 6–7 eV, which with B = 0.15 eV and C = 0.58 eV gives A ≈ 3–4 eV; however, Fig. 3 and Fig. 4 are computed at A = 6.0 eV, giving Udd ≈ 8.9 eV. Since A is one of the axes of the phase diagram and the quoted Δ range is only meaningful relative to Udd, this inconsistency must be resolved before the central claim can be assessed.","section":"Model (parameter values)"},{"comment":"The abstract's unqualified statement that the NiO2 layers 'fall inside' the critical region is not supported by the paper's own caveat in the Summary that the triplet 3B1 state is close in energy and small parameter changes could stabilize it, making superconductivity unlikely. Please make the abstract and the main conclusion reflect the near-degeneracy and the parameter sensitivity, rather than presenting the singlet ground state as an established fact.","section":"Abstract and Summary"}],"minor_comments":[{"comment":"The word 'stochiometric' should be 'stoichiometric'.","section":"Introduction"},{"comment":"The phrase 'the doped states of a a Mott insulator' contains a duplicated article 'a'.","section":"Results"},{"comment":"The text says 'The three lines show how the boundary shifts with tpd', but Fig. 2 contains four lines, for t_pd = 1.1, 1.3, 1.5, and 1.7 eV.","section":"Fig. 2 caption and text"},{"comment":"The shaded ellipse is described only as 'the area we believe to be relevant'; please specify the parameter range it represents and the source of those bounds.","section":"Fig. 2"},{"comment":"Reference [38] has an unmatched parenthesis in the journal/year field: 'Phys. Rev. B 101, 075107 (2020, URL ...'.","section":"References"},{"comment":"The word 'undertanding' should be 'understanding'.","section":"Summary and discussion"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is timely and the exact-diagonalization calculation is sound, but the central claim rests on parameter estimates with no quantified uncertainty, including an unpublished self-cited fit. The internal inconsistency between the stated Udd/A relation and the value A = 6.0 eV used in the main figures should be resolved before acceptance. The editor may wish to ask the authors to make the content of Ref. 21 available or to replace it with a published reference."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe genuinely new piece here is the full d8 multiplet exact-diagonalization phase diagram for a hole doped into an NiO2 layer, showing the crossover between a 1A1 Zhang-Rice-like singlet and a 3B1 Hund triplet. The spectral-weight analysis is clean, the ED is exact for the stated model, and the superexchange estimate — J roughly an order of magnitude smaller than in cuprates because Δ≈7–9 eV — is a useful negative constraint on magnon glue. That last point survives even if the ground state turns out to be the triplet, so it is the more robust part of the paper.\n\nThe soft spot is exactly what the stress-test flags: the placement of NdNiO2 inside the crossover region rests on t_pd and Δ values from ab initio fits, one of which is an unpublished self-cited manuscript, and no uncertainties are given. With t_pd=1.5 eV and A=6 eV the 1A1→3B1 boundary sits at Δ_c=8.1 eV, which is inside the quoted 7–9 eV range; at t_pd=1.3 eV the boundary moves to lower Δ, putting a substantial part of the quoted range in the triplet region. The paper's own summary concedes the triplet is close in energy. So the abstract's phrase \"we show that the NiO2 layers fall inside a critical region\" is a bit stronger than the evidence warrants; the robust statement is that the layers are close to a crossover, not that the singlet is the likely ground state. That is a moderate caveat, not a fatal one — the calculation is honest and the crossover is explicitly falsifiable through pressure or strain.\n\nThe citation pattern looks fine; the unpublished fit is worth checking but self-citation there is not unreasonable. This paper deserves a serious referee. A referee should ask for parameter sensitivity across the full quoted ranges and for a clearer separation between what is computed and what is estimated from external fits.\n\nI would bring it to reading group and would cite it for the phase diagram and for the magnon-glue caveat. Send it to peer review.","headline":"An honest exact-diagonalization phase diagram placing NdNiO2 near a singlet–triplet crossover, but the abstract oversells the singlet side and the input parameters carry no error bars.","tokens_in":9371,"tokens_out":2546,"would_cite":true,"duration_ms":27053,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Hole-doped NdNiO2 forms a copper-oxide-like singlet, so spin fluctuations cannot be the pairing glue.","keywords":["NdNiO2","nickelate superconductivity","Zhang-Rice singlet","spin-state crossover","Mott insulator","superexchange","exact diagonalization","charge transfer energy"],"falsifier":"Resolve the local symmetry and spin of the hole added to a NiO$_2$ layer--for example by Ni $L$-edge x-ray absorption or resonant inelastic x-ray scattering on a doped NdNiO$_2$ film, or by an independent embedded-cluster calculation with directly fitted $\\Delta$ and $t_{pd}$ values. Observation of a robust $S=1$ ($^3B_1$) ground state with no low-lying singlet would falsify the paper's central claim.","tokens_in":8150,"feed_emoji":"⚛️","tokens_out":12725,"duration_ms":122588,"temperature":0.7,"pith_summary":"Doped NdNiO$_2$ was expected to put an $S=1$ nickel spin on each site, which would be hard to reconcile with superconductivity. The authors show, using an impurity model with full nickel multiplet structure, that the NiO$_2$ layer sits in a critical region of parameter space: strong nickel-oxygen hybridization makes the added hole form a copper-oxide-like $^1A_1$ singlet (a Zhang-Rice-like state) instead of a triplet. The same parameters give a superexchange about an order of magnitude smaller than in cuprates, so antiferromagnetic spin fluctuations are very unlikely to glue the Cooper pairs. The triplet $^3B_1$ state remains close in energy, making the ground state tunable by strain or chemistry. If correct, this reframes what minimal model describes the nickelate layer and rules out one leading pairing mechanism.","feed_headline":"Doped nickelate holes form copper-like singlets","feed_subtitle":"If right, the nickelate's weak superexchange rules out magnon pairing and points to a different mechanism.","key_machinery":"The central object is an impurity-model Hamiltonian $H = U_{dd} + T_{pd} + T_{pp} + \\Delta + U_{pp}$ for a Ni ion embedded in a square lattice of O $2p$ orbitals, solved by exact diagonalization. What carries the argument is the competition between the full Coulomb multiplet on Ni ($U_{dd}$), the charge-transfer energy $\\Delta$, and the Ni-O ($t_{pd}$) and O-O ($t_{pp}$) hoppings; the phase diagram in the $A$-$\\Delta$ plane locates NiO$_2$ at the boundary between the $^1A_1$ singlet and $^3B_1$ triplet ground states of the two-hole bound state. The symmetry-labeled spectral functions and wavefunction weights show the switch from $d^9L_{b_1}$ to $d^8(a_1b_1)$ across the crossover.","core_discovery":"Using exact diagonalization of a one- or two-hole impurity Hamiltonian in which a Ni ion sits in an infinite square lattice of oxygen $2p$ orbitals, and including the full $3d^8$ multiplet through Racah parameters, the authors find that for $\\Delta\\approx 7$--$9$ eV and $t_{pd}\\approx 1.3$--$1.5$ eV the ground state of the added hole has $^1A_1$ symmetry. Its dominant wavefunction component is $d^9L_{b_1}$--one hole in $d_{x^2-y^2}$ and one in the in-plane oxygen $x^2-y^2$ ligand orbital--locked into a singlet, exactly the structure of a Zhang-Rice singlet. This is surprising because NiO$_2$ is nominally a Mott insulator, for which a Hund's-rule $^3B_1$ triplet ($d_{x^2-y^2}$ plus $d_{3z^2-r^2}$, total $S=1$) would be expected. The two states cross at $\\Delta\\approx 8.1$ eV for representative parameters, so NdNiO$_2$ is predicted to sit at a spin-state crossover. The paper further evaluates the superexchange $J_{dd} = 4t_{pd}^4/(\\Delta^2 U_{dd}) + 8t_{pd}^4/(\\Delta^2 (U_{pp}+2\\Delta))$ and finds it roughly an order of magnitude smaller than in cuprates, which makes magnon-mediated pairing implausible.","pith_inferences":["A testable extension not pursued in the paper: if the system sits at a spin-state crossover, the superconducting $T_c$ should respond strongly to strain; pushing the ground state into the $^3B_1$ region should suppress superconductivity.","The same impurity-model calculation could be repeated for LaNiO$_2$; the paper's logic predicts that a lowered $\\Delta$ or changed $t_{pd}$ tips its ground state toward the triplet, which would explain the absence of superconductivity there.","The mobile, strongly oxygen-centered hole in the $^1A_1$ state suggests pairing mechanisms based on charge-transfer fluctuations or interlayer coupling to Nd carriers, routes the paper leaves open."],"forward_implications":["The doped NiO$_2$ layer should be modeled by a Zhang-Rice-like singlet state, not by an $S=1$ local-moment system.","Because the $^3B_1$ triplet is nearly degenerate, modest changes in lattice parameters from strain, pressure, or chemical substitution can switch the ground state, explaining sample-to-sample and film-versus-bulk variation.","The large $\\Delta$ makes the superexchange $J_{dd}$ roughly an order of magnitude smaller than in cuprates, so antiferromagnetic correlations are weak and magnon exchange is unlikely to be the pairing glue.","If magnons are not the glue in this nickelate, then either nickelate superconductivity needs a different mechanism or cuprate superconductivity is not primarily magnon-mediated."],"supporting_citations":[{"why":"Reports the discovery of superconductivity in doped NdNiO$_2$, the experimental fact this work addresses.","marker":"[6]"},{"why":"Introduces the charge-transfer classification scheme used to distinguish Mott and charge-transfer insulators.","marker":"[14]"},{"why":"Supplies the large charge-transfer energy of NiO$_2$ compared with cuprates, placing the system near the crossover.","marker":"[15]"},{"why":"Defines the Zhang-Rice singlet, whose $^1A_1$ symmetry the hole-doped state is claimed to share.","marker":"[16]"},{"why":"Gives the superexchange expression used to conclude that $J_{dd}$ is an order of magnitude smaller than in cuprates.","marker":"[17]"},{"why":"Describes the impurity-model plus exact-diagonalization method, here extended from CuO$_2$ to NiO$_2$.","marker":"[19]"},{"why":"Fits ab initio results to extract the Ni-O hybridization and related parameters used in the model.","marker":"[21]"},{"why":"Provides additional ab initio estimates of the low-energy parameters for the nickelate layer.","marker":"[22]"},{"why":"Another ab initio source for the charge-transfer and hybridization parameters entering the phase diagram.","marker":"[23]"}],"fun_headline_variants":["Nickelate holes form copper-like singlets, kill magnon glue","Spin-state crossover in nickelates rules out magnon pairing","Cuprate-style singlets in nickelates, but weak superexchange","Nickelates: singlet holes, weak exchange, no magnon glue"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole conclusion depends on the guessed values of two energy scales: how much energy it costs to move an electron from oxygen to nickel, and how strongly nickel and oxygen orbitals mix; if either real value falls outside the assumed range, the doped hole would be a triplet rather than a singlet and the claim would collapse.","fun_headline_variants_meta":{"raw":{"variants":["Nickelate holes form copper-like singlets, kill magnon glue","Spin-state crossover in nickelates rules out magnon pairing","Cuprate-style singlets in nickelates, but weak superexchange","Nickelates: singlet holes, weak exchange, no magnon glue"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000231,"raw_usage":{"total_tokens":1513,"prompt_tokens":1003,"completion_tokens":510,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":619,"completion_tokens_details":{"reasoning_tokens":433}},"tokens_in":619,"tokens_out":510,"duration_ms":5260,"temperature":1.0,"reasoning_tokens":433,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T04:47:09.146159+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Resolve the local symmetry and spin of the hole added to a NiO$_2$ layer--for example by Ni $L$-edge x-ray absorption or resonant inelastic x-ray scattering on a doped NdNiO$_2$ film, or by an independent embedded-cluster calculation with directly fitted $\\Delta$ and $t_{pd}$ values. Observation of a robust $S=1$ ($^3B_1$) ground state with no low-lying singlet would falsify the paper's central claim.","supporting_citations":[{"cited_title":"The critical nature of the Ni spin state in doped NdNiO$_2$","cited_arxiv_id":"1909.02557","evidence_quote":"Reports the discovery of superconductivity in doped NdNiO$_2$, the experimental fact this work addresses."},{"cited_title":"Infinite-layer fluoro-nickelates as $d^9$ model materials","cited_arxiv_id":"1911.11788","evidence_quote":"Introduces the charge-transfer classification scheme used to distinguish Mott and charge-transfer insulators."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the superexchange expression used to conclude that $J_{dd}$ is an order of magnitude smaller than in cuprates."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Describes the impurity-model plus exact-diagonalization method, here extended from CuO$_2$ to NiO$_2$."},{"cited_title":"This is a considerable energy, which is why the proper modelling of the O band is important","cited_arxiv_id":null,"evidence_quote":"Fits ab initio results to extract the Ni-O hybridization and related parameters used in the model."},{"cited_title":"Foyevtsova, I","cited_arxiv_id":null,"evidence_quote":"Provides additional ab initio estimates of the low-energy parameters for the nickelate layer."},{"cited_title":"Si et al., Phys","cited_arxiv_id":null,"evidence_quote":"Another ab initio source for the charge-transfer and hybridization parameters entering the phase diagram."}],"review_version":1}