A Monte Carlo simulation of the physical vapor deposition of nickel

doi:10.1016/S1359-6454(96)00256-X

Acta Materialia

Volume 45, Issue 4, April 1997, Pages 1455-1468

https://doi.org/10.1016/S1359-6454(96)00256-X Get rights and content

Abstract

A two-step Monte Carlo method for atomistically simulating low energy physical vapor deposition processes is developed and used to model the two-dimensional physical vapor deposition of nickel. The method consists of an impact approximation for the initial adatom adsorption on a surface and a multipath diffusion analysis to simulate subsequent surface morphology and interior atomic structure evolution. An embedded atom method is used to determine the activation energies for each of the many available diffusional paths. The method has been used to predict the morphology/structure evolution of nickel films over the length and time scales encountered in practical deposition processes. The modeling approach has enabled determination of the effect of vapor processing variables such as flux orientation, deposition rate and substrate temperature on deposit morphology/microstructure as defined by packing density, surface roughness and growth column width (which appears closely related to grain size). Several aspects of the empirical Movchan-Demchishin structure zone model are well predicted by this approach.

References (67)

L.M. Hsiung et al.
Scripta Metall. et Mater.
(1993)
H.N.G. Wadley et al.
Composites. Engr.
(1995)
B.W. Dodson
Critical Reviews in Solid State and Materials Science
(1990)
X. W. Zhou, R. A. Johnson and H. N. G. Wadley, accepted for publication in Acta Metall. et...
K.H. Müller
Phys. Rev.
(1987)
K.H. Müller
J. Appl. Phys.
(1985)
R.D. Young et al.
J. Chem. Phys.
(1965)
Y. G. Yang, R. A. Johnson and H. N. G. Wadley, unpublished...
X. W. Zhou, R. A. Johnson and H. N. G. Wadley, to be...
A.G. Dirks et al.
Thin Solid Films
(1992)

R. Messier et al.

J. Appl. Phys.

(1985)

M.A. Cotta et al.

Phys. Rev. Lett.

(1993)

R.F. Bunshah

J. Vac. Sci. Technol.

(1974)

J.F. Butler

J. Vac. Sci. Technol.

(1970)

J.A. Thornton

SAE Trans.

(1974)

J. J. Schmitt, U.S. Patent No. 4, 788, 082, Nov. 29...

B.L. Halpern et al.

Metal Finishing

(1992)

B.L. Halpern et al.

J. F. Groves and H. N. G. Wadley, Composites. Engr. to be...

J.F. Groves et al.

J. Electrochem. Soc.

(1995)

E. Bauer

Z. Krist.

(1958)

F.C. Frank et al.F.C. Frank et al.

M. Volmer et al.

Z. Phys. Chem.

(1926)

J.N. Stranski et al.

Ber. Akad. Wiss. Wien

(1938)

B.A. Movchan et al.

Phys. Met. Metallogr.

(1969)

J.A. Thornton

Ann. Rev. Mater. Sci.

(1977)

R.W. Vook et al.

J. Vac. Sci. Technol.

(1965)

S.D. Dahlgren

J. Vac. Sci. Technol.

(1974)

K.L. Chopra

Thin Film Phenomena

(1969)

H.J. Leamy et al.

G.H. Gilmer

J.B. Adams et al.

J. Nuclear Materials

(1994)

M. Schneider et al.

Phys. Rev.

(1987)

Cited by (117)

Recent progress on the mesoscale modeling of architected thin-films via phase-field formulations of physical vapor deposition
2022, Computational Materials Science
Citation Excerpt :
As such, the development of reliable physics-based vapor deposition models, coupled with machine-learning-based surrogate modeling techniques, to predict realistic microstructures based upon material properties and processing parameters would greatly aid in overcoming this limitation and allow one to optimize or develop new vapor-deposition conditions for architected (i.e., designed with a purpose) thin-films. From the theoretical and computational materials perspective, different modeling tools from multiple length-scales have been utilized to investigate the various phenomena relevant to the vapor deposition and thin-film growth processes, namely vapor transport, kinetic surface roughening, surface adhesion and diffusion, bulk chemo-mechanical spinodal decomposition (phase separation), and bulk grain growth [13,16–25]. For example, the formative work of Keblinski et al. [19] laid out a phase-field-based model for a single-phase system capable of describing the dynamics of the depositing vapor phase(s) along with local and non-local growth mechanisms, which were shown to give rise to arbitrary surface growth and topology formation.
Thin-film coatings can be found everywhere in modern technological applications due to desirable electrical, mechanical, chemical, and optical properties. These properties directly depend upon the thin-film’s microstructural features, which are themselves influenced by the materials and vapor-deposition processing conditions used for fabrication. As such, understanding processing-microstructure relationships is essential to designing thin-films with optimized properties, and discovering new processing conditions that allow for novel thin-films with multifunctional microstructures. Here, a short review is presented on recent developments that utilize the phase-field method to simultaneously model the vapor-deposition process and corresponding microstructure formation at the mesoscale. Phase-field-based vapor-deposition models that simulate thin-film growth of immiscible alloy and polycrystalline systems are highlighted in addition to machine-learning-based surrogate models that can facilitate accelerated high-fidelity simulations along with materials design and exploration studies.
Investigating the structures and residual stress of Cu<inf>x</inf>(FeAlCr)<inf>100−x</inf> film on Ni substrate using molecular dynamics
2022, Materials Today Communications
This investigation explores the formation of a Cu_x(FeAlCr)_100−x film on a Ni(001) substrate performed by molecular dynamics simulation. The effects of compositions, incident angle, and temperature of the substrate on the intermixing, surface morphology, structure composition and residual stress are studied. The surface roughness of the film increases as the Cu ratio decreases from 55% to 10%. Besides, the roughness of the film surface achieves the maximum value at 30°. However, the surface roughness of the film decreases at a higher substrate temperature. The intermixing between Cu-Fe-Al-Cr atoms and Ni atoms happens at the interfacial zone. The intermixing rate is affected by the composition and incident angle in the coating process. The intermixing rate rises when rising the temperature of the substrate. Interestingly, during the coating process, the crystal structure of the high entropy alloys film is formed. The crystalline structure ratio of the Cu_x(FeAlCr)_100−x film slightly increases when the composition deviation decreases. In comparison, it reduces with the higher temperature of the Ni sample. The average biaxial stress and average normal stress of the Cu₂₅Fe₂₅Al₂₅Cr₂₅ film only have minor changes when increasing the incident angle. However, the average biaxial stress of the Cu₂₅Fe₂₅Al₂₅Cr₂₅ film increases as raising the temperature.
Effect of annealing and deposition of Cu atoms on Ni trench to interface formation and growth mechanisms of Cu coating
2020, Superlattices and Microstructures
Molecular dynamics (MD) simulation is used to simulate the deposition and annealing process of Cu atoms on the Ni (001) trench. The results indicate that the curvature of profile and surface roughness of the Cu layer film decrease as increasing the substrate temperature and/or incident energy. The surface roughness of the deposition Cu film is variability as increasing the incident angle from 0° to 30°. In the growth process, the intermixing occurs at the interface between Cu atoms and Ni trench, and the length increases as the increasing of substrate temperature and/or incident energy. These phenomena may be attributed to the enhanced atomic mobility for higher incident energy and thermal annealing. The peak-to-peak distances of g(r) clearly indicate that the crystal structure of Cu–Ni atoms after the annealing process is improved, and still FCC structure. The dislocation is formed during the deposition process and affected by the substrate temperature, as well as the incident energy. The annealing treatment reduces the residual stress in the deposited Cu layer. The formation of the vacancies appears during the deposition process. After the annealing process, the Cu–Ni atoms become more orderly in the crystal lattice, many vacancies and voids are destroyed. In addition, the surface of the Cu film layer becomes smoother.
The effect of incident energy, incident angle and substrate temperature on surface morphology and atomic distribution of NiTi films
2020, Materials and Design
Molecular dynamics simulation is employed to study the formation of NiTi thin films on Ni (001) substrates simultaneously deposited with Ni and Ti atoms in this work. The effects of different deposition parameters on surface morphology and atomic distribution of deposited NiTi films are investigated. The deposition parameters mainly include incident energy, incident angle and substrate temperature. The simulation results indicate that incident energy and substrate temperature have a significant effect on the surface roughness of deposited NiTi film. However, incident angles below 30° have little influence on surface roughness of NiTi film. The simulation results show that increasing substrate temperature and incident energy can enhance the intermixing between substrate and film. The content of Ni on the NiTi film surface increase significantly with the raise of incident energy. The defect analysis of the deposited film demonstrates that incident angle mainly affects the number of voids and vacancies in the film. For film samples of different deposition conditions, the films deposited at larger incident angles contain more voids and vacancies. While the films deposited at higher substrate temperatures contain fewer voids and vacancies.
Monte Carlo simulation of column growth in plasma spray physical vapor deposition process
2018, Surface and Coatings Technology
Citation Excerpt :
Monte Carlo methods start with random inputs from a probability distribution then perform a deterministic computation on the inputs to obtain numerical results. An overview on the Monte Carlo model for simulating microstructural evolution can be found in [13,19–25]. We use a direct Monte Carlo method for vapor-deposited films to simulate the impingement on a substrate in the PS-PVD process.
Plasma spray-physical vapor deposition is used to produce columnar microstructure coatings under particular operating parameters. Simulations of the growth of columns were carried out through a two-dimensional Monte Carlo model. The modeling was performed using inclined vapor flux impinging onto a substrate due to shadowing effects. An incoming particle travels along a straight line and attaches itself to already deposited particles. Furthermore, the newly deposited particle will relax to a stable surrounding position along the incoming velocity direction. The modeling results predicted the linking of an oblique vapor flux and column orientation. The numerical simulations were validated in three ways. Firstly, the porosity of simulated columns was predicted and compared to that obtained in the experimental columnar microstructure of coatings. Secondly, the morphology of simulated structures is compared to that of experimental coatings produced by plasma spray physical vapor deposition. Finally, the simulated orientation of columns is compared to the experimental one.
Physical vapor deposition of multiphase materials with phase nucleation via a coupled phase-field approach
2018, Computational Materials Science
A coupled phase-field model is presented for simulating physical vapor deposition (PVD) of multi-phase materials, including the effects of phase nucleation. This model is utilized to study the role of initial substrate phase and temperature distributions, which are important experimental deposition parameters, on PVD of a generic allotropic metal with two stable phases. PVD simulations are performed for the deposition of a high temperature phase below its phase transition temperature, and for bicrystal and amorphous substrates with a Gaussian temperature distribution. Several general observations are made from these simulations. During the initial stages of PVD, the substrate phase distribution acts as a template for the growing phases. As the thin film continues to grow, the release of latent heat due to deposition creates a temperature gradient within the film, i.e., regions near the film surface become hotter than near the substrate. Additionally, a substrate with a defined temperature distribution that encompasses temperatures above and below the phase transition temperature, allows for distinguishable regions within the thin film where different phases are stable and continue to grow. Three-dimensional simulations provide additional insight into the role of substrate temperature and phase distribution on the resulting microstructure with different flux rates.

View all citing articles on Scopus

View full text

A Monte Carlo simulation of the physical vapor deposition of nickel

Abstract

Scripta Metall. et Mater.

Composites. Engr.

Critical Reviews in Solid State and Materials Science

Phys. Rev.

J. Appl. Phys.

J. Chem. Phys.

Thin Solid Films

J. Appl. Phys.

Phys. Rev. Lett.

J. Vac. Sci. Technol.

J. Vac. Sci. Technol.

SAE Trans.

Metal Finishing

J. Electrochem. Soc.

Z. Krist.

Z. Phys. Chem.

Ber. Akad. Wiss. Wien

Phys. Met. Metallogr.

Ann. Rev. Mater. Sci.

J. Vac. Sci. Technol.

J. Vac. Sci. Technol.

Thin Film Phenomena

J. Nuclear Materials

Phys. Rev.