{"id":"7ca9f693-ffee-4899-a062-b0c95a59b589","arxiv_id":"1908.09412","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"A Poisson-Schrodinger band-bending model estimates the mobile carrier density in thin topological insulator films from surface Fermi level and gate-dependent transport, without knowing the film's dopant concentration.","lead":"This paper presents a way to estimate how much electrical current leaks through the interior of a thin topological insulator film, using surface measurements such as ARPES and gate-dependent four-point resistance. If reliable, it would let researchers separate the useful surface-state current from parasitic bulk and interface currents in future TI devices.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central dopant-independence claim rests on d << L, but for the BiSbTe3 example L is only ~2d at 4e17 cm^-3 doping, and E_bottom uncertainty changes n_film by an order of magnitude; the thin-film regime is asserted rather than demonstrated.","rationale":"The reader's CONDITIONAL verdict is appropriate. I independently identify the same load-bearing assumption: the d << L screening regime. The paper itself limits the claim in Sections II.B.4 and II.B.5, but does not verify L for the actual film; the numbers are suggestive (L ~ 20 nm at 4 x 10^17 cm^-3 vs d = 10 nm), not conclusive. The quantization treatment is a strength: the factor-2-2.5 reduction is plausible, and the subband energies make the 6 x 10^11 cm^-2 result consistent with the inputs. The main weakness is that the central quantitative prediction is not compared with any independent measurement, and the asymmetric input E_bottom^F comes from ref. [21] without derivation; Fig. 8 shows an order-of-magnitude spread. None of this proves the mechanism wrong, but it does mean the claim of a 'well-founded estimate' is conditional on parameters not established in this paper. Hence the verdict remains CONDITIONAL, with no change from the reader's assessment.","tokens_in":35353,"tokens_out":14459,"duration_ms":157533,"concrete_test":"Compute L from Eq. (3) with independently measured values: measure epsilon_r of the actual MBE BiSbTe3 film (THz or capacitance) and estimate n_b + p_b from transport on a separately grown thick film. Then re-evaluate n_film for d = 10 nm with Eq. (7) and Appendix B over dopant densities 10^16, 10^17, 4 x 10^17, and 10^18 cm^-3. If n_film varies by more than ~30% across this range, the dopant-independence claim is not valid for this material. Also, as a second check, report n_film from Fig. 8 as a function of E_bottom^F at E_top^F = 240 meV; if the spread exceeds factor 3, the asymmetric estimate needs an independent E_bottom^F measurement.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central claim is that in a thin TI film n_film is nearly independent of the unknown dopant density, so E_top^F and E_bottom^F from surface probes suffice. The mechanism is the screening-length condition: Eq. (3) defines L, and Eq. (2) shows band-bending curvature is proportional to 1/L^2; Sec. II.B.4 states that independence requires d << L. For the paper's own BiSbTe3 parameters (Table II, epsilon_r ~ 100, E_g = 0.26 eV, m* = 0.15 m_e), a non-degenerate dopant density of 4 x 10^17 cm^-3 gives L ~ 20 nm, only twice the d = 10 nm film; at 10x higher doping or a factor 2 lower epsilon_r, L <= d and Fig. 4's red curves would no longer be dopant-independent. No measurement of epsilon_r or of n_b + p_b for the actual MBE film is provided, so the key regime is asserted, not established. Separately, even granting d << L, Fig. 8 shows n_film changes by up to an order of magnitude as E_bottom^F varies over the plausible range; the symmetric and asymmetric results differ by a factor 3. Since E_bottom^F is imported from the quantum-capacitance analysis of ref. [21] without independent validation, the 'well-founded estimate' is not yet uniquely determined. Weighing these, the weakest load-bearing step is the unverified d << L regime.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript proposes a semi-classical calculation of the band bending and mobile carrier density in the interior of thin topological-insulator films, using experimentally determined surface and interface Fermi levels as input. The authors solve Poisson's equation with Boltzmann statistics and a Debye screening length, treat symmetric (E_top^F = E_bottom^F) and asymmetric boundary conditions, add quantization via square or triangular infinite wells, and apply the model to a 10 nm BiSbTe3 film using parameters from earlier ARPES and gate-dependent four-probe measurements. The central claim is that in the thin-film limit d << L the band bending and integrated film carrier density n_film are nearly independent of the unknown dopant concentration (represented by E_bulk^F), so that n_film, and with a known mobility the film conductivity, can be estimated from surface-sensitive measurements alone. For the BiSbTe3 example, the symmetric approximation gives n_film ~ 6e11 cm^-2 and the asymmetric approximation ~2e11 cm^-2 at zero gate voltage.","tokens_in":35756,"tokens_out":5996,"duration_ms":60777,"significance":"If the central claim holds, the paper offers a practically useful way to quantify the parasitic film-interior conduction channel that otherwise contaminates transport studies of topological surface states. The analytic Poisson-band-bending framework, the explicit distinction between dopant concentration and mobile carrier density, and the use of ARPES and gate-dependent transport data as inputs are valuable strengths, and the authors are unusually explicit about their approximations. However, the quantitative usefulness of the results is conditional on assumptions that are acknowledged but not fully validated: the d << L screening regime for the specific film, the Boltzmann approximation at a Fermi level only 20 meV below the conduction band, and the reliability of E_bottom^F and the mobility imported from the earlier transport analysis. Because these assumptions directly affect the quoted numbers, the paper is more convincing as a methodological proposal than as a validated quantitative estimate.","major_comments":[{"comment":"The central dopant-independence claim requires d << L, but for the paper's own BiSbTe3 parameters the condition is only marginal: with epsilon_r ≈ 100 and n_b + p_b = 4e17 cm^-3, Eq. (3) gives L ≈ 19 nm, only about twice the d = 10 nm film; a factor-two reduction in epsilon_r or a tenfold increase in doping makes L ≤ d. The paper provides no measured epsilon_r or n_b + p_b for the MBE film under consideration, so the regime in which Figs. 4 and 5 are claimed to be dopant-independent is asserted rather than demonstrated. Please quantify d/L over the plausible doping range and either restrict the claim to that range or provide direct evidence that the 10 nm BiSbTe3 film is in the d << L regime.","section":"II.B.4, Eq. (3)"},{"comment":"The asymmetric result depends strongly on E_bottom^F, which is imported from the quantum-capacitance transport analysis of ref. [21]. For fixed E_top^F = 240 meV, n_film varies by up to an order of magnitude with E_bottom^F across the plotted range, and the symmetric and asymmetric values (6e11 vs 2e11 cm^-2) differ by a factor of three. If the ref. [21] analysis already determines the film channel's carrier density or conductivity, the asymmetric calculation is not an independent estimate; if it does not, the uncertainty in E_bottom^F should be propagated into n_film. Please state explicitly what ref. [21] provides and give an uncertainty budget for the quoted n_film values.","section":"II.C, Fig. 8"},{"comment":"The Boltzmann approximation is knowingly used outside its validity for the example: E_top^F = 240 meV lies only 20 meV below the conduction-band edge, and the text concedes deviations of up to 50% from the Fermi-Dirac distribution. Since all quantitative values, including the factor 2 to 2.5 reduction from quantization and the final n_film estimates, are derived within this scheme, the paper should show at least one representative Fermi-Dirac calculation to demonstrate that the dopant-independence conclusion and the quoted n_film values survive.","section":"II.B.5"},{"comment":"The Schrödinger-Poisson self-consistency is truncated after one and a half iterations and wavefunction weighting is neglected. The authors state that these approximations may overestimate band bending near the surfaces and that quantization reduces n_film by a factor of 2 to 2.5 relative to the purely classical result. Because that factor materially changes the final estimates, a benchmark against a fully self-consistent solution for at least one parameter set is needed to establish the numerical accuracy of the method.","section":"II.B.2 and Supplemental Material B"}],"minor_comments":[{"comment":"The value v_Fermi = 5.6 × 10^-5 ms^-1 appears to be a units typo; the Dirac velocity in a topological insulator should be several orders of magnitude larger, likely 5.6 × 10^5 m/s.","section":"Table II"},{"comment":"The text says E_bottom^F = 156 meV at zero gate voltage while the Fig. 6 caption gives 155 meV; please make these values consistent.","section":"II.C and Fig. 6"},{"comment":"The typeset equation contains garbled square-root and brace symbols and should be checked carefully in production.","section":"Eq. (7)"},{"comment":"The distinction between dopant concentration and mobile carrier density is helpful and should be kept; consider stating the definition of the 'bulk Fermi energy of a corresponding extended crystal' more prominently at first use.","section":"II.B.7"}],"recommendation":"major_revision","confidential_remarks":"The main numerical inputs for the asymmetric branch come from the authors' own previous work (ref. [21]); the editor may wish to confirm the degree of overlap and whether the current paper adds a sufficiently distinct methodological contribution. The core Poisson-band-bending formalism is standard, but the surface CNL concept is clearly marked as non-measurable and ultimately not used as an input, so the presentation is honest. The stress-test concern about the d << L regime lands for this manuscript: the example film sits near the boundary of the claimed regime, and without a quantified screening length the central dopant-independence statement remains only partially supported."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe paper is worth a look if you work on transport in TI thin films. It offers a practical recipe: from an ARPES-measured top Fermi level and, in the asymmetric version, a bottom Fermi level extracted from gate-dependent four-probe data, you solve Poisson–Boltzmann with quantization corrections and get the mobile carrier density in the film interior. That fills a real need, because the film interior is otherwise hard to access experimentally. The extension beyond the Schottky approximation (refs 29,30) to weak bending, accumulation, and asymmetric boundary conditions is a genuine step forward. The interface conductivity measurements in Appendix A are also useful and concrete.\n\nThe paper is honest about its approximations: Boltzmann vs Fermi–Dirac can deviate by up to 50% (they note the bending is overestimated), the self-consistency loop is truncated after one and a half iterations, and the wave functions are neglected in the quantum corrections. That transparency is good.\n\nThe soft spots are real, though. The central claim — that the film carrier density is nearly independent of dopant concentration — rests on d << L. For their own BiSbTe3 parameters, L is about 20 nm, only twice the 10 nm film. That is not deep in the claimed regime, and the paper gives no measurement of epsilon_r or doping for the actual film. It says \"not too high dopant levels\" but doesn't say where the boundary is. Second, the asymmetric result depends on E_bottom^F and the mobility, both imported from the authors' earlier quantum-capacitance analysis (ref 21) without independent validation. Fig 8 shows n_film can change by an order of magnitude across the plausible range of E_bottom^F, and the symmetric vs asymmetric numbers differ by a factor of three. That is a wide spread for something called \"well-founded.\" Also, Table II lists v_Fermi = 5.6e-5 m/s, which looks like a units typo (surely 5.6e5 m/s?).\n\nThe math itself is standard semiconductor physics and seems internally consistent. The paper would be stronger with error bars on n_film, a quantitative statement of the d/L regime, and a comparison against an independent measurement or another group's data.\n\nVerdict: this deserves a serious referee. It is a useful modeling tool for a real experimental problem, with real data behind it, and the approximations are disclosed. I'd send it out, but I'd expect the reviewers to push for a sharper statement of the validity regime and for uncertainty propagation.","headline":"A practically useful band-bending recipe for TI thin-film transport, but the central dopant-independence claim is only marginally in its stated regime and the outputs are sensitive to inputs taken from the authors' own earlier analysis.","tokens_in":36260,"tokens_out":4086,"would_cite":false,"duration_ms":39194,"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":"The paper shows that the interior charge carrier density of a thin topological insulator film can be estimated from surface measurements alone, without knowing its dopant concentration.","keywords":["topological insulator","thin film","band bending","surface states","parasitic conduction channels","charge carrier density","ARPES","four-point transport"],"falsifier":"Grow thin TI films of the same material with intentionally varied dopant levels, for example by changing growth stoichiometry, while keeping the surface Fermi level fixed, and measure the interior sheet carrier density at zero gate voltage; if that density shifts by more than the calculation's predicted variation across the band gap, the screening-length assumption is violated and the dopant-independent estimate fails.","tokens_in":35166,"feed_emoji":"⚡","tokens_out":8642,"duration_ms":74571,"temperature":0.7,"pith_summary":"Thin films of topological insulators conduct current through several parallel channels: the protected surface states, the film interior, the interface layer, and the substrate. The interior contribution is hard to measure directly, so the authors calculate it from surface-sensitive data. Their central result is that, in the thin-film limit, the near-surface band bending is largely independent of the unknown dopant concentration inside the film, because the screening length is much longer than the film thickness. This makes it possible to estimate the total mobile charge carrier density, and hence the conductivity of the film interior, from ARPES measurements of the surface Fermi level, optionally combined with gate-dependent four-point transport. For BiSbTe$_3$, the symmetric calculation gives roughly $6\\times 10^{11}$ cm$^{-2}$ and the asymmetric, gate-informed calculation gives roughly $2\\times 10^{11}$ cm$^{-2}$ at zero gate voltage.","feed_headline":"Thin-film limit removes dopant guesswork from TI transport estimates","feed_subtitle":"Long screening length flattens band bending, so surface data alone yield the film's interior carrier density.","key_machinery":"The central mechanism is the comparison between the Debye screening length $L=\\sqrt{\\epsilon_0\\epsilon_r k_B T/(q^2(n_b+p_b))}$ and the film thickness $d$. When $d\\ll L$, Poisson's equation produces only weak band bending, so the boundary conditions at the top and bottom surfaces, represented by the surface Fermi levels $E_F^{\\rm top}$ and $E_F^{\\rm bottom}$, control the carrier distribution in the film rather than the unknown dopant concentration. Around this screening argument, the paper builds a semi-classical Poisson-Schrödinger calculation: it solves Poisson's equation with charge neutrality between the topological surface states and the space-charge layer, renormalizes the effective densities of states using quantized subbands in a square or triangular well, and iterates once to a second Poisson solution. A separate but supporting element is the charge-neutrality-level picture of the surface, which explains how trivial defect states shift the effective filling level of the Dirac cone and hence the measured surface Fermi energy.","core_discovery":"The paper establishes that in a thin topological insulator film, the top and bottom topological surface states exchange charge with the film interior until charge neutrality is reached, and because the Debye screening length $L$ far exceeds the film thickness $d$, the bands bend only weakly. As a result, a wide range of unintentional dopant concentrations leaves the band positions essentially fixed by the surface Fermi levels rather than by the bulk dopant density. The authors solve Poisson's equation for the band bending using symmetric boundary conditions (top and bottom surface Fermi levels equal) or asymmetric boundary conditions (top level from ARPES, bottom level extracted from gate-dependent four-point measurements), and they include confinement quantization through a square-well or triangular-well approximation of the band-bending potential. The integrated mobile carrier density $n_{\\rm film}$ is then nearly flat as a function of the bulk Fermi energy across the band gap. For a 10 nm BiSbTe$_3$ film they obtain $n_{\\rm film}\\approx 6\\times 10^{11}$ cm$^{-2}$ in the symmetric approximation and $\\approx 2\\times 10^{11}$ cm$^{-2}$ in the asymmetric approximation at zero gate voltage, with the latter being the more precise estimate. Combined with an interface conductivity measured on the bare substrate reconstruction before film growth, this allows the total measured 2D conductivity to be decomposed into surface-state, interior, interface, and substrate contributions.","pith_inferences":["A direct test of the claim would be to grow films with deliberately varied dopant concentrations and check that the interior sheet density stays on the flat plateau of the calculation's Fig. 5, rather than tracking the dopant level.","The screening-length criterion $L\\gg d$ offers a quick screening rule for other TI candidates: compute $L$ from the dielectric constant and intrinsic carrier density before investing in the full ARPES-plus-gate measurement program.","Low-temperature transport, where the film interior mobility rises relative to the surface channels, would provide a sharper falsifier, because the calculated interior conductivity could be compared with measured four-point data at temperatures where the mobility difference separates the channels.","The same logic could be extended to estimate the interior carrier density of other van-der-Waals thin films with large dielectric constants, not only topological insulators, as long as the weak-bending condition is met."],"forward_implications":["Given an ARPES measurement of the top surface Fermi level, the interior mobile carrier density of a thin TI film can be estimated without knowing its unintentional dopant concentration.","With additional gate-dependent four-point transport data, the asymmetric calculation gives a more precise value; for BiSbTe$_3$ at zero gate voltage this is about $2\\times 10^{11}$ cm$^{-2}$, one third of the symmetric estimate.","Because the interface conductivity can be measured on the bare substrate reconstruction before film growth, the total measured 2D conductivity can be decomposed into TSS, interior, interface, and substrate channels.","The approach is presented as general across TI material classes, and the same gate-dependent measurements also yield the carrier mobility, so the interior conductivity itself becomes accessible.","For device design, the results imply that low-conductivity substrate terminations such as Te/Si(111)-(1$\\times$1) keep parasitic interface currents negligible, whereas Bi/Si(111)-(√3$\\times$√3) would carry a substantial fraction of the current."],"supporting_citations":[{"why":"Supplies the measured top surface Fermi energy, gate-dependent four-point transport data, and material parameters for BiSbTe$_3$ that serve as input to the band bending calculations.","marker":"[21]"},{"why":"Reports the large dielectric constant ($\\epsilon_r\\approx 100$) for BiSbTe$_3$ that makes the Debye screening length exceed the film thickness.","marker":"[13]"},{"why":"Provides additional evidence for the large dielectric constant governing screening in TI thin films.","marker":"[14]"},{"why":"Adds further support for the large dielectric constant used in the screening-length argument.","marker":"[15]"},{"why":"Reports the Te/Si(111)-(1$\\times$1) interface conductivity measurement used to quantify the interface parasitic channel.","marker":"[12]"},{"why":"Provides the four-probe method and the Bi/Si(111)-(√3$\\times$√3) surface conductivity used to estimate interface channel contributions.","marker":"[6]"},{"why":"Supplies the semiconductor Poisson-equation band bending formalism on which the calculation is built.","marker":"[17]"},{"why":"Presents the Schottky approximation for band bending in thin TI films that this paper's Poisson-Schrödinger approach extends beyond.","marker":"[29]"}],"fun_headline_variants":["Bulk carrier density in thin TIs from surface data only","Thin-film limit negates need for dopant concentration","Surface measurements reveal interior conduction in thin films","New method separates surface and bulk currents in topological films"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculation stands on the assumption that the film is much thinner than the screening length of its mobile charges, so band bending stays weak; if a material screens more efficiently, through a smaller dielectric constant, higher carrier density, or degenerate doping, the unknown dopant concentration would control the result.","fun_headline_variants_meta":{"raw":{"variants":["Bulk carrier density in thin TIs from surface data only","Thin-film limit negates need for dopant concentration","Surface measurements reveal interior conduction in thin films","New method separates surface and bulk currents in topological films"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000298,"raw_usage":{"total_tokens":1807,"prompt_tokens":1109,"completion_tokens":698,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":725,"completion_tokens_details":{"reasoning_tokens":634}},"tokens_in":725,"tokens_out":698,"duration_ms":7603,"temperature":1.0,"reasoning_tokens":634,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:13:32.235429+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Grow thin TI films of the same material with intentionally varied dopant levels, for example by changing growth stoichiometry, while keeping the surface Fermi level fixed, and measure the interior sheet carrier density at zero gate voltage; if that density shifts by more than the calculation's predicted variation across the band gap, the screening-length assumption is violated and the dopant-independent estimate fails.","supporting_citations":[{"cited_title":"Zhang, C.-X","cited_arxiv_id":null,"evidence_quote":"Supplies the measured top surface Fermi energy, gate-dependent four-point transport data, and material parameters for BiSbTe$_3$ that serve as input to the band bending calculations."},{"cited_title":"It strongly depends both on the mate- rial system and the preparation parameters","cited_arxiv_id":null,"evidence_quote":"Reports the large dielectric constant ($\\epsilon_r\\approx 100$) for BiSbTe$_3$ that makes the Debye screening length exceed the film thickness."},{"cited_title":"10, the band bending corre- sponding to case A is shown","cited_arxiv_id":null,"evidence_quote":"Provides additional evidence for the large dielectric constant governing screening in TI thin films."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Adds further support for the large dielectric constant used in the screening-length argument."},{"cited_title":"However, exact values for the surface conductivities of the respective (1 ×1) and (7×7) reconstructions have not yet been reported in the literature","cited_arxiv_id":null,"evidence_quote":"Reports the Te/Si(111)-(1$\\times$1) interface conductivity measurement used to quantify the interface parasitic channel."},{"cited_title":"BiSbTe 3, which was studied extensively in ref","cited_arxiv_id":null,"evidence_quote":"Provides the four-probe method and the Bi/Si(111)-(√3$\\times$√3) surface conductivity used to estimate interface channel contributions."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the semiconductor Poisson-equation band bending formalism on which the calculation is built."},{"cited_title":"Durand, X.-G","cited_arxiv_id":null,"evidence_quote":"Presents the Schottky approximation for band bending in thin TI films that this paper's Poisson-Schrödinger approach extends beyond."}],"review_version":1}