{"id":"52e1acc6-a135-4af8-8172-78b95452dfbf","arxiv_id":"2501.04203","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Local Cd-P bonds in LnCd3P3 form short one-dimensional chains whose ordering is frustrated like kagome ice, coexisting with triangular-lattice magnetism.","lead":"This paper reports that LnCd3P3 compounds contain flat cadmium-phosphorus layers whose bonds distort into short chains that cannot settle into a regular pattern, a type of frustration usually seen in magnets. It matters because the same materials already host frustrated magnetic moments on a triangular lattice, offering a rare place to study two kinds of frustration in one dopable semiconductor.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Frustration claim rests on qualitative FMC fit; fitted Ising parameters actually order at low T, so quenched disorder/weak interactions are not ruled out.","rationale":"The paper presents a coherent multi-technique case that LnCd3P3 has local Cd–P bond-length disorder in the trigonal CdP3 layers: the Cmcm small-box PDF fits improve on the hexagonal model, the trigonal-layer-only model beats tetrahedral-only, the 3D-ΔPDF shows correlated displacive features, and the FMC/RMC models reproduce the main diffuse-scattering pattern. The deposited data and code are a genuine strength. My concern is not with the existence of local two-short/one-long bond distortions (that part is reasonably supported) but with the stronger interpretive step that this is a frustrated equilibrium state analogous to kagome ice. The FMC is run at a pseudo-temperature with J3/J2 tuned to 0.58, a point in the J1-to-negative-infinity phase diagram whose ground state is the ordered 'zig-zag' Herringbone phase; geometric frustration therefore does not forbid long-range order in the model. The experimental correlation length is temperature-independent and no superstructure is seen at 5 K, which could equally signal weak energy scales or quenched disorder. The paper reports no quantitative goodness-of-fit for the diffuse-scattering models, so a random-displacement null model is not excluded. This is exactly the kind of fitted-model ambiguity that a conditional verdict should flag; I do not think it falsifies the paper, so I keep the reader's CONDITIONAL verdict unchanged.","tokens_in":22176,"tokens_out":16243,"duration_ms":180873,"concrete_test":"Using the deposited PrCd3P3 data, compute a quantitative residual (e.g., normalized R-factor or reduced chi-square) between the experimental (H,K,5/2) diffuse plane and (i) the optimized FMC model, (ii) the RMC model, (iii) a control with the same Cdtrig/Ptrig displacement amplitudes but randomly permuted long-bond placements, and (iv) the Potts model. If the random control achieves a residual within about 2x the FMC residual or reproduces the half-integer maxima, the specific frustrated Ising/dimer interaction is not required by the data and the frustration claim is unsupported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim requires that the short-range bond correlations observed in LnCd3P3 are a frustrated equilibrium state (kagome-ice-like) rather than, say, a quenched distribution of local distortions. The modeling step that carries this weight is the FMC/RMC assignment of the diffuse scattering to a specific J1–J2–J3 Ising/dimer model with J3/J2 = 0.58. That assignment is only argued qualitatively ('well-captured', 'agrees well'); no R-factor, chi-square, or null-model comparison is reported. More importantly, the fitted parameter regime (J3/J2 > 0.5, J1 large) has an ordered 'zig-zag' (Herringbone) ground state, so geometric frustration does not actually prevent long-range order in the model. The observed absence of superstructure at 5 K and the near-temperature-independent correlation length between 80 K and 300 K are therefore equally compatible with (i) interactions so weak that the ordering temperature lies below 5 K, (ii) quenched disorder or glassiness inherited from synthesis, or (iii) the claimed frustration. Without a quantitative test that rules out (i)/(ii), the 'highly frustrated' qualifier and the kagome-ice analogy are not established.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports a combined X-ray and neutron scattering study of the LnCd3P3 family (Ln = La, Ce, Pr, Nd), arguing that the trigonal-planar CdP3 layers undergo a local Cmcm-type distortion in which two Cd–P bonds shorten and one lengthens, creating a dimer-like bond order on a honeycomb network. The authors map the dimer correlations to a frustrated J1–J2–J3 Ising model on the dual triangular lattice and, using forward and reverse Monte Carlo modeling, propose that a Herringbone-like pattern with J3/J2 = 0.58 reproduces the observed diffuse scattering and 3D-ΔPDF. They interpret the lack of long-range bond order as geometric frustration analogous to kagome ice, and support this with powder PDF fits, single-crystal diffuse scattering, and inelastic neutron scattering showing a split non-Kramers doublet that is attributed to the rare-earth ion feeling the local symmetry lowering.","tokens_in":22395,"tokens_out":5035,"duration_ms":54283,"significance":"If established, the paper reports a rare material family in which a frustrated bond-order instability coexists with a triangular-lattice antiferromagnet in a dopable semiconducting host. The independent structural observations—Cmcm small-box PDF fits, quasi-2D diffuse planes, and 3D-ΔPDF quadrupole features—are mutually consistent, and the mapping to a dimer/Ising/kagome-ice description is conceptually elegant. The public data and code availability are also strengths. However, the quantitative link between the model and the data is not yet established: the key parameter J3/J2 is selected by visual comparison, no null model is tested, and the temperature-independent short-range correlations are equally compatible with weak interactions or quenched disorder as with geometric frustration. The central claim is therefore plausible but not fully demonstrated.","major_comments":[{"comment":"The FMC model with J3/J2 = 0.58 is chosen on the basis of qualitative agreement ('well-captured', 'agrees well'), and no R-factor, chi-square, or residual metric is reported for either the FMC or RMC comparison against the experimental diffuse scattering or 3D-ΔPDF. Because the displacement amplitudes are fixed from the same powder PDF data and the pseudotemperature is adjustable, the agreement is not a parameter-free confirmation. Please provide a quantitative goodness-of-fit measure and compare explicitly with at least one null model (e.g., random distortions, the J3/J2 < 0.5 stripe limit, or the Potts model with optimized parameters) to show that the Herringbone Ising model uniquely reproduces the half-integer maxima.","section":"II.B / Fig. 4 / SI VI.A"},{"comment":"The paper's own phase diagram (ref. 66 and Fig. 1f-g) shows that the J3/J2 = 0.58 region has an ordered zig-zag/Herringbone ground state in the Ising model. The observed absence of superstructure at 5 K and the near-temperature-independent correlation length between 80 and 300 K are therefore equally consistent with (i) an ordering temperature below 5 K because the interactions are weak, (ii) quenched disorder or glassiness inherited from synthesis, or (iii) the claimed geometric frustration. The manuscript needs an explicit test that distinguishes these scenarios—for example, a temperature-dependent correlation-length/lineshape analysis, a comparison of diffuse intensity versus thermal population expectations, or an estimate of the interaction energy scale—before the 'highly frustrated' qualifier and the kagome-ice analogy are supported.","section":"II.B / Fig. 3(g) / Discussion"},{"comment":"The FMC and RMC cells retain only Cdtrig and Ptrig displacements, while the paper itself reports secondary Cdtet/Ptet correlations with a real-space feature at (1/3, 2/3, 1/8)-type vectors. The statement 'These do not impact the mapping to the 2D Ising model' is an assertion rather than a demonstrated result. Please show quantitatively that including or omitting tetrahedral-layer displacements leaves the simulated diffuse scattering and ΔPDF unchanged (or within statistical error), or otherwise delimit which features of the data require the tetrahedral contribution. Without this, the assignment of the full diffuse scattering to a purely in-plane dimer model remains incomplete.","section":"II.B / SI IV / Fig. S7"},{"comment":"The FMC procedure is described only partially: it uses 'occupational correlation vectors' with fixed energies J1, J2, J3, a pseudotemperature of 1.0, and it discards overlapping dimers by randomly choosing one and removing the others. Each of these choices can affect the simulated correlation length and the appearance of half-integer maxima. Please quantify their influence (e.g., by varying the pseudotemperature, supercell size, thermalization length, and the rule for resolving dimer overlaps) so that the 'thermalized ground state' result is robust.","section":"SI VI.A"}],"minor_comments":[{"comment":"Please state explicitly whether the reported correlation lengths are extracted only from the Lorentzian component of the pseudo-Voigt fits and how (or whether) instrumental resolution was deconvolved from the diffuse rods.","section":"Fig. 3(g)"},{"comment":"The manuscript alternates between '3D-ΔPDF' for the experimental maps and '2D-ΔPDF' for the simulated maps; a sentence defining this convention would help readers understand that the model is two-dimensional by construction.","section":"Fig. 4 / Fig. S12"},{"comment":"Several fitted Stevens parameters differ from the point-charge model by orders of magnitude (for example B4^0 and B6^0); a brief comment on the identifiability and conditioning of the CEF fit would strengthen the claim that the split doublet is a robust probe of local symmetry breaking.","section":"SI VII / Table S2"}],"recommendation":"major_revision","confidential_remarks":"I see no scope, citation, or data-availability problems; the Zenodo deposit is a clear strength. The main revision should focus on model selection and on ruling out weak-interaction or quenched-disorder explanations for the short-range bond order. I would be willing to look at a revised version."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"What you should know: this is a careful experimental paper that convincingly shows a local bond distortion in the CdP3 layers of LnCd3P3, with short-range correlations that look Herringbone-like. The \"frustrated bond order\" framing, however, is stronger than the modeling supports. The fitted Ising parameters actually have an ordered zig-zag ground state, so the absence of long-range order down to 5 K is not explained by frustration in their own model.\n\nThe genuinely new content is the local structure: powder PDF showing a Cmcm distortion with two short and one long Cd-P bond, quasi-2D diffuse scattering in single crystals, 3D-DeltaPDF features, and a plausible mapping of the long bonds to dimers on a honeycomb lattice. The evidence chain from PDF to diffuse scattering to Monte Carlo is internally consistent, and the data are on Zenodo. The supplement is honest, including the alternative Potts model and the hedge on CEF coupling.\n\nThe soft spots are real but not fatal to the empirical core. The FMC/RMC comparison to the data is qualitative—no R-factor, no chi-square, no null model. The J3/J2 = 0.58 value and the displacement amplitudes come from the same data they are used to reproduce. More importantly, the stress-test concern holds: in the dual triangular Ising model with large J1 and J3/J2 > 0.5, the ground state is the zig-zag (Herringbone) state. That is an ordered state. So using this model to argue that frustration prevents long-range order is unconvincing. The temperature-independent correlation length and the lack of superstructure at 5 K are equally compatible with quenched disorder or interactions too weak to order above 5 K. The authors don't rule those out.\n\nThat said, the local distortion itself and its short-range correlations are not a restatement of the model—those observations stand. The paper would be stronger if it presented a quantitative comparison against a disordered random-distortion model and addressed the ordered-ground-state issue.\n\nWho is this for? Experimentalists in frustrated magnetism and local structure. It deserves a serious referee, but the referee should push on the distinction between \"frustrated equilibrium state\" and \"quenched or weakly interacting distortions.\" I'd recommend sending it to peer review with that in mind.\n\nRecommendation: engage with it, but treat the frustration claim as provisional.","headline":"Solid local-structure study with a frustrated-bond narrative that outruns the evidence.","tokens_in":23060,"tokens_out":2437,"would_cite":true,"duration_ms":24529,"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":"The paper demonstrates that the LnCd3P3 family of triangular-lattice antiferromagnets (Ln = La, Ce, Pr, Nd) hosts a frustrated bond order in its CdP3 layers, where each trigonal-planar CdP3 unit forms one long and two short Cd–P bonds…","keywords":["frustrated bond order","kagome ice","honeycomb dimer model","triangular lattice antiferromagnet","diffuse X-ray scattering","pair distribution function","LnCd3P3","crystalline electric field"],"falsifier":"Measure the diffuse scattering at 5 K in a strain-free crystal: if the half-integer maxima sharpen into resolution-limited Bragg peaks, long-range bond order has formed and the frustrated kagome-ice description is wrong.","tokens_in":21958,"feed_emoji":"🧊","tokens_out":6014,"duration_ms":56130,"temperature":0.7,"pith_summary":"The paper sets out to show that the LnCd3P3 compounds are not merely triangular-lattice antiferromagnets but also host a frustrated bond order in their cadmium-phosphide layers. In each trigonal-planar CdP3 unit, two Cd–P bonds shorten and one lengthens, creating a local dimer that wants to order; the honeycomb network of such dimers is frustrated, so instead of long-range order the bonds form short one-dimensional chains that alternate direction in a Herringbone pattern over about 3–5 unit cells and persist to room temperature. The authors argue this maps onto a two-dimensional Ising antiferromagnet on the dual triangular lattice and, equivalently, onto kagome ice correlations. The same lattice that frustrates the rare-earth magnetic moments is a narrow-gap, dopable semiconductor, so charge and spin frustration coexist in one tunable platform, and the distortion is shown to reach the rare-earth crystal field.","feed_headline":"Bond order stays frustrated to 300 K in CdP3 layers","feed_subtitle":"Each CdP3 unit lengthens one Cd–P bond; the long bonds form short chains that map to kagome ice, coupling to the rare-earth moments.","key_machinery":"The central object is the dimer covering of the honeycomb Cd–P network, where each dimer is the single long Cd–P bond within a CdP3 unit. This covering maps one-to-one onto an antiferromagnetic Ising model on the dual triangular lattice, with spins on hexagon centers and dimers on frustrated up-up or down-down edges, and in the J1–J2–J3 regime it reduces to kagome ice correlations. The paper tunes J3/J2 to select the Herringbone (zig-zag) local order that reproduces the observed half-integer diffuse maxima, with the optimized ratio J3/J2 = 0.58.","core_discovery":"The paper claims that the LnCd3P3 family hosts a frustrated lattice instability rooted in the trigonal-planar CdP3 layers: each CdP3 unit distorts so that two of its three Cd–P bonds become shorter and one becomes longer, and this long bond acts as a dimer on the honeycomb Cd–P network. Because the honeycomb dimer tiling is frustrated, the dimers do not lock into long-range order; instead they organize into short one-dimensional CdP chains, about 3–5 unit cells long, with a weak Herringbone-like alternation of chain direction. The diffuse X-ray scattering, modeled by forward and reverse Monte Carlo, is consistent with a dual triangular-lattice Ising model with J3/J2 > 0.5, equivalent to emergent kagome ice correlations, and the distortion is shown to lower the crystal field at the Pr3+ site, demonstrating coupling between the bond order and the magnetic layer.","pith_inferences":["If the bond order couples to the crystal field in the insulating state and survives doping, carrier doping may screen or amplify the local distortion; measuring diffuse scattering in metallic Pr1−xSrxCd3P3 would test this directly.","The same dimer-to-Ising mapping should apply to the isostructural arsenide analogues LnCd3As3, where off-centering has already been reported; comparing their diffuse scattering would test whether the frustrated bond order is generic to the family.","The model assumes only in-plane dimer correlations; if tetrahedral-layer displacements contribute more substantially at lower temperature or under pressure, an extended three-dimensional model would be required.","A natural extension is to apply uniaxial strain and watch the half-integer diffuse maxima: sharpening or rotation of the pattern would confirm the degeneracy and the kagome-ice description."],"forward_implications":["No long-range bond order forms up to 300 K; the diffuse scattering is essentially temperature-independent, so the frustrated bond correlations persist over the entire measured range.","The bond distortion is not confined to the CdP layer: inelastic neutron scattering shows the rare-earth crystal-field ground state senses the local symmetry lowering, so the bond and magnetic channels are coupled.","The family is dopable: intentional Sr substitution into PrCd3P3 drives a transition toward metallic behavior with hole carriers, opening a route to tune carriers inside a frustrated lattice.","The mapping to a dual Ising model predicts two local orders, staggered/stripe versus Herringbone, and the observed half-integer diffuse maxima select the Herringbone branch with J3/J2 around 0.58.","Because the distortion pathways are three-fold degenerate, three 120°-rotated domains form, implying that external strain or field can bias among degenerate local configurations."],"supporting_citations":[{"why":"Establishes kagome ice correlations as the classical solution of frustrated Ising spins on the kagome lattice, the target pattern the bond order is compared to.","marker":"57"},{"why":"Supplies the exact mapping between dimer coverings of the honeycomb lattice and Ising degrees of freedom on the dual triangular lattice used to build the forward Monte Carlo model.","marker":"58,59"},{"why":"Provides the J1–J2–J3 triangular Ising phase diagram whose stripe and zig-zag regimes select the staggered versus Herringbone dimer configurations.","marker":"66"},{"why":"Define the three-dimensional difference pair distribution function method used to extract real-space displacement correlations from diffuse scattering.","marker":"63,64"},{"why":"The modified punch-and-fill procedure is used to isolate the diffuse scattering from the single-crystal data before ΔPDF analysis.","marker":"65"},{"why":"Reports analogous off-centering of Cd and As layers in LnCd3As3 compounds, supporting a generic in-plane instability across the family.","marker":"56"},{"why":"Provides the electronic-structure context locating the valence band maximum on phosphorus p states in the CdP layer, which motivates the bonding (Jahn-Teller) instability.","marker":"39"},{"why":"Reports the weakly metallic self-doped state of CeCd3P3, the baseline for the doping claims in this paper.","marker":"42"}],"fun_headline_variants":["Frustrated CdP3 bonds mimic kagome ice without long-range order","Long Cd-P bonds stay disordered down to 300 K and couple to magnetism","Triangular antiferromagnet hosts both frustrated charges and spins","CdP3 layers show emergent kagome ice from bond order frustration"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The entire interpretation hinges on the local distortion being exactly a dimer pattern—one long and two short Cd–P bonds per CdP3 unit—confined to the trigonal CdP plane, with the tetrahedral layers' motions too weak to matter.","fun_headline_variants_meta":{"raw":{"variants":["Frustrated CdP3 bonds mimic kagome ice without long-range order","Long Cd-P bonds stay disordered down to 300 K and couple to magnetism","Triangular antiferromagnet hosts both frustrated charges and spins","CdP3 layers show emergent kagome ice from bond order frustration"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000657,"raw_usage":{"total_tokens":3013,"prompt_tokens":960,"completion_tokens":2053,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":576,"completion_tokens_details":{"reasoning_tokens":1973}},"tokens_in":576,"tokens_out":2053,"duration_ms":13813,"temperature":1.0,"reasoning_tokens":1973,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T21:39:13.847150+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the diffuse scattering at 5 K in a strain-free crystal: if the half-integer maxima sharpen into resolution-limited Bragg peaks, long-range bond order has formed and the frustrated kagome-ice description is wrong.","supporting_citations":[{"cited_title":"& Sondhi, S","cited_arxiv_id":null,"evidence_quote":"Establishes kagome ice correlations as the classical solution of frustrated Ising spins on the kagome lattice, the target pattern the bond order is compared to."},{"cited_title":"& Mila, F","cited_arxiv_id":null,"evidence_quote":"Provides the J1–J2–J3 triangular Ising phase diagram whose stripe and zig-zag regimes select the staggered versus Herringbone dimer configurations."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"The modified punch-and-fill procedure is used to isolate the diffuse scattering from the single-crystal data before ΔPDF analysis."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Reports analogous off-centering of Cd and As layers in LnCd3As3 compounds, supporting a generic in-plane instability across the family."},{"cited_title":"& Kitagawa, J","cited_arxiv_id":null,"evidence_quote":"Provides the electronic-structure context locating the valence band maximum on phosphorus p states in the CdP layer, which motivates the bonding (Jahn-Teller) instability."},{"cited_title":"R., Broun, D","cited_arxiv_id":null,"evidence_quote":"Reports the weakly metallic self-doped state of CeCd3P3, the baseline for the doping claims in this paper."}],"review_version":1}