Pith. sign in

REVIEW 4 major objections 5 minor 67 references

Quantum effects in surface diffusion: application to diffusion of nitrogen adatoms over GaN(0001) surface

T0 review · 4 major / 5 minor · reviewed 2026-08-09 · deepseek-v4-flash

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

Pith's one-line read The paper claims that on GaN(0001), the Fermi-level position controls nitrogen adatom diffusion, so 1/6 ML gallium coverage lowers the barrier from 1.18 eV to 0.92 eV through electron redistribution.

desk verdict New DFT barriers and a Fermi-level mechanism for N diffusion on Ga-covered GaN(0001), with a plausible interpretation but several unsupported claims that need revision. read the letter →

arxiv 2501.19079 v2 pith:AWPWKU2M submitted 2025-01-31 cond-mat.mtrl-sci

Pawe{l} Strak , Cyprian Sobczak , Stanislaw Krukowski This is my paper · ORCID
classification cond-mat.mtrl-sci
keywords surfacediffusiongalliumnitridenitrogenadatomdensityfunctionaltheoryFermilevelquantumstatisticsresonantbondingnudgedelasticband
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 tries to establish that quantum statistics—specifically the position of the Fermi level and the resulting occupation of surface quantum states—controls the diffusion of nitrogen adatoms on GaN(0001), not just the direct bonding of the moving atom. It claims that on the clean gallium-terminated surface the N adatom moves from a threefold hollow site (H3) over a bridge saddle point with a 1.18 eV barrier, while adding two gallium adatoms, corresponding to 1/6 monolayer coverage, shifts the Fermi level and lowers the barrier to 0.92 eV. At full gallium coverage the stable site changes to on-top, the H3 site becomes the energy maximum, and the barrier is 1.23 eV. The decrease of 0.26 eV at fractional coverage is attributed to electron redistribution between resonant quantum states, not to a fast diffusion channel beneath a gallium overlayer. If this is right, models of Ga-rich molecular-beam epitaxial growth must include electronic-structure effects in the nitrogen incorporation kinetics, because the rate-limiting diffusion step depends on surface Fermi-level pinning.

What carries the argument

The central object is a set of resonant bonding states: four nonorthogonal, fractionally occupied quantum states formed from three nitrogen $|2p\rangle$ orbitals in overlap with four gallium $|4sp^3\rangle$ hybrids. These states carry occupation probabilities such as 3/4 in the bulk and 2/3 for the surface adatom, so their energy alone does not set the diffusion barrier—the Fermi level decides their occupation and the energy cost of losing a bond at the saddle point. The paper locates these states with projected densities of states and Crystal Orbital Hamilton Population (COHP) analysis, and maps diffusion paths with nudged elastic band calculations. The Fermi-level position relative to the resonant and broken-bond states is the mechanism that carries the argument.

What would settle it

A decisive check is to recompute the 1/6 ML N-adatom barrier with a denser k-point grid (e.g., 2×2×1 or 3×3×1) and a thicker slab: if the clean-to-covered barrier difference shrinks below about 0.1 eV, the Fermi-level mechanism is a numerical artifact. A complementary experimental check is to measure the nitrogen adatom diffusion activation energy on GaN(0001) as a function of gallium coverage; the paper's picture requires a minimum near 1/6 ML rather than a monotonic decrease with coverage.

Watch

Extended reading notes

Core claim

In the paper's own terms, the discovery is that GaN is bonded by four resonant states created from three $|N2p\rangle$ orbitals overlapping four $|Ga4sp^3\rangle$ hybrids, with fractional occupation probabilities, and that the same resonant-state physics dictates the surface diffusion barrier. The N adatom at the H3 minimum bonds through two $|N2p\rangle$ orbitals forming three partially occupied resonant states plus a $|N2p_z\rangle$ broken-bond state below the Fermi level. At the bridge saddle point one resonant bond is lost and the resulting broken-bond state's energy relative to the Fermi level sets the barrier. Adding two Ga adatoms at 1/6 ML changes the Fermi-level pinning from the $|Ga4sp^3\rangle$ state to the $|Ga4sp^3\text{-}4p_z\rangle$ bond states, placing that state closer to the Fermi level and reducing the barrier to 0.92 eV. At full 13/12 ML Ga coverage, the N adatom forms saturated bonds in the on-top position; moving it to the H3 site raises the bonding-state energies by about 0.35–0.43 eV, giving a 1.23 eV barrier, and no low-barrier diffusion channel beneath the Ga layer was found.

Load-bearing premise

The load-bearing premise is that a hydrogen-terminated (2√3×2√3) GaN slab with one nitrogen adatom, two gallium adatoms placed in T4 sites, a single k-point, and a 0.005 eV/Å force tolerance represents the real Ga-rich MBE surface closely enough that the 0.26 eV barrier difference is meaningful.

Editorial extensions

If this is right

  • The minimum N-adatom diffusion barrier under Ga-rich MBE conditions is 0.92 eV at 1/6 ML Ga coverage, not 1.18 eV as on the clean surface, so growth models should use the coverage-dependent value.
  • Because the barrier reduction comes from Fermi-level shifts, changing the Fermi level by doping, surface charge, or adsorbate coverage should alter diffusion kinetics even without direct adatom interactions.
  • At full Ga coverage the stable N site is on-top, so incorporation of nitrogen into a new wurtzite layer can proceed from that position; the H3 site is the saddle point instead of the minimum.
  • The earlier proposal of a low-barrier diffusion channel beneath a Ga overlayer is not supported; the full-coverage barrier is 1.23 eV, higher than the fractional-coverage barrier.
  • Surface electron counting rules extended to resonant states with fractional occupation are needed to predict which adsorption site and barrier apply at a given Ga chemical potential.

Reading between the lines

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

  • Editorial extension: the Fermi-level mechanism is general, so any adsorbate whose saddle-point state lies near the Fermi level should show a coverage- or doping-dependent diffusion barrier, and the same analysis could be applied to AlN(0001), InN(0001), and other polar semiconductor surfaces.
  • Editorial extension: because the reported clean-versus-covered barrier difference is only 0.26 eV and the calculations use a single k-point, repeating the NEB with denser k-point sampling and larger supercells would test whether the mechanism survives convergence; the paper does not report such checks.
  • Editorial extension: a direct atomistic prediction worth testing is that the nitrogen jump rate on GaN(0001) should have a non-monotonic dependence on gallium coverage, with fastest diffusion near 1/6 ML, which could be observed in island-density or step-flow growth experiments.
  • Editorial extension: the on-top stability at full coverage implies that N atoms arriving on the Ga adlayer can be incorporated directly beneath them, which would show up as a preference for Ga-polar step advancement during MBE; comparing predicted step-edge incorporation rates with measured growth morphologies could validate this.
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

4 major / 5 minor

Summary. The manuscript reports DFT/NEB calculations of nitrogen adatom diffusion on clean, partially Ga-covered, and fully Ga-covered GaN(0001), yielding barriers of 1.18 eV, 0.92 eV, and 1.23 eV, respectively. It attributes the 0.26 eV lowering at 1/6 ML coverage to Fermi-level shifts and to quantum-statistical redistribution of electrons among resonant surface states, and it reports that at full Ga coverage the N adatom's stable site is on-top rather than H3. The paper also introduces a resonant-bonding description of GaN based on two valence subbands and uses it to interpret the barrier changes.

Significance. If the central claim holds, the paper would establish that N-adatom diffusion on GaN(0001) under Ga-rich MBE conditions is controlled by electron redistribution and Fermi-level pinning rather than by direct adsorbate interactions, with a concrete falsifiable prediction: a minimum barrier of 0.92 eV at 1/6 ML and a stable on-top N site at full coverage. The computed barriers are ab initio total-energy differences rather than fitted values, and the NEB paths are explicit; these are genuine strengths. The resonant-state analysis is an interpretation layer, not the source of the barrier numbers, so the core numbers are not circular. However, the significance is conditional on numerical convergence and on the completeness of the diffusion path search, both of which are currently unresolved.

major comments (4)
  1. [Section II (calculation procedure)] The calculations use a single Monkhorst-Pack 1x1x1 k-point grid, and no convergence tests with respect to k-point sampling are reported. Since the partially covered surface is metallic (the text states the gap closes), the Fermi level and total-energy differences can be sensitive to k-point sampling; the 0.26 eV effect is the same order as typical errors from a 1x1x1 grid. Please report barriers and Fermi-level positions for at least 2x2x1 and 3x3x1 grids for the clean and two-Ga-adatom cases.
  2. [Sections II and III (model)] No convergence tests with respect to slab thickness or vacuum spacing are reported; only a single hydrogen-terminated (2√3 x 2√3) slab is used. The Fermi-level pinning mechanism invoked in Section IV.d depends on surface-state electrostatics, which is exactly what slab-thickness and vacuum tests would check. Please add tests with thicker slabs and larger vacuum, reporting changes in the 0.92 eV barrier and the Fermi-level difference.
  3. [Sections V(d) and V(m)] The negative conclusion that no diffusion channel exists beneath the Ga adlayer is not supported by any NEB calculation. Every reported path has the N adatom on top of the surface or adlayer (Figs. 5, 8, 10, 12); no initial or final configuration places N inside, below, or exchanging with the Ga adlayer. Since NEB connects only the chosen endpoints, it cannot rule out the subsurface/exchange channel proposed in Ref. [7]. The phrase in V(m) that the 0.26 eV reduction has 'no reference to surface diffusion barrier beneath Ga overlayer' is therefore a claim about an unsimulated channel. Please either simulate subsurface and exchange paths, or explicitly limit the 'overall barrier' claim to the on-surface paths.
  4. [Section IV.d] The causal attribution of the 0.26 eV barrier reduction to the Fermi-energy change is asserted from the same DFT output rather than demonstrated by controlled variation. The text states that at the bridge the destabilized state sits about 0.3 eV closer to the Fermi level and that 'this difference is responsible for the observed energy barrier,' but no decomposition of the NEB total-energy difference into state-energy and occupation contributions is provided. A quantitative test is needed: for example, vary the Fermi level by doping or by a background charge, or compute the separate energy contributions, and show the barrier tracks the Fermi level while the geometry change is held fixed.
minor comments (5)
  1. [Section IV.a] The printed resonant states (iii) and (iv) both contain two |N2p_x> terms, and one of them presumably should be |N2p_y>; as written, the normalization statement |<u_i|u_j>|^2 = (3/4) delta_ij cannot hold for i=j with normalized states. Please correct the definitions and state the intended overlap convention.
  2. [Section IV.e] The text says 'H2 site' where H3 is meant; the H3 label is used consistently everywhere else and the maximum-energy site is described as H3.
  3. [Section IV.d] The energy profile is referred to as 'Fig. 2' in several places, but the NEB barriers are plotted in Fig. 5; please correct the cross-references.
  4. [Section V] In the 'state of the art before publication' summary, the two valence-subband picture is cited to Refs. [49,50], but the experimental and ab initio evidence for that picture is presented in Ref. [53] (Magnuson et al.) and Ref. [54]; please update the citation.
  5. [Sections IV.e and V] The full-coverage calculation uses 13 Ga atoms on a 12-site cell, i.e., θ_Ga = 13/12 ML, while the text variously says 'full Ga coverage' and 'full Ga layer coverage'; please define the coverage precisely and state whether the conclusions depend on the extra Ga atom beyond one monolayer.

Circularity Check

0 steps flagged · score 1.0 of 10

No significant circularity: the diffusion barriers are ab initio NEB total-energy differences, and the quantum-statistics/Fermi-level explanation is a post-hoc interpretation rather than an input fitted to produce the barriers.

full rationale

The central numerical results (Delta E_bar = 1.18 eV clean, 0.92 eV at 1/6 ML Ga, 1.23 eV full Ga coverage) are obtained by DFT/NEB as total-energy differences between relaxed endpoints and saddle configurations (Section II and Figs. 5, 12). No parameter is fitted to the barrier values, and no equation reduces a reported barrier to an input definition. The Fermi-level/quantum-statistics discussion (Sections IV.c, IV.d, V(m)) is an explanatory overlay on those total energies: the authors identify which quantum states sit near the Fermi level and how their occupation changes, but they do not use that analysis to compute the barriers. The EECR framework is cited from prior work by the same group (Refs. [23,26,27]), but it is used to rationalize, not to generate, the diffusion energies, so it is not load-bearing circularity. The self-citation to Ref. [7] (Neugebauer et al.) correctly reports prior work and is used as a comparison, not as a uniqueness theorem. The paper contains a genuine evidentiary gap, which the manuscript itself exposes: Section V(d) states 'No diffusion channel under full Ga adlayer was identified,' but every NEB path shown is for an N adatom on or above the Ga adlayer, and no subsurface/exchange path is simulated. That is an unsupported negative conclusion and a correctness risk, but it is not a circular derivation. Similarly, the coarse 1x1x1 k-point grid and single-slab setup are convergence/representativeness concerns, not circularity. Thus the derivation chain is self-contained in the sense that the headline numbers do not reduce to their inputs; the quantum-mechanism claim is under-tested but not definitionally forced.

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

No numerical free parameters were fitted to data; the DFT barriers are ab initio outputs, and the coverage values are chosen physical scenarios, not fitted parameters. The central interpretive burden rests on the EECR framework and the resonant bonding model, both of which are introduced or self-cited by the authors. No new physical particles, forces, or conserved quantities are introduced; the resonant states are a basis-set interpretation of existing DFT wavefunctions.

assumptions (5)
  • domain assumption Born-Oppenheimer approximation and static NEB barriers ignore vibrational entropy and zero-point corrections.
    Section II states NEB is used to find minimum energy paths; no vibrational contributions are included, yet barriers are reported as absolute eV values.
  • domain assumption DFT with GGA-PBE (PBEJsJrLO), Troullier-Martins pseudopotentials, triple-zeta basis, and 1x1x1 k-point grid is sufficiently converged for claimed 0.01 eV barrier differences.
    Section II gives numerical settings but no convergence tests with respect to k-points, basis, or slab thickness; the 0.26 eV reduction is a small energy difference.
  • domain assumption The (2√3 x 2√3) slab with hydrogen termination and a compensating dipole layer represents the GaN(0001) surface under MBE Ga-rich conditions.
    Sections II and III; coverage configurations with 1, 2, 3, 4, and 13 Ga adatoms are chosen by hand. Deviations from experimental surface reconstructions would change barriers.
  • domain assumption The electron counting rule (ECR/EECR) from the authors' prior work correctly assigns occupation of surface quantum states.
    Sections IV.c and IV.d rely on ECR to conclude which states are filled and to attribute the barrier change to Fermi level position; references [23,26,27] are self-citations.
  • ad hoc to paper The resonant bonding model: four nonorthogonal states from three N2p orbitals with occupation probability 3/4 describe Ga-N bonding.
    Section IV.a introduces states |u1>-|u4> and the P=3/4 rule; this is the paper's new interpretive framework, not independently established.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Quantum effects in surface diffusion: application to diffusion of nitrogen adatoms over GaN(0001) surface." pith.science (2026). https://pith.science/paper/AWPWKU2M

@misc{pith2026250119079,
  author       = {Pith},
  title        = {Pith review of: Quantum effects in surface diffusion: application to diffusion of nitrogen adatoms over GaN(0001) surface},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/AWPWKU2M}},
  note         = {Machine review of arXiv:2501.19079}
}
read the original abstract

It is shown that quantum effects play determining role in nitrogen adatom diffusion due to several different factors. This could be related to the change of the energy of the quantum states and also due to the redistribution of electrons between the quantum states, both full and resonant, via quantum statistics partially governed by the Fermi energy level. These effects were studied in the case of nitrogen diffusion over clean and gallium covered Ga-terminated GaN(0001) surface. For the fractional coverage the density functional theory (DFT) calculations show that at the saddle point configuration the redistribution of electrons between different quantum states may affect the surface diffusion barrier significantly. The other quantum influence occurs via the change of the minimal energy configuration. Under fractional Ga coverage of GaN(0001) surface the nitrogen diffusion energy barrier proceeds from the resonant states governed energy minimal H3 site across the saddle point in the bridge configuration. At this path the barrier is affected the electron redistribution between surface quantum states both in the initial and the saddle point. In the case of the full GaN coverage the diffusion path is from on-top N adatom configuration via H3 site that corresponds to maximal energy. Therefore the diffusion barrier is Ebar= 1.18 eV for clean and Ebar= 0.92 eV for (1/6) ML to finally Ebar= 1.23 eV for full Ga coverage. Thus the overall barrier is reduced to Ebar= 0.92 eV due to quantum statistics effects. The identified stable N on-top configuration for the full coverage is essential for atomic mechanism of GaN growth in Ga-rich regime.

Figures

Figures reproduced from arXiv: 2501.19079 by the authors.

Figure 1
Figure 1. The factors affecting the effective barrier for diffusion via modification of the energy of the quantum states (red line) and the resulting total energy change (green line) in the course of the jump between two sites: a) standard case: the energy of quantum state remains far below Fermi energy, b) the energy and the occupation of quantum states, standard and [PITH_FULL_IMAGE:figures/full_fig_p009_1.png] view at source ↗
Figure 4
Figure 4. Energies of the quantum states of nitrogen adatom and neighboring gallium atoms states [PITH_FULL_IMAGE:figures/full_fig_p014_4.png] view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

67 extracted references · 46 canonical work pages

  1. [7]

    Neugebauer, T

    J. Neugebauer, T. K. Zywietz, M. Scheffler, J. E. Northrup, H. Chen, R.M. Feenstra, Adatom Kinetics On and Below the Surface: The Existence of a New Diffusion Channel. Phys. Rev. Lett. 90, 056101 (2003). https://doi.org/10.1103/PhysRevLett.90.056101

    work page doi:10.1103/physrevlett.90.056101 2003

  2. [1]

    Eyring, The activated complex in chemical reactions

    H. Eyring, The activated complex in chemical reactions. J. Chem. Phys. 3, 107 (1935). https://doi.org/10.1063/1.1749604

    doi:10.1063/1.1749604 1935

  3. [2]

    Ganeshan, L

    S. Ganeshan, L. G. Hector, Z. K. Liu, First -principles calculation of impurity diffusion coefficients in dilute Mg alloys using the 8 -frequency model. Acta Mater. 59,3214 (2011). https://doi.org/10.1016/j.actamat.2011.01.062

    work page doi:10.1016/j.actamat.2011.01.062 2011

  4. [3]

    Hartmann, Crystal Growth

    P. Hartmann, Crystal Growth. An Introduction, N. Holland Amsterdam 1973

    work page 1973

  5. [4]

    Rosenberger, Fundamentals of Crystal Growth, Springer Berlin 1979

    F. Rosenberger, Fundamentals of Crystal Growth, Springer Berlin 1979

    work page 1979

  6. [5]

    J. P. van der Eerden, Crystal growth mechanisms. In Handbook of Crystal Growth I: Fundamentals. Ed by D. T. J. Hurle N. Holland Amsterdam 1993

    work page 1993

  7. [6]

    Gomer, Diffusion of adsorbates on metal surfaces

    R. Gomer, Diffusion of adsorbates on metal surfaces. Rep. Prog. Phys. 53, 917 (1990). https://doi.org/10.1088/0034-4885/53/7/002

    work page doi:10.1088/0034-4885/53/7/002 1990

  8. [8]

    T. A. Witten, L. M. Sander, Diffusion-Limited Aggregation, A kinetic Critical Phenomenon. Phys. Rev. Lett. 47, 1400 (1981). https://doi.org/10.1103/PhysRevLett.47.1400

    doi:10.1103/physrevlett.47.1400 1981

Show all 67 references
  1. [9]

    B. B. Mandelbrot, The Fractal Geometry of Nature. Freeman, San Francisco 1982

    1982

  2. [10]

    Vicsek Fractal Growth Phenomena, World Scientific, Singapore 1992

    T. Vicsek Fractal Growth Phenomena, World Scientific, Singapore 1992

    1992

  3. [11]

    Meakin, Formation of fractal clusters and networks by irreversible Diffusion- Limited Aggregation

    P. Meakin, Formation of fractal clusters and networks by irreversible Diffusion- Limited Aggregation. Phys. Rev. Lett. 51, 1119 (1983). https://doi.org/10.1103/PhysRevLett.51.1119

    1983 doi

  4. [12]

    J. Krug, P. Meakin, Columnar Growth in Oblique-Incidence Ballistic Deposition – Faceting, Noise-Reduction and Mean-Field Theory. Phys. Rev. A. 43, 900 (1991). https://doi.org/10.1103/PhysRevA.43.900

    1991 doi

  5. [14]

    Krukowski, J

    S. Krukowski, J. C. Tedenac, Fractal to compact transition during growth of 2D Kossel crystal in vapor diffusion field. J. Cryst. Growth 160, 167 (1996). https://doi.org/10.1016/0022-0248(95)00892-6

    1996 doi

  6. [15]

    T. Taki, M. Strassburg, Visible LEDs: More than Efficient Light. ECS J. Solid State Sci. Technol. 9, 015017 (2020). https://www.doi.org/10.1149/2.0402001JSS

    2020 doi

  7. [16]

    Kneissl, T

    M. Kneissl, T. Y. Seong, J. Han, H. Amano, The emergence and prospects of deep- ultraviolet light-emitting diode technologies. Nature Photonics 13, 233-244 (2019). https://www.doi.org/10.1038/s41566-019-0359-9

    2019 doi

  8. [17]

    I. M. Watson, Metal organic vapour phase epitaxy of AlN, GaN, InN and their alloys: A key chemical technology for advanced device applications. Coord. Chem. Rev. 257 2120 – 2141 (2013). https://doi.org/10.1016/j.ccr.2012.10.020

    2013 doi

  9. [18]

    T. D. Moustakas, B. Paiella, Optoelectronic device physics and technology of nitride semiconductors from the UV to terahertz. Rep. Prog. Phys. 80, 106501 (2017). https://doi.org/10.1088/1361-6633/aa7bb2

    2017 doi

  10. [19]

    Skierbiszewski, Z

    C. Skierbiszewski, Z. R. Wasilewski, I. Grzegory, S. Porowski, Nitride-based laser diodes by plasma-assisted MBE- from violet to green emission. J. Cryst. Growth 311, 1632- 1639 (2008). https://doi.org/10.1016/j.jcrysgro.2008.12.040

    2008 doi

  11. [20]

    Koukitu, Y

    A. Koukitu, Y. Kumagai, Hydride Vapor Phase Epitaxy of GaN in Technology of Gallium nitride crystal growth, ed by D. Ehrentraut, E. Meissner, M. Bockowski, Springer Verlag 2010, pp 31 – 60

    2010

  12. [21]

    Grzegory, J

    I. Grzegory, J. Jun, M. Bockowski, S. Krukowski, M. Wroblewski, B. Lucznik, S. Porowski, III-V nitrides – thermodynamics and crystal growth at high N2 pressure. J. Phys. Chem. Solids, 56, 639 – 647 (1995). http://dx.doi.org/10.1016/0022- 3697(94)00257-6

    1995 doi

  13. [22]

    Krukowski, I

    S. Krukowski, I. Grzegory, M. Bockowski, B. Lucznik, G. Nowak, J. Borysiuk, M. Wroblewski, M. Leszczynski, P. Perlin, S. Porowski, J. L. Weyher, Growth of AlN, GaN and InN from the solution. Int. J. Mater. Prod. Technol. 22, 226 – 261 (2004). https://doi.org/10.1504/IJMPT.2005.005766

    2004 arXiv

  14. [23]

    Strak, W

    P. Strak, W. Miller, S. Krukowski, Surface charge induced 𝑠𝑝3 − 𝑠𝑝2 transformation: new energy optimization universal mechanism of (4 × 4) reconstruction of stoichiometric GaN(0001) surface. Materials 17, 4581 (2024). https://doi.org/10.3390/ma17184581 38

    2024 doi

  15. [24]

    D. J. Chadi, Atomic structure of GaAs(100)-(2x1) and (2x4) reconstructed surfaces. J. Vac. Sci. Technol. A 5, 834 (1987). https://doi.org/10.1116/1.574366

    1987 doi

  16. [25]

    M. D. Pa shley, Electron counting model and its application to island structures on molecular-beam epitaxy grown GaAs(001) and ZnSe(001) . Phys. Rev. B 40, 10481 (1989). https://doi.org/10.1103/PhysRevB.40.10481

    1989 doi

  17. [26]

    Kempisty, S

    P. Kempisty, S. Krukowski, Adsorption of ammonia at GaN(0001) surface in the mixed ammonia/hydrogen ambient - a summary of ab initio data, AIP. Adv. 4, 117109 (2014). http://dx.doi.org/10.1063/1.4901922

    2014 doi

  18. [27]

    Kempisty, P

    P. Kempisty, P. Strąk, K. Sakowski, S. Krukowski, Thermodynamics of GaN(s)- NH3(v)+N2(v)+H2(v) system - electronic aspects of the processes at GaN(0001) Surface, Surf. Sci. 63, 12-33 (2017). https://doi.org/10.1016/j.susc.2017.03.007

    2017 doi

  19. [28]

    Brocks, P

    G. Brocks, P. J. Kelly, R. Car, Energetics of adatoms on the Si(100) surface, Surf. Sci. 269, 860 (1992). https://doi.org/10.1016/0039-6028(92)91362-F

    1992 doi

  20. [29]

    K. N. He, X. S. Kong, C. S. Liu. Robust activation energy predictions of solute diffusion in machine learning method. Comp. Mater. Sci. 184, 109948 (2020). https://doi.org/10.1016/j.commatsci.2020.109948

    2020

  21. [30]

    Garcia et al

    A. Garcia et al. SIESTA: Recent developments and applications. J. Chem. Phys. 152, 204108 (2020). https://doi.org/10.1063/5.0005077

    2020 doi

  22. [31]

    Junquera, O

    J. Junquera, O. Paz, D. Sanchez-Portal, E. Artacho, Numerical atomic orbitals for linear-scaling calculations. Phys. Rev. B 64, 235111 (2001). https://doi.org/10.1103/PhysRevB.64.235111

    2001 doi

  23. [32]

    Anglada, J

    E. Anglada, J. M. Soler, J. Junquera, E. Artacho, Systematic generation of finite-range atomic basis sets for linear-scaling calculations. Phys. Rev. B 66, 205101 (2002). https://doi.org/10.1103/PhysRevB.66.205101

    2002 doi

  24. [33]

    Troullier, J

    N. Troullier, J. L. Martins, Efficient pseudopotentials for plane-wave calculations. Phys. Rev. B, 43, 1993-2006. (1991). https://doi.org/10.1103/PhysRevB.43.1993 39

    1991 doi

  25. [34]

    Troullier, J

    N. Troullier, J. L. Martins, Efficient pseudopotentials for plane-wave calculations. II. Operators for fast iterative diagonalization, Phys. Rev. B, 43, 8861-8869 (1991) . https://doi.org/10.1103/PhysRevB.43.8861

    1991 doi

  26. [35]

    H.J Monkhorst, J. D. Pack, Special points for Brillouin-zone integrations. Phys. Rev. B, 13, 5188 - 5192 (1976). https://doi.org/10.1103/PhysRevB.13.5188

    1976 doi

  27. [36]

    L. S. Pedroza, A. J. R. da Silva, K. Capelle, Gradient-dependent density functionals of the Perdew-Burke-Ernzerhof type for atoms, molecules, and solids, Phys. Rev. B 79, 201106(R) (2009). https://doi.org/10.1103/PhysRevB.79.201106

    2009 doi

  28. [37]

    M. M. Odashima, K. Capelle, S. B. Trickey, Tightened Lieb-Oxford Bound for Systems of Fixed Particle Number, J. Chem. Theory Comput. 5, 798 - 807 (2009) . https://doi.org/10.1021/ct8005634

    2009 doi

  29. [38]

    Leszczynski, , H

    M. Leszczynski, , H. Teisseyre, T. Suski, I. Grzegory, M. Bockowski, J. Jun, S. Porowski, K. Pakula, J.M. Baranowski, C.T. Foxon, T.S. Cheng, Lattice parameters of gallium nitride Appl. Phys. Lett. 69, 73 (1996). https://doi.org/10.1063/1.118123

    1996 doi

  30. [39]

    Ferreira, M

    L.G. Ferreira, M. Marques, and L.K. Teles, Approximation to density functional theory for the calculation of band gaps of semiconductors. Phys. Rev. B 78, 125116, (2008). https://doi.org/10.1103/PhysRevB.78.125116

    2008 doi

  31. [40]

    Ribeiro, L.R.C

    M. Ribeiro, L.R.C. Fonseca, and L.G. Ferreira, Accurate prediction of the Si/SiO2 interface band offset using the self-consistent ab initio DFT/LDA-1/2 method. Phys. Rev. B, 79, 241312 (2009). https://doi.org/10.1103/PhysRevB.79.241312

    2009 doi

  32. [41]

    Monemar, J.P

    B. Monemar, J.P. Bergman, I. A. Buyanova, H. Amano, I. Akasaki, T. Detchprohm, K. Hiramatsu, and N. Sawaki, The excitonic bandgap of GaN: Dependence on substrate. Solid State Electron. 41, 239 (1997). https://doi.org/10.1016/S0038- 1101(96)00208-0

    1997 doi

  33. [42]

    Yeo, T.C

    Y.C. Yeo, T.C. Chong, M.F. Li, Electronic band structures and effective-mass parameters of wurtzite GaN and InN. J. Appl. Phys. 83, 1429 (1997). https://doi.org/10.1063/1.366847

    1997 doi

  34. [43]

    Bengtsson, Dipole correction for surface supercell calculations

    L. Bengtsson, Dipole correction for surface supercell calculations. Phys. Rev. B 59 (1999) 12301 – 12304. https://doi.org/10.1103/PhysRevB.59.12301

    1999 doi

  35. [44]

    Neugebauer and M

    J. Neugebauer and M. Scheffler, Adsorbate-substrate and adsorbate-adsorbate interactions of Na and K adlayers on Al(111). Phys. Rev. B 46, 16067 (1992). https://doi.org/10.1103/PhysRevB.46.16067 40

    1992 doi

  36. [45]

    Krukowski, P

    S. Krukowski, P. Kempisty, P. Strak, Electrostatic condition for the termination of the opposite face of the slab in DFT simulations of semiconductor surfaces. J. Appl. Phys. 105, 113701 (2009). http://dx.doi.org/10.1063/1.3130156

    2009 doi

  37. [46]

    Henkelman, B.P

    G. Henkelman, B.P. Uberuaga and H. Jonsson, A climbing image nudged elastic band method for finding saddle points and minimum energy paths. J. Chem. Phys. 113, 9901 – 9904 (2000). https://doi.org/10.1063/1.1329672

    2000 doi

  38. [47]

    Henkelman, H

    G. Henkelman, H. Jonsson, Improved tangent estimate in the nudged elastic band method for finding minimum energy paths and saddle points. J. Chem. Phys. 111, 9978 - 9985 (2000). https://doi.org/10.1063/1.1323224

    2000 doi

  39. [48]

    Sheppard, R

    D. Sheppard, R. Terrell, G. Henkelman, Optimization method for finding minimum energy path. J. Chem. Phys. 128, 134106 (2008). http://dx.doi.org/10.1063/1.2841941

    2008 doi

  40. [49]

    Hrytsak, P

    R. Hrytsak, P. Kempisty, E. Grzanka, M. Leszczynski, M. Sznajder, DFT study on point defects migration through the pseudomorphic and lattice-matched InN/GaN interfaces. Comp. Mater. Sci. 186, 110039 (2021). https://doi.org/10.1016/j.commatsci.2020.110039

    2021

  41. [50]

    Smalc-Koziorowska, E

    J. Smalc-Koziorowska, E. Grzanka, A. Lachowski, R. Hrytsak, M. Grabowski, S. Grzanka, S. Kret, R. Czernecki, H. Turski, L. Marona, T. Markut, T. Schulz, M. Albrecht, M. Leszczynski, Role of Metal Vacancies in the Mechanism of Thermal Degradation of InGaN Quantum Wells, ACS App...

    2021 doi

  42. [51]

    Hrytsak, P

    R. Hrytsak, P. Kempisty, M. Leszczynski, M. Sznajder, Modeling of the Point Defect Migration across the AlN/GaN Interfaces-Ab Initio Study. Materials. 15, 478 (2022). https://doi.org/10.3390/ma15020478

    2022 doi

  43. [52]

    Hrytsak, P

    R. Hrytsak, P. Kempisty, E. Grzanka, M. Leszczynski, M. Sznajder, Identification of the Kirkendall effect as a mechanism responsible for thermal decomposition of the InGaN/GaN MQWs system. New J. Phys. 24, 123007 (2022). https://doi.org/10.1088/1367-2630/aca698 41

    2022 doi

  44. [53]

    Magnuson, F

    M. Magnuson, F. Mattesini, C. Hoglund, J. Birch, L. Hultman, Electronic structure of GaN and Ga investigated by soft x-ray spectroscopy and first-principles methods. Phys. Rev. B 82, 085125 (2010). https://dor.org/10.1103/PhysRevB.81.085125

    2010 doi

  45. [54]

    Vacuum 99, 166 (2013)

    Maria Ptasinska, Jakub Sołtys, Jacek Piechota, Stanisław Krukowski, Electronic properties on GaN(0001) surface – ab initio investigation. Vacuum 99, 166 (2013). http://dx.doi.org/10.1016/j.vacuum.2013.05.006

    2013 doi

  46. [55]

    F. A. Kekule, Untersuchungen uber aromatische Verbindungen. Liebigs Annalen der Chemie und Pharmacie 137, 129 (1866). https://doi.org/10.1002%2Fjlac.18661370202 (in German)

  47. [56]

    P. A. Schultz, R. P. Messmer, Are There  Bonds in Benzene? Phys. Rev. Lett. 58, 2416 (1987). https://doi.org/10.1103/PhysRevLett.58.2416

    1987 doi

  48. [57]

    P. W. Anderson, The Resonating Bond State in LA2CuO4 and Superconductivity. Science 235, 1196 (1987). https://doi.org/10.1126/science.235.4793.1196

    1987

  49. [58]

    Weber, Variational Study of Strongly Correlated Electronic Models

    C. Weber, Variational Study of Strongly Correlated Electronic Models. PhD thesis. EPFL Lausanne 2007

    2007

  50. [59]

    Dronskowski, P

    R. Dronskowski, P. E. Bloechl, Crystal orbital Hamilton populations (COHP): energy- resolved visualization of chemical bonding in solids based on density-functional calculations. J. Phys. Chem. 97, 8617 (1993). https://doi.org/10.1021/j100135a014

    1993 doi

  51. [60]

    V. L. Deringer, A. L. Tchougreeff, R. Dronskowski, Crystal Orbital Hamilton Population (COHP) Analysis as Projected from Plane-Wave Basis Sets. J. Phys. Chem. 115, 5461 (2011). https://doi.org/10.1021/jp202489s

    2011 doi

  52. [62]

    Bernardini, V

    F. Bernardini, V. Fiorentini, and D. V anderbilt, Spontaneous polarization -and piezoelectric constants of III-V nitrides, Phys. Rev. B 56 (2001) R10024

    2001

  53. [63]

    Bernardini, V

    F. Bernardini, V. Fiorentini, D. Vanderbilt, Accurate Calculation of Polarization-Related Quantities in Semiconductors, Phys. Rev. B 63 (2001) 193201

    2001

  54. [64]

    C. E. Dreyer, A. Janotti, C. G. Van de Walle, D. Vanderbilt, Correct Implementation of Polarization Constant in Wurtzite Materials and Impact on III -Nitrides, Phys. Rev. X 6 (2016) 021038. 42

    2016

  55. [65]

    Strak, P

    P. Strak, P. Kempisty, K. Sakowski, A. Kaminska, D. Jankowski, K. P. Korona, K. Sobczak, J. Borysiuk, M. Beeler, E. Grzanka, E. Monroy, S. Krukowski, Ab initio and experimental studies of polarization and polarization related fields in nitrides and nitride structures, AIP Adv....

    2017

  56. [66]

    Krukowski, Z

    S. Krukowski, Z. Romanowski, I. Grzegory, S. Porowski, Interaction of N2 molecule with liquid Ga surface-quantum mechanical calculations (DFT). J. Cryst. Growth 189-190, 159 (1998). http://dx.doi.org/10.1016/S0022-0248(98)00195-X

    1998 doi

  57. [67]

    Krukowski, Z

    S. Krukowski, Z. Romanowski, I. Grzegory, S. Porowski, Surface reaction of nitrogen with liquid group III metals. J. Chem Phys. 114, 6353 ( 2001). http://dx.doi.org/10.1063/1.1355984

    2001 doi

  58. [68]

    Kempisty, K

    P. Kempisty, K. Kawka, A. Kusaba, Y. Kangawa, Polar gaN surfaces under Gallium Rich Conditions: Revised Thermodynamic Insights from Ab Initio Calculations. Materials 16, 5982 (2023). https://doi.org/10.390/ma16175982

    2023

  59. [69]

    Kempisty, P

    P. Kempisty, P. Strak, K. Sakowski, S. Krukowski, DFT study of ammonia desorption from the GaN(0001) surface covered with a NH3 /NH2 mixture. J.Cryst. Growth 403, 105 (2014). http://dx.doi.org/10.1016/j.jcrysgro.2014.06.016

    2014 doi

Pith tools

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