{"id":"ac6f62a1-8f3a-4997-bcee-b47cd31d311f","arxiv_id":"1908.02431","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"First-principles simulations predict that a free-standing hexagonal 2D sodium sheet is a stable 2D electron gas with a doping-independent Bloch-Grüneisen temperature near 50 K and Wiedemann-Franz behavior.","lead":"This paper uses computer simulations to predict the electrical and thermal behavior of a hypothetical two-dimensional sheet of sodium atoms. It finds the sheet's resistance can be tuned by doping and that it might outperform graphene at high temperatures.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Claimed doping-independent Θ_BG ≈ 50 K conflicts with the paper's own formula Θ_BG ∝ k_F ∝ √n; the reported pinning may be a fitting artifact rather than a physical prediction.","rationale":"The reader's weakest assumption was the rigid-band shift, and that is a real limitation. However, for a nearly-free-electron half-filled band, a rigid-band shift is a standard first approximation, so it is not by itself the most damaging risk. The sharper issue is internal: even under the rigid-band model, the claimed doping-independence of Θ_BG contradicts the paper's own definition Θ_BG = 2ℏk_F v_s/k_B, since k_F varies with carrier density. The authors present the 50 K crossover as a visual reading of Fig. 3(b) and give fitted A and B parameters only for one (pure) case, so the 'independence' claim lacks the numerical support needed for a headline. The proposed check directly settles whether the physical Θ_BG is actually pinned. The computational machinery (EPW+BTE, phonon stability, MD stability check) is appropriate and the remaining transport numbers (resistivity crossover with graphene, WF law) are plausible, so the paper warrants a conditional rather than outright rejection. I therefore keep the reader's CONDITIONAL verdict unchanged, with the explicit condition that the Θ_BG extraction be validated by the phonon-based test.","tokens_in":8855,"tokens_out":17202,"duration_ms":190965,"concrete_test":"Using the DOS-based n(E_F) relation from Section III, compute k_F for each E_F in Fig. 3, read the LA/TA sound velocity v_s from the phonon dispersion in Fig. 1(a), and evaluate Θ_BG^phys = 2ℏk_F v_s/k_B. Compare this analytic sequence with the fitted 50 K crossover. If Θ_BG^phys varies by more than about 20% across the doping range (especially at E_F = ±1.5 eV) while the fitted crossover does not, then the density-independent Θ_BG headline is unsupported and should be revised.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central transport claim includes a Θ_BG that is 'almost constant at 50 K, independent of the type or density of the charge carriers' (Abstract, Section III). This is the same Θ_BG the authors define as 2ℏk_F v_s/k_B (Section I). Under their own rigid-band model, k_F ∝ √n, so the definition demands that Θ_BG vary with doping; the density range in Fig. 3 spans roughly an order of magnitude from hole doping E_F = -1.5 eV to electron doping +1.5 eV, which would shift Θ_BG by a factor of about 3 unless the sound velocity compensates, and no mechanism is given. The paper also notes that BG theory is not satisfied at E_F = 1.5 eV, yet the abstract asserts density/type independence over that same range. The 50 K value is extracted from power-law fits that are only shown for one (pure) case, with no convergence test or error estimate for the crossover. Thus the pinned Θ_BG is not robustly supported, and this matters because it is one of the three headline predictions.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript uses DFT, DFPT, and Boltzmann transport-equation calculations to argue that a free-standing hexagonal monolayer of sodium is thermodynamically stable, behaves as a 2D electron gas, and has phonon-limited electrical and thermal transport. It reports that the intrinsic resistivity of electron-doped 2D Na (E_F = 0.5 eV) is about 1.4 times that of graphene at room temperature but drops below graphene beyond 450 K, that the Bloch-Grüneisen temperature is pinned near 50 K independently of carrier type and density, that the electronic thermal conductivity at 300 K is about 1.24 times that of bulk Na, and that the Wiedemann-Franz law holds with a Lorenz number close to 2.41e-8 V^2/deg^2.","tokens_in":9100,"tokens_out":5495,"duration_ms":57775,"significance":"If the results are correct, 2D Na would be a new stable elemental 2D metal whose transport can be gate-tuned and is comparable to graphene's phonon-limited resistivity, which is of interest for nanoscale electronics and thermoelectrics. Strengths of the paper include the use of a standard first-principles workflow, the explicit comparison between Allen's transport model and an exact BTE solution, the stability checks via phonons and molecular dynamics, and the presentation of the Lorenz-number derivation in the supplementary material. The main risk is that the quantitative predictions rest on rigid-band doping and on resistivity fits whose uncertainty is not quantified; the claimed doping-independence of Theta_BG is not demonstrated and is in tension with the definition of Theta_BG.","major_comments":[{"comment":"The claim that Theta_BG is pinned at about 50 K independent of carrier type and density is not supported by the data shown. The fits to A T^4 and B T are presented only for the undoped case, with A = 2.9e-8 K^-4 and B = 3.3e-3 K, and the crossover is then read off as about 50 K. For the doped cases in Fig. 3(b), only visual inspection is offered; no per-doping fits or confidence intervals are given. More fundamentally, the definition Theta_BG = 2 hbar k_F v_s / k_B in Section I implies Theta_BG proportional to k_F. In a 2DEG with constant DOS, n is proportional to E_F, so k_F is proportional to sqrt(n); across the plotted range E_F = -1.5 to +1.5 eV, one would expect Theta_BG to vary by roughly sqrt(1.5/0.5) ~ 1.7 unless v_s changes to cancel, and no such mechanism is identified. The text itself states that at E_F = 1.5 eV the Fermi surface is no longer simple and Bloch-Grüneisen theory is not satisfied, so the abstract's 'independent of the type or density' is overbroad. Please either restrict the claim to the single-band regime |E_F| <= 0.79 eV with fits for each doping level, or give a quantitative explanation for the near-constancy.","section":"Section III, Fig. 4(a) and Abstract"},{"comment":"The entire carrier-density dependence, including the resistivity reduction at E_F = 0.5 eV and the crossover with graphene, rests on the rigid-band approximation: the Fermi level is shifted while electronic bands, phonons, and electron-phonon matrix elements are kept fixed. The manuscript does not test this approximation, and the text does not acknowledge its limitations directly in the doping discussion. Please validate the rigid-band assumption at least for a few representative densities, for example by recomputing band structures and electron-phonon matrix elements self-consistently, or clearly state the expected error and the consequences for the predicted 1/E_F and 1/E_F^2 power laws and for the claimed Theta_BG pinning.","section":"Section III, first paragraph of Results"},{"comment":"The comparison with graphene is not reproducible as presented. The text says 'we have presented the phonon limited temperature dependence electrical resistivity of graphene' but does not state whether the graphene curve was computed in this work or taken from the literature. To support the quantitative claims that the doped 2D Na resistivity is 'about 1.4 times larger' and 'falls below the latter 450 K onwards', please specify the provenance of the graphene resistivity data, including the computational method, parameters, and any references, or compute it with the same workflow for a controlled comparison.","section":"Section III, Fig. 5(a)"},{"comment":"No convergence tests or numerical error estimates are reported. The manuscript mentions a fine grid of 200 x 200 x 1 and Fig. 4(a) labels 'different grids', but the differences between grids are not quantified. Since the central numbers are ratios of resistivities and a fitted crossover temperature, please report the convergence of the resistivity and thermal conductivity with respect to k-grid, q-grid, smearing, and vacuum thickness, and provide error bars on the fitted coefficients A and B and on the resulting Theta_BG.","section":"Section II and Figs. 3-4"}],"minor_comments":[{"comment":"The word 'satbility' should be 'stability'.","section":"Fig. 1 caption"},{"comment":"The phrase 'The gird of 200 x 200 x 1' should be 'The grid of 200 x 200 x 1'.","section":"Section III, Fig. 4 discussion"},{"comment":"The sentence 'Since we have enhanced conductivity for E_F = 0.5 eV' would be clearer as 'Since 2D Na has enhanced conductivity at E_F = 0.5 eV'.","section":"Section III, discussion of E_F = 0.5 eV"},{"comment":"The statement that Theta_BG is independent of the type and density of carriers is inconsistent with the body-text caveat that at E_F = 1.5 eV the Fermi surface is no longer simple and Bloch-Grüneisen theory does not apply; please reconcile the abstract and conclusions with the single-band limitation.","section":"Abstract and Conclusions"},{"comment":"The sentence 'Theta_D in bulk systems is equivalent to Theta_BG in 2D, systems with low electron density' is confusing; please clarify the intended analogy between Debye temperature and Bloch-Grüneisen temperature.","section":"Section I"}],"recommendation":"major_revision","confidential_remarks":"The paper is suitable in scope for cond-mat.mtrl-sci, but the current abstract overstates the doping independence of Theta_BG relative to the body text. The authors should also strengthen the convergence and rigid-band tests before resubmission; otherwise the quantitative headline claims will remain difficult for readers to assess."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a solid DFPT/BTE transport paper with one genuinely new set of results and one overreaching claim. The transport quantities—resistivity, thermal conductivity, Lorenz number—are not in the earlier stability paper by Nevalaita and Koskinen, so that part is new. The methodology is standard: DFT/DFPT, EPW Wannier interpolation, BTE solved with Allen's model and an exact solver, and two grids tested. Good that they compare the two approaches and show the T⁴ and T regimes.\n\nThe soft spot is the Θ_BG claim. They define Θ_BG = 2ħk_F v_s/k_B, and in their own 2DEG picture k_F ∝ √n. Over the E_F range shown (roughly −1.5 to +1.5 eV, densities spanning about an order of magnitude), that definition would move Θ_BG by a factor of about three unless v_s compensates, and no mechanism is given. The observed 'pinning' at 50 K is read from power-law fits to resistivity, shown explicitly for only one case, with no error estimate or convergence test for the crossover. So the abstract's 'almost constant at 50 K, independent of type or density' does not follow from the calculation. The paper itself notes that Bloch–Grüneisen theory is not satisfied at E_F = 1.5 eV, yet the abstract includes that density. That is an internal inconsistency worth fixing.\n\nOther limitations are addressable rather than fatal: the doping is treated with a rigid-band shift; there is no code/data release, no convergence data, and no experimental validation. The comparison with graphene is reasonable for phonon-limited resistivity but should not be read as a device-level claim.\n\nIn short, this is a useful computational contribution for people working on 2D metals and electron–phonon transport. It deserves a serious referee, but the Θ_BG headline needs to be either derived properly or removed.","headline":"Competent first-principles transport study of a predicted 2D sodium sheet, but the headline 'doping-independent Θ_BG ≈ 50 K' contradicts the paper's own Θ_BG ∝ k_F definition and is likely a fitting artifact.","tokens_in":9642,"tokens_out":2683,"would_cite":false,"duration_ms":30816,"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":"This paper predicts that a free-standing single layer of sodium is a thermodynamically stable 2D metal whose electrical resistivity is tunable by doping and drops below graphene above 450 K.","keywords":["2D sodium","two-dimensional electron gas","Bloch-Grüneisen temperature","electron-phonon coupling","intrinsic electrical resistivity","Lorenz number","Boltzmann transport equation","first-principles calculations"],"falsifier":"A gate-tuned resistivity measurement on a free-standing monolayer sodium sheet would settle the claims: if the crossover below graphene's resistivity near 450 K does not appear, or if the Bloch-Grüneisen temperature shifts by more than a few kelvin with carrier density, the rigid-band clean-Fermi-surface picture fails. A first-principles calculation that re-optimizes the doped band structure self-consistently, instead of rigidly shifting $E_F$, would test the same assumption without experimental access.","tokens_in":8646,"feed_emoji":"⚡","tokens_out":9581,"duration_ms":81428,"temperature":0.7,"pith_summary":"This paper argues that a free-standing single layer of sodium is thermodynamically stable and behaves as a two-dimensional electron gas (2DEG) with a clean, circular Fermi surface. Using first-principles density-functional theory and an accurate solution of the Boltzmann transport equation, it predicts that the intrinsic electrical resistivity is dominated by electron-phonon scattering and can be tuned by gate doping. For electron doping at a Fermi energy of $E_F = 0.5$ eV, the resistivity is about 1.4 times that of graphene at low temperature but falls below graphene above 450 K. The Bloch-Grüneisen temperature, the crossover from $T^4$ to linear-$T$ resistivity, stays pinned at about 50 K regardless of carrier type or density. The electronic thermal conductivity of the pure sheet is about 1.24 times that of bulk sodium at 300 K, and the Wiedemann-Franz law holds with a Lorenz number of $2.41 \\times 10^{-8}\\ \\mathrm{V}^2/\\mathrm{deg}^2$.","feed_headline":"2D sodium resistivity dips below graphene above 450 K","feed_subtitle":"First-principles study predicts a stable single-layer metal whose resistivity is gate-tunable and beats graphene at high temperature.","key_machinery":"The load-bearing objects are a clean circular Fermi surface and the Bloch-Grüneisen temperature $\\Theta_{BG} = 2\\hbar k_F v_s / \\kappa_B$, which marks where resistivity changes from $T^4$ to linear-$T$ behavior. The calculation chain uses density-functional theory to obtain bands, phonons, and electron-phonon matrix elements on a coarse grid, Wannier interpolation to reach a fine grid, and an accurate Boltzmann transport equation to compute resistivity and thermal conductivity. Doping is treated by rigidly shifting the Fermi energy of the pristine band structure, with carrier densities read from the density of states, which is nearly constant in this 2DEG. This clean-parabolic-band description is what produces the doping-independent $\\Theta_{BG}$ and the resistivity scalings $\\rho \\propto 1/E_F$ and $1/E_F^2$ in different doping regimes.","core_discovery":"The central claim is that 2D Na is a new thermodynamically stable elemental 2D metal whose transport is governed by a clean 2DEG Fermi surface originating from the half-filled $3s$ orbital. First-principles DFT/DFPT calculations combined with an accurate Boltzmann-transport solution show that the electron-phonon limited resistivity $\\rho_{e-ph}(T)$ obeys the Bloch-Grüneisen picture, with two regimes crossing at $\\Theta_{BG} \\approx 50$ K that is independent of the carrier type or density. At experimentally accessible electron doping ($E_F = 0.5$ eV, carrier density $2.23 \\times 10^{14}$ cm$^{-2}$), $\\rho_{e-ph}$ is about 1.4 times that of graphene at low temperature but becomes smaller than graphene's above 450 K. The electronic thermal conductivity of pure 2D Na is about 1.24 times that of bulk Na at 300 K, and the calculated Lorenz number is $2.41 \\times 10^{-8}\\ \\mathrm{V}^2/\\mathrm{deg}^2$, so the Wiedemann-Franz law is satisfied. The paper concludes that the same transport mechanism should appear in all Na-like systems with a clean Fermi surface, including bulk compounds with planes of Na atoms.","pith_inferences":["The rigid-band approximation is likely the fragile step: self-consistent doping could renormalize the $3s$ band and alter the velocity enhancement that produces the graphene crossover, so a direct test would be valuable.","The pinned 50 K Bloch-Grüneisen temperature suggests a universal low-temperature $T^4$ regime in low-density 2D metals with soft flexural phonons, potentially applying to other alkali or alkaline-earth monolayers.","If the free-standing limit is difficult to realize, bulk compounds with Na planes, such as sodium cobalt oxides, could serve as a proxy to look for the clean-Fermi-surface transport signature.","Because the Lorenz number is constant, the thermopower of 2D Na should be directly tied to the resistivity ratio, hinting at a use in thermoelectric or bolometric devices that gate-tune the Fermi level."],"forward_implications":["2D Na becomes a new elemental 2D metal whose carrier density can be tuned by a gate voltage, enabling electronic-device applications.","At electron doping with $E_F = 0.5$ eV, the intrinsic resistivity lies below graphene's above 450 K, suggesting competitive high-temperature electrical transport.","The doping-independent $\\Theta_{BG} \\approx 50$ K shows that transport is controlled by soft phonon modes rather than Fermi-surface geometry, as long as the Fermi surface remains clean.","The Wiedemann-Franz law holds with $L \\approx 2.41 \\times 10^{-8}$ V$^2$/deg$^2$, so the electronic thermal conductivity tracks the electrical conductivity.","The transport mechanism is expected to generalize to other Na-like systems with clean Fermi surfaces, including bulk materials containing Na planes."],"supporting_citations":[{"why":"Survey that identifies hexagonal 2D Na as mechanically stable, providing the structural starting point.","marker":"[23]"},{"why":"Exact Boltzmann transport equation solver for electron-phonon resistivity and thermal conductivity.","marker":"[33]"},{"why":"Wannier interpolation of electron-phonon matrix elements to fine Brillouin-zone grids.","marker":"[34]"},{"why":"DFT/DFPT implementation used for band structures, phonons, and electron-phonon matrix elements.","marker":"[35]"},{"why":"Introduces the Bloch-Grüneisen temperature and the T^4 and linear-T resistivity regimes.","marker":"[24]"},{"why":"Demonstrates gate-tunable carrier density in graphene, establishing the experimental feasibility of the doping levels used.","marker":"[25]"},{"why":"Borophene study showing a constant Bloch-Grüneisen temperature under doping, the direct comparison case for 2D Na.","marker":"[27]"},{"why":"Supplies the resistivity formula and the 1/E_F and 1/E_F^2 scaling for parabolic-band 2D systems.","marker":"[44]"},{"why":"Reference value of the Lorenz number used for the Wiedemann-Franz comparison.","marker":"[46]"}],"fun_headline_variants":["2D sodium beats graphene resistivity above 450 K","Wiedemann-Franz law holds in atomically thin sodium","Gate-tunable resistivity in stable 2D sodium","2D sodium: new stable metal with tunable carrier density"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The predictions depend on the rigid-band approximation for doping, where shifting the Fermi energy leaves the bands, phonons, and electron-phonon coupling unchanged, and on ignoring all scattering besides electron-phonon coupling, such as impurities, electron-electron interactions, and substrate effects.","fun_headline_variants_meta":{"raw":{"variants":["2D sodium beats graphene resistivity above 450 K","Wiedemann-Franz law holds in atomically thin sodium","Gate-tunable resistivity in stable 2D sodium","2D sodium: new stable metal with tunable carrier density"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000963,"raw_usage":{"total_tokens":4158,"prompt_tokens":1060,"completion_tokens":3098,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":676,"completion_tokens_details":{"reasoning_tokens":3030}},"tokens_in":676,"tokens_out":3098,"duration_ms":26478,"temperature":1.0,"reasoning_tokens":3030,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T14:45:19.366541+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A gate-tuned resistivity measurement on a free-standing monolayer sodium sheet would settle the claims: if the crossover below graphene's resistivity near 450 K does not appear, or if the Bloch-Grüneisen temperature shifts by more than a few kelvin with carrier density, the rigid-band clean-Fermi-surface picture fails. A first-principles calculation that re-optimizes the doped band structure self-consistently, instead of rigidly shifting $E_F$, would test the same assumption without experimental access.","supporting_citations":[{"cited_title":"Nevalaita and P","cited_arxiv_id":null,"evidence_quote":"Survey that identifies hexagonal 2D Na as mechanically stable, providing the structural starting point."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Introduces the Bloch-Grüneisen temperature and the T^4 and linear-T resistivity regimes."},{"cited_title":"Zhang, J","cited_arxiv_id":null,"evidence_quote":"Borophene study showing a constant Bloch-Grüneisen temperature under doping, the direct comparison case for 2D Na."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the resistivity formula and the 1/E_F and 1/E_F^2 scaling for parabolic-band 2D systems."},{"cited_title":"Kittel, Wiley, New York , 7th edition, p","cited_arxiv_id":null,"evidence_quote":"Reference value of the Lorenz number used for the Wiedemann-Franz comparison."}],"review_version":1}