Pith. sign in

REVIEW 3 major objections 6 minor 7 references

Emergence of a Bandgap in Nano-Scale Graphite: A Computational and Experimental Study

T0 review · 3 major / 6 minor · reviewed 2026-08-12 · deepseek-v4-flash

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

Pith's one-line read Nano-patterned graphite opens a ~112 meV bandgap, and the paper attributes it to tensile strain from patterning.

desk verdict A reproducible ARPES gap in patterned HOPG is worth taking seriously, but the strain explanation is a post hoc fit, not a prediction. read the letter →

arxiv 2411.14244 v4 pith:CJXJANN7 submitted 2024-11-21 cond-mat.mtrl-sci

Sujinda Chaiyachad , Trung-Phuc Vo , Warakorn Jindata , Sirisak Singsen , Tanachat Eknapakul , Chutchawan Jaisuk , Patrick Le Fevre , Francois Bertran
show 7 more authors
Donghui Lu Yaobo Huang Hideki Nakajima Watchara Liewrian Ittipon Fongkaew Jan Minar Worawat Meevasana
This is my paper · ORCID
classification cond-mat.mtrl-sci
keywords graphitebandgapengineeringangle-resolvedphotoemissionspectroscopytensilestrainHOPGdensityfunctionaltheoryterahertzoptoelectronicsfocusedionbeampatterning
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

Conventionally, graphite is a gapless semimetal, which blocks it from semiconductor uses such as terahertz photodetectors. This paper reports that nanoscale square patterns cut into highly oriented pyrolytic graphite (HOPG) with a focused ion beam open a bandgap of about 112 meV at the Dirac point, seen with angle-resolved photoemission spectroscopy. The authors argue that the gap comes from tensile strain left behind by the patterning rather than from quantum confinement, disorder, or layer twist: the Raman 2D band shifts $4\,\mathrm{cm}^{-1}$ to lower frequency, a sign of stretched carbon rings. Two independent calculations reproduce the effect, one with a stretched top layer and one with biaxially strained trilayer graphene. If the strain picture holds, structural patterning becomes a route to terahertz-scale gaps in thick graphitic carbon.

What carries the argument

The mechanism is strain-induced bandgap opening at the graphite $K$ point, captured by two independent calculations. In the first, a semi-infinite HOPG crystal is treated with a Green's-function multiple-scattering formalism combined with the one-step model of photoemission, which computes the measured intensity including matrix-element and surface effects; a displacement of the top layer along the in-plane $b$ direction by $0.080$ Å opens a gap in the simulated ARPES spectrum. In the second, a trilayer-graphene unit cell is stretched biaxially in a plane-wave DFT calculation, and $1.1\%$ tensile strain produces a gap near the measured value. On the experimental side, the gap is read from symmetrized energy distribution curves, in which the ARPES intensity at $k_F$ is mirrored about the Fermi level to expose the missing spectral weight.

What would settle it

Measure the strain field of the patterned squares directly, for example by nanoscale X-ray diffraction or a calibrated Raman strain map. If the measured strain is much smaller than the 0.080 Å shift or the 1.1% biaxial values used in the calculations while the 112 meV gap persists, the strain origin is falsified; if the strain matches, the claim is supported.

Watch

Extended reading notes

Core claim

The central claim is that a bandgap of $112 \pm 15$ meV opens at the $K$ point of nano-patterned HOPG, a material normally considered gapless, and that the gap is caused by tensile strain introduced during focused-ion-beam patterning. The evidence is a density-of-states peak about 59 meV below the Fermi level in the energy distribution curve and a clear gap in the symmetrized ARPES spectrum of the patterned squares, while the unpatterned substrate and nearby regions remain gapless; the same gap was reproduced at a second synchrotron. Raman spectroscopy shows a $4\,\mathrm{cm}^{-1}$ redshift of the 2D band, which the authors read as a sign of stretched carbon rings. Calculations support the strain interpretation: stretching the topmost layer by $0.080$ Å in a semi-infinite one-step photoemission calculation opens a 137 meV gap, and a trilayer-graphene DFT calculation under $1.1\%$ biaxial tensile strain gives 93 meV. The authors therefore conclude that strain, rather than quantum confinement, disorder, or twist, is the origin of the observed gap.

Load-bearing premise

The load-bearing premise is that the patterned sample really is stretched by roughly the amount used in the calculations—0.080 Å in the top layer or about 1.1% biaxial strain—yet the paper only infers this from a $4\,\mathrm{cm}^{-1}$ Raman redshift and never measures the strain magnitude.

Editorial extensions

If this is right

  • If the strain picture is right, the bandgap is not tied to the few-nanometer quantum-confinement regime, so comparatively large patterned graphitic structures can behave as semiconductors.
  • A roughly 100 meV gap corresponds to terahertz and mid-infrared photon energies, making nano-patterned graphite a candidate platform for terahertz detectors and emitters.
  • The two calculated strain configurations (0.080 Å top-layer shift and 1.1% biaxial strain) suggest that strain magnitude is a tuning parameter for the gap size.
  • Because the one-step photoemission calculations reproduce the measured ARPES spectra only when matrix-element and surface effects are included, the method offers a way to predict photoemission signatures of strained layered materials.

Reading between the lines

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

  • The paper leaves implicit that the actual strain state of the patterned sample is not measured; the agreement between 137 meV, 93 meV, and 112 meV depends on strain values chosen after the fact, so a direct strain measurement is the decisive test.
  • If strain is the mechanism, the gap should be systematically tunable by changing pattern geometry, ion dose, or pattern size, and it should disappear when strain relaxes; the paper does not report such a control experiment.
  • The symmetrized-EDC analysis assumes a symmetric background, so a tunneling or transport measurement of the same samples would independently confirm whether the 112 meV feature is a true gap rather than a matrix-element or lineshape effect.
  • Because the generalized-gradient DFT approximation tends to underestimate gaps, the closeness of the 93 and 137 meV values to experiment may be partly fortuitous; hybrid or GW calculations on the same strained geometries would give cleaner quantitative predictions.
Share X Bluesky LinkedIn Reddit HN

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. This manuscript reports the observation of a ~112 meV bandgap in focused-ion-beam patterned nano-scale HOPG using ARPES at two synchrotron beamlines, alongside Raman spectroscopy and two DFT-based theoretical approaches. The one-step KKR photoemission calculation with a distorted top layer (0.080 Å) yields a gap of 137 meV, and a VASP trilayer calculation under 1.1% biaxial tensile strain yields 93 meV; both are compared with the measured 112 meV gap. The paper attributes the gap to tensile strain induced by patterning, citing a 4 cm^-1 redshift of the Raman 2D band as evidence of strain.

Significance. The experimental demonstration of a bandgap in a quasi-bulk graphitic system via nanopatterning is interesting and potentially relevant for THz and optoelectronic devices, and the observation is strengthened by reproducibility at two beamlines. The theoretical machinery—one-step photoemission with matrix-element effects and TR-LEED final states—is sophisticated and appropriate, and the paper explicitly offers machine-checkable computational results. However, the central quantitative claim that the calculations 'predict' the observed gap is not supported, because the strain parameters are not measured but chosen post hoc to match the experiment. The only strain indicator is a small Raman shift that is not calibrated, and standard calibrations suggest strain values far below those used in the calculations. The paper's value lies mainly in the experimental gap observation and the plausibility of the strain mechanism, not in a parameter-free quantitative prediction.

major comments (3)
  1. [§3, Figs. 4 and 5] The strain values used in the two calculations—0.080 Å top-layer distortion and 1.1% biaxial strain—are not determined from any measurement on the patterned sample, and the text itself states only that a 'small amount of tensile strain may be present' based on the 4 cm^-1 redshift of the 2D Raman band. Because the computed gaps (137 and 93 meV) are then presented as agreeing with the measured 112 meV, the quantitative agreement is a post hoc fit over the strain parameters, not a prediction. The abstract's claim of 'predict this bandgap with remarkable accuracy' is therefore overstated. Using standard biaxial strain calibrations for the graphene 2D band (approximately -60 to -140 cm^-1 per 1% strain), the observed 4 cm^-1 shift corresponds to roughly 0.03–0.07% strain, one to two orders of magnitude smaller than the 1.1% used in the VASP calculation; the manuscript should either provide a mechanism for local strain enhancement within the 300 nm squares or explicitly treat the strain values as illustrative rather than measured.
  2. [§3, Fermi-level offset] The Fermi level in the KKR/DFT calculation is shifted by an ad hoc 0.08 eV offset to align the calculated valence-band top with the ARPES spectrum (page 14). This offset is an additional free parameter in the comparison, so the agreement between the calculated 137 meV and the measured 112 meV is obtained by tuning both the distortion and the offset. The sensitivity of the computed gap to this offset should be reported or the offset should be justified independently.
  3. [§3, kz dependence] The statement that the bandgap does not close at any k point (page 20) is based on a three-layer slab at a single strain value (1.1%), and Fig. 6 shows only the dispersion along kz, without a quantitative definition of the gap criterion or a demonstration that the gap remains open over the entire Brillouin zone. The claim should be made precise, and the same check should be reported for the 0.080 Å distortion model, since that model is the one used in the one-step KKR comparison with the ARPES data.
minor comments (6)
  1. [References] Reference [63] appears three times with different papers, reference [18] is omitted from the list, and the sentence 'energy gap of graphite is found in order of ~30 – 40 eV' (page 15) should read meV.
  2. [Fig. 3 inset and Raman analysis] The fitting procedure for the 2D Raman band is described only in a caption; the number of Lorentzian components, the fitting range, and the uncertainty of the 4 cm^-1 redshift should be stated.
  3. [§2, FIB parameters] The source voltage for the focused ion beam is stated as 300 kV (page 4); this is likely a typo for 30 kV and should be corrected.
  4. [§3, symmetrized EDC analysis] The extraction of the 112 ± 15 meV gap from the symmetrized EDC in Fig. 2m should be described explicitly, including the functional form used for the fit and how the uncertainty is derived; the fitted curve should be overlaid on the data.
  5. [Fig. 4b] The notation for the distortion values is inconsistent: the text uses negative values for compression and positive for stretching, but Fig. 4b labels are given as (0.155, -0.105, -0.006) and (0.042); please clarify the sign convention in the figure and caption.
  6. [Fig. 5(c)] The comparative plot of energy gaps for monolayer and trilayer graphene should include error bars or a sensitivity estimate for the strain values, and the statement that trilayer gap 'increases roughly' should be quantified.

Circularity Check

2 steps flagged · score 6.0 of 10

The claimed quantitative prediction of the 112 meV gap is obtained by scanning unmeasured strain parameters and retaining the values (0.080 Å, 1.1%) whose calculated gaps land near 112 meV; the qualitative strain-gap mechanism is independent, but the numerical agreement is a selection effect.

  1. fitted input called prediction [Section 3, discussion of Fig. 4 (one-step KKR/SPRKKR calculation)]
    "To investigate the distortion dependence of HOPG electronic structure, one C atom is shifted close to the other along b direction by distorted values Δ(Å) ∈ [−0.155,−0.006] with step 0.025 Å ... From computational point of view, stretching the topmost layer by 0.080 Å along b direction creates a bandgap Eg = 137 meV meanwhile Eg = 112 meV is observed in ARPES measurements."

    The strain parameter Δ is scanned over a range, and the specific value 0.080 Å is then used for the comparison with experiment. The paper does not measure or derive this distortion on the patterned sample; the only experimental strain indicator is a 4 cm−1 Raman 2D redshift, which the paper leaves unquantified. The calculated gap of 137 meV is therefore not a prediction from an independently determined strain state: the distortion magnitude was selected so that the resulting gap would be in the neighborhood of the measured 112 meV. The quantitative agreement is a selection from the scan, not an independent confirmation.

  2. fitted input called prediction [Section 3, discussion of Fig. 5 (VASP trilayer graphene calculation)]
    "In order to study band-gap opening, we considered the electronic band structure of trilayer graphene compared with monolayer graphene in the presence of biaxial tensile strain varying from 0 to 2% ... computational investigation reveals that the highest gap opening of trilayer graphene is under 1.1% strain with an energy gap of 93.0 meV, which is close to the measured energy gap of HOPG (112 meV)."

    The biaxial strain is scanned from 0% to 2%, and 1.1% is singled out because it gives the largest calculated gap (93.0 meV), which the paper then states is close to the measured 112 meV. Because the strain actually present in the FIB-patterned nano-HOPG is never quantified — the 4 cm−1 redshift of the Raman 2D band is only interpreted as a possibility of tensile strain — the 1.1% value is a fitting parameter chosen to make the calculation match the experiment. The qualitative strain-gap trend is real, but the reported numerical agreement is fixed by parameter choice rather than by an ab initio or experimentally constrained prediction.

full rationale

The ARPES observation of a 112 ± 15 meV gap in nano-HOPG is an independent experimental result and is not itself circular. The qualitative conclusion that tensile strain can open a gap in graphitic systems is also supported by the two calculations and by prior literature. However, the abstract's claim that the calculations 'predict this bandgap with remarkable accuracy' is not supported by the derivation chain. In the KKR one-step calculation, the distortion Δ is scanned and the value 0.080 Å is then used because it produces Eg = 137 meV, close to the measured 112 meV. In the VASP calculation, the biaxial strain is scanned from 0% to 2% and 1.1% is selected because it produces the largest gap, 93.0 meV, again close to the measured value. The paper never measures or derives either strain value; the only strain evidence is a 4 cm−1 redshift of the Raman 2D band, which the paper treats only as a qualitative possibility of tensile strain. Thus the quantitative agreement between theory and experiment for the gap magnitude reduces to a post hoc choice of free parameters. This is the fitted-input-called-prediction pattern: the strain is the fitted input, and the gap value is the target. I found no load-bearing self-citation or uniqueness-imported-from-authors circularity; the cited one-step model and VASP methods are standard independent computational frameworks. The experimental gap observation and the qualitative strain mechanism retain independent content, so the partial circularity score is 6 rather than higher.

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

The central quantitative agreement rests on two strain parameters (0.080 Å and 1.1%) that are tuned to match the measured gap, plus an ad hoc Fermi-level offset. The idealized model geometries (one distorted layer, trilayer slab) are not validated against measured strain or stacking. No new physical entities are introduced.

free parameters (3)
  • KKR distortion Delta = 0.080 Å (stretching along b)
    Chosen to produce a 137 meV gap matching the measured 112 meV (Fig. 4c).
  • Biaxial tensile strain = 1.1%
    Selected from a 0-2% scan because it gives a 93 meV trilayer gap near the measured 112 meV (Fig. 5c).
  • Fermi-level offset = +0.08 eV
    Added to align the calculated valence band top with ARPES; does not affect gap size but is an ad hoc energy alignment.
assumptions (4)
  • domain assumption DFT at LDA/PBE level adequately describes the band structure of graphite/graphene
    Used in both SPRKKR and VASP; known to underestimate gaps, and the paper invokes this to explain discrepancies.
  • domain assumption The one-step model of photoemission with TR-LEED final states reproduces ARPES intensities for this system
    Needed to compare calculated and measured ARPES maps (Section 2).
  • ad hoc to paper A single distorted top layer (or a trilayer slab) represents the patterned HOPG surface
    The actual disorder and strain distribution from FIB is not measured; these are idealized geometries.
  • domain assumption A 4 cm^-1 redshift of the Raman 2D band indicates tensile strain
    Based on prior Raman-strain studies (Ref 47); the magnitude of strain is not derived from the shift.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Emergence of a Bandgap in Nano-Scale Graphite: A Computational and Experimental Study." pith.science (2026). https://pith.science/paper/CJXJANN7

@misc{pith2026241114244,
  author       = {Pith},
  title        = {Pith review of: Emergence of a Bandgap in Nano-Scale Graphite: A Computational and Experimental Study},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/CJXJANN7}},
  note         = {Machine review of arXiv:2411.14244}
}
read the original abstract

Bandgaps in layered materials are critical for enabling functionalities such as tunable photodetection, efficient energy conversion, and nonlinear optical responses, which are essential for next-generation photonic and quantum devices. Gap engineering could form heterostructures with complementary materials like transition metal dichalcogenides or perovskites for multi-functional devices. Graphite, conventionally regarded as a gapless material, exhibits a bandgap of ~100 meV in nano-scale patterned highly oriented pyrolytic graphite (HOPG), as revealed by angle-resolved photoemission spectroscopy (ARPES) and Raman measurements. Our state-of-the-art calculations, incorporating photoemission matrix element effects, predict this bandgap with remarkable accuracy and attribute it to mechanical distortions introduced during patterning. This work bridges theory and experiment, providing the direct evidence of a tunable bandgap in HOPG. Beyond its fundamental significance, this finding opens new possibilities for designing materials with tailored electronic properties, enabling advancements in terahertz devices and optoelectronics.

Figures

Figures reproduced from arXiv: 2411.14244 by the authors.

Figure 4
Figure 4. For illustrative purposes, the unit cell of the distored and pristine HOPG are superimposed (Fig. 4a), implying the structural difference between a pristine and distorted case with AB stacking sequence. In this configuration, the second C sub-layer of the first layer (from top) contains distortion. To investigate the distortion dependence of HOPG electronic structure, one C atom is shifted close to the other along 𝑏… view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

7 extracted references · 7 canonical work pages

  1. [1]

    Introduction Two-dimensional graphene exhibits remarkable properties with a wide range of applications, including high-mobility metallic transistors [1], magneto-electronic devices [2], quantum computing devices based on quantum electrodynamics [3], and terahertz (THz) optoelectronics [4,5]. For the THz devices, THz radiation can pass through media withou...

    work page 2004

  2. [10]

    R. R. Nair, P. Blake, A. N. Girgorrenko, K. S. Novoselov, T. J. Booth, T. Stauber, N. M. R. Peres, A. K. Geim, Fine structure constant defines visual transparency of graphene, Science 320 (2008) 1308-1308. [11] T.D. Nguyen, S. Liu, G. Kumar, A. Nahata, Z. V. Vardeny, Terahertz plasmonic properties of highly oriented pyrolytic graphite. Appl. Phys. Lett. 1...

    work page 2008

  3. [19]

    K. A. Ritter, L. W. Lyding, The influence of edge structure on the electronic properties of graphene quantum dots and nanoribbons, Nat. Mater. 8 (2009) 235-242. [20] H. Ebert, D. Ködderitzsch, J. Minár, Calculating condensed matter properties using the KKR-Green's function method—recent developments and applications, Rep. Prog. Phys. 74 (2011) 096501. [21...

    work page Pith review arXiv 2009

  4. [36]

    Grimme, S Ehrlich, L

    S. Grimme, S Ehrlich, L. Goerigk, Effect of the damping function in dispersion corrected density functional theory. J. Comput. Chem. 32 (2011) 1456-1465. [37] H. J. Monkhorst, J. D. Pack, Special points for Brillouin-zone integrations, Phys. Rev. B 13 (1976) 5188. [38] V. V. Karasiev, Communication: Simple and accurate uniform electron gas correlation ene...

    work page 2011

  5. [44]

    Thiyagarajan, A

    K. Thiyagarajan, A. Ananth, B. Saravanakumar, Y. S. Mok, S. J. Kim, Defect-induced metallic-to-semiconducting transition in multilayer graphene, RSC Adv. 5 (2015). 16821-16827. [45] M. Hashimoto, T. Yoshida, A. Fujimori, D. H. Lu, Z.-X. Shen, M. Kubota, K. Ono, M. Ishikado1, K. Fujita, S. Uchida, Effects of out-of-plane disorder on the nodal quasiparticle...

    work page 2015

  6. [53]

    García, P

    N. García, P. Esquinazi, J. Barzola-Quiquia, S. Dusari, Evidence for semiconducting behavior with a narrow band gap of Bernal graphite. New J. Phys. 14 (2012) 053015. [54] A. Zhamu, B. Z. Jang, U.S. Patent No. 10,822,725. Washington, DC: U.S. Patent and Trademark Office (2020). [55] V. Pantin, J. Avila, M. A. Valbuena, P. Esquinazi, M. E. Dávila, M. C. As...

    work page 2012

  7. [61]

    Mori-Sánchez, A

    P. Mori-Sánchez, A. J. Cohen, W. Yang, Localization and delocalization errors in density functional theory and implications for band-gap prediction, Phys. Rev. Lett 100 (2008) 146401. [62] B. S. Kim, J. W. Rhim, B. Kim, C. Kim, S. R. Park, Determination of the band parameters of bulk 2H-MX2 (M= Mo, W; X= S, Se) by angle-resolved photoemission spectroscopy...

    work page 2008

Pith tools

Reviewed August 12, 2026 · model on record in the stance chip above.