Kinetic Monte Carlo simulation of the nitridation of the GaAs (100) surfaces

A. P. Castro; H. W. Leite Alves

REVIEW

Kinetic Monte Carlo simulation of the nitridation of the GaAs (100) surfaces

A. P. Castro, H. W. Leite Alves

Other Condensed Matter · Submitted 5 Mar 2006 · 16 references

Pith's verdict: Openly preliminary but genuinely new site-resolved adsorption and KMC trends for N on GaAs(100); the central mechanism is a well-labeled guess and the KMC model is too thin to carry the roughness conclusion.

T0 · AI review

How this review was made: An AI referee (deepseek-v4-flash) read the paper on August 28, 2026 and answered a fixed list of questions about it. T0 is the starting level every paper gets. The levels describe how deeply Pith has checked the paper, never the authors.

T0 AI review An AI read the paper and answered a fixed list of questions.
T1 Audit Automatic checks that the references and files the paper relies on exist and match.
T2 Compute The paper's code or calculations re-run with the same results.
T3 Formal statement The claims written out in a formal math language a computer can read.
T4 Proof Those claims proved, and the proof checked by computer.

How Pith checks papers →

· 4 major objections · 4 minor · Ask about this paper · The compiled record · Read the full review →

Challenge this review Comment (0) Read on arXiv ↗

Cite ▾

arXiv IDcond-mat/0603111 (version 1)

Pith IDMJI5A7HB, Pith's permanent record of this paper and its review

CitationA. P. Castro, H. W. Leite Alves. "Kinetic Monte Carlo simulation of the nitridation of the GaAs (100) surfaces." arXiv:cond-mat/0603111 (2006). https://pith.science/paper/MJI5A7HB

Cite this reviewBibTeX for Pith's review of this paper

Claim paper

Why Pith says this

Highlight any text to comment on it.

Editor's note

read the letter

Quick take: this is a small, openly preliminary computational paper on N adsorption at the As-terminated GaAs(100) (2x1) surface. The specific site-resolved adsorption energies (C site most stable at 3.17 eV, D and E about 1.1 eV higher) and the KMC island-density trends are new, and the authors are unusually candid about what they have and haven't done. The soft spots are real, though: no convergence checks, a KMC model with only three processes and a fixed prefactor, and a central mechanism that is labeled as a guess.

The strongest part is the DFT adsorption energetics. The authors identify that N prefers the unoccupied C, D and E sites, that adsorption over the As dimer is 1.48 eV less stable, and that, unlike on GaAs(110), there is no adsorption barrier. The structural analysis (As dimer tilt at the C site, bond lengths) is plausible and useful as a starting point. The KMC results—high island density and roughness at low flux—are a natural consequence of the very high diffusion barriers, but that's exactly where the problem lies.

The KMC model includes only N adsorption at site C and diffusion along two directions, with a common 1e13 s^-1 prefactor. The authors explicitly exclude diffusion over As dimers because dimers act as traps. What's missing is any test of sensitivity to the process list. No As2 desorption, N desorption, dimer reformation, exchange, or defect-mediated pathways are included. If even one low-barrier process exists, the effective mobility and island density could change by orders of magnitude. The barriers themselves (5.5 and 7.34 eV) come from a single set of static calculations with no k-point/cutoff/slab convergence checks, and one of them (60 eV) is acknowledged to be physically meaningless. More importantly, the KMC uses those same DFT barriers to produce the island densities, so the claimed agreement between KMC and DFT is not independent confirmation.

The anion-exchange mechanism—N pushes the As dimer out, leading to As2 desorption—is presented as a supposition, with calculations "currently underway." That's fine for a preliminary report, but it's not a result.

What's genuinely good: the authors are upfront. They call the 60 eV barrier nonsense, state the model limitations in footnotes, and don't oversell the mechanism. That honesty is rare and should be respected. For a reader working on III-nitride growth or GaAs surface chemistry, these numbers are a reasonable starting point. But I wouldn't cite them as established, and the roughness conclusion isn't a general result.

My recommendation: send this to peer review, but the referee should require convergence tests, a sensitivity analysis of the KMC process list, and either a defended prefactor or an argument that the conclusions are insensitive to it. If the authors can do that, this could become a solid contribution. As is, it's a promising but incomplete preliminary study.

Referee report

Major objections (4)

[Section II, footnote [15]] The manuscript reports a 60.0 eV diffusion barrier for N over As dimers and simultaneously states that "We cannot affirm that this energy barrier has such value, which it does not make any sense." This is an explicit admission that the numerical result is an artifact, yet the barrier is retained and used as the basis for treating As dimers as traps and for excluding dimer-crossing diffusion from the KMC model (footnote [16]). Because the completeness of the KMC process set is load-bearing for the island-density and roughness conclusions, the unphysical result must not be kept as evidence. The authors should either compute this barrier with a proper transition-state search that allows substrate relaxation, or explicitly treat the omitted process as an assumption and test how its inclusion changes the conclusions.
[Section III, Figs. 6-7] The KMC model contains only N adsorption at site C and N diffusion along the [-1 1 0] and [1 1 0] directions, with a common prefactor of 10^13 s^-1. No As2 desorption, N desorption, dimer reformation, or exchange pathways are included, even though the paper's own anion-exchange mechanism (Section II) involves displacement of the As dimer and eventual As2 desorption. Without a sensitivity analysis showing that adding any of these processes leaves the predicted island density and roughness unchanged, the central claim that low N mobility produces high island density and rough surfaces is underdetermined.
[Section II, Fig. 5 and III] The diffusion barriers (7.34 eV and 5.50 eV) are extracted from a single set of static DFT-LDA total-energy calculations with the surface atoms kept fixed at their relaxed positions, and no convergence tests are reported for the k-point mesh, plane-wave cutoff, or slab thickness. Since these barriers enter the KMC rates exponentially, the absence of any uncertainty estimate leaves the absolute island densities unquantified. At minimum, the numerical setup should be documented and convergence checks should be provided.
[Section III, last paragraph] The statement that the KMC results are "in agreement with our obtained DFT-LDA results" is not informative, because the KMC input rates are derived from exactly those DFT results. The agreement is therefore an internal-consistency check rather than an independent confirmation. The text should be revised to avoid implying independent validation, which it cannot provide.

Minor objections (4)

[Abstract and Section II] The abstract and conclusion state that N atoms prefer the unoccupied sites "close to the As dimer," but the distinction among sites C, D and E and their energy ordering (C most stable, D and E less stable by about 1.1 eV) should be stated explicitly in both places.
[Section II, after Fig. 4] The anion-exchange scenario is introduced with "we suppose" and "Calculations to check this model are currently underway." This is a hypothesis and should be clearly labeled as such in the abstract, where it currently reads as a deduced conclusion.
[Section III] The use of a common prefactor of 10^13 s^-1 (Ref. [8]) for all processes should be justified for N on GaAs(100), or at least stated as a standard-order-of-magnitude assumption.
[General] The text repeatedly describes the results as "preliminary" (abstract, Sections I and III). If this is intended as a short communication, the scope should be stated clearly in the introduction so that the reader understands the level of support expected for each claim.

How Pith checks papers → · Full audit of this paper →

The paper

The paper claims

On the As-terminated GaAs(100) $(2\times1)$ surface, nitrogen binds most strongly at the unoccupied C site, pushes the arsenic dimer out of the surface, and, because N diffusion is extremely slow, drives anion exchange and rough island…

The reading

This paper uses first-principles total-energy calculations and kinetic Monte Carlo simulations to establish how nitrogen atoms begin to nitridize the arsenic-terminated GaAs(100) $(2\times1)$ surface, the surface present during high-temperature growth of cubic III-nitrides. It finds that N strongly prefers the unoccupied C site near the arsenic dimer, with an adsorption energy of $3.17$ eV, and that adsorption distorts and pushes the dimer outward, setting up the anion-exchange mechanism in which N replaces As and eventually As$_2$ desorbs. The same calculations give very high diffusion barriers along the available surface directions, so N is nearly immobile compared with Ga or As$_2$. In the Monte Carlo model built from those rates, low N mobility produces high island densities and rough surfaces at low flux and around $0.3$ ML coverage, with efficient nitridation at low flux only below 700 K.

What carries the argument

The central machinery is a two-stage pipeline: parameter-free density functional theory total-energy calculations for N adsorbed at six candidate sites on a five-layer As-terminated GaAs(100) $(2\times1)$ slab, giving relaxed bond lengths, bond angles, adsorption energies, and diffusion barriers; and a solid-on-solid kinetic Monte Carlo simulation that takes the C-site adsorption energy and the two allowed diffusion paths as rate inputs, with a common $10^{13}\,\mathrm{s}^{-1}$ prefactor, to evolve island statistics over coverage and flux. The DFT energies carry the chemical claim (site preference, dimer push-out), and the barriers carry the kinetic claim (low mobility, rough growth).

Extended reading notes

Core claim

On the As-terminated GaAs(100) $(2\times1)$ surface, the paper finds that N adsorption is site-selective: the C site is the most stable ($3.17$ eV adsorption energy), followed by D and E (about $1.09$ and $1.12$ eV less stable), all unoccupied positions adjacent to the As dimer. The adsorbate distorts the dimer, and at C the dimer tilts and is pushed outward, so a second arriving N can complete the anion exchange and expel As$_2$. Diffusion is nearly frozen: $7.34$ eV along $[-1\,1\,0]$, $5.50$ eV along $[1\,1\,0]$, with the As dimer acting as an immobile trap and only the C-E-C-B-C and E-D-E paths available. Feeding these energetics into a kinetic Monte Carlo model with a common $10^{13}\,\mathrm{s}^{-1}$ prefactor yields high island density, tall islands, and rough surfaces at low flux and around $0.3$ ML coverage, matching the claim that slow N mobility controls nitridation morphology.

Load-bearing premise

The kinetic Monte Carlo conclusions rest on assuming that nitridation is fully described by just three events, namely N adsorption at the C site and N hops along the [-1 1 0] and [1 1 0] paths, with rates from static DFT energies and a single $10^{13}\,\mathrm{s}^{-1}$ prefactor, and with no desorption, dimer reformation, or defect-mediated processes.

Editorial extensions

If this is right

Nitridation should nucleate at the unoccupied C sites, making the As-dimer rows a template for N islands rather than a barrier.
When a second N arrives at a C site whose dimer has already been pushed out, As is expelled as As$_2$, providing a microscopic route for the experimentally observed surfactant effect of As in cubic III-nitride growth.
Because N diffusion is so slow, island density at fixed coverage will be much higher than for Ga or As$_2$ diffusion, and surfaces will roughen instead of smoothing during nitridation.
At low N flux, efficient nitridation should only be achieved below about 700 K; above that, the balance of adsorption and immobility still gives dense, rough islands but the process becomes inefficient.
The C-site adsorption energy of $3.17$ eV provides a reference value for further atomistic studies of N/GaAs(100) reactivity.

Reading between the lines

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

If the C-site preference is real, a natural testable extension is to compute the full potential-energy surface for N diffusion with the surface relaxed at each step; the frozen-surface barriers reported here (5.50 and 7.34 eV) may lower substantially, which would soften the predicted immobility.
The model's restriction to three events omits As$_2$ desorption, N desorption, and dimer reformation; including those channels should suppress island density at higher coverage, so the reported high-density and rough regime is best read as an upper bound until those processes are added.
The paper itself flags the 60 eV barrier over the As dimer as a numerical artifact; a properly relaxed calculation of that path would decide whether dimer sites are true traps or merely high-barrier sites, changing the KMC diffusive network.
The asymmetry between the [1 1 0] barrier (5.50 eV) and the [-1 1 0] barrier (7.34 eV) suggests anisotropic diffusion and therefore elongated islands; this prediction is not drawn in the paper but could be checked by Monte Carlo with anisotropic hop rates or by STM.

Original abstract

We present, in this work, our preliminary results of a systematic theoretical study of the adsorption of N over As-terminated GaAs (100) (2$\times$1) surfaces. We analyzed the changes in the bond-lenghts, bond-angles and the energetics involved before and after deposition. Our results show that the N-atoms will prefer the unoccupied sites of the surface, close to the As dimer. The presence of the N pushes the As dimer out of the surface, leading to the anion exchange between the N and As atoms. Based on our results, we discussed about the kinetics of the N islands formation during epitaxial growth of the III-Nitrides.

The compiled record

The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.

Compiled by classifier-1.13-20260917 via OpenRouter, question bank abstract-p1-20260923, 2026-09-23: 1 request, 4,269 input tokens, $0.0002 reported.

What the abstract itself states, as probabilities: claims first 0.03 · states limitation 0.65 · human participants 0.01 · controlled comparison 0.19 · randomized assignment 0.01 · states quantitative result 0.05 · compares with prior results 0.06 · mentions code or data availability 0.02.

The abstract, sentence by sentence

arXiv listing abstract v1. Results and the main claim are underlined; show every role.

We present, in this work, our preliminary results of a systematic theoretical study of the adsorption of N over As-terminated GaAs (100) (2$\times$1) surfaces. We analyzed the changes in the bond-lenghts, bond-angles and the energetics involved before and after deposition. Our results show that the N-atoms will prefer the unoccupied sites of the surface, close to the As dimer. The presence of the N pushes the As dimer out of the surface, leading to the anion exchange between the N and As atoms. Based on our results, we discussed about the kinetics of the N islands formation during epitaxial growth of the III-Nitrides.

background contribution method result interpretation main claim

What the classifier reports for each sentence. The number is the probability it gave the label.
#	sentence	role	the authors' own work	states the main conclusion	strength
s0	We present, in this work, our preliminary results of a systematic theoretical study of the adsorption of N over As-terminated GaAs (100) (2$\times$1) surfaces.	contribution 1.00	0.87	0.04	asserted 0.74
s1	We analyzed the changes in the bond-lenghts, bond-angles and the energetics involved before and after deposition.	method 0.65	0.94	0.04	asserted 0.92
s2	Our results show that the N-atoms will prefer the unoccupied sites of the surface, close to the As dimer.	result 1.00	0.94	0.32	asserted 0.87
s3	The presence of the N pushes the As dimer out of the surface, leading to the anion exchange between the N and As atoms.	result 0.93	0.88	0.39	asserted 1.00
s4	Based on our results, we discussed about the kinetics of the N islands formation during epitaxial growth of the III-Nitrides.	interpretation 0.76	0.78	0.12	unclear 0.59

What kind of paper the classifier reports this is: simulation study 1.00 · evidence from simulation or computation 1.00.

What was examined, and what was not

sources read

arXiv listing abstract v1

not examined

full text
figures and tables
references and citing sentences
relations to other papers

levels run

P abstract

not run

S sections
T body sentences
C citations
R cross-paper relations

Every quoted sentence on this record is the authors'. Every label is from Pith's closed vocabulary, with the probability shown beside it. Every number in the measurements table is the authors', checked by code. Record as JSON

Figures

**Figure 1.** Top view of the possible N adsorption sites over

**Figure 2.** Adsorption energy variation as a function of the

**Figure 4.** Relaxed atomic positions after the N adsorption at

**Figure 5.** Total energy variation as a function of the N

Showing 4 of the paper's 5 figures. See all figures in the PDF ↗ Figures reproduced from arXiv:cond-mat/0603111 by the authors.

References

Works cited (16)

[16]

We did not include the N diffusion over the As dimers in the kinetic Monte Carlo simulations, once that the As dimers act as a N traps

View on Pith
[1]

S. C. Jain et al., J. Appl. Phys. 87, 965 (2000)

View on Pith 2000
[2]

C. C. Silva, H. W. Leite Alves, L. M. R. Scolfaro, and J. R. Leite, phys. stat. sol. (c) 2, 2468 (2005)

View on Pith 2005
[3]

P. Hill et al., J. Vac. Sci. Technol. B15, 113 (1997)

View on Pith 1997
[4]

H. Okumura et al., Appl. Phys. Lett. 72, 3056 (1998)

View on Pith 1998
[5]

D. Schikora et al., Phys. Rev. B54, R8381 (1996); D. J. As, D. Schikora and K. Lischka, phys. stat sol. (c) 0, 1607 (2003)

View on Pith 1996
[6]

F. Widmann et al. Appl. Phys. Lett. 73, 2642 (1998)

View on Pith 1998
[7]

X. Gonze et al., Comput. Mater. Sci. 25, 478 (2002), and references therein

View on Pith 2002
[8]

P. Kratzer, E. Penev and M. Scheffler, Appl. Phys. A 75, 79 (2002), and references therein

View on Pith 2002
[9]

M. Fuchs and M. Scheffler, Comput. Phys. Commun. 119, 67 (1999); C. Stampfl and C. G. Van de Walle, Phys. Rev. B59, 5521 (1999)

View on Pith 1999
[10]

A. de Paiva, J. L. A. Alves and H. W. Leite Alves, Braz. J. Phys. 34, 608 (2004)

View on Pith 2004
[11]

W. G. Schmidt and F. Bechstedt, Phys. Rev. B54, 16742 (1996)

View on Pith 1996
[12]

V. P. LaBella et al ., Int. J. Mod. Phys. B15, 2301 (2001)

View on Pith 2001
[13]

A. P. Castro and H. W. Leite Alves, Microelectron. J. 36, 1045 (2005)

View on Pith 2005
[14]

M. Itoh et al., Phys. Rev. Lett. 81, 633 (1998)

View on Pith 1998
[15]

The obtained energy barrier for the N diffusion over the As dimers was, numerically, 60.0 eV, a very high value. We cannot affirm that this energy barrier has such value, which it does not make any sense. However, we can only infer, from our numerical results, that this barrier is very high so that, if the N is adsorbed over As dimers, it stays immobile o...

View on Pith

Discussion

Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.

[1] [16]

We did not include the N diffusion over the As dimers in the kinetic Monte Carlo simulations, once that the As dimers act as a N traps

View on Pith

[2] [1]

S. C. Jain et al., J. Appl. Phys. 87, 965 (2000)

View on Pith 2000

[3] [2]

C. C. Silva, H. W. Leite Alves, L. M. R. Scolfaro, and J. R. Leite, phys. stat. sol. (c) 2, 2468 (2005)

View on Pith 2005

[4] [3]

P. Hill et al., J. Vac. Sci. Technol. B15, 113 (1997)

View on Pith 1997

[5] [4]

H. Okumura et al., Appl. Phys. Lett. 72, 3056 (1998)

View on Pith 1998

[6] [5]

D. Schikora et al., Phys. Rev. B54, R8381 (1996); D. J. As, D. Schikora and K. Lischka, phys. stat sol. (c) 0, 1607 (2003)

View on Pith 1996

[7] [6]

F. Widmann et al. Appl. Phys. Lett. 73, 2642 (1998)

View on Pith 1998

[8] [7]

X. Gonze et al., Comput. Mater. Sci. 25, 478 (2002), and references therein

View on Pith 2002

[9] [8]

P. Kratzer, E. Penev and M. Scheffler, Appl. Phys. A 75, 79 (2002), and references therein

View on Pith 2002

[10] [9]

M. Fuchs and M. Scheffler, Comput. Phys. Commun. 119, 67 (1999); C. Stampfl and C. G. Van de Walle, Phys. Rev. B59, 5521 (1999)

View on Pith 1999

[11] [10]

A. de Paiva, J. L. A. Alves and H. W. Leite Alves, Braz. J. Phys. 34, 608 (2004)

View on Pith 2004

[12] [11]

W. G. Schmidt and F. Bechstedt, Phys. Rev. B54, 16742 (1996)

View on Pith 1996

[13] [12]

V. P. LaBella et al ., Int. J. Mod. Phys. B15, 2301 (2001)

View on Pith 2001

[14] [13]

A. P. Castro and H. W. Leite Alves, Microelectron. J. 36, 1045 (2005)

View on Pith 2005

[15] [14]

M. Itoh et al., Phys. Rev. Lett. 81, 633 (1998)

View on Pith 1998

[16] [15]

The obtained energy barrier for the N diffusion over the As dimers was, numerically, 60.0 eV, a very high value. We cannot affirm that this energy barrier has such value, which it does not make any sense. However, we can only infer, from our numerical results, that this barrier is very high so that, if the N is adsorbed over As dimers, it stays immobile o...

View on Pith