Pith. sign in

REVIEW 3 major objections 4 minor 2 references

Boron Arsenide Heterostructures: Lattice-Matched Heterointerfaces, and Strain Effects on Band Alignments and Mobility

T0 review · 3 major / 4 minor · reviewed 2026-08-14 · deepseek-v4-flash

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

Pith's one-line read Boron arsenide, already prized for its thermal conductivity, is predicted to get faster and more useful under tensile strain: in-plane electron and hole mobilities rise to 2417 and 3550 cm²/Vs, and the material becomes a…

desk verdict Useful predictive study of strained BAs and two heterojunctions, but the headline hole-mobility gain mixes a no-SOC strained calculation with an SOC-included unstrained baseline and should be reframed before publication. read the letter →

arxiv 1909.01449 v1 pith:E5BA6RNP submitted 2019-09-03 cond-mat.mtrl-sci

Kyle Bushick , Sieun Chae , Zihao Deng , John Heron , Emmanouil Kioupakis This is my paper · ORCID
classification cond-mat.mtrl-sci
keywords boronarsenidebiaxialstraincarriermobilityelectron-phononcouplingbandalignmenttype-IIheterojunctionInGaNZnSnN2
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

Boron arsenide (BAs) is usually studied for its record thermal conductivity; this paper asks what it can do electronically and as a heterostructure partner. Using hybrid-functional density functional theory and first-principles electron-phonon calculations, the authors argue that 1% biaxial tensile strain raises the in-plane electron mobility from 1341 to 2417 cm²/Vs and the hole mobility from 1387 to 3550 cm²/Vs at 300 K. The strain works by splitting the conduction-band valleys, so electrons move with a lighter transverse mass, and by splitting the valence-band maximum, so holes lose scattering channels. The same calculations place BAs as a nearly lattice-matched substrate for In-rich InGaN and ZnSnN2, with type-II band alignments that would separate photogenerated carriers. If these predictions hold, strained BAs films and BAs-based heterojunctions become concrete options for fast transistors, full-visible-spectrum LEDs, and solar cells.

What carries the argument

The central object is the strain-split band edge of cubic BAs. In the unstrained crystal the six conduction minima near X are equivalent, so transport averages a light transverse mass and a heavy longitudinal mass; 1% biaxial strain lowers the out-of-plane pair of valleys by about 200 meV, so in-plane electrons carry only the lighter transverse mass ($m^* = 0.24$), while out-of-plane electrons see the heavier mass ($m^* = 1.09$). At the valence edge, strain removes the degeneracy of the top two bands by about 115 meV, cutting the hole interband scattering channels from two to one. For heterojunctions, the key instrument is the vacuum-level alignment of slab electrostatic potentials, applied directly to BAs and ZnSnN2 and taken from published hybrid-functional results for InGaN; that alignment produces the type-II offsets.

What would settle it

Grow a BAs film on a substrate imposing about 1% biaxial tensile strain, measure the 300 K in-plane and out-of-plane electron and hole mobilities (predicted 2417/482 and 3550/2956 cm²/Vs), and check the predicted valley splitting (roughly 200 meV) and valence-band splitting (roughly 115 meV) by optical or angle-resolved photoemission measurements. If the mobilities do not increase and the band splittings are absent, the central claim fails. For the alignments, measure the valence-band offset at a BAs/InGaN interface by X-ray photoelectron spectroscopy; a value far from 0.30 eV would overturn the type-II prediction.

Watch

Extended reading notes

Core claim

The paper's central discovery is that tensile strain is a powerful lever on boron arsenide's electronic performance, not just its phonons. At 1% biaxial tensile strain, the in-plane electron mobility rises from 1341 to 2417 cm²/Vs and the in-plane hole mobility from 1387 to 3550 cm²/Vs at 300 K, an 80% and 68% increase, because strain splits the conduction-band valleys, so electrons move with a lighter transverse effective mass, and splits the valence-band maximum, so holes lose an interband scattering channel. The paper also finds that BAs is nearly lattice-matched to In-rich InGaN and to ZnSnN2, with type-II (staggered) band alignments: conduction- and valence-band offsets of 0.42/0.30 eV against In0.5Ga0.5N and 0.503/0.516 eV against ordered ZnSnN2. These numbers make p-BAs/n-InGaN and p-BAs/n-ZnSnN2 junctions plausible for carrier-separating optoelectronic devices on a substrate that also removes heat efficiently.

Load-bearing premise

The reported offsets assume that the way the paper lines up each material's energy bands to the vacuum level is consistent across BAs, ZnSnN2 (computed with a fitted mixing parameter), and InGaN (taken from an earlier calculation); if those references disagree, the offsets and even the type-II classification could shift.

Editorial extensions

If this is right

  • At 1% tensile strain, in-plane electron and hole mobilities reach 2417 and 3550 cm²/Vs, an 80% and 68% improvement over unstrained BAs, making strained BAs a candidate channel for fast, low-power transistors.
  • The out-of-plane electron mobility falls to 482 cm²/Vs, so strained BAs is strongly anisotropic; this matters for any vertical device geometry and is a direct fingerprint of the valley-repopulation mechanism.
  • BAs substrates introduce less than 1% misfit strain for InGaN near 47% In, which would let InGaN devices reach red wavelengths without the dislocation problem that limits GaN substrates.
  • Type-II alignments with In0.5Ga0.5N and ZnSnN2 mean photogenerated electrons and holes separate across the junction, supporting photovoltaic and photodetector designs on a high-thermal-conductivity substrate.
  • Strain of 4% can lower the BAs band gap from 1.78 eV to roughly 0.89 to 1.15 eV, a range that could supply the 1.1 eV subcell in a multijunction solar cell.

Reading between the lines

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

  • A direct experimental test is to measure the mobility anisotropy and the roughly 200 meV valley splitting of a strained BAs film; if the anisotropy and splitting do not appear together, the proposed mechanism, rather than the mobility numbers, would be the point of failure.
  • The same valley-splitting logic that works for BAs and strained silicon should apply to other cubic multi-valley semiconductors, so the calculation could be repeated for boron phosphide or boron antimonide to see whether strain is a general mobility lever in the boron pnictides.
  • The paper stops at intrinsic junctions; a natural extension is to model the actual space-charge region and doping levels at BAs/ZnSnN2 interfaces, where charged defects in ZnSnN2 may bend bands and shift the effective barrier from the ideal intrinsic offsets.
  • The InGaN offsets come from an external calculation rather than a fresh alignment in this paper's framework, so a direct photoemission measurement of the BAs/InGaN valence-band offset is the cleanest way to confirm or revise the type-II assignment.
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 / 4 minor

Summary. The manuscript uses hybrid-functional DFT and EPW electron-phonon calculations to study biaxially strained boron arsenide (BAs), reporting that tensile strain reduces the band gap, increases in-plane electron and hole mobilities (to 2417 and 3550 cm2/Vs at 1% strain), and that BAs forms nearly lattice-matched, type-II heterojunctions with InGaN and ZnSnN2. The authors propose BAs as both a strain-engineered thin-film channel material and a high-thermal-conductivity substrate for optoelectronic junctions.

Significance. If the quantitative claims hold, this work would provide a useful predictive framework for BAs heterostructures and identify specific device-relevant interfaces. The methodological strengths are the converged EPW mobility calculations with Brillouin-zone extrapolation, the direct comparison to prior electron-mobility results from Liu et al., and the physically clear mechanism attributed to valley and valence-band degeneracy lifting. The band-alignment section also addresses an important gap in the literature. However, the central hole-mobility enhancement claim rests on a cross-framework comparison (no-SOC strained vs. SOC-included unstrained), and the band offsets are assembled from three different computational frameworks, so the quantitative conclusions need revision or careful qualification.

major comments (3)
  1. [Section 3.2, Figure S5, and Conclusion] The claimed 68% hole-mobility increase is not a same-framework comparison. The paper's own unstrained hole mobility, computed without spin-orbit coupling, is 1387 cm2/Vs, and the strained in-plane value is 3550 cm2/Vs, which is a 156% increase. The 68% figure instead compares 3550 cm2/Vs to the SOC-included unstrained literature value 2110 cm2/Vs from Liu et al. The argument in Section 3.2 that SOC should not strongly affect the strained hole mobility because the strained VBM is non-degenerate with or without SOC (Figure 3(e)-(f)) is an inference, not a quantitative calculation; SOC can renormalize effective masses and electron-phonon matrix elements even when degeneracy is already lifted. Because the abstract and conclusion advertise the 'more than 60%' improvement, this load-bearing number must either be computed with SOC in the strained case or be restated with explicit caveats about the different baselines.
  2. [Section 2 and Section 3.3] The reported band offsets (0.42/0.30 eV for InGaN, 0.503/0.516 eV for ZnSnN2) are assembled from three different reference schemes: BAs is aligned to vacuum using an unrelaxed HSE slab, InGaN data are taken from ref. [47] without recomputation, and ZnSnN2 is aligned using an HSE mixing parameter of 31% fitted to the experimental gap. If the vacuum-level references are not mutually consistent, the offsets and the type-II classification could shift. The manuscript should either recompute InGaN in the same slab framework, test the sensitivity of the ZnSnN2 offsets to the mixing parameter, or explicitly quantify the expected uncertainty from these procedural differences. This point is central because the proposed heterojunction applications depend on the type-II alignment and the absolute offset magnitudes.
  3. [Section 3.2, temperature dependence] The temperature-dependent mobility analysis compares strained no-SOC results with the SOC-included literature data of Liu et al. in Figure 4(a)-(b), and the text states that the power laws in the 200-500 K range are 'similar to the values found by Liu et al.' This is an additional instance of the same cross-framework issue: the unstrained no-SOC hole mobility differs from the SOC-included value by a factor of about 1.5, so any quantitative comparison of scattering mechanisms or power-law exponents should be made on a consistent footing. The qualitative conclusions about acoustic versus optical deformation potential may survive, but the current presentation conflates two different computational settings.
minor comments (4)
  1. [Abstract and Section 1] The opening sentence 'BAs is III-V semiconductor' is missing an article; it should read 'BAs is a III-V semiconductor.'
  2. [Supporting Information, Figure S4/S5 references] The Supporting Information text contains apparent figure-reference mismatches: the ordered ZnSnN2 slab potential is said to be shown in 'Figure S4,' but Figure S4 is a band-structure plot, and the disorder-alignment discussion refers to 'Figure S5' for electrostatic potentials while Figure S5 is the mobility convergence plot. The numbering should be corrected.
  3. [Figure 5(b)] The abscissa of Figure 5(b) mixes two different variables (In content x for InGaN and long-range order parameter S for ZnSnN2) on a single axis, which makes the plot difficult to read; a two-panel layout or separate axes would improve clarity.
  4. [Section 2] The mobility calculations use LDA while the band-structure and band-alignment calculations use HSE; the manuscript should state explicitly whether the strained geometries used in the two parts are consistent, since the strain relaxation procedure is described separately for each code.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: the strain-mobility and band-alignment results are obtained from direct DFT/EPW calculations with independent external benchmarks, not from parameters that encode the predicted quantities.

full rationale

The central claims are computed rather than fitted. The unstrained and 1%-strained electron and hole mobilities come from EPW calculations on relaxed LDA geometries, with no parameter adjusted to reproduce any mobility target; the unstrained electron value is benchmarked against the independent Liu et al. calculation (1400 cm2/Vs). The strain dependence of the band gap and absolute band positions is obtained directly from HSE slab calculations, and the InGaN alignment is imported from an external calculation (ref. 47) rather than being derived from this paper's inputs. The 31% HSE mixing parameter used for ZnSnN2 is calibrated to the experimental gap following refs. 28 and 35, but the reported band offsets (0.503/0.516 eV) are not equal to that gap by construction; they also depend on the vacuum-level electrostatic alignment and on the independently calculated BAs band positions. This is a functional-calibration caveat, not a circular reduction. Self-citations (refs. 4, 11, 28) supply supporting inputs such as p-type doping, effective-mass values, and the mixing convention, but none of these is the load-bearing derivation of the claimed predictions. Therefore no step in the derivation chain reduces by construction to its own input. The 68% hole-mobility enhancement should be read with the caveat that the strained no-SOC value is ratioed to the SOC-included unstrained literature baseline (3550 vs 2110 cm2/Vs), while the same-framework unstrained no-SOC value is 1387 cm2/Vs; this is a comparison-consistency issue, not circularity.

Assumptions & free parameters 1 free parameters · 6 assumptions · 0 invented entities

The BAs strain and mobility claims are parameter-free DFT and EPW calculations benchmarked against Liu et al. for unstrained mobilities. The main fitted element is the 31% HSE mixing parameter for ZnSnN2, which affects only the heterojunction alignment portion. The lattice-match argument assumes Vegard behavior in InGaN, and the band alignment procedure assumes a shared, method-insensitive vacuum reference.

free parameters (1)
  • HSE mixing parameter for ZnSnN2 = 0.31
    Chosen to match the experimental band gap of ZnSnN2 (Section 2, Methodology); this enters the calculated ZnSnN2 band edges and hence the reported ZnSnN2/BAs band offsets. It is a fitted number, not derived.
assumptions (6)
  • domain assumption Vegard's law accurately interpolates InGaN lattice constants between GaN and InN.
    Used in Section 3.3.1 to identify In compositions near 0.47 as lattice-matched to BAs with less than 1% misfit; any nonlinearity in the alloy lattice constant would shift the matched composition.
  • domain assumption Unrelaxed slab calculations of the plane-averaged electrostatic potential give the bulk-to-vacuum level offset needed for absolute band alignment.
    Used in Section 2 to align BAs (and ordered ZnSnN2) bands to vacuum; surface relaxation, reconstruction, and quantum confinement in the 8 to 12 layer slabs are neglected.
  • domain assumption Spin-orbit coupling does not significantly change the strained hole mobility because the strained valence band maximum is non-degenerate.
    Invoked in Section 3.2 to justify comparing a no-SOC strained hole mobility (3550 cm2/Vs) with an SOC-included unstrained baseline (2110 cm2/Vs from ref. 13); the size of the SOC correction in the strained case is not computed.
  • domain assumption LDA accurately describes the band structure and electron-phonon coupling needed for mobility calculations in BAs.
    All mobility values in Section 3.2 are computed with LDA and no SOC; LDA underestimates band gaps, and the accuracy of the electron-phonon matrix elements is assumed.
  • ad hoc to paper HSE06 with a mixing parameter of 31% gives accurate band edges for ordered ZnSnN2.
    The mixing parameter is tuned to reproduce the experimental gap (as stated in Section 2), and the resulting band edges are then used to compute the ZnSnN2/BAs offsets; this is an ad hoc choice specific to this study.
  • domain assumption Cation disorder in ZnSnN2 changes the bulk-to-vacuum electrostatic potential by less than 150 meV, so ordered and disordered structures can be aligned via DFT eigenvalues.
    Stated in the Supporting Information (Figure S3 discussion); the 150 meV variation is treated as negligible for the reported offsets of about 0.5 eV.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Boron Arsenide Heterostructures: Lattice-Matched Heterointerfaces, and Strain Effects on Band Alignments and Mobility." pith.science (2026). https://pith.science/paper/E5BA6RNP

@misc{pith2026190901449,
  author       = {Pith},
  title        = {Pith review of: Boron Arsenide Heterostructures: Lattice-Matched Heterointerfaces, and Strain Effects on Band Alignments and Mobility},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/E5BA6RNP}},
  note         = {Machine review of arXiv:1909.01449}
}
read the original abstract

BAs is III-V semiconductor with ultra-high thermal conductivity, but many of its electronic properties are unknown. This work applies predictive atomistic calculations to investigate the properties of BAs heterostructures, such as strain effects on band alignments and carrier mobility, considering BAs as both a thin film and a substrate for lattice-matched materials. The results show that strain decreases the band gap independent of sign or direction. In addition, biaxial tensile strain increases the in-plane electron and hole mobilities by more than 60% compared to the unstrained values due to a reduction of the electron effective mass and of hole interband scattering. Moreover, BAs is shown to be nearly lattice-matched with InGaN and ZnSnN2, two important optoelectronic semiconductors with tunable band gaps by alloying and cation disorder, respectively. The results predict type-II band alignments and determine the absolute band offsets of these two materials with BAs. The combination of the ultra-high thermal conductivity and intrinsic p-type character of BAs, with its high electron and hole mobilities that can be further increased by tensile strain, as well as the lattice-match and the type-II band alignment with intrinsically n-type InGaN and ZnSnN2 demonstrate the potential of BAs heterostructures for electronic and optoelectronic devices.

Figures

Figures reproduced from arXiv: 1909.01449 by the authors.

Figure 1
Figure 1. Potential configurations of BAs in semiconductor devices. (a) Epitaxially grown BAs thin film as part of a transistor. Epitaxially straining BAs (inset) increases both the electron and the hole mobility. (b) Schematic of an optoelectronic device utilizing the junction between a thermally conducting BAs substrate for efficient heat extraction in conjunction with nearly-lattice-matched direct-band-gap semiconductor fi… view at source ↗
Figure 2
Figure 2. (a) Calculated band structure of unstrained BAs along X to G and G to Z using HSE hybrid functional considering spin-orbit effects. The X and Z directions are equivalent due to the cubic symmetry. (b) The effect of biaxial strain within the (100) plane on the band structure of BAs. The upper panel shows the splitting of the degeneracy of the CBM valleys near X and Z under strain due to the broken symmetry. Tensile s… view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

2 extracted references · 2 canonical work pages

  1. [3]

    Results and Discussion 3.1. Strain Effects on the Band Structure and Absolute Band Positions We first examine the effects of strain on the band structure along the G–X and G–Z directions, which include the conduction band minimum (CBM) and valence band maximum (VBM). These results, calculated with the HSE hybrid functional and spin-orbit effects, are show...

    work page Pith review arXiv 2018

  2. [30]

    Joubert, Phys

    D. Joubert, Phys. Rev. B 1999, 59, 1758. [31] J. Heyd, G. E. Scuseria, M. Ernzerhof, J. Chem. Phys. 2003, 118, 8207. [32] G. Kresse, J. Hafner, Phys. Rev. B 1993, 47, 558. [33] G. Kresse, J. Furthmüller, Comput. Mater. Sci. 1996, 6, 15. [34] G. Kresse, J. Furthmüller, Phys. Rev. B - Condens. Matter Mater. Phys. 1996, 54, 11169. [35] S. Chen, P. Narang, H....

    work page 1999

Pith tools

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