Pith. sign in

REVIEW 3 major objections 6 minor 26 references

STM study of single phosphorus incorporation into silicon by heating PBr3 on Si(100)

T0 review · 3 major / 6 minor · reviewed 2026-07-31 · grok-4.5

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read A single phosphorus atom from dissociated PBr3 on Si(100) exchanges with nearby silicon on heating, forming a stable P–Si–Br heterodimer with Br atop the Si, with incorporation starting as low as 175 °C.

desk verdict Clean same-area STM tracking of single-P incorporation from PBr3, with a coherent P–Si–Br assignment and a 1.44 eV barrier that matches the low onset temperature. read the letter →

arxiv 2607.27939 v1 pith:RWCKY6PC submitted 2026-07-30 cond-mat.mtrl-sci cond-mat.mes-hall

Tatiana V. Pavlova , Vladimir M. Shevlyuga This is my paper · ORCID
classification cond-mat.mtrl-scicond-mat.mes-hall
keywords STMphosphorusdopingSi(100)PBr3atomic-precisionP-SiheterodimerDFTincorporationbarrier
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper establishes the atomic pathway by which one phosphorus atom enters the Si(100) surface when delivered as PBr3. The molecule fully dissociates at room temperature, leaving phosphorus in a preferred end-bridge site. By annealing the identical surface region inside the STM, the authors track that same atom before and after heating and show it swaps with a neighboring silicon atom, ending as a P–Si heterodimer with a bromine atom sitting on the silicon. DFT gives a 1.44 eV barrier for the lowest path, matching the observed onset near 175 °C. The result supplies a concrete, site-resolved picture of single-atom doping that can guide more precise placement of donors in silicon.

What carries the argument

The P–Si–Br heterodimer—the substituted surface dimer with Br bonded on top of the remaining Si atom—together with same-area in-situ STM annealing and NEB barriers that map the exchange from the initial end-bridge Seb geometry to this product.

What would settle it

A filled-state bias series or local spectroscopy on the same tracked sites that disagreed with the calculated electronic structure of P–Si–Br, or repeated appearance of that feature at locations that never held an end-bridge phosphorus before annealing, would refute the assignment.

Watch

Extended reading notes

Core claim

Upon annealing, phosphorus that starts in the end-bridge site after room-temperature dissociation of PBr3 on Si(100) exchanges with a nearby silicon atom and forms a stable P–Si–Br complex in which bromine sits atop the silicon atom of the heterodimer. Same-area in-situ STM before and after heating, combined with DFT, identifies this structure and yields a minimum activation barrier of 1.44 eV, consistent with incorporation beginning at temperatures as low as 175 °C.

Load-bearing premise

The post-anneal depression-plus-protrusion feature is identified as the brominated P–Si heterodimer mainly by visual match to empty-state STM simulations and relative DFT energies of a limited set of candidate geometries.

Editorial extensions

If this is right

  • Atomic-precision doping with PBr3 can be designed around the end-bridge starting geometry and the 1.44 eV exchange path.
  • Silicon-island formation is a secondary consequence of ejected Si, not the primary signature needed to locate incorporated P.
  • Bromine pairs diffuse with a lower barrier than P–Si exchange, so Br can leave the incorporation site before or during doping.
  • Phosphorus dimers form readily but require extra energy to break before incorporation, raising the temperature needed if pairing occurs.
  • In-situ same-area STM separates true incorporation sites from water-related defects that look similar after annealing.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • Residual Br that preferentially caps the heterodimer may need to be accounted for in halogen-based lithography masks before clean substitutional P is obtained.
  • The lower thermal onset relative to typical PH3 routes could allow gentler budgets when building single-donor devices.
  • Tracking one precursor molecule through the anneal could become a routine check for other dopant precursors on Si(100).
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 6 minor

Summary. The manuscript reports a combined STM/DFT study of single-P incorporation into Si(100) from room-temperature-dissociated PBr3. Using in situ annealing inside the STM, the authors track the same surface regions before and after heating and assign the post-anneal “object B” (depression adjacent to a protrusion on one dimer) to a P–Si heterodimer with Br atop the Si atom. Sequential anneals show Si-island formation beginning at 175 °C. NEB calculations give a minimum barrier of 1.44 eV from the end-bridge Seb geometry to this P–Si–Br structure, stated to be consistent with the observed onset temperature. Additional discussion covers P2 dimer formation, Br-pair diffusion, and water-related features that can mimic brominated heterodimers.

Significance. Atomic-precision donor placement in Si is a central goal for single-impurity devices; a site-resolved incorporation pathway from a halogenated precursor relevant to halogen resist lithography is therefore of clear technological and surface-science interest. The decisive experimental strength is same-area STM before and after in situ anneal (Fig. 2), which directly links pre-anneal Seb loci to post-anneal features and avoids reliance on island morphology alone. DFT STM simulations and NEB barriers are used in a conventional supporting role. If the P–Si–Br assignment holds, the work supplies a concrete atomic pathway and a lower onset temperature than typical PH3 reports, with practical implications for process windows in atomic-precision doping.

major comments (3)
  1. [§III.B, Fig. 3] §III.B and Fig. 3: The central structural assignment of object B as P–Si–Br rests on (i) visual agreement of one empty-state Tersoff–Hamann image with experiment and (ii) relative DFT energies among a small candidate set (Br-on-P higher by 1.26 eV; distant Br higher by 0.35 eV). Bare P–Si heterodimer depressions are already known from PH3 work; the bright lobe at Us > +2 V is ascribed to Br, but other adsorbates or tip states can produce similar contrast. No filled-state images, bias series, or STS on the same tracked sites are shown, and the simulated library is not exhaustive. Because the 1.44 eV pathway and the claim of a stable P–Si–Br product both depend on this ID, the manuscript should either add orthogonal contrast (filled states / bias dependence on tracked sites) or substantially expand the candidate set and state the assignment more cautiously.
  2. [§III.A, Fig. 2] §III.A and Fig. 2: Same-area images show object B at former Seb sites “in most cases,” but no quantitative yield is reported (fraction of Seb → object B vs bare heterodimer, islands, or other products; number of independent tracked molecules). Without counts or error bars, the claim that the dominant pathway is Seb → P–Si–Br cannot be weighed against minority channels. A simple tally over the imaged areas would make the pathway claim falsifiable and proportionate to the data.
  3. [§II, §III.A] §II and §III.A: The 175 °C onset is load-bearing for consistency with the 1.44 eV barrier, yet temperature is inferred from a heating-current calibration on a separate test wafer (thermocouple 150–400 °C; pyrometer check only at 400 °C). Contact thermal gradients and emissivity differences can shift the low-T end by tens of degrees. The paper should quantify calibration uncertainty at ~175–220 °C and discuss how a ±25–50 °C error would affect the barrier–temperature comparison (and the contrast with ~250 °C PH3 reports).
minor comments (6)
  1. [Title, passim] Title and running text contain PDF-extraction spacing artifacts (“b y heating”, “incorpo ration”, “functiona l”, etc.). Clean for production.
  2. [Fig. 2] Fig. 2 panel label “Afret heating” should be “After heating.”
  3. [Fig. 1] Fig. 1 caption and panel letters: panel (c) is written as “( ) с” (Cyrillic) in the source; unify Latin labels.
  4. [§II, Fig. 4] NEB details (five images, 0.03 eV/Å) are given; a brief check that the 1.44 eV barrier is stable to image count or a climbing-image result would strengthen §III.B without new experiments.
  5. [§III.A] Clarify early that all displayed STM frames are empty-state only, so “protrusion/depression” language is polarity-specific (mentioned once but easy to miss when comparing to PH3 heterodimer literature).
  6. [§III.D] §III.D water-related diagonal features: the tentative Br/OH assignment is useful for avoiding misidentification; a single sentence on how often such sites appear relative to true object B would help readers.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: incorporation pathway is fixed by same-area STM and independent DFT, not by construction from inputs or self-citation.

full rationale

The load-bearing chain is experimental and non-tautological. PBr3 dissociates at RT into identifiable P and Br sites (Seb/Sbr); in-situ annealing inside the STM tracks the same surface region before and after heating (Fig. 2), showing new depression+protrusion features (object B) at former Seb loci plus Si islands. Object B is assigned as the P–Si–Br heterodimer by matching empty-state STM to Tersoff–Hamann simulations of candidate geometries and by relative DFT energies (Fig. 3, §III.B)—standard structure identification, not a fit of a free parameter that is then re-predicted. The NEB barrier from Seb to Sin1 (1.44 eV) is computed independently and only afterward compared to the observed onset near 175 °C; the temperature is not used to set the barrier. Prior author papers supply the RT adsorption structures and Br visibility, which are separate measurements of the initial state, not algebraic identities that force the incorporation result. No uniqueness theorem, ansatz smuggled via self-citation, fitted-input-as-prediction, or renaming of a known empirical law appears. Correctness risk in the limited STM–DFT fingerprint library for object B is real but is not circularity.

Assumptions & free parameters 3 free parameters · 5 assumptions · 0 invented entities

The claim rests on standard UHV-STM practice, PBE-DFT total energies and NEB barriers, Tersoff–Hamann STM simulation, and a lab-specific current-to-temperature calibration. No new physical entities are postulated. Free parameters are ordinary computational/experimental choices (cutoff, slab, NEB image count, anneal times) rather than values fitted to force the incorporation story.

free parameters (3)
  • NEB intermediate image count (five) and force criterion 0.03 eV/Å = 5 images; 0.03 eV/Å
    Controls resolution of the reported 1.44 eV and 1.17 eV barriers; not converged systematically in the text.
  • Sample temperature vs heating-current calibration on a test wafer = calibration curve 150–400 °C
    Maps anneal currents to 175–400 °C; thermocouple plus pyrometer check at 400 °C. Absolute T for the in-STM 300 °C and 175 °C onset inherit this calibration uncertainty.
  • DFT setup (PBE, 350 eV cutoff, 4×4/8×8 slabs, 3×3×1 k-mesh) = PBE, 350 eV, 8-layer slab
    Standard but choice-dependent absolute energies and barriers; relative ordering of Sin structures could shift with functional/dispersion.
assumptions (5)
  • domain assumption PBE-DFT total energies and NEB barriers adequately rank Si(100)–P–Br configurations and approximate experimental activation energies.
    Used throughout §II and §III.B–C to identify P–Si–Br as most stable and to claim consistency with 175 °C onset.
  • domain assumption Tersoff–Hamann empty-state simulations, averaged over two Si-dimer buckling configurations, are sufficient to fingerprint embedded P and Br-terminated heterodimers in experimental STM.
    Central to assigning objects A and B in Fig. 3 and §III.B.
  • domain assumption Appearance of Si islands after anneal indicates ejection of Si by substitutional P incorporation.
    Inherited from PH3 literature and used as supporting indicator in §III.A; authors correctly note they do not rely on it alone.
  • domain assumption In-STM direct-current anneal with tip retracted leaves the imaged region representative and does not create spurious P/Br rearrangements from tip or field effects.
    Required for the before/after same-area claim in Fig. 2 and §II.
  • domain assumption Standard UHV surface preparation and dissociative adsorption assignments for PBr3 on Si(100) from the authors’ prior work hold.
    Initial Seb/Sbr structures (Fig. 1f) are taken from Refs. 13–14 as the starting point for incorporation.

how reviews work

0 comments
Cite this review

Pith. "Pith review of STM study of single phosphorus incorporation into silicon by heating PBr3 on Si(100)." pith.science (2026). https://pith.science/paper/RWCKY6PC

@misc{pith2026260727939,
  author       = {Pith},
  title        = {Pith review of: STM study of single phosphorus incorporation into silicon by heating PBr3 on Si(100)},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/RWCKY6PC}},
  note         = {Machine review of arXiv:2607.27939}
}
read the original abstract

The objective of miniaturizing doped areas in silicon, with the ultimate goal of achieving atomic-precision doping, requires a fundamental understanding of the dopant incorporation process at the atomic level. We present a combined scanning tunneling microscopy (STM) and density functional theory (DFT) investigation of single phosphorus atom incorporation into the Si(100) surface. Phosphorus was supplied via PBr3 molecules, which completely dissociate on Si(100) at room temperature. By performing in situ annealing within the STM, we directly tracked the same phosphorus atom before and after heating. Upon annealing, the P atom undergoes an exchange with a nearby Si atom, forming a stable P-Si-Br complex with a Br atom located atop the Si atom of the heterodimer. The activation barrier calculated using DFT is consistent with our observation of doping starting at temperatures as low as 175 C. These results provide detailed atomic-scale insight into the phosphorus incorporation pathway and offer a foundation for improving methods of precise, single-atom doping in silicon.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

26 extracted references · 19 canonical work pages

  1. [1]

    Warschkow , author H

    author author O. Warschkow , author H. F. \ Wilson , author N. A. \ Marks , author S. R. \ Schofield , author N. J. \ Curson , author P. V. \ Smith , author M. W. \ Radny , author D. R. \ McKenzie , \ and\ author M. Y. \ Simmons ,\ 10.1103/PhysRevB.72.125328 journal journal Phys. Rev. B \ volume 72 ,\ pages 125328 ( year 2005 ) NoStop

    work page doi:10.1103/physrevb.72.125328 2005

  2. [2]

    author author H. F. \ Wilson , author O. Warschkow , author N. A. \ Marks , author N. J. \ Curson , author S. R. \ Schofield , author T. C. G. \ Reusch , author M. W. \ Radny , author P. V. \ Smith , author D. R. \ McKenzie , \ and\ author M. Y. \ Simmons ,\ 10.1103/PhysRevB.74.195310 journal journal Phys. Rev. B \ volume 74 ,\ pages 195310 ( year 2006 ) NoStop

    work page doi:10.1103/physrevb.74.195310 2006

  3. [3]

    Warschkow , author N

    author author O. Warschkow , author N. J. \ Curson , author S. R. \ Schofield , author N. A. \ Marks , author H. F. \ Wilson , author M. W. \ Radny , author P. V. \ Smith , author T. C. G. \ Reusch , author D. R. \ McKenzie , \ and\ author M. Y. \ Simmons ,\ 10.1063/1.4939124 journal journal J. Chem. Phys. \ volume 144 ,\ pages 014705 ( year 2016 ) NoStop

    work page doi:10.1063/1.4939124 2016

  4. [4]

    author author S. R. \ Schofield , author N. J. \ Curson , author M. Y. \ Simmons , author F. J. \ Rue , author T. Hallam , author L. Oberbeck , \ and\ author R. G. \ Clark ,\ 10.1103/PhysRevLett.91.136104 journal journal Phys. Rev. Lett. \ volume 91 ,\ pages 136104 ( year 2003 ) NoStop

    work page doi:10.1103/physrevlett.91.136104 2003

  5. [5]

    author author S. R. \ McKibbin , author W. R. \ Clarke , author A. Fuhrer , author T. C. G. \ Reusch , \ and\ author M. Y. \ Simmons ,\ 10.1063/1.3269924 journal journal Appl. Phys. Lett. \ volume 95 ,\ pages 233111 ( year 2009 ) NoStop

    work page doi:10.1063/1.3269924 2009

  6. [6]

    author author J. M. \ Bennett , author O. Warschkow , author N. A. \ Marks , \ and\ author D. R. \ McKenzie ,\ 10.1103/PhysRevB.79.165311 journal journal Phys. Rev. B \ volume 79 ,\ pages 165311 ( year 2009 ) NoStop

    work page doi:10.1103/physrevb.79.165311 2009

  7. [7]

    author author N. J. \ Curson , author S. R. \ Schofield , author M. Y. \ Simmons , author L. Oberbeck , author J. L. \ O'Brien , \ and\ author R. G. \ Clark ,\ 10.1103/PhysRevB.69.195303 journal journal Phys. Rev. B \ volume 69 ,\ pages 195303 ( year 2004 ) NoStop

    work page doi:10.1103/physrevb.69.195303 2004

  8. [8]

    Reusch , author N

    author author T. Reusch , author N. Curson , author S. Schofield , author T. Hallam , \ and\ author M. Simmons ,\ https://doi.org/10.1016/j.susc.2005.10.031 journal journal Surf. Sci. \ volume 600 ,\ pages 318 ( year 2006 ) NoStop

    work page doi:10.1016/j.susc.2005.10.031 2005

Show all 26 references
  1. [9]

    author author T. C. G. \ Reusch , author M. W. \ Radny , author P. V. \ Smith , author O. Warschkow , author N. A. \ Marks , author N. J. \ Curson , author D. R. \ McKenzie , \ and\ author M. Y. \ Simmons ,\ 10.1021/jp068834c journal journal J. Phys. Chem. C \ volume 111 ,\ pa...

    2007 doi

  2. [10]

    \ Hsueh , author D

    author author Y.-L. \ Hsueh , author D. Keith , author Y. Chung , author S. K. \ Gorman , author L. Kranz , author S. Monir , author Z. Kembrey , author J. G. \ Keizer , author R. Rahman , \ and\ author M. Y. \ Simmons ,\ https://doi.org/10.1002/adma.202312736 journal journal ...

    2024 doi

  3. [11]

    author author T. V. \ Pavlova \ and\ author V. M. \ Shevlyuga ,\ 10.1063/5.0185671 journal journal J. Chem. Phys. \ volume 160 ,\ pages 054701 ( year 2024 ) NoStop

    2024 doi

  4. [12]

    author author T. V. \ Pavlova ,\ 10.1063/5.0266110 journal journal J. Chem. Phys. \ volume 162 ,\ pages 194701 ( year 2025 ) NoStop

    2025 doi

  5. [13]

    author author V. M. \ Shevlyuga , author Y. A. \ Vorontsova , \ and\ author T. V. \ Pavlova ,\ 10.1021/acs.jpcc.3c00421 journal journal J. Phys. Chem. C \ volume 127 ,\ pages 8978 ( year 2023 ) NoStop

    2023 doi

  6. [14]

    author author T. V. \ Pavlova \ and\ author V. M. \ Shevlyuga ,\ 10.1039/D5CP03294F journal journal Phys. Chem. Chem. Phys. \ volume 27 ,\ pages 23421 ( year 2025 ) NoStop

    2025 doi

  7. [15]

    author author G. W. \ Brown , author B. P. \ Uberuaga , author H. Grube , author M. E. \ Hawley , author S. R. \ Schofield , author N. J. \ Curson , author M. Y. \ Simmons , \ and\ author R. G. \ Clark ,\ 10.1103/PhysRevB.72.195323 journal journal Phys. Rev. B \ volume 72 ,\ p...

    2005 doi

  8. [16]

    Horcas , author R

    author author I. Horcas , author R. Fernández , author J. M. \ Gómez-Rodríguez , author J. Colchero , author J. Gómez-Herrero , \ and\ author A. M. \ Baro ,\ @noop journal journal Rev. Sci. Instrum. \ volume 78 ,\ pages 013705 ( year 2007 ) NoStop

    2007

  9. [17]

    author author J. P. \ Perdew , author K. Burke , \ and\ author M. Ernzerhof ,\ 10.1103/PhysRevLett.77.3865 journal journal Phys. Rev. Lett. \ volume 77 ,\ pages 3865 ( year 1996 ) NoStop

    1996 doi

  10. [18]

    Kresse \ and\ author J

    author author G. Kresse \ and\ author J. Furthm\"uller ,\ 10.1103/PhysRevB.54.11169 journal journal Phys. Rev. B \ volume 54 ,\ pages 11169 ( year 1996 ) NoStop

    1996 doi

  11. [19]

    Kresse \ and\ author D

    author author G. Kresse \ and\ author D. Joubert ,\ 10.1103/PhysRevB.59.1758 journal journal Phys. Rev. B \ volume 59 ,\ pages 1758 ( year 1999 ) NoStop

    1999 doi

  12. [20]

    Tersoff \ and\ author D

    author author J. Tersoff \ and\ author D. R. \ Hamann ,\ @noop journal journal Phys. Rev. B \ volume 31 ,\ pages 805 ( year 1985 ) NoStop

    1985

  13. [21]

    J \'o nsson , author G

    author author H. J \'o nsson , author G. Mills , \ and\ author K. W. \ Jacobsen ,\ title Nudged elastic band method for finding minimum energy paths of transitions , \ in\ 10.1142/9789812839664_0016 booktitle Classical and quantum dynamics in condensed phase simulations ,\ edi...

    1998 doi

  14. [22]

    Rioux , author M

    author author D. Rioux , author M. Chander , author Y. Z. \ Li , \ and\ author J. H. \ Weaver ,\ 10.1103/PhysRevB.49.11071 journal journal Phys. Rev. B \ volume 49 ,\ pages 11071 ( year 1994 ) NoStop

    1994 doi

  15. [23]

    author author J. Y. \ Lee \ and\ author M.-H. \ Kang ,\ @noop journal journal Phys. Rev. B \ volume 69 ,\ pages 113307 ( year 2004 ) NoStop

    2004

  16. [24]

    author author M. W. \ Radny , author P. V. \ Smith , author T. C. G. \ Reusch , author O. Warschkow , author N. A. \ Marks , author H. F. \ Wilson , author N. J. \ Curson , author S. R. \ Schofield , author D. R. \ McKenzie , \ and\ author M. Y. \ Simmons ,\ 10.1103/PhysRevB.7...

    2006 doi

  17. [25]

    Sagisaka , author M

    author author K. Sagisaka , author M. Marz , author D. Fujita , \ and\ author D. Bowler ,\ 10.1103/PhysRevB.87.155316 journal journal Phys. Rev. B \ volume 87 ,\ pages 155316 ( year 2013 ) NoStop

    2013 doi

  18. [26]

    author author J. M. \ Bennett , author O. Warschkow , author N. A. \ Marks , \ and\ author D. R. \ McKenzie ,\ 10.1103/PhysRevB.82.235417 journal journal Phys. Rev. B \ volume 82 ,\ pages 235417 ( year 2010 ) NoStop

    2010 doi

Pith tools

Reviewed July 31, 2026 · model on record in the stance chip above.