{"id":"9c4608c5-8fe8-4bcc-8c98-d39cf344dbd0","arxiv_id":"2505.07999","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"In a stochastic Gross-Pitaevskii model with non-Markovian reservoir memory, polariton condensate lattices show a temperature-driven crossover from ferromagnetic to antiferromagnetic phase order and retain long-range coherence at 50 K, where a single spot does not condense.","lead":"This paper simulates lattices of exciton-polariton condensates and shows that the phase order switches from in-phase to checkerboard as temperature rises, while the lattice stays ordered at temperatures where a single spot fails to condense. It matters because it suggests polariton lattices can remain phase-coherent at much higher temperatures than single condensates, which is relevant for future polaritonic devices.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"FM/AFM crossover rests on single realizations, not an order parameter; stochastic selection could explain the apparent temperature dependence.","rationale":"The reader's weakest_assumption is the exponential-memory approximation; that is a legitimate model-fidelity concern. However, the single most load-bearing issue for the central claim is internal: even if the model and the exponential kernel are accepted, the FM/AFM crossover is evidenced by only two single snapshots, with no order parameter or realization statistics. Because the system is stochastic and FM/AFM states can be metastable, the headline claim is not falsifiable from the presented data. A modest ensemble-statistics check would settle it. I therefore leave the reader's CONDITIONAL verdict unchanged but shift the emphasis from the memory kernel to the missing phase-order statistics and the underived steady-state assumption in Eqs. (40)-(42).","tokens_in":9262,"tokens_out":12593,"duration_ms":146033,"concrete_test":"Rerun the 3x3 lattice with the same code and parameters for T = 5, 20, and 50 K using N independent noise seeds from zero initial conditions (e.g., N = 100). Evaluate the late-time order parameter P = (1/N_b) sum over nearest-neighbor spot pairs of <cos(theta_i - theta_j)>, together with a histogram of final FM vs AFM configurations. If P is consistently positive at 5 K and negative at 20 K with non-overlapping seed-to-seed distributions, the crossover stands; if both phases occur at both temperatures with comparable probability, the claimed temperature control is a selection artifact.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that temperature controls phase order is inferred from one \"typical\" realization at T = 5 K and one at T = 20 K (Figs. 1 and 2). The model is stochastic by construction (noise terms in Eqs. (4)-(8), (13)-(14)), and the text itself describes path-dependent vortex expulsion; in polariton lattices FM and AFM configurations are known to be near-degenerate and selected by initial conditions (Refs. [9,26]). No lattice order parameter, no realization histogram, and no statement of how many independent runs were used is given. The supporting stability argument in Sec. V does not repair this: Eq. (40) assumes phi = C psi with \"C real and small\" without derivation, and Eq. (42)'s alpha_eff < 0 is asserted by \"using this relation\" without the algebra. Thus the temperature-driven crossover, the paper's headline phase-ordering result, is not established even within the model; it may be a single-run sampling artifact.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies a two-dimensional lattice of incoherently pumped exciton-polariton condensates using a stochastic Gross-Pitaevskii equation with a non-Markovian coupling to the excitonic reservoir. After reducing the memory integral to a time-local system via a Markov-embedding auxiliary field φ, the authors numerically simulate 2×2 and 3×3 lattices at T = 5, 20, and 50 K. They report (i) an in-phase (“ferromagnetic”) phase configuration at 5 K and a checkerboard (“antiferromagnetic”) configuration at 20 K, accompanied by a change in steady-state condensate density; (ii) that the lattice retains spatial coherence at 50 K where a single spot does not condense; and (iii) a stability argument in Sec. V that in the steady state φ = Cψ with small real C leads to an effective attractive nonlinearity, α_eff < 0, which they associate with suppression of modulational instability. The abstract and conclusions frame the lattice stabilization as the main result.","tokens_in":9431,"tokens_out":7625,"duration_ms":72058,"significance":"If established, the claims would be significant: they suggest that the coupling between condensate spots via outflowing matter waves can protect phase order against thermal fluctuations well above the single-spot condensation temperature, and that reservoir temperature can switch the lattice between FM and AFM order. This would be a useful theoretical prediction for experiments using multi-spot incoherent pumping. The manuscript also has methodological strengths: Markov embedding turns the non-Markovian problem into a solvable time-local stochastic system, and the use of the scintillation index and g(1)(Δr,t) gives concrete coherence diagnostics. However, the present evidence does not yet establish the headline results.","major_comments":[{"comment":"The simulations are not reproducible because the parameter values are never stated. The model contains α_c, α_r, γ_cav, γ_exR, ρ_0, m_eff, and the temperature-dependent memory rate γ_eff(T), but none of these is given numerically; only the pump-spot radius (2.5 μm) and spacing (10 μm) are specified. Without this information the reader cannot check whether the reported temperature scales (5–50 K) are physically reasonable, whether the FM/AFM density change is within the model's parameter regime, or whether Eq. (42) indeed gives α_eff < 0. A parameter table (or a precise pointer to values used in Refs. [15,16]) is required.","section":"Section II (Eqs. (4)–(19)) and Section IV"},{"comment":"The FM-to-AFM crossover is inferred from one “typical” realization at T = 5 K and one at T = 20 K. Since the model is stochastic (noise terms in Eqs. (4), (13), (23)) and the text itself notes that phase configurations are path-dependent, a crossover claim requires an order parameter (e.g., staggered nearest-neighbor phase correlation) evaluated over an ensemble of independent runs. No number of runs, no histogram, and no error bars are given; the temperature dependence could therefore be a sampling artifact of near-degenerate FM/AFM configurations (cf. Refs. [9,26]).","section":"Section IV, Figs. 1 and 2"},{"comment":"The stability mechanism rests on an ansatz φ = Cψ with “C purely real and small” that is justified only by an appeal to numerical simulation, and the conclusion α_eff = dμ/dρ_1 < 0 is then stated without the algebra. The sign of α_eff depends on α_r, α_c, and dρ_2/dρ_1 obtained from Eq. (41); none of these is quantified. As it stands, the derivation does not establish the claimed effective attractive nonlinearity, and the connection to lattice stabilization is not made explicit.","section":"Section V, Eqs. (40)–(42)"},{"comment":"The paper attributes lattice stability to “suppression of modulational instability,” but the preceding calculation finds α_eff < 0, which is the regime where modulational instability is usually expected (Ref. [30]). No analysis of the lattice modes, no growth-rate computation, and no comparison of single-spot versus lattice MI is provided. The causal explanation is therefore speculative (as the abstract itself acknowledges with “probably”), and it should be either substantiated or clearly separated from the numerical observation.","section":"Section V, last paragraph"}],"minor_comments":[{"comment":"The word “checkboard” should be “checkerboard”; the spelling error appears in the abstract and in the text.","section":"Abstract and Section IV"},{"comment":"The first sentence contains a duplicated “the” (“using the the step-wise scheme”), and “abcense” should be “absence.”","section":"Section III"},{"comment":"The stochastic increment is written with ΔzΔy in the denominator, whereas the grid spacing elsewhere is introduced as Δx and Δy; this looks like a typo.","section":"Eq. (24)"},{"comment":"The autocorrelation ⟨η*_cav(r,t′)η_cav(r,t′)⟩ uses the same time argument on both fields; it should involve two different spacetime arguments or an explicitly corrected convention.","section":"Eq. (5)"},{"comment":"The captions say “at different time instants” but do not state which times are shown; this makes it impossible to identify the early and late stages from the figures alone.","section":"Figs. 1 and 2"},{"comment":"The ensemble average appears to be taken inside the ratio rather than on numerator and denominator separately; as written, the normalization is ambiguous. Please clarify the definition.","section":"Eq. (38)"}],"recommendation":"major_revision","confidential_remarks":"To the editor: the paper is built on the authors' prior work (Refs. [14–16]) for both the memory-kernel approximation and the numerical method, and it seems likely that the missing parameters and γ_eff(T) are given there; the authors should be asked to make the Letter self-contained. The main risk is not the model itself but the statistical support for the crossover. I recommend major revision rather than rejection because the missing evidence (ensemble statistics, a parameter table, and the explicit algebra for Eq. (42)) is well within the scope of a revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Let me cut to the chase: this paper reports two genuinely interesting numerical claims—a temperature-driven FM-to-AFM crossover in the phase pattern of a polariton condensate lattice, and lattice-coherence survival at 50 K where a single spot fails. The model is a non-Markovian stochastic Gross-Pitaevskii equation with an exponential memory kernel and Markov embedding, which the authors have used before. Applied to 2x2 and 3x3 lattices, it's a sensible extension.\n\nWhat it does well: the numerical setup is described in enough detail to follow the scheme (split-step with OU noise), the use of g(1) and the scintillation index is appropriate, and the single-spot comparison in Fig. 3 makes the lattice-stabilization point concrete. The authors are not overclaiming in the abstract, but the phrase 'probably due to suppression of the modulational instability' is honest about the tentative mechanism.\n\nThe soft spots are real and load-bearing. First, no material parameters are given: no alpha_c, alpha_r, gamma_cav, gamma_exR, rho0, or gamma_eff(T). That alone prevents any reproduction. Second, the FM/AFM crossover rests on one 'typical' realization at T=5 K and one at T=20 K, with no order parameter, no histogram over runs, and no statement of how many runs were done. Given that FM and AFM configurations are near-degenerate and selected by initial conditions/vortex expulsion (Refs. 9, 26), the temperature dependence could easily be a sampling artifact. The stress-test note is right: this is the paper's headline result and it is not established even within the model as reported. Third, the stability analysis in Sec. V is sketchy: Eq. (40) assumes phi = C psi with C real and small, without derivation, and the sign of alpha_eff is asserted via 'using this relation' without algebra. The self-referential extraction of C from the same simulations makes the argument circular-ish, though not formally circular.\n\nThe exponential memory kernel approximation is also load-bearing; if the true non-Markovian kernel deviates from exponential at 50 K, the lattice stabilization could be an artifact. That is worth flagging but not fatal; the authors do cite prior work justifying it, and the model is at least explicit.\n\nBottom line: the paper deserves serious refereeing, not desk rejection. The question is important for polaritonic devices and the model is not silly. But as it stands, the evidence is insufficient. I'd send it to peer review with a strong request for a parameter table, ensemble statistics (at minimum 20-50 runs per temperature), a proper order parameter (e.g., average staggered phase), and a complete derivation of alpha_eff. If those come back, the result could be credible.","headline":"Plausible but under-supported numerical claims about lattice-protected polariton coherence; needs parameters and ensemble statistics before it can be trusted.","tokens_in":9986,"tokens_out":2867,"would_cite":false,"duration_ms":28187,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A lattice of exciton-polariton condensates retains coherent phase order at temperatures that destroy a single condensate, with an order switch from in-phase to checkerboard as temperature rises.","keywords":["exciton-polariton condensates","condensate lattices","non-Markovian reservoir memory","stochastic Gross-Pitaevskii equation","phase ordering crossover","ferromagnetic and antiferromagnetic order","modulational instability","Markov embedding"],"falsifier":"A homodyne interferometry measurement on an incoherently pumped 2x2 or 3x3 polariton lattice over 5-50 K: observing no long-range phase pattern above the single-spot condensation temperature, or no density drop coinciding with the in-phase-to-checkerboard switch, would refute the central claim. Alternatively, a two-time pump-probe measurement showing a non-exponential reservoir memory at 50 K would undercut the model's foundation.","tokens_in":9034,"feed_emoji":"⚛️","tokens_out":9099,"duration_ms":76143,"temperature":0.7,"pith_summary":"This paper asks whether coupling between separate exciton-polariton condensate spots protects their phase coherence against thermal noise. Using a stochastic Gross-Pitaevskii equation (the standard condensate wave-function model) with non-Markovian reservoir memory, it argues that a square lattice of incoherently pumped spots keeps a macroscopically ordered phase pattern at temperatures up to 50 K, where a single spot fails to condense at all. The model also predicts a temperature-controlled crossover: low temperatures favor an in-phase (ferromagnetic) ordering of the spots, while higher temperatures favor a checkerboard (antiferromagnetic) ordering, and the switch comes with a drop in steady-state condensate density. If right, this makes polariton lattices a tunable platform for studying order in open quantum systems and for building stable phase-coherent devices.","feed_headline":"Polariton lattices keep phase order when single spots fail","feed_subtitle":"Simulations show 2x2 and 3x3 lattices stay coherent at 50 K, flipping from in-phase to checkerboard as temperature rises.","key_machinery":"The load-bearing mechanism is the Markov embedding of non-Markovian reservoir memory. The reservoir's memory kernel (the retarded self-energy) is approximated by an exponential $e^{-\\gamma_{\\rm eff}(t-t')}$ whose decay rate $\\gamma_{\\rm eff}$ grows linearly with temperature; introducing an auxiliary field $\\phi(r,t)$ converts the integro-differential equations into a closed set of time-local differential equations. This makes large-scale stochastic simulations feasible and lets the authors sweep temperature. A second piece is the coupling of lattice spots through outflowing matter waves, which locks the phases of neighboring condensates and, near the steady state, yields a negative effective nonlinearity $\\alpha_{\\rm eff}<0$; the authors argue that this suppresses the modulational instability responsible for single-spot fragmentation.","core_discovery":"The central discovery, on the paper's own terms, is that the phase order of a polariton condensate lattice is temperature-governed and collectively stabilized. In numerical solutions of the non-Markovian stochastic Gross-Pitaevskii equations, 2x2 and 3x3 lattices reach steady states with well-defined relative phases at 50 K, a temperature at which an isolated pumping spot shows no condensate. The steady-state ordering changes from in-phase (ferromagnetic) at 5 K to alternating-phase (antiferromagnetic) at 20 K, accompanied by a lower condensate density. The authors trace the protection to suppression of the modulational instability that would otherwise destroy a single condensate: near the steady state the effective nonlinearity becomes attractive, and the lattice geometry stabilizes the resulting pattern.","pith_inferences":["Beyond the paper: if the effective attractive nonlinearity is the stabilizer, varying lattice geometry or spot spacing should move the temperature window for order, a prediction that can be checked without new physics.","Beyond the paper: the ferromagnetic/antiferromagnetic crossover resembles order selection in XY-type models, so the same two-temperature story may appear in other driven-dissipative condensate lattices even where the microscopic reservoir is Markovian.","Beyond the paper: the density jump accompanying the ordering switch could serve as an experimental readout that is easier to measure than phase interferometry.","Beyond the paper: the exponential-memory approximation should be stress-tested by repeating the simulations with a realistic non-exponential kernel at 50 K; if order survives, the mechanism is genuine, and if not, the prediction is an artifact."],"forward_implications":["Incoherently pumped polariton lattices of only a few spots (2x2 and 3x3) can serve as coherent many-spot sources at temperatures where single-spot condensation is absent.","Temperature becomes a control knob for choosing between ferromagnetic and antiferromagnetic phase patterns, with a measurable density change marking the crossover.","The lattice protection mechanism operates with as few as four spots, suggesting that collective stabilization does not require an extended lattice.","The model's suppression of modulational instability in finite lattices offers a route to stable multispot condensates for polariton simulators and optical devices."],"supporting_citations":[{"why":"introduces the Markov embedding technique that turns the non-Markovian evolution into a time-local system.","marker":"[14]"},{"why":"supplies the non-Markovian reservoir model and the numerical evidence of memory-driven coherence.","marker":"[15]"},{"why":"provides the low-temperature exponential approximation and the linear temperature dependence of the memory decay rate.","marker":"[16]"},{"why":"establishes the in-phase and checkerboard configurations for polariton condensate lattices and their energy ordering.","marker":"[9]"},{"why":"supplies the outflowing-wave approximation and the phase-locking picture for coupled condensate spots.","marker":"[26]"},{"why":"showed the temperature-dependent steady-state density variation in earlier numerical simulations.","marker":"[25]"},{"why":"also demonstrated the density variation with temperature that the crossover here builds on.","marker":"[28]"},{"why":"is the basis for identifying modulational instability as the destructive mechanism that the lattice suppresses.","marker":"[30]"},{"why":"reported lattice stabilization via negative effective mass in micropillar lattices, the phenomenon extended here to incoherent pumping.","marker":"[31]"}],"fun_headline_variants":["Lattices shield polariton condensates from heat-induced disorder","Temperature flips polariton lattice from in-phase to checkerboard","Polariton arrays stay ordered at 50 K, single spots don't","Collective stability lets polariton lattices beat the heat","How polariton lattices keep phase coherence at high temps"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The model assumes the reservoir's memory of the past decays as a simple exponential whose rate rises linearly with temperature; if the real memory at tens of kelvin is not exponential, the predicted lattice order and its temperature flip could be artifacts of that assumption.","fun_headline_variants_meta":{"raw":{"variants":["Lattices shield polariton condensates from heat-induced disorder","Temperature flips polariton lattice from in-phase to checkerboard","Polariton arrays stay ordered at 50 K, single spots don't","Collective stability lets polariton lattices beat the heat","How polariton lattices keep phase coherence at high temps"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00038,"raw_usage":{"total_tokens":1968,"prompt_tokens":844,"completion_tokens":1124,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":460,"completion_tokens_details":{"reasoning_tokens":1037}},"tokens_in":460,"tokens_out":1124,"duration_ms":11035,"temperature":1.0,"reasoning_tokens":1037,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T22:05:56.522261+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A homodyne interferometry measurement on an incoherently pumped 2x2 or 3x3 polariton lattice over 5-50 K: observing no long-range phase pattern above the single-spot condensation temperature, or no density drop coinciding with the in-phase-to-checkerboard switch, would refute the central claim. Alternatively, a two-time pump-probe measurement showing a non-exponential reservoir memory at 50 K would undercut the model's foundation.","supporting_citations":[{"cited_title":"Makarov, A.A","cited_arxiv_id":null,"evidence_quote":"introduces the Markov embedding technique that turns the non-Markovian evolution into a time-local system."},{"cited_title":"Alliluev, D.V","cited_arxiv_id":null,"evidence_quote":"supplies the non-Markovian reservoir model and the numerical evidence of memory-driven coherence."},{"cited_title":"Alliluev, D.V","cited_arxiv_id":null,"evidence_quote":"provides the low-temperature exponential approximation and the linear temperature dependence of the memory decay rate."},{"cited_title":"Berloff, M","cited_arxiv_id":null,"evidence_quote":"establishes the in-phase and checkerboard configurations for polariton condensate lattices and their energy ordering."},{"cited_title":"Ohadi, R.L","cited_arxiv_id":null,"evidence_quote":"supplies the outflowing-wave approximation and the phase-locking picture for coupled condensate spots."},{"cited_title":"Alliluev, D.V","cited_arxiv_id":null,"evidence_quote":"showed the temperature-dependent steady-state density variation in earlier numerical simulations."},{"cited_title":"Helluin, L","cited_arxiv_id":null,"evidence_quote":"also demonstrated the density variation with temperature that the crossover here builds on."},{"cited_title":"Kuznetsova, D.V","cited_arxiv_id":null,"evidence_quote":"is the basis for identifying modulational instability as the destructive mechanism that the lattice suppresses."},{"cited_title":"Baboux, D","cited_arxiv_id":null,"evidence_quote":"reported lattice stabilization via negative effective mass in micropillar lattices, the phenomenon extended here to incoherent pumping."}],"review_version":1}