Stability of multi-component epilayers and nanopattern formation

A uniform multi-component epilayer may lose stability under the combined action of spinodal decomposition and epilayer–substrate interaction, separating into multiple phases. The phases may further self-organize into regular patterns. This paper investigates the compositional stability of a ternary epliayer and the subsequent emergence of nanoscale patterns. Multiple energetic forces and kinetic processes involving phase separation, phase coarsening and phase refining are incorporated into a continuous phase field model. Linear stability analysis is performed by perturbing a uniform concentration field into a sinusoidal field with small amplitude and arbitrary wavelength. The analysis shows that the epilayer–substrate interaction counteracts the coarsening effect of phase boundary energy and may lead to the formation of steady nanoscale patterns. Detailed analysis also reveals the interaction of multi-phases and its effect on the stability condition. Numerical simulation of evolving concentration field is discussed at the end of the paper. The simulations show that the pattern formation process of multi-component epilayers involves remarkably rich dynamics.     

Ripening of surface phases coupled with oscillatory dynamics and self-induced spatial chaos through surface roughening

Some pattern formation processes on single-crystal catalytic surfaces involve transitions between alternative surface phases coupled with oscillatory reaction dynamics. We describe a two-tier symmetry-breaking model of this process, based on nanoscale boundary dynamics interacting with oscillations of adsorbate coverage on microscale. The surface phase distribution oscillates together with adsorbate coverage, and, in addition, undergoes a slow coarsening process due to the curvature dependence of the drift velocity of interphase boundaries. The coarsening is studied both statistically, assuming a circular shape of islands of the minority phase, and through detailed Lagrangian modeling of boundary dynamics. Direct simulation of boundary dynamics allows us to take into account processes of surface reconstruction, leading to self-induced surface roughening. As a result, the surface becomes inhomogeneous, and the coarsening process is arrested way before the thermodynamic limit is reached, leaving a chaotic distribution of surface phases. (c) 1999 American Institute of Physics.     

纳米尺度下Sm-Co合金体系中相组成与相稳定性的研究

建立了纳米晶合金相的热力学模型,可定量描述纳米尺度下合金体系中化合物相的热力学性质,并预测合金相的稳定性及其转变规律.利用该模型全面计算了纳米晶Sm-Co合金体系中各化合物相在不同晶粒尺寸下的摩尔吉布斯自由能随温度的变化关系,预测了纳米尺度下Sm-Co合金体系中各物相的相对稳定性及转变规律.模型预测结果示出,在室温附近,随着纳米晶粒尺寸的减小,某些纳米晶合金相的摩尔吉布斯自由能将由负值变为正值,预示着将向其他更稳定的纳米晶合金相转变,这是与传统粗晶材料中合金相的稳定性仅依赖于温度条件而完全不同的纳米晶合金 关键词： 纳米晶材料热力学 Sm-Co合金 相稳定性 相变     

Simulation of fluid-solid coexistence via thermodynamic integration using a modified cell model

Despite recent advances, precise simulation of fluid-solid transitions still remains a challenging task. Thermodynamic integration techniques are the simplest methods to study fluid-solid coexistence. These methods are based on the calculation of the free energies of the fluid and the solid phases, starting from a state of known free energy which is usually an ideal-gas state. Despite their simplicity, the main drawback of thermodynamic integration techniques is the large number of states that must be simulated. In the present work, a thermodynamic integration technique, which reduces the number of simulated states, is proposed and tested on a system of particles interacting via an inverse twelfth-power potential energy function. The simulations are implemented at constant pressure and the solid phase is modeled according to the constrained cell model of Hoover and Ree. The fluid and the solid phases are linked together by performing constant-pressure simulations of a modified cell model. The modified cell model, which was originally proposed by Hoover and Ree, facilitates transitions between the fluid and the solid phase by tuning a homogeneous external field. This model is simulated on a constant-pressure path for a series of progressively increasing values of the field, thus allowing for direct determination of the free energy difference between the fluid and the solid phase via histogram reweighting. The size-dependent results are analyzed using histogram reweighting and finite-size scaling techniques. The scaling analysis is based on studying the size-dependent behavior of the second- and higher-order derivatives of the free energy as well as the dimensionless moment ratios of the order parameter. The results clearly demonstrate the importance of accounting for size effects in simulation studies of fluid-solid transitions.     

Study of the solid-liquid-vapour phase equilibria of flexible chain molecules using Wertheim's thermodynamic perturbation theory

The phase diagram of flexible molecules formed by freely-jointed tangent spheres is studied using the first-order thermodynamic perturbation theory of Wertheim for both fluid and solid phases. A mean-field term is added to the free energy of the fluid and solid phase in order to account for attractive dispersion forces. The approach is used to determine the global (solid-liquid-vapour) phase diagrams and triple points of chain molecules of increasing chain length. It is found that the triple point temperature is not affected strongly by the length of the chain, whereas the gas-liquid critical temperature increases dramatically. The asymptotic limits of the phase diagram for infinitely long chains are discussed. The reduced critical temperature of infinitely long chains as given by the mean-field theory is 2/3, and the reduced triple point temperature is 0.048 56, so that an asymptotic value of T _t/T _c = 0.07284 for the ratio of the triple to critical point temperatures is obtained. This indicates that fully-flexible tangent chains present an enormous liquid range. The proposed theory, while being extremely simple, provides a useful insight into the phase behaviour of chain molecules, showing the existence of finite asymptotic limits for the triple and critical point temperatures. However, since n-alkanes present an asymptotic limit of about T _t/T _c, = 0.40, the agreement With experiment is not quantitative. This suggests that fully flexible models may not be appropriate to model the solid phases of real chain molecules.     

界面润湿性及固相体积分数对颗粒粗化动力学影响的相场法研究

利用多相场模型模拟了液-固两相体系中固相颗粒的粗化过程, 分析了界面润湿性及固相体积分数对粗化指数、粗化速率及颗粒尺寸分布的影响.结果表明, 不同固相体积分数下粗化指数基本不变, 但粗化速率常数及尺寸分布与固相体积分数及界面润湿性密切相关.在完全润湿条件下, 随着固相体积分数的增加, 粗化速率常数逐渐增大; 而非完全润湿条件下, 随着固相体积分数的增加, 粗化速率常数增大速度变缓, 且当润湿性较低、 固相分数较大时, 粗化速率常数还将随体积分数的增加而下降. 此外, 模拟结果表明各种润湿条件下颗粒的尺寸分布均随着固相分数增加而变宽, 分布峰值降低, 但非完全润湿条件下峰值下降变缓.模拟结果为理解不同实验观测结果之间的分歧提供了依据. 关键词： 粗化 相转变 相场法 润湿性     

Fluid-fluid phase separation in colloid-polymer mixtures studied with small angle light scattering and light microscopy

Mixtures of colloidal silica spheres and polydimethylsiloxane in cyclohexane with a colloid-polymer size ratio of about one were found to phase separate into two fluid phases, one which is colloid-rich and one which is colloid-poor. In this work the phase separation kinetics of this fluid-fluid phase separation is studied for different compositions of the colloid-polymer mixtures, and at several degrees of supersaturation, with small angle light scattering and with light microscopy. The small angle light scattering curve exhibits a peak that grows in intensity and that shifts to smaller wave vector with time. The characteristic length scale that is obtained from the scattering peak is of the order of a few μm, in agreement with observations by light microscopy. The domain size increases with time as , which might be an indication of coarsening by diffusion and coalescence, like in the case of binary liquid mixtures and polymer blends. For sufficiently low degrees of supersaturation the angular scattering intensity curves satisfy dynamical scaling behavior.     

From molecular flexibility to shape, curvature and thermotropic transitions in pure surfaces

A material surface of pure constituents with a flexible molecular chain (amphiphilics) is considered; thermodynamic behaviour is studied in the chain length-temperature plane. The Hamiltonian of the system is modelled as the sum of a formation term which refers to the polymer nature of the chain, and of a fluctuation term with a specific elastic form. For closed systems the model exhibits phases with uniform curvature and conformational order/disorder or, alternatively, modulated phases; a critical chain length is found for the existence of modulated phases; the dependence of transition temperature on energy parameters is determined. A critical region is found for open systems, where conformational disorder drives spontaneous generation of curvature; this lies above a characteristic chain length and around the shape transition temperature. Received: 13 November 1996 / Revised: 9 May 1997 / Received in final form: 4 November 1997 / Accepted: 10 November 1997     

Modeling the thermodynamic properties of bimetallic nanosolids

A model has been developed to account for size, shape, surface segregation, composition and dimension dependent cohesive energy of bimetallic nanosolids, and further been extended to predict the size dependent thermodynamic properties, such as melting temperature, Curie temperatures, ordering temperature and phase diagram. The cohesive energy, melting temperature, Curie temperatures and ordering temperature of bimetallic nanosolids decrease with decreasing the particle size. The depression is dramatic in the lower range of size, while it becomes smoothly in large size. For nano phase diagram, the solidus and liquidus curves drop and the two-phase zones become small, as the size of the nanosolids decreases. The two-phase zones of the nano phase are always lower than the regions indicated in the bulk Ag-Pd alloy phase diagram, and they may deteriorate into a curve at a critical size. It is also found that the thermodynamic properties of nanosolids not only depend on the compositions, the atomic diameter and the cohesive energy of each component, but also depend on the size and the shape. The model predictions are consistent with the corresponding simulation, semi-empirical model and experimental data.     

Free energy studies of freezing in slit pores: an order-parameter approach using Monte Carlo simulation

We report a molecular simulation study of freezing transitions for simple fluids in narrow slit pores. A major stumbling block in previous studies of freezing in pores has been the lack of any method for calculating the free energy difference between the confined solid and liquid phases. Conventional thermodynamic integration methods often fail for confined systems, due to the difficulty in choosing a suitable path of integration. We use a different approach that involves calculating the Landau free energy as a function of a suitable order parameter, using the grand canonical Monte Carlo simulation method. The grand free energy for each phase can then be obtained by one-dimensional integration of the Landau free energy over the order parameter. These calculations are carried out for two types of wall—fluid interaction, a hard wall and a strongly attractive wall modelled on carbon. The grand free energy results for both cases clearly indicate a first order fluid to solid transition. In the case of the attractive carbon wall, there are three phases. Phase A corresponds to all layers having a liquid-like structure; phase B corresponds to the contact layers (the layers adjacent to the two pore walls) being frozen and the rest of the layers being fluid-like; phase C corresponds to all the layers being frozen. Our results for the angular structure function in the individual molecular layers show strong evidence of a transition from a two-dimensional liquid phase to a hexatic phase. This is followed by a transition from the hexatic to a crystal phase.     

Thermodynamic modeling of oxide phases in the Fe–Mn–O system

A critical evaluation and thermodynamic modeling for thermodynamic properties of all oxide phases and phase diagrams in the Fe–Mn–O system are presented. Optimized Gibbs energy parameters for the thermodynamic models of the oxide phases were obtained which reproduce all available and reliable experimental data within error limits from 298 K to above the liquidus temperatures at all compositions covering from known oxide phases, and oxygen partial pressure from metal saturation to 0.21 bar. The optimized thermodynamic properties and phase diagrams are believed to be the best estimates presently available. Two spinel phases (cubic and tetragonal) were modeled using Compound Energy Formalism (CEF) with the use of physically meaningful parameters. The present Fe–Mn spinel solutions can be integrated into a larger spinel solution database, which has been already developed. The database of the model parameters can be used along with a software for Gibbs energy minimization in order to calculate any type of phase diagram section and thermodynamic properties.     

The role of the “Casimir force analogue” at the microscopic processes of crystallization and melting

Melting (crystallization), a phase transition from a crystalline solid to a liquid state, is a common phenomenon in nature. We suggest a new factor, “the Casimir force analogue”, to describe mechanisms of melting and crystallization. The Casimir force analogue is a force occurring between the surfaces of solid and liquid phases of metals caused by different energy density of phonons of these phases. It explains abrupt changes in geometry and thermodynamic parameters at a melting point. “The Casimir force analogue” helps to estimate latent melting heat and to gain an insight into a solid–liquid transition problem.     

Trees at an Interface

A lattice tree at an interface between two solvents of different quality is examined as a model of a branched polymer at an interface. Existence of the free energy is shown, and the existence of critical lines in its phase diagram is proven. In particular, there is a line of first order transitions separating a positive phase from a negative phase (the tree being predominantly on either side of the interface in these phases), and a curve of localization–delocalization transitions which separate the delocalized positive and negative phases from a phase where the tree is localized at the interface. This model is generalized to a branched copolymer which is examined in a certain averaged quenched ensemble. Existence of a thermodynamic limit is shown for this model, and it is also shown that the model is self-averaging. Lastly, a model of branched polymers interacting with one another through a membrane is considered. The existence of a limiting free energy is shown, and it is demonstrated that if the interaction is strong enough, then the two branched polymers will adsorb on one another.     

Coarsening of cracks in a uniaxially strained solid

We discuss the stress relaxation in a uniaxially strained solid due to the coarsening of a system of parallel cracks. We emphasize similarities and differences of this process to Ostwald ripening in a first order phase transition. A conventional mean-field approximation breaks down and several independent length scales have to be taken into account. Strong elastic interactions between the cracks determine the growth behavior. We derive scaling laws for the coarsening of the different length scales involved and the time evolution of stress relaxation, finally leading to the equilibrium state of a fractured body. The characteristic size of the cracks grows linearly in time which is much faster than in usual Ostwald ripening.     

Structural state and phase composition of nickel surface layers modified by high-intensity implantation of titanium ions

This paper reports on the results of experimental investigations into the microstructure and the elemental and phase compositions of ion-alloyed nickel surface layers modified by high-intensity implantation of titanium ions. It is established that the implantation of titanium ions into nickel surface layers up to 1600 nm thick leads to the formation of intermetallic phases (namely, NiTi, Ni₃Ti, and NiTi₂), a solid solution of titanium in nickel, titanium oxides of different stoichiometries, and titanium carbide TiC. It is demonstrated that the phases formed in the ion-alloyed nickel surface layers have a nanocrystalline structure. The average size of the nanocrystal grains is equal to 40 nm.     

On the size dependence of the melting temperature of nanoparticles

The size dependence of the melting temperature of nanocrystals has been investigated within the thermodynamic approach. A formula is obtained, which, in contrast to the classical Thomson formula, takes into account the metastable character of equilibrium between the crystal core and melt shell. Comparative investigation of the size dependence of the melting temperature, disregarding and taking into account the size dependences of the surface tension of the solid and liquid phases and the interface tension, has been performed by the example of aluminum, tin, and copper nanoparticles.     

The hard ellipsoid-of-revolution fluid I. Monte Carlo simulations

We present the results of Monte Carlo simulations on a system of hard ellipsoids of revolution with length-to-breadth ratios a/b = 3, 2.75, 2, 1.25 and b/a = 3, 2.75, 2, 1.25. We identify four distinct phases, viz. isotropic fluid, nematic fluid, ordered solid and plastic solid. The coexistence points of all first order phase transitions are located by performing absolute free energy computations for all coexisting phases. We find nematic phases only for a/b ≥ 2.75 and a/b ≤ 1/2.75. A plastic solid is only observed for 1.25 ≥ a/b ≥ 0.8. It is found that the phase diagram is surprisingly symmetric under interchange of the major and minor axes of the ellipsoids.     

Critical thermodynamic evaluation and optimization of the Fe-Mg-O system

A complete literature review, critical evaluation and thermodynamic modeling of the phase diagrams and thermodynamic properties at 1 bar total pressure of all oxide phases in the Fe-Mg-O system are presented. Optimized model equations for the thermodynamic properties of all phases are obtained which reproduce all available thermodynamic and phase equilibrium data within experimental error limits from 25 °C to above the liquidus temperatures at all compositions and oxygen partial pressures. The complex phase relationships in the system have been elucidated and discrepancies among the data have been resolved. The database of the model parameters can be used along with software for Gibbs energy minimization in order to calculate any type of phase diagram section. Sublattice models, based upon the compound energy formalism, were used for the spinel, pyroxene, olivine and monoxide phases. The use of physically reasonable models means that the models can be used to predict properties, phase equilibria, and cation site distributions in composition and temperature regions where data are not available.     

The phase coexistence method to obtain surface free energies and nucleation barriers: a brief review

ABSTRACT

A recently developed method where one analyses the finite size effects associated with liquid–solid phase equilibria including vapour–crystal coexistence is briefly reviewed. It is shown that the estimation of the chemical potential of the vapour surrounding the crystal as function of the crystal volume yields information on the bulk coexistence conditions, when an extrapolation to the thermodynamic limit is performed. Estimating the pressure of the fluid surrounding the crystal nucleus in the finite simulation box and the volume of this nucleus that coexists with the fluid in thermal equilibrium, an estimate for the total surface excess free energy can be obtained, which to a very good approximation is independent of the size of the simulation box. The free energy barrier against homogeneous nucleation of crystals thus can be estimated as a function of the nucleus volume. Monte Carlo simulations for the soft effective Asakura–Oosawa model of colloid-polymer mixtures which form face-centered cubic colloidal crystals are used to exemplify this method, computing the surface excess free energy of these crystals over a wide range of crystal volumes, without the need to characterise the non-spherical crystal shape. A possible extension of these concepts to heterogeneous nucleation is also briefly discussed.     

The influence of size on the composition of nano-precipitates in coherent precipitation

Various experiments show that the solute content in small coherent precipitates in their early stages of growth-coarsening is significantly lower than its equilibrium value and gradually increases with the particle size until equilibrium composition is reached. In this paper, we investigate the thermodynamic stability of a mono-dispersed assembly of coherent precipitates in a finite matrix by minimizing the total free energy of the system to account for this size effect observed for nano-phases. Both interfacial energy and misfit elastic energy were taken into account and simple regular solid solutions were considered for both phases. It is found that for small sizes, precipitates with a low solute content are favoured energetically. The solute content in precipitates decreases when increasing either interfacial energy or coherency misfit. The solute concentration in nano-particles is observed to increase with size. The asymptotic composition remains under the equilibrium concentration even for small misfits. The model was confronted to atom probe tomography experiments performed in FeCr system. Predictions exhibit a composition trend with precipitate size that is in good agreement with experiments.