{"id":"e82851d8-6d1b-4a41-be53-0ff8101fcf39","arxiv_id":"2411.16957","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Classical Monte Carlo simulations confirm that the triplet boson model of orthonickelates phase-separates into macroscopic domains, matching the mean-field Maxwell construction.","lead":"This paper simulates a theoretical model of nickel oxide materials using two different methods: an approximate analytical calculation and a computer Monte Carlo simulation. The computer simulation matches the analytical prediction that the material separates into distinct regions with different electronic phases.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"MC algorithm omits Metropolis-Hastings proposal-ratio correction; non-symmetric pair proposal biases sampled distribution, so the phase-separation evidence in Figs. 4–7 is not tied to the stated Boltzmann weight.","rationale":"The reader's weakest assumption correctly notes that the classical product-state ansatz is unvalidated against exact quantum methods. That is a general concern about the physical relevance of the model. However, the more immediate and more load-bearing problem is that the MC algorithm, as described, does not even sample the intended Boltzmann distribution for that classical ansatz. The proposal distribution in Sec. 3 is manifestly non-symmetric, and the paper gives no indication of the required Hastings factor. This is an internal inconsistency, not merely a lack of external validation. It directly affects the numerical evidence for the central claim: the phase-separated states and the apparent agreement with MFA could be artifacts of the biased sampling, especially at the temperatures where T_c,max is determined. Because the numerical simulation is the only evidence connecting the MFA Maxwell construction to the 'observed state of the system,' an invalid sampling algorithm voids the central claim as presented. The concern is concrete and testable: a simple modification of the acceptance probability or proposal scheme would settle it. If the corrected simulation reproduces the same phase-separation regions, the paper could be rehabilitated; but as written, the central numerical result is unreliable. This justifies a REJECT verdict rather than CONDITIONAL, since the main evidence is currently invalid, although the analytic MFA part may survive.","tokens_in":7599,"tokens_out":7026,"duration_ms":73466,"concrete_test":"Re-run the same N=96×96 simulations at V/J=4, t/J=1.5 using the same pair-proposal scheme but with the correct Metropolis-Hastings acceptance min(1, exp(-ΔE/T) · (n1,old n2,old / n1,new n2,new)^2), or equivalently replace the proposal with a symmetric one (e.g., small random rotations of the on-site 8D unit vectors). Then recompute the phase-separation regions and the order-parameter curves in Figs. 4–7. If the BS/CO and CO/AFM coexistence regions and the representative snapshots at points A, B, C remain qualitatively unchanged, the concern is moot; if the boundaries shift substantially or the phase separation disappears, the paper's central claim fails.","verdict_should_be":"REJECT","load_bearing_attack":"The central numerical claim rests entirely on the classical MC simulations of Sec. 3. The update rule selects a pair of sites, fixes their mean density n̄, and draws new densities from the conditional prior derived from uniform sampling on the 8D unit sphere. This proposal is not symmetric: q(new|old) depends on the old state only through n̄, but the joint density of (n1, n2) at fixed n̄ is proportional to n1^2 n2^2. Consequently q(old)/q(new) = (n1,old n2,old / n1,new n2,new)^2, which is not unity. The paper states that 'standard Metropolis algorithm' is used, but standard Metropolis requires symmetric proposals. Without the Hastings correction factor min(1, exp(-ΔE/T) · q(old)/q(new)), the stationary distribution is π(state) ∝ q(state) exp(-E/T) instead of the intended exp(-E/T). Since q(state) includes a factor ∏_i n_i^2, the simulation effectively samples a modified ensemble with an extra density-dependent weight. This bias is temperature-independent in the log and of order 2T ln n, which at the simulated T/J up to 0.63 is comparable to t/J = 1.5 and V/J = 4. It can shift coexistence boundaries, alter order-parameter onsets, and even stabilize or suppress phase separation. The reported T_c,max ratio of ~6.3 between MFA and MC may partly be an artifact of this sampling bias rather than a physical quantum-vs-classical difference. Thus the conclusion that 'numerical simulation by classical MC method demonstrated ... phase separation' is not supported by the algorithm as described.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies a spin-triplet composite-boson model for orthonickelates (Hamiltonian (1)). In Sec. 2 the authors derive MFA stability conditions for homogeneous charge-ordered (CO), boson superfluid (BS), and antiferromagnetic (AFM) phases, and construct phase-separation binodals via the Maxwell construction. In Sec. 3 they introduce a classical Monte Carlo algorithm that samples product-state wave functions (Eq. (3)) with a kinematically fixed total boson density. In Sec. 4 they compare MC order-parameter onsets with the MFA phase diagram for z=4, V/J=4, t/J=1.5, after rescaling the temperature scales, and report phase separation into macroscopic BS/CO and CO/AFM domains. The central conclusion is that the classical MC simulation confirms the MFA Maxwell-construction picture.","tokens_in":7947,"tokens_out":15188,"duration_ms":139591,"significance":"The analytical stability conditions and the Maxwell-construction analysis (Sec. 2) are clear and internally consistent within the mean-field framework. The MC simulation is a genuinely independent numerical procedure that does not reduce by construction to the MFA result, so the qualitative agreement is a meaningful test of the Maxwell construction inside the product-state manifold. The snapshots in Fig. 7 directly show macroscopic domain coexistence. The main limitations are that the MC algorithm as described is not a correct Metropolis update for the stated Boltzmann weight, and the quasi-classical product-state ansatz is not validated against quantum benchmarks; the paper itself flags that scaling-theory estimates are 'planned to be made in the future' (Sec. 4).","major_comments":[{"comment":"The pair-update proposal is not symmetric. The conditional density f1(n1; n̄) in Eq. (9) is proportional to n1^2 (2n̄ - n1)^2 = n1^2 n2^2, so for a move from an old state x to a new state x' on the same pair, q(x'|x) ∝ (n1' n2')^2 while q(x|x') ∝ (n1 n2)^2. The ratio q(x|x')/q(x'|x) is generally not unity. The statement in Sec. 3 that 'we use standard Metropolis algorithm' is therefore insufficient: standard Metropolis requires symmetric proposals. The acceptance probability must include the Metropolis-Hastings factor min(1, exp(-ΔE/T) (n1 n2 / n1' n2')^2). Without this correction, the stationary distribution is ∝ exp(-E/T) ∏_i n_i^2 instead of exp(-E/T), and since 2T ln n is of order unity over the simulated T/J range, the bias can shift phase boundaries and affect the order-parameter onsets. The authors should either implement the correction or demonstrate numerically that this bias does not alter the reported phase diagram.","section":"Sec. 3, after Eqs. (4)-(12)"},{"comment":"The quasi-classical product-state ansatz in Eq. (3) is an uncontrolled approximation to the quantum Hamiltonian (1). The MC simulation samples configurations in the manifold of product states with a uniform measure, and the MFA of Ref. [9] operates on the same manifold. Consequently, agreement between the MC and MFA results does not validate the ansatz itself; it only confirms the Maxwell-construction picture within the classical manifold. To support the claim that phase separation is a property of the quantum model, the paper should benchmark the classical MC against exact diagonalization or quantum Monte Carlo for small lattices, or explicitly restrict its conclusion to the quasi-classical version of the model.","section":"Sec. 3, Eqs. (3)-(7); Sec. 4"},{"comment":"The critical temperatures are extracted from an arbitrary 1% threshold of the order parameter, with no error bars, finite-size scaling, or thermalization/blocking analysis. The paper itself states that 'more extensive analyses with the estimate of critical temperatures as per scaling theory are planned to be made in the future.' Since the temperature rescaling T_c,max(MFA)/T_c,max(MC) ≈ 6.3 is a central element of the comparison, the choice of threshold directly affects the reported scale ratio. Please provide sensitivity tests (e.g., different thresholds, system sizes, and numbers of MC steps) and statistical error estimates, or explicitly present the phase diagram as preliminary.","section":"Sec. 4, Figs. 4-6"}],"minor_comments":[{"comment":"The phrase 'observed state of the system in numerical simulations' should be clarified as 'the state observed in the numerical simulations', since no experimental observation is made.","section":"Abstract"},{"comment":"The formulas for n1,min and n1,max could be derived explicitly; as written they are correct but not transparent. A short explanation that they enforce n1, n2 ∈ [0,1] would help the reader.","section":"Eq. (12)"},{"comment":"The formula for T_BS has ambiguous parentheses: T_BS = 4t(4n-3)[3 ln(n/(3(1-n)))]^{-1} would be clearer than the printed expression.","section":"Fig. 2 caption and Sec. 2"},{"comment":"The text says 'uniformly distributed in section [0,1]' where 'interval' or 'segment' is meant; this is a language issue that should be corrected throughout.","section":"Sec. 3"},{"comment":"The MC implementation does not specify how many initial steps are discarded for thermalization before averaging, nor whether the 4×10^6 steps include thermalization. Please specify the averaging protocol.","section":"Sec. 3, step 7"},{"comment":"The statement that 'AFM phase is unstable at small inter-center charge-to-charge interaction, V/J < 1' is an MFA result from Sec. 2; the MC simulations were only performed at V/J = 4. The conclusion should distinguish between MFA-derived and MC-verified claims.","section":"Sec. 5, Conclusion"}],"recommendation":"major_revision","confidential_remarks":"The paper has a useful analytical MFA analysis, but the central numerical claim is undermined by the incorrect Metropolis proposal. The fix is in principle within scope: implement the Metropolis-Hastings correction and rerun. The lack of any quantum benchmark and the arbitrary 1% threshold are further concerns. The paper would be better framed as a study of the quasi-classical variant of the model, not as a Monte Carlo method for the full quantum Hamiltonian. The topic fits cond-mat.stat-mech, and the snapshots are compelling, but the current version requires substantial revision before it can be considered for publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a modest but honest MFA plus classical Monte Carlo study of a triplet-boson model with fixed density. The new part is the MC implementation and the MC phase diagram, and the main qualitative conclusion—that the system phase separates instead of forming homogeneous mixed phases—is plausible. But there is a real algorithmic flaw in the MC proposal that the authors need to fix or explain: the pair update in Sec. 3 is not symmetric, and the paper gives no Hastings correction, so the sampled ensemble carries an extra n_i^2 weighting. That is a genuine problem, not a nitpick.\n\nCredit where due: the MFA stability analysis is clean, and the Maxwell construction is clearly laid out. The comparison between MFA and MC is a real comparison, not a fit, since the MC is an independent numerical procedure. Fig. 7 shows nice snapshots of coexisting domains, and the authors are upfront that the 1% threshold is ad hoc and that finite-size scaling has not been done.\n\nThe soft spots, in order. First, the Metropolis-Hastings issue. The proposal density for new densities n1,n2 at fixed n̄ is proportional to n1^2 n2^2, which depends on the new values. Accepting with the Boltzmann factor alone targets a modified distribution. The bias is of order 2T ln n, comparable to the couplings at the simulated T. It could shift T_c and coexistence boundaries, including the reported MFA/MC temperature-scale ratio of about 6.3. The authors should add the ratio, or show numerically that it doesn't matter. Second, there are no error bars, no multiple seeds, and only one lattice size. For a phase-separation claim, you need at least some finite-size evidence. Third, the quasi-classical product-state ansatz is not benchmarked against exact quantum results on small clusters, so the quantitative agreement with MFA could partly come from the shared classical approximation. The authors acknowledge this indirectly by saying scaling analyses are planned.\n\nWho is this for? People working on nickelate models and on constrained boson Monte Carlo. It is a valid piece of work, not a breakthrough, and the central message is probably right. I would send it to peer review: it deserves a serious referee. But I would ask the authors to fix the Hastings issue and add error bars or finite-size scaling before acceptance.","headline":"Plausible MFA/MC phase-separation story, but the MC algorithm has a non-symmetric proposal without Hastings correction that needs to be fixed.","tokens_in":8505,"tokens_out":4792,"would_cite":false,"duration_ms":47467,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["82B26","82B80","82B20"],"pacs":[],"model":"deepseek-v4-flash","headline":"The triplet-boson model for orthonickelates, simulated with a classical Monte Carlo method that keeps boson density fixed, forms macroscopic phase-separated domains rather than homogeneous states with several nonzero order parameters.","keywords":["orthonickelates","triplet boson model","phase separation","mean-field approximation","classical Monte Carlo","Maxwell construction","charge ordering","antiferromagnetism"],"falsifier":"Run an unbiased quantum calculation, such as exact diagonalization on a small cluster or a projective quantum Monte Carlo for Hamiltonian (1) at $z=4$, $V/J=4$, $t/J=1.5$ and $n=0.65$; if the ground state exhibits simultaneous nonzero CO and AFM order parameters on the same sites, or fails to show spatial separation into CO and AFM domains, the paper's central claim would be refuted.","tokens_in":7380,"feed_emoji":"🧲","tokens_out":7973,"duration_ms":67954,"temperature":0.7,"pith_summary":"This paper argues that the low-temperature states of the triplet-boson model for orthonickelates are not homogeneous phases with several simultaneous order parameters but spatial mixtures of pure phases. Using a classical Monte Carlo algorithm with kinematically fixed boson concentration, the authors find that for $z=4$, $V/J=4$, $t/J=1.5$ the system separates into boson-superfluid/charge-ordered domains for $0.15<n<0.50$ and into charge-ordered/antiferromagnetic domains for $0.5<n<1.0$. The result matters because it gives a concrete prediction for phase coexistence in nickelates and because it supports the Maxwell construction used in mean-field treatments as a qualitatively correct description.","feed_headline":"Simulation: nickelate model prefers phase separation","feed_subtitle":"Classical Monte Carlo and mean-field agree: no mixed phases, only CO/BS and CO/AFM domains.","key_machinery":"The central object is the quasi-classical on-site wave function $|\\psi_i\\rangle = c_{1,11}|1,11\\rangle_i + c_{1,10}|1,10\\rangle_i + c_{1,1-1}|1,1-1\\rangle_i + c_{0,00}|0,00\\rangle_i$, with four amplitudes parametrized by angles $\\theta_i,\\psi_i,\\phi_i$ and uniformly sampled on the 8-dimensional sphere of coefficients. The Metropolis updates act on pairs of sites: a new density on one site is drawn from the inverse distribution function $F_1(n_1;\\bar n)$ constructed from the density-weighted measure $f(n)=3n^2$, and the second site density is fixed by $\\bar n$, which keeps the total boson concentration constant at every step. This machine lets the simulation reveal whether mixed phases persist or decompose, and it is compared against the MFA Maxwell construction built from chemical-potential equality $\\mu_i(n,T)=\\mu^*$ and free energy $f=m_1f_1+m_2f_2$.","core_discovery":"On the paper's own terms, the central discovery is that the mixed-order-parameter solutions of the mean-field approximation are not realized as homogeneous states when the full lattice is simulated. The classical Monte Carlo method with fixed boson density shows the same phase separation that the MFA Maxwell construction predicts: BS/CO coexistence for $0.15<n<0.50$ and CO/AFM coexistence for $0.5<n<1.0$ at $z=4$, $V/J=4$, $t/J=1.5$. The temperature scales of MFA and MC differ by a factor of about $6.3$, yet the ordering of critical temperatures of CO and AFM phases persists, while the BS phase occupies a smaller region in the MC diagram than in MFA.","pith_inferences":["If the product-state ansatz is what drives the qualitative agreement, an entangled quantum calculation could find that quantum fluctuations stabilize some homogeneous mixed phases; this is an inference, not a result of the paper.","The fixed-density pair-update algorithm could be adapted to other conserved-density hard-core boson models, since it avoids the flat-chemical-potential problem that makes grand-canonical simulations slow at low temperature.","The apparent shift of the CO critical-temperature maximum from $n=0.5$ to about $n=0.55$ in the MC data may be a threshold or finite-size artifact; a scaling-theory analysis, which the authors say is planned, would settle whether it is physical."],"forward_implications":["Homogeneous 'supersolid-like' mixed phases of the triplet-boson model are thermodynamically unstable at low temperature and decompose into macroscopic phase-separated regions.","The Maxwell construction in mean-field theory gives qualitatively correct binodals for the model, so MFA phase diagrams can be trusted for locating coexistence regions at least at the qualitative level.","The boson-superfluid phase is unstable at high concentrations and separates with the non-ordered phase, unlike the singlet local-boson model.","The antiferromagnetic phase is unstable at weak charge-charge interaction ($V/J<1$) and gives way to AFM/NO phase separation, while at $V/J>1$ homogeneous AFM is stable at all concentrations below its critical temperature."],"supporting_citations":[{"why":"Supplies the mean-field phase diagrams, free energies, and the claim that mixed phases have higher free energy than phase separation, which the MC results are compared against.","marker":"[9]"},{"why":"Defines the spinless hard-core boson model whose thermodynamic behavior this triplet-boson model generalizes.","marker":"[10]"},{"why":"Provides the independent quantum Monte Carlo demonstration for local singlet bosons that the homogeneous supersolid-like phase is unstable and phase separation occurs.","marker":"[11]"},{"why":"Introduces the Maxwell phenomenological construction used to set phase-separation boundaries in mean-field theory.","marker":"[12]"},{"why":"Documents that free energy of mixed phases in the local-boson model is higher than the free energy of phase separation, supporting the same conclusion used here.","marker":"[13]"}],"fun_headline_variants":["Monte Carlo and mean-field agree: no mixed nickelate phases","Simulation: nickelate model splits into coexisting domains","No homogeneous mixed phases in orthonickelate model","Mean-field and Monte Carlo both show phase separation","Nickelate model: no mixed phases, only coexisting domains"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The classical Monte Carlo method replaces the quantum state by a product of single-site quasi-classical wave functions, so if quantum entanglement changes which phases are stable, the agreement between Monte Carlo and mean-field could reflect the shared classical approximation rather than the true quantum model.","fun_headline_variants_meta":{"raw":{"variants":["Monte Carlo and mean-field agree: no mixed nickelate phases","Simulation: nickelate model splits into coexisting domains","No homogeneous mixed phases in orthonickelate model","Mean-field and Monte Carlo both show phase separation","Nickelate model: no mixed phases, only coexisting domains"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000413,"raw_usage":{"total_tokens":2036,"prompt_tokens":743,"completion_tokens":1293,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":359,"completion_tokens_details":{"reasoning_tokens":1210}},"tokens_in":359,"tokens_out":1293,"duration_ms":9455,"temperature":1.0,"reasoning_tokens":1210,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T12:42:27.632051+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run an unbiased quantum calculation, such as exact diagonalization on a small cluster or a projective quantum Monte Carlo for Hamiltonian (1) at $z=4$, $V/J=4$, $t/J=1.5$ and $n=0.65$; if the ground state exhibits simultaneous nonzero CO and AFM order parameters on the same sites, or fails to show spatial separation into CO and AFM domains, the paper's central claim would be refuted.","supporting_citations":[{"cited_title":"Ryumshin, S.V","cited_arxiv_id":null,"evidence_quote":"Supplies the mean-field phase diagrams, free energies, and the claim that mixed phases have higher free energy than phase separation, which the MC results are compared against."},{"cited_title":"Micnas, J","cited_arxiv_id":null,"evidence_quote":"Defines the spinless hard-core boson model whose thermodynamic behavior this triplet-boson model generalizes."},{"cited_title":"Batrouni, R.T","cited_arxiv_id":null,"evidence_quote":"Provides the independent quantum Monte Carlo demonstration for local singlet bosons that the homogeneous supersolid-like phase is unstable and phase separation occurs."},{"cited_title":"Kapcia, S","cited_arxiv_id":null,"evidence_quote":"Introduces the Maxwell phenomenological construction used to set phase-separation boundaries in mean-field theory."},{"cited_title":"Kapcia, J Supercond Nov Magn 26, 8, 2647 (2013)","cited_arxiv_id":null,"evidence_quote":"Documents that free energy of mixed phases in the local-boson model is higher than the free energy of phase separation, supporting the same conclusion used here."}],"review_version":1}