{"id":"b02ea33a-c821-498c-a5ed-aee47ed0e3e2","arxiv_id":"2608.08178","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Quantum Monte Carlo predicts a Fermi supersolid with coexisting checkerboard order and pair superfluidity in the half-filled square-lattice three-color Hubbard model with weak third-color coupling.","lead":"This paper predicts a supersolid state, where crystalline density order coexists with superfluid pairing order, in a simple three-color Fermi-Hubbard model on a square optical lattice. The predicted state could be tested in existing ultracold lithium experiments, since it needs only on-site attractions and no long-range interactions.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The supersolid claim relies on long-range U(1) CSF order in 2D at T>0, which Mermin–Wagner forbids; the paper never addresses quasi-long-range order.","rationale":"The paper's central claim is that the attractive three-color Fermi-Hubbard model on a 2D square lattice hosts a Fermi supersolid at T=1/12, with coexisting CDW and CSF orders. The CDW breaks a discrete Z2 symmetry and can order at finite T in 2D, so that part is in principle allowed. The CSF, however, is a U(1) superfluid order parameter; in two dimensions at finite temperature, continuous symmetry breaking is forbidden by Mermin–Wagner–Hohenberg for short-range interactions. The paper does not mention this theorem, nor does it discuss the possibility of algebraic (quasi-long-range) superfluid order. Instead, P_CSF in Eq. (4) is extrapolated to a nonzero thermodynamic limit, which is inconsistent with the theorem. The finite-size extrapolations in the SM rely on four system sizes and linear/quadratic fits that cannot reliably distinguish algebraic from exponential decay. The correlation-ratio crossing method, while standard for locating quantum phase transitions, requires a BKT-specific finite-size scaling analysis in this context; the simple aL^{-b} fit in Eq. (S1) is not justified for a Kosterlitz–Thouless transition. The reader's weakest assumption correctly identified this same issue, and our reading confirms it is the most load-bearing concern. The supersolid claim as literally stated ('coexistence of CSF and CDW orders' with nonzero order parameters) is therefore unsupported. However, a revised interpretation in terms of quasi-long-range CSF order coexisting with true CDW order could still be valuable and is experimentally relevant, so a conditional verdict remains appropriate. The paper should be revised to address Mermin–Wagner, to perform a proper BKT analysis, and to report correlation functions rather than extrapolated order parameters for the CSF channel.","tokens_in":13219,"tokens_out":4601,"duration_ms":45459,"concrete_test":"Compute the raw CSF correlation function C(r)=⟨Δ†(r)Δ(0)⟩ at the largest lattice L=12, T=1/12, at the claimed supersolid point (|U|=2.0, |U'|=0.05); plot ln C(r) versus r and versus ln r. If the decay is exponential, the CSF order is absent and the supersolid claim fails. If it is algebraic, the state is quasi-long-range ordered and the paper must be revised to describe a BKT-type supersolid, abandoning the P_CSF thermodynamic-limit extrapolation as an order parameter.","verdict_should_be":"UNCHANGED","load_bearing_attack":"In 2D at finite temperature, continuous U(1) symmetry cannot be spontaneously broken for short-range interactions (Mermin–Wagner–Hohenberg). The order parameter P_CSF defined in Eq. (4) as sqrt(S_CSF(L,Γ)/N) therefore must tend to zero in the thermodynamic limit. Yet the paper extrapolates P_CSF to nonzero values using 1/L fits on L=6, 8, 10, 12 (SM §S4) and uses these to claim a supersolid phase with coexisting CDW and CSF order. No discussion of algebraic (BKT) order appears; the correlation-ratio crossing in Eq. (3) is used as if it located a conventional second-order transition, but for a BKT transition the finite-size scaling is different and a nonzero P_CSF extrapolation is invalid. The finite-size data cannot distinguish exponential from power-law decay of the CSF correlation function, so the central claim of a true finite-temperature 2D supersolid is unsupported as stated. The mean-field treatment in Eqs. (6)–(8) spontaneously breaks U(1) at T=1/12, but mean-field ignores the thermal fluctuations that restore the symmetry in 2D, so the mean-field coexistence does not rescue the interpretation.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies the half-filled attractive three-color Fermi-Hubbard model on a square lattice at T=1/12 using sign-problem-free determinant quantum Monte Carlo. It defines CDW and CSF structure factors, extrapolates the corresponding order parameters to the thermodynamic limit, and constructs a phase diagram with liquid-like, CDW, and 'supersolid' regions. The central claim is that weak coupling of a third color to a strongly paired color-1/color-2 subsystem produces a lattice Fermi supersolid with coexisting CDW and color-superfluid order. A mean-field argument is offered as a mechanistic explanation, and triple-occupancy data are used to argue experimental feasibility.","tokens_in":13481,"tokens_out":4868,"duration_ms":54122,"significance":"If correct, the result would identify a new finite-temperature supersolid phase in a simple fermionic lattice model that is directly accessible to current ultracold-atom experiments, going beyond the bosonic and Rydberg-based proposals. The manuscript has notable strengths: the DQMC simulations are sign-free, the order parameters are measured directly rather than inferred from a fitting target, finite-size data are reported for several lattice sizes, and the triple-occupancy analysis addresses a real experimental constraint. However, the central thermodynamic interpretation is undermined by the fact that the putative CSF order is a continuous U(1) order in two dimensions at finite temperature, where true long-range order is forbidden and only quasi-long-range (BKT-type) order is possible. The paper neither acknowledges this nor performs the scaling analysis needed to distinguish the two scenarios.","major_comments":[{"comment":"The nonzero thermodynamic-limit values of P_CSF obtained from 1/L polynomial extrapolations are incompatible with the Mermin-Wagner-Hohenberg theorem: at T=1/12>0 in two dimensions with short-range interactions, the continuous U(1) symmetry associated with the 12-pair superfluid cannot be spontaneously broken, so P_CSF, defined in Eq. (4) as sqrt(S_CSF(L,Γ)/N), must tend to zero in the thermodynamic limit. The finite-size data for L=6, 8, 10, and 12 cannot distinguish exponential decay from algebraic decay of the CSF correlation function. The authors should either perform a BKT-style analysis of the CSF correlations (e.g., extracting the power-law exponent at the transition) or explicitly reframe the claim as quasi-long-range order, and revise the abstract, the phase diagram in Fig. 1, and the word 'supersolid' accordingly.","section":"Three-color Fermi Supersolidity, Eqs. (3)-(4); SM §S4"},{"comment":"The correlation-ratio crossing of R_CSF is interpreted as locating a conventional continuous transition into a CSF-ordered phase. For a two-dimensional U(1) system at finite temperature, the expected transition is of BKT type, where the correlation ratio exhibits different finite-size behavior and the ordered phase is characterized by algebraic correlations rather than a nonzero order parameter. Without an analysis of the real-space or momentum-space decay of the CSF correlations, the extrapolated critical couplings |U'_c|=0.14(3), 0.10(2), and 0.09(3) do not establish a phase boundary to a state with true off-diagonal long-range order.","section":"SM §S1, Eq. (S1), and Fig. S1"},{"comment":"The mean-field Hamiltonian H_MF_12 is constructed by hand with nonzero order parameters Δ and D inserted as inputs, so the coexistence found in the self-consistent solution is imposed rather than emergent. Moreover, the mean-field treatment spontaneously breaks the U(1) symmetry at finite temperature in two dimensions, so it cannot rescue the DQMC extrapolation or serve as evidence for a true supersolid. The section is useful as a qualitative explanation of why pairing survives in the presence of CDW order, but it should be repositioned accordingly and should not be used to justify the thermodynamic-limit interpretation of P_CSF.","section":"Origin of robust supersolidity, Eq. (6)"}],"minor_comments":[{"comment":"There are several typographical errors, including 'relaization' in the Introduction and 'clor 3' in the Model section; these should be corrected.","section":"Introduction"},{"comment":"The caption contains 'fuction' instead of 'function', and the notation Δ and D is easy to confuse with the density of states D(ε) in Eq. (8) and Fig. 4; a clearer notation would improve readability.","section":"Fig. 3 caption"},{"comment":"The liquid-like phase is defined by the arbitrary threshold D123<0.03, and the crossover criterion in Eq. (S5) also uses an ad hoc threshold. These thresholds should be justified or shown not to affect the qualitative phase boundaries.","section":"Fig. 1 and SM §S4"},{"comment":"The Hamiltonian is written with U_αβ<0, but the text and figures use |U| and |U'|; please make the sign convention consistent throughout, especially in Eq. (6) where |U|D multiplies the density term.","section":"Model and Eq. (1)"},{"comment":"The phase diagram is generated by interpolating discrete DQMC data points, but the interpolation procedure and the absence of error bars on the interpolated color map are not described; this should be clarified in the Supplemental Material.","section":"Fig. 1"}],"recommendation":"major_revision","confidential_remarks":"The manuscript contains potentially interesting finite-temperature coexistence physics, but the central claim of a true supersolid with long-range CSF order is not supported by the present analysis. The necessary fixes—a BKT-style scaling analysis of the CSF correlations and a careful reframing of the phase as quasi-long-range ordered—are substantial but within the manuscript's scope, so I recommend major revision rather than rejection. I also note that the paper leans heavily on the authors' prior work (Refs. [39-43]); while the citations are appropriate, the novelty relative to those works should be stated more explicitly."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This paper reports DQMC simulations of the three-color attractive Hubbard model on a square lattice at half-filling, with U12=U and U13=U23=U'<U. The new physics is the claim of a supersolid phase where a charge-density wave coexists with color-superfluid order for weak U'. That coexistence region, mapped in the (|U|,|U'|) plane, is not present in earlier work on SU(3)-symmetric or pi-flux lattices. The numerics are standard, the sign-free argument is plausible, and the finite-size data show clear correlation-ratio crossings and consistent order-parameter behavior. Credit is due for the concrete experimental mapping to 6Li and for showing that triple occupancy remains low in the proposed regime.\n\nThe problem is that the central claim, as stated, cannot be true in the thermodynamic limit at T=1/12 on a 2D lattice. A continuous U(1) symmetry cannot be spontaneously broken at finite temperature in two dimensions (Mermin–Wagner–Hohenberg). The CSF order parameter P_CSF defined in Eq. (4) is an ODLRO-like quantity, and it must vanish in the thermodynamic limit. The authors extrapolate it to a nonzero value using 1/L fits and treat the correlation-ratio crossing as a conventional second-order transition. They never address algebraic/BKT order or the distinction between true long-range order and quasi-long-range order. The finite-size data cannot distinguish exponential from power-law decay. So the phrase \"lattice Fermi supersolid state\" with coexisting CSF and CDW orders overstates what the data show. The mean-field section does not help: it imposes both order parameters by hand and ignores the thermal fluctuations that restore U(1) symmetry.\n\nThis is a load-bearing flaw, not a minor omission. That said, the work still has value: it identifies a parameter regime where strong CDW order coexists with robust superfluid correlations at finite size, and that may well be a finite-temperature precursor or a BKT-type superfluid. The paper would need a serious revision: either extrapolate to T=0 to establish true ODLRO, or reframe the claim as quasi-long-range superfluid order coexisting with an Ising CDW. The experimental proposal is reasonable and the interaction regime is accessible.\n\nWho is this for? Ultracold-atom theorists and experimentalists working on three-component Fermi gases. It deserves a serious referee because the question is important and the work is mostly careful, but the conclusion needs major surgery before it can be taken at face value.\n\nMy recommendation: send it to review, with a request for major revision addressing the finite-temperature symmetry issue.","headline":"A promising finite-size observation of coexisting CDW and superfluid correlations, but the thermodynamic-limit supersolid claim violates Mermin–Wagner and needs reframing.","tokens_in":14036,"tokens_out":3908,"would_cite":false,"duration_ms":38539,"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 half-filled attractive three-color Fermi-Hubbard model on a square lattice hosts a Fermi supersolid when a weakly coupled third color lifts the degeneracy between charge-density-wave and color-superfluid order.","keywords":["three-color Fermi-Hubbard model","fermionic supersolid","color superfluid","charge density wave","determinant quantum Monte Carlo","trion CDW","square optical lattice","attractive Hubbard model"],"falsifier":"Recompute the CSF correlation ratio and the CDW order parameter on lattices of linear size 16 and 18 at $T=1/12$ for $|U|=2.0$ and $|U'|=0.05$ to $0.15$; if the correlation-ratio crossing disappears or the extrapolated order parameters tend to zero, the claimed thermodynamic supersolid is not there. A complementary check is the superfluid stiffness: a genuine two-dimensional finite-temperature superfluid coexisting with CDW should exhibit quasi-long-range algebraic order rather than a conventional nonzero order parameter.","tokens_in":12922,"feed_emoji":"⚛️","tokens_out":8906,"duration_ms":85790,"temperature":0.7,"pith_summary":"The paper claims that a Fermi supersolid—simultaneous charge-density-wave (CDW) and color-superfluid (CSF) order—exists in the half-filled attractive three-color Fermi-Hubbard model on a square lattice. The mechanism is asymmetric attraction: colors 1 and 2 form moderately to strongly bound pairs, and a weakly coupled third color acts as a sublattice pinning field that lifts the degeneracy between CDW and CSF. Using determinant quantum Monte Carlo at $T=1/12$ on lattices $L=6$ through $12$, the authors find a finite interaction window in which both order parameters survive extrapolation to the thermodynamic limit. If correct, this gives a concrete cold-atom route to fermionic lattice supersolidity using only on-site attractive interactions and a simple square optical lattice.","feed_headline":"Three-color Fermi gas predicted to host supersolid phase","feed_subtitle":"Weakly coupled third color turns a strongly paired two-color subsystem into coexisting crystalline and superfluid order.","key_machinery":"The load-bearing objects are the color-superfluid correlation ratio $R_{\\mathrm{CSF}}$ of Eq. (3), whose finite-size crossings locate the CSF phase boundary, and the CDW order parameter $D_{123}$ of Eq. (2), whose $1/L$ extrapolation maps the density order. The microscopic mechanism is captured by the mean-field quasiparticle spectrum $E_k=\\sqrt{\\epsilon_k^2+(|U|\\Delta)^2+(|U|D)^2}$: the CDW field $D$ couples the $\\vec k$ and $\\vec k+\\vec Q$ sectors, and although increasing $D$ opens a gap, the density of states stays peaked at the gap edges, so the pairing kernel $I(\\epsilon)$ remains large and CSF order survives. The third color's weak attraction acts as the sublattice potential that lifts the SU(2) pseudospin degeneracy and selects the coexistence direction.","core_discovery":"The central discovery, stated as the authors would state it, is that weak coupling $U'$ between color 3 and the paired color-1/color-2 subsystem converts the mutually exclusive superfluid and density-wave orders of the isolated two-color model into coexisting orders. The CDW background of 12-pairs produces a sublattice potential that slightly modulates the color-3 density, which in turn deepens the potential for the pairs; this enhances CDW while the localized wells weaken but do not destroy the color superfluid. The quantitative evidence is the coexistence of an extrapolated CDW order parameter $D_{123}$ with a finite CSF correlation ratio over a window around $|U|=2.0$--$3.5$ and $|U'|$ up to roughly $0.15$, with a crossover from 12-CDW to trion 123-CDW inside the supersolid. The authors support the mechanism with a mean-field theory in which the CDW field opens a gap while the density of states near the gap edges keeps the pairing kernel substantial.","pith_inferences":["The authors leave implicit that the same 'third color as pinning bath' idea could be tested in other geometries or fillings; a natural next calculation is whether the gap-edge density-of-states effect survives on lattices with different Van Hove singularities.","Because the numerical case rests on finite-size extrapolation at one temperature, a decisive extension would be a superfluid-stiffness or algebraic-order analysis of the CSF channel, which the paper does not provide.","An immediate experimental discriminator is simultaneous measurement of the $(\\pi,\\pi)$ density structure factor and pair coherence in one cloud; both appearing together in the predicted interaction window would confirm the coexistence, and either alone would falsify the supersolid interpretation.","The mean-field picture suggests a spectral prediction: the single-particle gap should open with CDW order while pairing remains, which could be probed by momentum-resolved spectroscopy as a sharper test of the proposed mechanism."],"forward_implications":["A fermionic supersolid can be reached in a half-filled square lattice with only on-site attractions, without the off-site interactions required by bosonic and Rydberg supersolid proposals.","The supersolid occupies a finite pocket of the phase diagram, bounded by the loss of CSF order at $|U'_c|\\approx0.09$--$0.14$ for $|U|=2.0$--$3.5$ and by thermal suppression at large $|U|$.","Within the supersolid, increasing $|U'|$ continuously transfers spectral weight from superfluid to density order, and eventually the density order crosses over from 12-CDW to trion 123-CDW.","The predicted regime is experimentally plausible: at $|U|=2.0$ the triple occupancy stays below 0.20, within the stability window of current three-color lattice experiments.","If the coexistence is genuine, the same pairing plus pinning mechanism should be visible in both the CDW structure factor and the pair coherence of the same experimental cloud, giving two independent detection channels."],"supporting_citations":[{"why":"Supplies the experimentally realized half-filled square-lattice three-color Fermi-Hubbard system and the low-three-body-loss stability window the prediction targets.","marker":"[34]"},{"why":"Establishes the doublon-versus-trion competition and the CSF/CDW ordering tendencies that the supersolid proposal builds on.","marker":"[35–43]"},{"why":"Provides the sign-problem-free determinant quantum Monte Carlo formulation used to generate all numerical data.","marker":"[39–42, 47]"},{"why":"Documents three-body loss rates and 6Li Feshbach properties used to argue the predicted supersolid regime is experimentally stable.","marker":"[44–46]"},{"why":"Supplies the finite-size scaling precedent behind the linear and quadratic $1/L$ extrapolations of order parameters.","marker":"[52]"}],"fun_headline_variants":["Three-color Fermi gas shows supersolid order","Weak third color flips two-color gas to supersolid","Supersolid predicted in simple three-color lattice","Attractive three-color model yields lattice supersolid"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the determinant quantum Monte Carlo data at $T=1/12$ on lattices $L=6$, $8$, $10$, and $12$, extrapolated with linear or quadratic curves in $1/L$, correctly identify a thermodynamic supersolid in two dimensions at finite temperature rather than a quasi-long-range algebraic state.","fun_headline_variants_meta":{"raw":{"variants":["Three-color Fermi gas shows supersolid order","Weak third color flips two-color gas to supersolid","Supersolid predicted in simple three-color lattice","Attractive three-color model yields lattice supersolid"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000245,"raw_usage":{"total_tokens":1523,"prompt_tokens":918,"completion_tokens":605,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":534,"completion_tokens_details":{"reasoning_tokens":545}},"tokens_in":534,"tokens_out":605,"duration_ms":6635,"temperature":1.0,"reasoning_tokens":545,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T00:19:28.686365+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the CSF correlation ratio and the CDW order parameter on lattices of linear size 16 and 18 at $T=1/12$ for $|U|=2.0$ and $|U'|=0.05$ to $0.15$; if the correlation-ratio crossing disappears or the extrapolated order parameters tend to zero, the claimed thermodynamic supersolid is not there. A complementary check is the superfluid stiffness: a genuine two-dimensional finite-temperature superfluid coexisting with CDW should exhibit quasi-long-range algebraic order rather than a conventional nonzero order parameter.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the experimentally realized half-filled square-lattice three-color Fermi-Hubbard system and the low-three-body-loss stability window the prediction targets."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the finite-size scaling precedent behind the linear and quadratic $1/L$ extrapolations of order parameters."}],"review_version":1}