Abstract
Gas-surface reactions of the type that contribute to growth during the chemical vapor deposition (CVD) of diamond films are generally completed in picoseconds, well within timescales accessible by molecular dynamics (MD) simulations. For low-pressure deposition, however, the time between collisions for a surface site can be microseconds, which makes direct modeling of CVD crystal growth impossible using standard MD methods. To effectively bridge this discrepancy in timescales, the gas-surface reactions can be modeled using MD trajectories, and then this data can be used to define probabilities in a Monte Carlo algorithm where each step represents a gas-surface collision. We illustrate this approach using the reaction of atomic hydrogen with a diamond (111) surface as an example, where we use abstraction and sticking probabilities generated using classical trajectories in a simple Monte Carlo algorithm to determine the number of open sites as a function of temperature. We also include models for the thermal desorption of hydrogen that predict that growth temperatures are not restricted by the thermal loss of chemisorbed hydrogen.
Similar content being viewed by others
References
J.C. Angus and C.C. Hayman, Science 241, 913 (1988); F.G. Celii and J.E. Butler, Annu. Rev. Phys. Chem. 42, 643 (1991).
A.V. Hamza, G.D. Kubiak and R.H. Stulen, Surf. Science Lett. 206, L833 (1988).
B.B. Pate, Surf. Science 165, 83 (1986).
K.A. Snail and L.M. Hanssen, J. Cryst. Growth 112, 651 (1991).
H.J.C. Berendsen, J.P.M. Postma, W.F. van Gunsteren, A. DiNola and J.R. Haak, J. Chem. Phys. 81, 3684 (1984).
C.W. Gear, Numerical Initial Value Problems in Ordinaryj Differential Equations, (Prentice-Hall, Englewood Cliffs, 1971).
D.W. Brenner, Phys. Rev. B 42, 9458 (1990).
C.K. Westbrook, J. Warnatz, W.J. Pitz Twenty Second Symposium (International) on Combustion, The Combustion Institute, 1988, p.893.
M. Page and D.W. Brenner, J. Am. Chem. Soc. 113, 3271 (1991).
M. Frenklach, Phys. Rev. B, in press.
G.C. Canavos, Applied Probability and Statistical Methods, (Little, Brown and Company, Boston, 1984) p. 84.
M. Schulberg (private communication). We also tried other related rates and secondorder desorption with little change in the results.
Y. Sato, Proceedings of the First International Conference on New Diamond Science and Technology, p. 4.
Author information
Authors and Affiliations
Rights and permissions
About this article
Cite this article
Brenner, D., Robertson, D., Carty, R. et al. Combining Molecular Dynamics and Monte Carlo Simulations to Model Chemical Vapor Deposition: Application to Diamond. MRS Online Proceedings Library 278, 255–260 (1992). https://doi.org/10.1557/PROC-278-255
Published:
Issue Date:
DOI: https://doi.org/10.1557/PROC-278-255