Molecular-dynamics simulations using embedded atom method potentials were carried out to study the growth and subsequent annealing of a 6.3 monolayer thick Cu film on a bcc-Ta ͑100͒ substrate of a very large area ͑100ϫ 100 nm 2 ͒. The purpose is to obtain, for this typical example of a crystallographically incompatible system, atomic-level insight into the microstructural evolution of the film/substrate interface and the film itself, with the limiting influence of in-plane boundary conditions

more »

... pt to a minimum. It is found that the first Cu plane grows heteroepitaxially on the bcc-Ta ͑100͒ surface. The second Cu plane is the most interesting plane; it grows in the form of a pattern of 26-atom misfit supercells which after prolonged film deposition relax energetically, by breaking up in groups and forming 30 Å wide in-plane island strips separated by fcc ͑1 / 6͒͗211͘ vectors. The distorted monolayer represents a way of enabling epitaxy between very different crystal structures without introducing misfit dislocations. The third and higher Cu planes are fcc ͑111͒ planes. The Cu film is polycrystalline, starting from the second plane up, with two different in-plane crystal orientations. The influence of a one-monolayer Ta terrace on the substrate is minimal. Upon annealing the film at elevated temperatures, a tendency toward island coarsening and agglomeration is observed, but a full threedimensional island formation such as found experimentally is not observed. Compared to an earlier simulation on a smaller substrate area ͑18ϫ 18 nm 2 ͒, the current simulations have produced a much richer microstructure. The present results therefore warn against the use of too small simulation domains when complex strain fields can be expected.