{"id":"0b6dce00-7b86-44e8-bd89-78facd6636da","arxiv_id":"2606.07199","paper_version":1,"verdict":"UNVERDICTED","confidence":"LOW","novelty_score":4.0,"correctness_risk":"unknown","formal_verification":"none","parameter_count":1,"one_line_summary":"Dynamical singlets from 3z²-r² orbitals hybridized with x²-y² orbitals control the distinct pressure versus strain responses in bilayer nickelates.","lead":"Calculations within density functional theory plus cluster dynamical mean-field theory indicate that interlayer dynamical singlets formed by singly occupied 3z²-r² orbitals hybridized with planar x²-y² orbitals control much of the physics in bilayer Ruddlesden-Popper nickelates. The hybridization responds differently to hydrostatic pressure than to epitaxial strain, reproducing the observed difference between bulk crystals and thin films along with ARPES and transport data.","discovery_kind":"extension","skeptic_critique":{"model":"grok-4.3","headline":"Cluster DMFT setup may miss longer-range or multi-orbital effects that alter the pressure-strain response of the hybridization","rationale":"The reader's weakest assumption directly identifies the same methodological vulnerability; the full-text calculations would need explicit robustness checks (absent from the abstract) to remove it, so the UNVERDICTED status is appropriate.","tokens_in":1634,"tokens_out":278,"duration_ms":14678,"concrete_test":"Recompute the hybridization strength and spectral functions under the same pressure and strain conditions but with an enlarged cluster (e.g., 4-site bilayer instead of minimal) or extended orbital basis (including apical p_z); if the differential response reverses or vanishes, the load-bearing assumption fails.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim requires that interlayer dynamical singlets (singly occupied 3z²-r²) hybridized to x²-y² control the physics and produce opposite responses to hydrostatic pressure versus epitaxial strain. This depends on the specific cluster DMFT implementation and orbital projection accurately isolating that hybridization without significant contamination from apical oxygen, longer-range hoppings, or in-plane momentum dependence outside the cluster. If those contributions modify the effective interlayer coupling under the two strain protocols, the claimed dichotomy (and match to bulk vs film experiments) would not follow from the reported calculations.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.3","summary":"The manuscript uses density-functional theory plus cluster dynamical mean-field theory to argue that the electronic physics of bilayer Ruddlesden-Popper nickelates is dominated by interlayer dynamical singlets formed from singly occupied Ni 3z²-r² orbitals that hybridize with itinerant x²-y² states. This hybridization is shown to respond oppositely to hydrostatic pressure versus epitaxial strain, thereby explaining the experimental dichotomy between bulk crystals and thin films while also reproducing ARPES spectra and transport measurements.","tokens_in":1801,"tokens_out":474,"duration_ms":13184,"significance":"If the central claim is substantiated, the work supplies a concrete orbital-resolved mechanism for the pressure-strain dichotomy in these materials and demonstrates how cluster DMFT can isolate dynamical singlet physics. The explicit comparison to both ARPES and transport data, together with the parameter-free character of the reported strain response once U and J are fixed, would constitute a useful benchmark for future studies of nickelate heterostructures.","major_comments":[{"comment":"§3.2 and the computational-methods paragraph: the hybridization strength is extracted from a two-site bilayer cluster; the manuscript does not demonstrate that enlarging the cluster or adding longer-range in-plane hoppings leaves the opposite pressure-versus-strain response intact, which is load-bearing for the claimed dichotomy.","section":"§3.2"},{"comment":"Table 1 and Fig. 4: the reported ARPES and resistivity comparisons are shown only for the final parameter set; no sensitivity plot versus U or J is provided, leaving open whether the hybridization response is an independent prediction or an artifact of fitting to the same data.","section":"Table 1, Fig. 4"}],"minor_comments":[{"comment":"The notation for the singlet operator in Eq. (7) is introduced without an explicit definition of the orbital projectors used to define the 3z²-r² subspace.","section":"Eq. (7)"},{"comment":"Figure 2 caption does not state the temperature at which the DMFT spectra are evaluated, complicating direct comparison with the cited ARPES data.","section":"Fig. 2"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":{"model":"grok-4.3","summary":"We thank the referee for the careful review and constructive comments. We appreciate the positive assessment of the work's significance. We address each major comment below and will revise the manuscript accordingly.","responses":[{"response":"The two-site bilayer cluster is the minimal setup required to capture the interlayer dynamical singlet formation between singly occupied 3z²-r² orbitals, which constitutes the central physics. In-plane x²-y² states are treated as itinerant within the DMFT self-consistency loop, which incorporates longer-range effects through the lattice Green's function. While a larger cluster would be desirable for full convergence, such calculations are computationally prohibitive and the qualitative pressure-versus-strain dichotomy originates from c-axis compression modulating the interlayer hybridization, a feature robustly encoded in the bilayer geometry. We will add an explicit discussion of this cluster choice and its limitations to the methods section in the revision.","revision_made":"partial","referee_comment":"[§3.2] §3.2 and the computational-methods paragraph: the hybridization strength is extracted from a two-site bilayer cluster; the manuscript does not demonstrate that enlarging the cluster or adding longer-range in-plane hoppings leaves the opposite pressure-versus-strain response intact, which is load-bearing for the claimed dichotomy."},{"response":"The interaction parameters U = 8 eV and J = 1 eV were fixed from literature constrained-RPA estimates for nickelates before any comparison to ARPES or transport data; no fitting to those observables was performed. The opposite pressure-strain response of the hybridization is therefore a direct prediction. To address the concern we will add a supplementary figure showing the hybridization under pressure and strain across a range of U and J values, confirming that the qualitative dichotomy persists.","revision_made":"yes","referee_comment":"[Table 1, Fig. 4] Table 1 and Fig. 4: the reported ARPES and resistivity comparisons are shown only for the final parameter set; no sensitivity plot versus U or J is provided, leaving open whether the hybridization response is an independent prediction or an artifact of fitting to the same data."}],"tokens_in":1238,"tokens_out":455,"duration_ms":18530,"standing_objections":[]},"desk_editor":{"model":"grok-4.3","letter":"The paper's main claim is that interlayer dynamical singlets from singly occupied 3z²-r² orbitals, hybridized to the itinerant x²-y² states, control how hydrostatic pressure and epitaxial strain produce opposite effects in these materials. The calculations are said to reproduce the bulk-versus-film experimental difference plus some ARPES and transport features.\n\nThis is a direct extension of prior DFT+DMFT studies on nickelates, with the new element being the explicit tracking of how that hybridization responds differently under the two strain protocols. The approach is straightforward and stays within an established computational framework.\n\nThe abstract states that the results match experiment, which is the useful part for the subfield. At the same time, no information appears on the precise cluster size, the values chosen for U and J, sensitivity tests, or direct comparisons to alternative orbital projections or longer-range hoppings. The stress-test concern about missing contributions from apical oxygen or in-plane momentum dependence outside the cluster therefore lands as a real open question rather than a minor quibble.\n\nThe work is aimed at researchers already following nickelate superconductivity and layered transition-metal oxides. It supplies a mechanistic story for a documented experimental pattern, so it is worth sending to referees even though the current write-up leaves the robustness of the hybridization response under-specified.","headline":"Cluster DMFT work ties the pressure-strain split in bilayer nickelates to hybridization of dynamical singlets in 3z²-r² orbitals with x²-y² bands, but parameter details and validation are thin.","tokens_in":2299,"tokens_out":349,"would_cite":false,"duration_ms":11459,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.3","headline":"Interlayer dynamical singlets from 3z²-r² orbitals control the physics of bilayer nickelates and respond differently to pressure and strain.","keywords":["bilayer nickelates","dynamical singlets","Ruddlesden-Popper","cluster dynamical mean-field","epitaxial strain","hydrostatic pressure","3z2-r2 orbitals","x2-y2 orbitals"],"falsifier":"A direct measurement showing that the hybridization between 3z²-r² and x²-y² orbitals changes identically under hydrostatic pressure and epitaxial strain would contradict the calculated dichotomy.","tokens_in":2542,"feed_emoji":"","tokens_out":678,"duration_ms":21983,"temperature":0.7,"pith_summary":"Calculations within density functional theory plus cluster dynamical mean-field theory show that the electronic properties of bilayer Ruddlesden-Popper nickelates are largely governed by interlayer dynamical singlets. These singlets arise from electrons in the 3z²-r² orbitals and hybridize with the more itinerant x²-y² orbitals in the planes. The strength of this hybridization changes in distinct ways when the material is subjected to hydrostatic pressure versus epitaxial strain. This distinction explains why bulk single crystals and epitaxial thin films exhibit different behaviors, and the calculations reproduce key features seen in angle-resolved photoemission spectroscopy and transport experiments.","feed_headline":"Dynamical singlets explain nickelate pressure-strain split","feed_subtitle":"Interlayer singlets in 3z²-r² orbitals hybridize differently with x²-y² under hydrostatic pressure than epitaxial strain, matching bulk vers","key_machinery":"Interlayer dynamical singlets formed from singly occupied 3z²-r² orbitals hybridized with itinerant planar x²-y² orbitals.","core_discovery":"The physics of the bilayer Ruddlesden-Popper nickelates is to a significant extent controlled by interlayer dynamical singlets which are formed from the 3z²-r² orbitals singly occupied by electrons and are hybridized with itinerant planar x²-y² orbitals. The hybridization responds differently to hydrostatic pressure and to epitaxial strain, capturing the experimentally observed dichotomy between bulk single crystals and epitaxial thin films and reproducing several experimental results including angle-resolved photoemission and transport measurements.","pith_inferences":["The singlet picture may extend to explain similar tuning dichotomies in other layered nickelates or related transition-metal compounds.","Targeting orbital-selective hybridization in material design could be tested by combining pressure and strain in a single experiment.","If the singlets dominate, changes in their formation energy under doping would offer a testable prediction for the evolution of metallic or ordered states."],"forward_implications":["The distinct responses to pressure and strain account for the different properties seen in bulk crystals versus thin films.","The calculations reproduce angle-resolved photoemission spectra and transport data under these tuning conditions.","Tuning the occupancy and hybridization of the 3z²-r² orbitals provides a mechanism to alter the overall electronic behavior."],"fun_headline_variants":["Singlets from 3z2-r2 control bilayer nickelate pressure-strain split","Interlayer dynamical singlets respond differently to pressure and strain","Dynamical singlets in nickelates explain bulk versus thin film behavior","Hybridization of singlets differs under hydrostatic pressure and strain"],"cache_read_input_tokens":2112,"weakest_assumption_plain":"The cluster dynamical mean-field setup and orbital choice accurately capture the interlayer hybridization physics without significant contributions from other degrees of freedom or longer-range effects.","fun_headline_variants_meta":{"raw":{"variants":["Singlets from 3z2-r2 control bilayer nickelate pressure-strain split","Interlayer dynamical singlets respond differently to pressure and strain","Dynamical singlets in nickelates explain bulk versus thin film behavior","Hybridization of singlets differs under hydrostatic pressure and strain"]},"model":"grok-4.3","cost_usd":0.006051,"raw_usage":{"total_tokens":2812,"prompt_tokens":568,"num_sources_used":0,"completion_tokens":73,"cost_in_usd_ticks":60512000,"prompt_tokens_details":{"text_tokens":568,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":2171,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":568,"tokens_out":73,"duration_ms":16546,"temperature":1.0,"reasoning_tokens":2171,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-06-27T20:41:07.653837+00:00","model_set":{"reader":"grok-4.3"},"falsifier":"A direct measurement showing that the hybridization between 3z²-r² and x²-y² orbitals changes identically under hydrostatic pressure and epitaxial strain would contradict the calculated dichotomy.","supporting_citations":[],"review_version":1}