{"id":"955fa866-6ee0-48c0-82f9-ac7ead407310","arxiv_id":"1908.07513","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Specific lattice-matched 2D semiconductor heterobilayers are predicted to host equilibrium excitonic superfluidity with estimated critical temperatures up to about 31 K.","lead":"This paper predicts that certain stacked pairs of atom-thin semiconductor layers can host a superfluid made of electron-hole pairs, with no voltage applied. If correct, it hands experimenters a concrete list of materials for testing a long-sought quantum state at temperatures of tens of kelvin.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"LDA-only band alignments are the load-bearing input for the three candidate pairs; hybrid/GW checks are needed to confirm the type-III overlaps that set n and Tc.","rationale":"I agree with the reader's weakest_assumption: LDA band alignment is the critical step. I weighed the alternative concern that Sb2Te2Se/BiTeCl and LiAlTe2/BiTeI may suffer interband hybridization (Discussion, end of first paragraph), which the authors acknowledge and propose to mitigate with hBN; this is secondary because it leaves the Hf pair intact and is not the basis of the screening. The LDA uncertainty is more fundamental because all three candidates depend on it and no hybrid/GW verification is presented. The paper is a credible screening proposal with a clear model and phonon evidence for one pair, but the central claim 'should occur in true equilibrium' is conditional on the electronic-structure check and on synthesis. Thus my concern does not alter the reader's CONDITIONAL verdict.","tokens_in":11126,"tokens_out":7544,"duration_ms":80857,"concrete_test":"Recompute the heterostructure band alignments and band structures for Sb2Te2Se/BiTeCl, Hf2N2I2/Zr2N2Cl2, and LiAlTe2/BiTeI with HSE06 or G0W0 using the same slab geometry and vacuum; extract W and n with the charge-transfer model in Methods. If W becomes negative or falls below ~1 meV, or n drops below ~10^10 cm^-2, that pair no longer lies in the condensate region of Fig. 2b, and the central claim for that pair fails. For the Hf pair, also re-run the BCS gap equation with the updated W to check whether Tc ~31 K survives.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Results 'Identifying optimal 2D bilayer heterostructures' and Fig. 3 select the three headline pairs from LDA Kohn-Sham band alignments relative to vacuum, with final overlaps W = 67, 47, and 38 meV. These W values determine the equilibrium carrier densities n > 10^12 cm^-2 and therefore the Tc marks in Fig. 2b (Hf2N2I2/Zr2N2Cl2, Tc ~31 K). LDA is known to underestimate band gaps and can misplace band edges; a hybrid or GW calculation could change W by hundreds of meV or reverse its sign, moving a pair to the staggered-gap DEBG or high-density plasma side of Fig. 2b. The authors only state in the Discussion that the functional-dependent spread of n would change Tc estimates but not the conclusion; no hybrid/GW data for the selected pairs are shown. Because the claim is specifically that these three equilibrium platforms have optimal carrier density, this unverified electronic-structure input is the load-bearing assumption. (Separately, the Discussion admits hybridization in Sb2Te2Se/BiTeCl and LiAlTe2/BiTeI can destroy superfluidity; that is a real limitation for two pairs, though the proposed hBN spacer would need quantitative testing.)","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes that intrinsically stable 2D semiconductor heterobilayers with doubly indirect, type-III broken-gap band alignments can host spontaneous excitonic condensation in true equilibrium, without external gating. The authors combine a model RPA/BCS phase diagram for interlayer excitons with a database screening of 258 exfoliable 2D materials, identifying three lattice-matched pairs—Sb2Te2Se/BiTeCl, Hf2N2I2/Zr2N2Cl2, and LiAlTe2/BiTeI—as the most promising. For Hf2N2I2/Zr2N2Cl2 they present projected band structures and phonon spectra showing the absence of a Peierls-like instability, and they estimate critical temperatures up to about 31 K. The central claim is that these chemically specific, lattice-matched, broken-gap heterobilayers provide an optimal material platform for excitonic superfluidity.","tokens_in":11387,"tokens_out":4468,"duration_ms":47111,"significance":"If the central claim holds, this would be a valuable step toward equilibrium excitonic condensates in two-dimensional materials, with concrete, falsifiable predictions of specific material pairs and estimated critical temperatures. The paper is strongest in its systematic database screening, its explicit attention to lattice matching and thermal stability, and its phonon calculation for the primary candidate showing no detrimental CDW-like instability. The conceptual distinction from voltage-doped nonequilibrium bilayers and from monolayer Peierls-prone systems is well motivated. However, the quantitative predictions are not yet fully load-bearing because they depend on electronic-structure approximations and model parameters whose robustness is not established for the final candidates; the current evidence is suggestive rather than definitive.","major_comments":[{"comment":"The selection of Sb2Te2Se/BiTeCl, Hf2N2I2/Zr2N2Cl2, and LiAlTe2/BiTeI as optimal rests on LDA Kohn-Sham band alignments relative to vacuum, with final overlaps W = 67, 47, and 38 meV that set the carrier density n > 10^12 cm^-2 and hence the Tc marks in Fig. 2b. LDA is known to underestimate gaps and can misplace band edges by hundreds of meV; a hybrid or GW calculation could change W by enough to move a pair out of the condensate region or even reverse the sign of the overlap. The Discussion only states that a functional-dependent spread in n would change Tc estimates without changing the conclusion, but no hybrid/GW data or error bar is presented for the selected pairs. Because the headline quantitative claim (Tc up to ~31 K) is n-sensitive, this is a load-bearing uncertainty rather than a presentation issue.","section":"Results, 'Identifying optimal 2D bilayer heterostructures' and Discussion"},{"comment":"The phase diagram in Fig. 2 and the quantitative Tc values are obtained with fixed model choices d = 3 Å, κ = 1, T = 300 K for screening, and isotropic parabolic bands with equal carrier masses. The RPA potential of Eq. (2) and the BCS equations (3)-(5) are sensitive to these parameters; for real heterobilayers d is determined by vdW relaxation, κ is affected by the environment and other layers, and the bands in Fig. 4c are neither isotropic nor equal-mass. The manuscript reports no sensitivity analysis for d or κ in the main text, so the quoted Tc values and the placement of the three materials in Fig. 2b should be read as semi-quantitative; the claim that these particular materials sit in the optimal window needs robustness checks.","section":"Methods, Model Hamiltonian and BCS mean-field gap equations"},{"comment":"The paper itself states that in Sb2Te2Se/BiTeCl and LiAlTe2/BiTeI the zero modulation vector leads to band crossing and hybridization, and that single-particle hybridization can fix the phase of the order parameter and destroy superfluidity, with hBN insertion proposed as a remedy. No calculation of the interlayer tunneling matrix element or of the hBN-spaced geometry is shown. Since two of the three headline candidates are affected, the statement that all three 'emerge promising' is not yet fully supported; quantitative treatment of hybridization (or explicit demonstration that the hBN spacer suppresses it) is needed.","section":"Discussion"}],"minor_comments":[{"comment":"The statement that LDA 'fortuitously captures well' interlayer van der Waals interactions is supported by a single bilayer-graphene reference; a sentence acknowledging the limited benchmark or citing broader validation would be helpful.","section":"Methods, First-principles DFT calculation"},{"comment":"The caption refers to materials labelled (1, 2, and 3) in Fig. 3, but the mapping of these labels to the three named pairs is left implicit in the main text; a sentence listing the pairs and their final overlaps would improve readability.","section":"Fig. 2 caption"},{"comment":"The distinction between the initial overlap W0 and the final overlap W, and how each is obtained from the heterostructure calculation, could be stated more explicitly; the numerical formulas are clear but the connection to Fig. 3's ΔE values is not spelled out.","section":"Methods, Charge transfer model"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a technically competent and timely proposal, and the database-screening approach is a strength. The main risk is whether the authors can validate the LDA band alignments for the final candidates; without that, the quantitative Tc claims remain conditional. This is a standard concern for first-principles predictions and does not by itself warrant rejection, but it does require substantive additional work."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Gupta, Kutana, and Yakobson report a high-throughput search for 2D heterobilayers that could host excitonic condensation without gates. The genuinely new content is the material prediction: three lattice-matched, broken-gap pairs, with Hf2N2I2/Zr2N2Cl2 at Tc ~31 K. If the band alignments hold, this is a concrete route to equilibrium exciton condensate experiments, a real step past gated double-layer devices.\n\nThe paper does the necessary homework. The screening is systematic, the phase diagram is clearly laid out, and the phonon spectrum for Hf2N2I2/Zr2N2Cl2 directly addresses the Peierls-instability problem that plagues monolayer TiSe2. The authors also volunteer that two of the three pairs hybridize and would need an hBN spacer—more candor than usual.\n\nThe soft spot is exactly where the reader's report puts it: LDA band edges relative to vacuum are the load-bearing input. Overlaps W of 67, 47, and 38 meV set the carrier densities and hence the Tc marks. LDA is known to underestimate gaps and can shift edges by tens of meV. A hybrid or GW calculation could shrink the overlap or even flip the alignment. The authors mention functional-dependent spread in n but do not show the checks for the selected pairs. This makes the quantitative Tc numbers provisional, not the design idea. The RPA potential with d=3 Å and κ=1, and the parabolic isotropic BCS equations, are also first-pass choices; fine for a screening paper, less so for a precision prediction.\n\nI would not call this a takedown. The central claim—that these specific pairs are promising equilibrium platforms—is plausible and worth testing. The paper deserves referee time and likely publication after revision. A referee should request hybrid or GW alignments for the three candidates, or at least a quantified error bar, and ask the authors to be explicit that two of the three pairs rely on hBN spacers whose thickness and dielectric response have not been calculated.\n\nI would bring this to reading group and would cite it if I worked on 2D exciton condensates. It is not the last word, but it is a legitimate, honest one.","headline":"A credible materials-screening proposal for equilibrium excitonic condensation in 2D heterobilayers; the specific pairs are new, but LDA-only band alignments are the main uncertainty.","tokens_in":11909,"tokens_out":2502,"would_cite":true,"duration_ms":25021,"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":"Three lattice-matched 2D heterobilayers can host equilibrium excitonic superfluidity, with estimated critical temperatures up to about 31 K.","keywords":["excitonic condensation","excitonic insulator","heterobilayers","type-III band alignment","two-dimensional materials","superfluidity","BEC-BCS crossover","lattice-matched heterostructures"],"falsifier":"A clean falsifier would be a many-body calculation (hybrid functional or GW) of the stacked pairs: if the conduction-band minimum no longer lies below the valence-band maximum, or the overlap falls below roughly 1 meV, the carrier density leaves the condensate window. The experimental counterpart is angle-resolved photoemission and inverse photoemission on exfoliated stacks: a staggered (type-II) alignment for any of the three pairs would rule out the equilibrium excitonic condensate.","tokens_in":10928,"feed_emoji":"⚛️","tokens_out":11109,"duration_ms":96624,"temperature":0.7,"pith_summary":"Excitonic superfluidity—dissipationless flow of bound electron-hole pairs—has long been pursued in bulk semimetals, quantum wells, and 2D layers, but existing platforms require external doping or are destabilized by lattice distortions. This paper proposes that lattice-matched 2D semiconductor heterobilayers with a broken (type-III) gap can host the condensate in true equilibrium: electrons and holes settle on different layers and different valleys, so no gate voltage is needed and phonon-driven Peierls instability is suppressed. Screening hundreds of exfoliable 2D materials, the paper identifies three chemically specific pairs—Sb2Te2Se/BiTeCl, Hf2N2I2/Zr2N2Cl2, and LiAlTe2/BiTeI—with the right band overlap and carrier density. Solving a model Hamiltonian and BCS-type gap equations gives critical temperatures up to about 31 K. If correct, the result converts a long-sought quantum state into a materials-design problem rather than a device-engineering one.","feed_headline":"Three 2D bilayers host excitonic superfluid without gating","feed_subtitle":"Lattice-matched semiconductor pairs overlap their bands to form stable condensates up to about 31 K.","key_machinery":"The load-bearing object is the doubly-indirect broken-gap heterobilayer: two lattice-matched monolayer semiconductors stacked with van der Waals spacing, with the valence-band maximum on one layer at one valley and the conduction-band minimum on the other layer at another valley. The quantitative estimates ride on two pieces of machinery: an RPA-screened interlayer Coulomb potential used in a two-band effective-mass Hamiltonian to compute exciton binding energies and radii, and a self-consistent BCS-type gap equation whose order parameter satisfies $k_B T_c = 0.57 \\Delta$. The resulting phase diagram organizes the BKT superfluid phase at low density and the BCS regime at high density, with the optimal carrier-density window around $10^{11}$ to $10^{12}$ cm$^{-2}$.","core_discovery":"The paper's central claim is that the excitonic condensate can be the intrinsic ground state of certain 2D semiconductor heterobilayers, requiring no optical pumping and no electrostatic gating. The essential electronic configuration is a doubly-indirect band overlap: the valence-band maximum lies on one layer at one valley, while the conduction-band minimum lies on the other layer at a different valley, separating the electron and hole in both real space and momentum space. In that geometry the interlayer Coulomb attraction drives spontaneous pairing, while the spatial separation keeps interlayer electron-phonon coupling weak enough that no Peierls distortion or charge-density-wave reconstruction preempts the electronic condensate. From a high-throughput search of exfoliable 2D materials, the paper singles out Sb2Te2Se/BiTeCl, Hf2N2I2/Zr2N2Cl2, and LiAlTe2/BiTeI, with calculated band overlaps of 67, 47, and 38 meV after stacking, carrier densities above $10^{12}$ cm$^{-2}$, and a highest estimated critical temperature of about 31 K for Hf2N2I2/Zr2N2Cl2. The same phase diagram places these materials across the BEC-BCS crossover, which the paper argues can be traversed by strain or electric-field fine-tuning.","pith_inferences":["The same screening logic could be rerun with hybrid functionals or GW band structures; better band-edge accuracy might move the three named pairs but still leave other lattice-matched combinations inside the condensate window.","A natural control experiment would be to twist one of the layers: the commensurability loss should destroy the momentum-matched condensate, confirming that the doubly-indirect overlap is the operative mechanism.","If the equilibrium condensate exists, it should reveal itself in single-particle spectra as a gap $2\\Delta$ at the Fermi level and in response functions as a softened collective mode, distinguishing it from a conventional charge-density-wave state.","The same type-III design might extend to bent or cylindrical heterobilayers, where flexoelectric charge separation would supply an additional tuning knob beyond the flat-layer electric-field control the paper quantifies."],"forward_implications":["If the prediction is right, these three heterobilayers are equilibrium excitonic superfluids: the condensate forms without applied voltage, avoids leakage currents, and is mediated by purely electronic interactions rather than lattice coupling.","The estimated carrier densities and critical temperatures place the materials near the BEC-BCS crossover, so they could be used to study strong- and weak-coupling exciton condensation in a single tunable system.","Because electrons and holes occupy opposite layers, the condensate should show enhanced interlayer Josephson-like tunneling and dissipationless charge counterflow, both measurable in transport experiments.","External tuning by an electric field of order 1 V/nm shifts the carrier density by about $10^{12}$ cm$^{-2}$, which would let experiments move along the phase diagram and raise $T_c$.","The design rule—lattice-matched type-III heterobilayers with doubly-indirect band overlap—extends beyond the three specific pairs to other combinations from the 2D materials family."],"supporting_citations":[{"why":"Supplies the database of exfoliable 2D materials from which lattice-matched candidate pairs are screened.","marker":"[45]"},{"why":"Provide the RPA-screened interlayer Coulomb potential used to compute exciton binding energies and radii.","marker":"[19,25]"},{"why":"Introduces the idea that Coulomb attraction between overlapping conduction and valence bands can spontaneously form excitons, delaying metallicity.","marker":"[14]"},{"why":"Gives the BCS-like excitonic-insulator gap equations and the treatment of a finite modulation vector used in the self-consistent solution.","marker":"[15]"},{"why":"Supply the self-consistent gap equations for spatially indirect equilibrium exciton condensates used to compute $\\Delta$ and $T_c$.","marker":"[21,44]"},{"why":"Provide the BKT transition temperature relation for the two-dimensional superfluid phase.","marker":"[39,40]"},{"why":"Gives the BCS relation $k_B T_c = 0.57 \\Delta$ used to estimate critical temperatures from the order parameter.","marker":"[42]"},{"why":"Identify 1T-TiSe2 as the competing monolayer candidate whose excitonic phase the heterobilayers are designed to improve on.","marker":"[32,33]"},{"why":"Demonstrates exciton condensation in voltage-tuned 2D double layers and shows twisting destroys it, motivating the lattice-matched equilibrium design.","marker":"[31]"}],"fun_headline_variants":["Intrinsic excitonic superfluid in three 2D heterobilayers","Excitonic condensate up to 31 K without gating in 2D stacks","Lattice-matched 2D pairs condense excitons at equilibrium","No voltage needed: 2D bilayers host excitonic superfluid"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The prediction rests on LDA band alignments relative to vacuum giving the right type-III overlaps for the three pairs; if a more accurate functional moves the bands apart or into a staggered lineup, the proposed materials fall outside the condensate region.","fun_headline_variants_meta":{"raw":{"variants":["Intrinsic excitonic superfluid in three 2D heterobilayers","Excitonic condensate up to 31 K without gating in 2D stacks","Lattice-matched 2D pairs condense excitons at equilibrium","No voltage needed: 2D bilayers host excitonic superfluid"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000232,"raw_usage":{"total_tokens":1538,"prompt_tokens":1040,"completion_tokens":498,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":656,"completion_tokens_details":{"reasoning_tokens":416}},"tokens_in":656,"tokens_out":498,"duration_ms":5241,"temperature":1.0,"reasoning_tokens":416,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T12:54:09.624056+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A clean falsifier would be a many-body calculation (hybrid functional or GW) of the stacked pairs: if the conduction-band minimum no longer lies below the valence-band maximum, or the overlap falls below roughly 1 meV, the carrier density leaves the condensate window. The experimental counterpart is angle-resolved photoemission and inverse photoemission on exfoliated stacks: a staggered (type-II) alignment for any of the three pairs would rule out the equilibrium excitonic condensate.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the database of exfoliable 2D materials from which lattice-matched candidate pairs are screened."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces the idea that Coulomb attraction between overlapping conduction and valence bands can spontaneously form excitons, delaying metallicity."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the BCS-like excitonic-insulator gap equations and the treatment of a finite modulation vector used in the self-consistent solution."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the BCS relation $k_B T_c = 0.57 \\Delta$ used to estimate critical temperatures from the order parameter."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Demonstrates exciton condensation in voltage-tuned 2D double layers and shows twisting destroys it, motivating the lattice-matched equilibrium design."}],"review_version":1}