Microstructure and surface morphology of YSZ thin films deposited by e-beam technique

In present study yttrium-stabilized zirconia (YSZ) thin films were deposited on optical quartz (amorphous SiO₂), porous Ni-YSZ and crystalline Alloy 600 (Fe-Ni-Cr) substrates using e-beam deposition technique and controlling technological parameters: substrate temperature and electron gun power which influence thin-film deposition mechanism. X-ray diffraction, scanning electron microscopy (SEM), and atomic force microscopy (AFM) were used to investigate how thin-film structure and surface morphology depend on these parameters. It was found that the crystallite size, roughness and growth mechanism of YSZ thin films are influenced by electron gun power. To clarify the experimental results, YSZ thin-film formation as well evolution of surface roughness at its initial growing stages were analyzed. The evolution of surface roughness could be explained by the processes of surface mobility of adatoms and coalescence of islands. The analysis of these experimental results explain that surface roughness dependence on substrate temperature and electron gun power non-monotonous which could result from diffusivity of adatoms and the amount of atomic clusters in the gas stream of evaporated material.     

动力学晶格蒙特卡洛方法模拟Cu薄膜生长

利用动力学晶格蒙特卡洛方法模拟了Cu薄膜在Cu(100)面上的三维生长过程。模型中考虑了四个动力学过程:原子沉积、增原子迁移、双原子迁移和台阶边缘原子迁移,各动力学过程发生的概率由多体势函数确定。讨论了基底温度、沉积速率及原子覆盖率对Cu原子迁移、成核和表面岛生长等微观生长机制的影响;获得了Cu薄膜的表面形貌图并计算了表面粗糙度。模拟结果表明,随基底温度升高或沉积速率下降,岛的平均尺寸增大,数目减少,形状更加规则。低温时,Cu薄膜表现为分形的离散生长,高温时,Cu原子迁移能力增强形成密集的岛。Cu薄膜表面粗糙度随着基底温度的升高而迅速减小;当基底温度低于某一临界温度时,表面粗糙度随原子覆盖率或沉积速率的增大而增大;当基底温度超过临界温度时,表面粗糙度随原子覆盖率或沉积速率的变化很小,基本趋于稳定。     

Optimal conditions for submonolayer nanoisland growth in ion beam assisted deposition

We study the optimal conditions for nanoisland growth in ion beam assisted deposition (IBAD). This situation occurs when adatom islands remain small enough to prevent the onset of three-dimensional growth, while at the same time preventing ion-induced surface erosion. To this end, we develop a rate equation model of IBAD, which embodies continuous deposition of adatoms and creation of vacancies, recombination of vacancies at adatom island edges, as well as recombination of adatoms at vacancy island edges. These rate equations are solved by numerical simulations based on the particle coalescence method. To determine the optimal growth condition, we find the largest mean size of the vacancy islands leading to their survival. We show that at this onset between the rough and smooth layer-by-layer growth regimes there is a simple exponential relation between the largest size of the vacancy islands and the external control parameters of the growth.     

Growth mechanism and morphology of Ge on Pb covered Si(111) surfaces

We study the mechanism and surface morphology in epitaxial growth of Ge on Pb covered Si(111) using a scanning tunneling microscope (STM). We find that Ge adatoms have a very large diffusion length at room temperature. The growth is close to perfect layer-by-layer for the first two bilayers. Surface roughness increases gradually with the film thickness, but no 3D islands are found at room temperature. For growth at 200°C, 3D Ge islands appear after completion of the second bilayer. At room temperature, we believe, the Pb layer enhances surface diffusion and the descending-step motion of Ge adatoms, but the ascending-step motion is hindered and thus 3D island growth is suppressed.     

非均匀基底上三维薄膜生长的模拟研究

考虑原子在基底表面的扩散、沿岛周界的扩散和不同层间的扩散以及非均匀基底上表面吸附能分布的各向异性，建立起非均匀基底表面上原子扩散和三维薄膜生长的动力学蒙特卡罗模型.模拟得到在不同生长条件下出现的层状生长、岛状生长和混合生长三种生长模式和相应的多层薄膜生长形貌图.通过统计三维薄膜中原子在各层的分布，计算薄膜的表面粗糙度，得到薄膜生长模式与生长条件之间的关系. 关键词： 薄膜生长 非均匀基底 动力学蒙特卡罗模拟     

沉积粒子能量对薄膜早期生长过程的影响

利用Monte Carlo(MC)模型研究了能量粒子对薄膜生长的初始阶段岛膜的形貌和岛的尺寸的影响,沉积粒子的能量范围为:0—0.7eV.在模型中考虑了原子沉积、吸附原子扩散和蒸发等过程,并详细考虑了临近和次临近原子的影响.结果表明,在所采用的参量范围内不同的基底温度情况下,能量粒子的影响有很大的区别.低基底温度情况下,沉积粒子强烈地影响着薄膜的生长过程中,岛膜的形貌、数量和尺寸随能量粒子的能量增加而有很大的变化.分析表明,这些变化都是由于能量粒子的介入使得表面吸附粒子的扩散能力增强所致 关键词： 薄膜生长 Monte Carlo方法 扩散     

薄膜外延生长的计算机模拟

以Cu膜为例，用Monte-Carlo算法模拟了薄膜生长的随机过程，并提出了更加完善的模型.在合理选择原子间相互作用计算方法的基础上，考虑了原子的吸附、在生长表面的迁移及迁移所引起的近邻原子连带效应、从生长表面的脱附等过程.模拟计算了薄膜的早期成核情况以及表面粗糙度和相对密度.结果表明，随着衬底温度的升高或入射率的降低，沉积在衬底上的原子逐步由离散型分布向聚集状态过渡形成一些岛核，并且逐步由二维岛核向三维岛核过渡.在一定的原子入射率下，存在三个优化温度，成核率最高时的最大成核温度Tn、薄膜的表面粗糙度最低 关键词： Monte-Carlo算法 计算机模拟 薄膜生长     

EFFECT OF SMALL CLUSTER DIFFUSION DURING TWO-DIMENSIONAL THIN FILM GROWTH ON METAL SURFACE

The diffusion of small clusters such as dimers and trimers on metal surface and the growth of two-dimensional thin films are studied by Monte Carlo simulation, using realistic growth model and physical parameters. It is found that small cluster diffusion plays an important role in the process of thin film growth at not very low temperature. It affects not only the island density and the size of islands but also the critical value of saturation occurring during growth of thin films. The effect of small cluster diffusion depends on both the size of critical nucleus and the growth temperature. The simulation results also show that the larger the cluster allowed to diffuse, the easier the saturation that takes place, giving rise to the lower critical coverage of saturation occurring. It is suggested that the effect of small cluster diffusion should be included in establishing the growth models of thin films.     

薄膜生长中的表面动力学(Ⅱ)

本文较全面地从理论上研究了薄膜生长过程中原子在表面上的各种动力学表现，涉及的内容包括亚单层生长时，原子在表面上的扩散，粘接，成核，以及已经形成的原子岛之间的相互作用，兼并，失稳，退化等一系列过程。在第一部分(即0—6章，刊登在《物理学进展》23卷1期上)介绍了薄膜生长动力学的基础之后，从第7章开始我们侧重研究一个在向同性和各向异性表面都普遍成立的原子岛尺寸大小及密度分布的标度定律；建立了一套研究在各向同性和各向异性表面上二维原子岛退化过程的广义动力学标度理论。基于对层间质量扩散通道的研究，本文提出了两种不同的原子下跃机制，既任意跃迁机制和选择跃迁机制，并做了进一步的实验验证。利用这一新的层间质量扩散机制，我们成功地解释了实验上观测到的三维原子岛的退化规律。更加有趣的是，本文讨论了应力对岛的形状和各种动力学规律的影响。在最后我们还提出了一个利用凝聚能转化来控制二维原子岛生长的方法，其目的是希望能够找到一种人为有效地在表面上制备低维量子结构的方法。     

薄膜生长中的表面动力学(Ⅱ)

本文较全面地从理论上研究了薄膜生长过程中原子在表面上的各种动力学表现 ,涉及的内容包括亚单层生长时 ,原子在表面上的扩散 ,粘接 ,成核 ,以及已经形成的原子岛之间的相互作用 ,兼并 ,失稳 ,退化等一系列过程。在第一部分 (即 0～ 6章 ,刊登在《物理学进展》2 3卷 1期上 )介绍了薄膜生长动力学的基础之后 ,从第 7章开始我们侧重研究一个在向同性和各向异性表面都普遍成立的原子岛尺寸大小及密度分布的标度定律 ;建立了一套研究在各向同性和各向异性表面上二维原子岛退化过程的广义动力学标度理论。基于对层间质量扩散通道的研究 ,本文提出了两种不同的原子下跃机制 ,既任意跃迁机制和选择跃迁机制 ,并做了进一步的实验验证。利用这一新的层间质量扩散机制 ,我们成功地解释了实验上观测到的三维原子岛的退化规律。更加有趣的是 ,本文讨论了应力对岛的形状和各种动力学规律的影响。在最后我们还提出了一个利用凝聚能转化来控制二维原子岛生长的方法 ,其目的是希望能够找到一种人为有效地在表面上制备低维量子结构的方法。     

二氧化锆薄膜表面粗糙度的研究

采用电子束蒸发工艺,利用泰勒霍普森相关相干表面轮廓粗糙度仪,研究了不同基底粗糙度、不同二氧化锆薄膜厚度以及不同的离子束辅助能量下所沉积的二氧化锆薄膜的表面粗糙度。结果表明：随着基底表面粗糙度的增加,二氧化锆薄膜表面粗糙度呈现出先缓慢增加,当基底的粗糙度大于10nm后呈现快速增加的趋势;随着二氧化锆薄膜厚度的增加,其表面均方根粗糙度（RMS）先减小后增大;随着辅助沉积离子能量的增加,其表面粗糙度呈现出先减小后增加的趋势。     

Multilayer growth of BaTiO3 thin films via pulsed laser deposition: An energy-dependent kinetic Monte Carlo simulation

An energy-dependent kinetic Monte Carlo approach was proposed to simulate the multilayer growth of BaTiO₃ thin films via pulsed laser deposition, in which the four steps, such as the deposition of atoms, the diffusion of adatoms, the bonding of adatoms, and the surface migration of adatoms, were considered. Distinguishing with the traditional solid-on-solid (SOS) model, the adatom bonding and the overhanging of atoms, according to the perovskite structure, were specially adopted to describe the ferroelectric thin film growth. The activation energy was considered from the interactions between the ions, which were calculated by Born-Mayer-Huggins (BMH) potential. From the simulation the relative curves of the each layer coverage and roughness vs total coverage were obtained by varying the parameter values of the incident kinetic energy, laser repetition rate and mean deposition rate. The relationship between growth modes and the different parameters was also acquired.     

Studies on morphology and molecular arrangement of pentacene on different substrates

Pentacene (C₂₂H₁₄) thin films with different thicknesses were fabricated to study the dynamic growth process and morphology of pentacene on different substrates. A discontinuous monomolecular layer was observed when a pentacene thin film is about 0.5 nm thick on native oxide silicon wafer. The terraced islands and dendritic structure gradually formed with increasing pentacene thin film thickness. The height of each layer is about 1.4 nm which corresponds well with the length of the long axis of the pentacene molecule at 1.45 nm. Experimental results show that the pentacene molecule is perpendicular to the silicon wafer surface with a slight tilted angle. However, the pentacene molecular orientation on a polymer pre-covered indium tin oxide coated substrate could not give any indication on the scale of nanometers. The surface roughness of substrates strongly influences pentacene molecular diffusion and the morphology of pentacene thin films.     

Investigation of interface roughness cross-correlation properties of optical thin films from total scattering losses

The interface roughness and interface roughness cross-correlation properties affect the scattering losses of high-quality optical thin films. In this paper, the theoretical models of light scattering induced by surface and interface roughness of optical thin films are concisely presented. Furthermore, influence of interface roughness cross-correlation properties to light scattering is analyzed by total scattering losses. Moreover, single-layer TiO₂ thin film thickness, substrate roughness of K9 glass and ion beam assisted deposition (IBAD) technique effect on interface roughness cross-correlation properties are studied by experiments, respectively. A 17-layer dielectric quarter-wave high reflection multilayer is analyzed by total scattering losses. The results show that the interface roughness cross-correlation properties depend on TiO₂ thin film thickness, substrate roughness and deposition technique. The interface roughness cross-correlation properties decrease with the increase of film thickness or the decrease of substrates roughness. Furthermore, ion beam assisted deposition technique can increase the interface roughness cross-correlation properties of optical thin films. The measured total scattering losses of 17-layer dielectric quarter-wave high reflection multilayer deposited with IBAD indicate that completely correlated interface model can be observed, when substrate roughness is about 2.84 nm.     

沉积工艺对二氧化锆薄膜生长特性影响的研究

利用反应离子束溅射、反应磁控溅射和电子束蒸发在K9基底上沉积ZrO2薄膜，并用原子力显微镜对薄膜表面形貌进行测量。通过数值相关运算，对不同工艺条件下薄膜生长界面进行定量描述，得到了薄膜表面的粗糙度指数、横向相关长度、标准偏差粗糙度等参量。由于沉积条件的不同，薄膜生长具有不同的动力学过程。在反应离子束溅射和反应磁控溅射沉积薄膜过程中，薄膜生长动力学行为均可用Kuramoto-Sivashinsky方程来描述，电子束蒸发制备薄膜的过程可以用Mullins扩散模型来描述，并发现在沉积薄膜过程中基底温度和沉积过程的稳定性对薄膜表面特征影响很大。     

Effect of Incident Intensity on Films Growth in Pulsed Laser Deposition

Incident intensity, defined by the amount of particles deposited per pulse, is an important parameter in the film growth process of pulsed laser deposition （PLD）. Different from previous models, we investigate the irreversible and reversible growth processes by using a kinetic Monte Carlo method and find that island density and film morphology strongly depend on pulse intensity. At higher pulse intensities, lots of adatoms instantaneously diffuse on the substrate surface, and then nucleation easily occurs between the moving adatoms resulting in more smaller-size islands. In contrast, at the lower pulse intensities, nucleation event occurs preferentially between the single adatom and existing islands rather than forming new islands, and therefore the average island size becomes larger in this case. Additionally, our results show that substrate temperature plays an important role in film growth. In particular, it can determine the films shape and weaken the effect of pulse intensity on film growth at the lower temperatures by controlling the mobility rate of atoms. Our results can match the related theoretical and experimental results.     

Atomic-level investigation of Al and Ni thin film growth on Ni(1 1 1) surface: Molecular dynamics simulation

Using molecular dynamics (MD) simulation, the structural characteristics of Al and Ni thin film growth on Ni(1 1 1) substrate according to the incident energy of adatoms were investigated. In case of Al on Ni(1 1 1), Al adatoms were grown basically through the layer-by-layer growth mode. On the other hand, Ni thin films on Ni(1 1 1) surface at low incident energy were shown to favor island growth. The steering effect due to atomic attraction, which results in rougher surface, was significantly observed at low incident energy. The growth mode of Ni film was, however, changed to follow layer-by-layer growth mode for the incident energy of 6 eV. The different aspects of surface morphology between Al and Ni deposition on Ni(1 1 1) surface could be successfully explained by the surface diffusion and impact cascade diffusion.     

Brownian motion of 2D vacancy islands by adatom terrace diffusion.

We have studied the Brownian motion of two-dimensional (2D) vacancy islands on Ag(110) at temperatures between 175 and 215 K. While the detachment of adatoms from the island and their diffusion on the terrace are permitted in this temperature range, the periphery diffusion of single adatoms is prohibited. The present scanning tunneling microscopy results provide the first direct experimental proof that the Brownian motion of the islands follows a simple scaling law with terrace diffusion being the rate limiting process. The activation energy of the vacancy island motion is determined to 0.41 eV.     

Electron and atom dynamics at solid surfaces and relation to epitaxy

At the surface, the three dimensional symmetry of a solid is broken. Electrons and atoms near the surface may rearrange to lower the free energy of the system. Scattering by defects and confinement by boundaries of electrons may produce long-range charge density oscillations. Adatoms interact with each other via mutual perturbation of the surface, known as indirect electronic and elastic interactions. These interactions are very weak and are also oscillatory. For some systems, formation of adsorption layer superstructures can be directly correlated to adatom–adatom interactions. When the temperature is raised, adatoms and admolecules can start to diffuse, interact, or react. They may aggregate into clusters and islands, and grown into a thin film. The stability of clusters may exhibit magic numbers in size and thickness. When the temperature is changed, island shape transitions may occur. The growth of islands and ultra-thin films can also be influenced by electronic effects as well as by the addition of a surfactant layer. All these growth behaviors in epitaxy can be understood from the mechanisms and energetics of elementary surface atomic processes, and atom and electron dynamics. They, in turn, can be studied in details using atomic resolution microscopy.     

Neural network approach to diffusion of B and N adatoms on the Pt(111) surface

Artificial neural network boasts an outstanding class in terms of its performance and efficiency and is being spread widely in most fields of physics. We report our investigation of the diffusions of boron and nitrogen adatoms on the Pt(111) surface by performing molecular dynamics simulations equipped with machine-learned potentials. Platinum is commonly used as a substrate for the growth of hexagonal boron nitride (h-BN) thin films, and the diffusion of B and N atoms on the substrate, which are decomposed from the precursor molecules, plays important roles in the initial stages of h-BN growth. The two-dimensional potential energy surfaces and the trajectories of the B and N adatoms are consistent with the DFT calculation. The Arrhenius plots of the diffusion coefficients produce the diffusion barriers of the B and N adatoms on Pt(111), which agree well with the DFT barriers.