{"id":"70c49d2d-8c8c-4140-9745-15f576d2bf89","arxiv_id":"2602.03788","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"In spin-imbalanced 1D Fermi gases, spatial jumps in pairing or polarization create BCS/FFLO/normal coexistence, with FFLO correlations penetrating normal regions and a buffer FFLO zone forming at BCS-normal junctions.","lead":"Numerical Bogoliubov–de Gennes solutions show that spatially varying the pairing strength or spin polarization in a one-dimensional imbalanced Fermi gas makes BCS, FFLO, and normal phases coexist inside one box. The FFLO pair correlations leak into normal regions, and a buffer FFLO zone appears when a BCS region is forced to meet a polarized normal gas.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"1D BdG proximity claims need beyond-mean-field check; FFLO buffer may be a mean-field artifact.","rationale":"The reader identified the load-bearing assumption as the validity of BdG mean-field in 1D, and I agree this is the single most important concern. Without a beyond-mean-field calculation, the central predictions—especially the buffer FFLO phase—remain unverified. I considered other potential weaknesses: (i) the grand-canonical treatment with fixed μ and h instead of fixed total N and polarization is conventional and unlikely to change qualitative phase coexistence; (ii) the robustness to ramp width is claimed but not shown in detail, yet that is secondary because even a step-function artifact would not be central if the effect disappears with a gradual ramp; (iii) the lack of code/data is a reproducibility issue, not a correctness one. The DMRG test directly targets the gap in the argument: whether the predicted proximity and buffer phenomena are real in 1D or mean-field artifacts. Therefore the reader's CONDITIONAL verdict remains appropriate, pending such a benchmark.","tokens_in":17586,"tokens_out":5804,"duration_ms":70030,"concrete_test":"Perform ground-state DMRG (or MPS) on a 1D Fermi-Hubbard lattice with open boundary conditions, using a spatially dependent attractive interaction U(x) (U=0 on the 'normal' half) and a uniform spin-dependent potential h to realize the BCS–normal junction of Fig. 5 and the FFLO–normal junction of Fig. 2. For the same effective filling and polarization, compute the real-space pair correlation ⟨c_{i↑}c_{i↓}⟩ and its Fourier transform separately on the normal side. If a clear peak at the FFLO momentum q=k_F↑−k_F↓ and a buffer region with oscillatory pair correlations persist, the BdG claims are robust; if the peak is absent or heavily suppressed, the central claims are mean-field artifacts.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim—that an inhomogeneous population-imbalanced 1D Fermi gas exhibits a stable buffer FFLO phase at BCS–normal interfaces and FFLO-momentum proximity leakage into normal regions—rests entirely on self-consistent BdG mean-field theory with a local order parameter Δ(x). In one dimension, quantum fluctuations destroy true long-range order: the exact Gaudin–Yang solution and Luttinger-liquid analyses show only power-law correlations, and the FFLO state is not a true phase. The paper acknowledges this in Sec. IV but does not quantify how fluctuations modify the interfacial correlations. The predicted buffer FFLO could be an artifact of the mean-field treatment: BdG artificially stabilizes a nonzero Δ(x) and suppresses phase fluctuations, so the 'phase' boundaries and proximity penetration lengths extracted from Δ(x) and F(x) may not survive beyond mean-field. Since the novelty of the paper is precisely these inhomogeneous coexistence and proximity effects, the lack of any beyond-mean-field benchmark (e.g., DMRG or MPS) is the most load-bearing weakness. If fluctuations significantly suppress the finite-momentum peak in the normal region or eliminate the buffer zone, the central conclusion fails.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript studies quasi-one-dimensional two-component Fermi gases with population imbalance confined in box potentials, using self-consistent Bogoliubov–de Gennes (BdG) theory. It first maps the homogeneous phase diagram of the BCS, FFLO, and normal phases, then constructs four inhomogeneous configurations by step-like spatial profiles of either the pairing strength or the spin-polarization field: FFLO–normal (case a), FFLO–BCS by a spin-field step (case b), FFLO–BCS by a pairing-strength step (case c), and BCS–normal (case d). The central claims are that FFLO pair correlations penetrate into adjacent normal regions while retaining the finite-momentum signature q = k_F↑ − k_F↓; that BCS regions expel population imbalance and maintain a pairing plateau; and that a 'buffer FFLO phase' emerges when a BCS phase is joined to a polarized normal phase. The appendix provides an internal consistency check between the real-space FFLO modulation and the density-derived q (q_real ≈ 0.238 vs q_theory ≈ 0.234) and documents hard-wall boundary effects.","tokens_in":17838,"tokens_out":6246,"duration_ms":76656,"significance":"If the results are viewed as BdG mean-field predictions, the paper is a useful extension of superconducting proximity studies to spin-imbalanced quasi-1D atomic Fermi gases. The self-consistent calculations are carefully executed: the authors quote a convergence tolerance of 1e-5, check grid-size sensitivity, and verify robustness against smooth interface profiles. The momentum-space signature of finite-momentum pairing leaking into a normal region is a concrete, falsifiable prediction that could be probed by momentum-resolved spectroscopy. However, the significance is limited by the mean-field nature of the calculation in exactly the regime where 1D quantum fluctuations are known to be strong; the authors acknowledge this in Sec. IV but do not quantify it. No code or data is shipped, which reduces verifiability, though the numerical setup is described in sufficient detail for replication.","major_comments":[{"comment":"The phrase 'buffer FFLO phase' overstates what is computed. In the normal region of Fig. 5 the pairing interaction is set to zero, so the local order parameter Δ(x) vanishes identically by the self-consistency condition Δ = -gF. The finite-momentum peak in F(x) on the normal side is a proximity-induced pair correlation, not a thermodynamic FFLO phase. The text in Sec. III.D itself uses the more careful wording 'buffer zone supporting the FFLO behavior,' but the abstract and conclusion retain 'buffer FFLO phase.' This distinction matters because the paper's headline claim is about a new phase; revise the terminology throughout, or provide a clear operational definition of 'phase' that applies when Δ=0.","section":"Abstract and Sec. III.D"},{"comment":"The central claims are ground-state structures in 1D, where exact solutions (Gaudin–Yang), Luttinger-liquid theory, and DMRG/QMC studies show only power-law correlations and no genuine long-range order. The manuscript acknowledges this in Sec. IV but does not provide any beyond-mean-field benchmark. In particular, the 'buffer FFLO' region and the FFLO-momentum proximity peak could be BdG artifacts caused by the mean-field suppression of phase fluctuations. Because these are the paper's main new results, the authors should either (i) add a DMRG/MPS calculation for at least the BCS–normal case of Fig. 5, or (ii) explicitly and consistently frame the paper as a mean-field study and discuss, with estimates, how Luttinger-liquid fluctuations are expected to modify the proximity length and the spectral peak. As written, the load-bearing claim is left unsupported outside the mean-field approxim","section":"Sec. IV"}],"minor_comments":[{"comment":"The 'bulk' Fourier transform is not precisely defined. Please specify the spatial intervals over which F(x) is windowed for each half of the box, since the finite decay length in the normal region could affect the peak shape and position.","section":"Sec. II.C"},{"comment":"The FFLO momentum q is defined from the computed densities via k_Fσ = πρσ. Observing the Fourier peak at this q is therefore partly a self-consistency check rather than an independent prediction. The text should state this explicitly, for example by noting that q_theory is evaluated from the self-consistent density profiles rather than from input parameters.","section":"Sec. II.C and Appendix A"},{"comment":"The dashed phase boundary labeled 'best estimations' is vague. Please specify the criterion used to assign the BCS/FFLO/normal regions in the low-g, low-h corner, or remove the dashed line if it is only a guide.","section":"Fig. 1"},{"comment":"The model is strictly one-dimensional, while the abstract and introduction refer to 'quasi-one-dimensional' gases. Clarify whether transverse confinement is assumed to freeze out only the transverse modes, or whether the BdG calculation includes any transverse width.","section":"Sec. II.A"}],"recommendation":"major_revision","confidential_remarks":"The paper is a competent and clearly written BdG study, but the novelty rests on the 'buffer FFLO phase' and proximity predictions, which are not benchmarked beyond mean field. The most straightforward path to acceptance is a small DMRG/MPS calculation for one representative case (e.g., the BCS–normal configuration of Fig. 5) or a consistent reframing of all claims as mean-field predictions. The circularity in the q definition is not fatal but should be acknowledged. The absence of code or data is a verifiability issue but not, by itself, a blocker."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's my read. The paper does something concrete: it takes a quasi-1D population-imbalanced Fermi gas in a box, solves self-consistent BdG with step profiles of pairing strength or spin polarization, and maps out what happens at the interfaces. The genuinely new results are the finite-momentum FFLO pair-correlation peak that leaks into a noninteracting normal region, and the buffer FFLO zone that appears when you try to join BCS directly to a polarized normal phase. That buffer is a nice qualitative surprise, and it is not something I saw in the earlier equal-population proximity work.\n\nWhat the paper does well: the numerics are careful at the level of the mean-field calculation. They check convergence with a 1e-5 tolerance, verify grid insensitivity, compare open and periodic boundary conditions, test gradual ramps of the parameters, and in Appendix A they get good agreement between the real-space FFLO wavelength and the momentum-space peak from density-defined Fermi wavevectors. That internal consistency matters. The paper also honestly flags the 1D absence of true long-range order and cites the Gaudin–Yang and DMRG literature, so the authors are not overclaiming the status of 'phases' in one dimension.\n\nThe soft spots are real but not disqualifying. The biggest one is the one the stress test names: everything rests on BdG mean-field, and in 1D quantum fluctuations can change correlation exponents and possibly suppress or shift the finite-momentum peak in the normal region. The buffer FFLO region in particular might not survive a DMRG or MPS calculation. The authors acknowledge the limitation qualitatively but don't benchmark any inhomogeneous case against a beyond-mean-field method. That's the natural next step, not a fatal flaw. A second, smaller issue is the definition of q from the computed densities; that makes the peak position partly a self-consistency check. The reader's 'circularity burden' of 2 is about right. Third, no code or data is shipped, which makes it harder to reproduce the exact figures, though the method is standard.\n\nI'd send this to peer review. The results are new, the numerics look solid at the mean-field level, and the buffer FFLO claim is exactly the kind of thing a referee should probe. A serious referee will ask for a beyond-mean-field benchmark or at least a clear statement of why it is not feasible. If the authors add that, the paper becomes much more valuable. For me, the paper is worth citing for the proximity calculations even with the caveat.","headline":"A careful BdG study of 1D imbalanced Fermi gases with spatial parameter steps; the new proximity features are plausible but need a beyond-mean-field check before betting on them.","tokens_in":18356,"tokens_out":2199,"would_cite":true,"duration_ms":22759,"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":"Inhomogeneous pairing or polarization in a quasi-1D Fermi gas stabilizes coexisting BCS, FFLO, and normal phases, with finite-momentum FFLO correlations leaking across interfaces and a buffer FFLO phase forming between BCS and normal region","keywords":["inhomogeneous Fermi gas","population imbalance","FFLO phase","proximity effect","Bogoliubov-de Gennes","superfluid phases","cold atoms","box potential"],"falsifier":"Compute the momentum-resolved pair correlation on the normal side of a BCS–normal junction with a small spin imbalance using a numerically exact method (e.g., density-matrix renormalization group for the same box potential); the paper's central claim predicts a clear peak at q = kF↑ − kF↓ in that region. Absence of such a peak would indicate that the mean-field proximity effect and buffer FFLO phase are artifacts of the approximation.","tokens_in":17448,"feed_emoji":"⚛️","tokens_out":2703,"duration_ms":31127,"temperature":0.7,"pith_summary":"The paper asks what happens when a quasi-1D attractive Fermi gas is forced to host different phases side by side, by making the pairing strength or the spin polarization vary in space. Using self-consistent Bogoliubov–de Gennes theory, it shows that FFLO pairing correlations, characterized by a finite center-of-mass momentum, do not stop at an interface but penetrate into adjacent normal or BCS regions. When a BCS superfluid is joined to a polarized normal gas, a buffer FFLO phase appears in between, even though no bulk FFLO phase would exist in either uniform region. These predictions give concrete momentum-space signatures—peaks at the FFLO momentum in the pair-correlation spectrum—that could be seen in cold-atom experiments with box potentials.","feed_headline":"Spin-polarized Fermi gas spawns a buffer FFLO phase at BCS-normal interfaces","feed_subtitle":"Self-consistent BdG solutions show finite-momentum pair correlations penetrating normal regions, giving measurable signatures in cold-atom b","key_machinery":"The machinery is the self-consistent Bogoliubov–de Gennes (BdG) formalism for a two-component Fermi gas with spin-dependent chemical potentials and a spatially varying order parameter Δ(x). The key object is the pair correlation function F(x) = ⟨ψ↓(x)ψ↑(x)⟩, which can be nonzero even where the local pairing interaction vanishes, and its discrete Fourier transform identifies the momentum content of pairing. The load-bearing scales are the FFLO momentum q = kF↑ − kF↓ (from the mismatch of spin-resolved Fermi wavevectors) and the BCS coherence length ξ_BCS ≈ ℏ²kF/(mΔ), whose product ξ_BCS q controls whether finite-momentum correlations can penetrate an interface. The iterative solution of the B","core_discovery":"The central claim is that spatially varying pairing interaction or spin polarization in a quasi-1D population-imbalanced Fermi gas produces stable real-space coexistence of BCS, FFLO, and normal phases, with the FFLO pair correlations penetrating neighboring regions at the characteristic momentum q = kF↑ − kF↓. In the FFLO–normal interface, the pair correlation leaks into the noninteracting side while retaining its oscillatory structure and a spectral peak at q. In the FFLO–BCS interface, the order parameter transitions smoothly from oscillations to a plateau, and the BCS side stays essentially unpolarized. Most strikingly, when a BCS phase is directly joined to a polarized normal phase, the","pith_inferences":["The appearance of a buffer phase at an interface between two phases that are not adjacent on the uniform phase diagram may be a general principle: when a continuous order-parameter field cannot connect two regions without a sign change or large gradient, a third phase with intermediate momentum structure emerges to soften the mismatch.","Momentum-resolved radio-frequency spectroscopy on the normal side of such a junction could directly test the predicted finite-momentum peak in |F(k)|, which would be a sharper experimental signature than density profiles or the real-space oscillatory tail.","In strictly one-dimensional systems, true long-range order is absent and Luttinger-liquid corrections may replace the mean-field buffer region with a crossover of enhanced FFLO-type correlations; a numerically exact calculation for the same box setup would clarify whether the buffer remains a distinct phase.","The same inhomogeneous-BdG approach could be transferred to two-dimensional geometries or to Bose–Fermi mixtures, where interface-induced finite-momentum pairing might appear even when the uniform phase diagram has no such phase."],"forward_implications":["In an FFLO–normal junction, the pair correlation penetrates the normal region and displays a Fourier peak at the FFLO momentum q, giving a clear finite-momentum proximity signature absent in equal-population systems.","In an FFLO–BCS junction, the BCS side remains essentially featureless at finite momentum; the penetration of FFLO correlations is suppressed when ξ_BCS q is appreciably less than 1.","Joining a BCS superfluid to a polarized normal gas under a small but uniform spin-polarization field produces a buffer FFLO phase at the interface, so the final structure is BCS–FFLO–normal.","The spatial profiles are qualitatively robust when the step-like parameter changes are replaced by gradual ramps, as long as the transition width stays below ξ_BCS or 1/q.","The BCS phase resists population imbalance in its bulk, maintaining a flat order parameter and balanced spin densities, while FFLO and normal regions accommodate the imbalance."],"fun_headline_variants":["Buffer FFLO phase emerges when BCS meets polarized normal gas","Spatial inhomogeneity stabilizes coexisting BCS, FFLO, and normal phases","FFLO order penetrates normal gas at finite momentum","Pairing leaks into normal regions at FFLO momentum","Buffer FFLO phase at BCS-normal interface from polarization gradients"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The load-bearing premise is that BdG mean-field theory with a local order parameter Δ(x) faithfully captures ground-state pairing correlations across interfaces in a quasi-1D Fermi gas, despite the fact that one-dimensional systems have no true long-range order and are strongly affected by fluctuations.","fun_headline_variants_meta":{"raw":{"variants":["Buffer FFLO phase emerges when BCS meets polarized normal gas","Spatial inhomogeneity stabilizes coexisting BCS, FFLO, and normal phases","FFLO order penetrates normal gas at finite momentum","Pairing leaks into normal regions at FFLO momentum","Buffer FFLO phase at BCS-normal interface from polarization gradients"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00081,"raw_usage":{"total_tokens":3409,"prompt_tokens":781,"completion_tokens":2628,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":525,"completion_tokens_details":{"reasoning_tokens":2550}},"tokens_in":525,"tokens_out":2628,"duration_ms":17365,"temperature":1.0,"reasoning_tokens":2550,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T04:48:58.060811+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the momentum-resolved pair correlation on the normal side of a BCS–normal junction with a small spin imbalance using a numerically exact method (e.g., density-matrix renormalization group for the same box potential); the paper's central claim predicts a clear peak at q = kF↑ − kF↓ in that region. Absence of such a peak would indicate that the mean-field proximity effect and buffer FFLO phase are artifacts of the approximation.","supporting_citations":[],"review_version":1}