Electromigration revisited: competition between Kirkendall shift and backstress in pure metals and two-phase alloys

The simple phenomenological model and analytical approach of electromigration in the two-phase alloy (solder) under combined influence of the Kirkendall effect, backstress and sedimentation is presented. It is compared with electromigration in pure metal under condition of quasi-equilibrium vacancies (unlimited power of vacancy sinks-sources) and electromigration in pure metal with account of nonequilibrium vacancies.     

Monte Carlo simulation of phase separation and clustering in the ABV model

As a model for a binary alloy undergoing an unmixing phase transition, we consider a square lattice where each site can be either taken by an A atom, a B atom, or a vacancy (V), and there exists a repulsive interaction between AB nearest neighbor pairs. Starting from a random initial configuration, unmixing proceeds via random jumps of A atoms or B atoms to nearest neighbor vacant sites. In the absence of any interaction, these jumps occur at jump rates _A and _B, respectively. For a small concentration of vacancies (c _v=0.04) the dynamics of the structure factorS(k,t) and its first two momentsk ₁(t),k ₂ ² (t) is studied during the early stages of phase separation, for several choices of concentrationc _B of B atoms. Forc _B=0.18 also the time evolution of the cluster size distribution is studied. Apart from very early times, the mean cluster sizel(t) as well as the moments of the structure function depend on timet and the ratio of the jump rates (= _B/ _A) only via a scaled timet/(). Qualitatively, the behavior is very similar to the direct exchange model containing no vacancies. Consequences for phase separation of real alloys are briefly discussed.     

Dynamic percolation transition induced by phase separation: A Monte Carlo analysis

The percolation transition of geometric clusters in the three-dimensional, simple cubic, nearest neighbor Ising lattice gas model is investigated in the temperature and concentration region inside the coexistence curve. We consider quenching experiments, where the system starts from an initially completely random configuration (corresponding to equilibrium at infinite temperature), letting the system evolve at the considered temperature according to the Kawasaki spinexchange dynamics. Analyzing the distributionn _l(t) of clusters of sizel at timet, we find that after a time of the order of about 100 Monte Carlo steps per site a percolation transition occurs at a concentration distinctly lower than the percolation concentration of the initial random state. This dynamic percolation transition is analyzed with finite-size scaling methods. While at zero temperature, where the system settles down at a frozen-in cluster distribution and further phase separation stops, the critical exponents associated with this percolation transition are consistent with the universality class of random percolation, the critical behavior of the transient time-dependent percolation occurring at nonzero temperature possibly belongs to a different, new universality class.     

Polymer mixtures in confined geometries: Model systems to explore phase transitions

While binary (A,B) symmetric polymer mixtures ind = 3 dimensions have an unmixing critical point that belongs to the 3d Ising universality class and crosses over to mean field behavior for very long chains, the critical behavior of mixtures confined into thin film geometry falls in the 2d Ising class irrespective of chain length. The critical temperature always scales linearly with chain length, except for strictly two-dimensional chains confined to a plane, for whichT _c ∝N ^5/8 (this unusual exponent describes the fractal contact line between segregated chains in dense melts in two spatial dimensions,d = 2). When the walls of the thin film are not neutral, but preferentially attract one species, complex phase diagrams occur due to the interplay between capillary condensation and wetting phenomena. For ‘competing walls’ (one wall prefers A, the other prefers B) particularly interesting interface localization-delocalization transitions occur, while analogous phenomena in wedges are related to the ‘filling transition’.     

The effect of second neighbour repulsion on fcc binary alloy phase diagrams: A Monte Carlo study

We report the results of the Monte Carlo simulation of the phase diagram of fcc binary alloys using a 3-D Ising model with nearest and next-nearest neighbour repulsive interactions. Calculations are carried out for a ratio of second- to first-neighbour interaction energies of 0.4. The resulting phase diagram contains three superstructures A₂B₂, A₂B and A₅B, each separated by a disordered fcc phase. There was no evidence for the formation of an A₃B phase. These results are in qualitative agreement with CVM results.     

Influence of atomic size mismatch on binary alloy phase diagrams

Abstract

The aim of this paper is to investigate the consequences of atomic size mismatch on the thermodynamics and the topology of binary phase diagrams of face centred cubic alloys. Simple pairwise interatomic potentials with few controlling parameters are used to identify general tendencies. Thermodynamic states are computed by Monte Carlo simulations on a non-rigid lattice. A special attention has been paid to the comparison between calculations in the canonical ensemble, where composition–temperature phase diagrams are determined through van der Waals loops, and in the grand canonical ensemble, where phase diagrams are computed using an interface migration technique. It is shown that these two procedures lead essentially to the same incoherent phase diagram. In the case of phase separating systems, we argue that the introduction of a size mismatch leads to a shrinkage of the solid solution domain and that the asymmetry of the miscibility gap is essentially controlled by the anharmonicity of the heteroatomic potential. Finally, in the case of ordering systems, we show that the asymmetry of the phase diagram may be due to the anharmonicity of the pair potentials or to the differences between their curvatures, the former effect being dominant if the atomic size mismatch is large.     

三维X-Y模型的滞后标度和动态相变行为

采用Monte Carlo方法对3DX-Y模型进行数值模拟计算，研究了在非线性外场驱动下3DX-Y模 型的滞后标度和动态相变.得出了滞后标度关系为Area～h^α₀ωβ(1-T/T_c)^γ， 其中α=0.57，β=0.34，γ=0.90.发现其动态相变行为在一定的临界参数条件下，初始短周 期(周期数PN≤10)内的结果具有与Ising模型类似的对称性破缺；但在长周期内(PN≥200)的 结果则明显区别于Ising模型而与Heisenberg模型相近，也即无稳定动态有序铁磁相的存在. 关键词： X-Y模型 滞后标度 动态相变 Monte Carlo方法     

两相钛合金再结晶退火组织与织构演变的蒙特卡罗模拟

蒙特卡罗(MC)方法被广泛应用于模拟金属材料在退火过程中的静态再结晶行为. 在已有两相材料晶粒长大MC模型基础上, 引入形核阶段, 综合考虑再结晶晶粒吞并形变晶粒和再结晶晶粒竞争长大两种情况, 建立了退火时两相合金再结晶MC模型.结合电子背散射衍射所测 初始晶粒形貌、相成分、晶体学取向及应变储能相对值, 该模型被应用于TC11钛合金退火过程中的微观组织及织构演变模拟.结果表明, 所建模型能够较好体现退火过程中两相晶粒的形核及晶粒长大行为. 与β相相比较, α相具有较低的再结晶速率和较高的晶粒长大速率, 前者主要归结于α相较低的初始应变储能, 后者则体现了该条件下初始组织形貌、分布及两相比例对晶粒长大具有重要影响; 由于非均匀形核的影响, 模拟得到的再结晶速率变化与 假设均匀形核的Johnson-Mehl-Avrami-Kolmogorov 再结晶方程存在明显差异.同时, 两相的基本织构特征在退火过程中无明显变化, 但织构强度增加. 关键词： 两相钛合金 再结晶 蒙特卡罗方法 织构     

Resistivity and critical temperature of (La1-xPrx)0.7Ca0.3MnO3

An empirical formula of the critical temperature that is concentration dependent for polycrystalline (La_1-xPr_x)_0.7Ca_0.3MnO₃ is presented in this paper. With this formula, the temperature dependence of resistance is simulated for various values of x by using the random resistor network model and the Monte Carlo method. The hysteresis effect in \rho -- T curves is reasonably explained. The simulation results are in good agreement with the relevant experimental measurements.     

Effect of Polymer-Substrate Interactions on the Surface Morphology of Polymer Blend Thin Films: Comparison between Simulation and Experiment

Microphase and macrophase separation phenomena can simultaneously appear in ABA/C copolymer blend systems due to the immiscibility among monomers A, B, and C. In this work, the surface morphologies and compositions of ABA/C blend thin films confined between two walls, which were used to mimic SEBS/PMMA films, have been simulated by a lattice Monte Carlo (MC) method. The effect of the polymer-wall interaction on the surface morphologies and compositions of thin films was investigated as a function of blend composition and film thickness. It is shown that the simulated surface morphologies of thin films resulting from the macrophase separation between copolymer ABA and homopolymer C and the microphase separation between block A and block B in ABA copolymer are similar to the experimental surface morphology of SEBS/PMMA polymer blend films observed by atomic force microscope (AFM). The effect of substrate on the surface morphologies by MC simulation is qualitatively consistent with the experimental results. The composition profiles of thin films are given to characterize the micro- and macrophase separation in thin films. It is indicated that the surface energy of the substrate (substrate/air) plays a crucial role on the surface composition. For a fixed surface, the adsorptions of polymer on the substrate and film thickness are also important.     

高压作用下相分离液体玻璃转变的分子动力学研究

采用分子动力学模拟方法,研究了二元混合液体在不同外压作用下的相分离与玻璃转变过程,计算了相分离液体在玻璃转变过程中的结构和动力学特征.研究发现,外压会促进相分离的产生,并提高玻璃转变温度,会使β弛豫出现的温度更高、存在的时间更长,导致系统扩散性降低.同时还发现,相分离液体的玻璃转变过程存在微观不均匀现象. 关键词： 相分离 玻璃转变 分子动力学模拟 外压影响     

二元Lennard-Jones液体的相分离过程及其扩散性质的分子动力学研究

采用分子动力学(MD)模拟方法,研究了二元体系中相分离过程、粒子的扩散系数以及相分离域尺寸大小随温度的变化规律.发现,相分离域随温度的生长过程可以分为两个阶段,分别是温度比较高的快速生长阶段和低温时的稳定生长阶段;相分离体系中系统的扩散激活能不是常数,而是一个随温度变化的函数,并且当温度高于60 K时,满足关系式E(T)=a+bT^c.讨论了组元尺寸的变化对相分离过程的影响.结果表明,随两组元中某一组元 关键词： 相分离 扩散 分子动力学模拟     

The deconfining phase transition in lattice quantum chromodynamics

We present a large-scale Monte Carlo calculation of the deconfining phase transition temperature in lattice quantum chromodynamics without fermions. Using the Wilson action, the transition temperature as a function of the lattice couplingg is consistent with scaling behavior dictated by the perturbativeα function for 6/g²>6.15. Speaker at the conference; on leave from CRIP, Budapest.     

Monte Carlo study of dynamic phase transition in Ising metamagnet driven by oscillating magnetic field

The dynamical responses of Ising metamagnet (layered antiferromagnet) in the presence of a sinusoidally oscillating magnetic field are studied by Monte Carlo simulation. The time average staggered magnetisation plays the role of dynamic order parameter. A dynamical phase transition was observed and a phase diagram was plotted in the plane formed by field amplitude and temperature. The dynamical phase boundary is observed to shrink inward as the relative antiferromagnetic strength decreases. The results are compared with that obtained from pure ferromagnetic system. The shape of dynamic phase boundary observed to be qualitatively similar to that obtained from previous meanfield calculations.     

Phase separation of symmetrical polymer mixtures in thin-film geometry

Monte Carlo simulations of the bond fluctuation model of symmetrical polymer blends confined between two neutral repulsive walls are presented for chain lengthN _A=N _B=32 and a wide range of film thicknessD (fromD=8 toD=48 in units of the lattice spacing). The critical temperaturesT _c (D) of unmixing are located by finite-size scaling methods, and it is shown that , wherev ₃0.63 is the correlation length exponent of the three-dimensional Ising model universality class. Contrary to this result, it is argued that the critical behavior of the films is ruled by two-dimensional exponents, e.g., the coexistence curve (difference in volume fraction of A-rich and A-poor phases) scales as , where ₂ is the critical exponent of the two-dimensional Ising universality class ( ₂=1/8). Since for largeD this asymptotic critical behavior is confined to an extremely narrow vicinity ofT _c (D), one observes in practice effective exponents which gradually cross over from ₂ to ₃ with increasing film thickness. This anomalous flattening of the coexistence curve should be observable experimentally.     

Metastable phase separation and rapid solidification of undercooled Co-Cu alloy under different conditions

The metastable liquid phase separation and rapid solidification behaviours of Co_61.8Cu_38.2 alloy were investigated by using differential thermal analysis (DTA) in combination with glass fluxing, electromagnetic levitation (EML) and drop tube techniques. It is found that the liquid phase separation process and the solidification microstructures intensively depend on the experimental processing parameters, such as undercooling level, cooling rate, gravity level, liquid surface tension and the wetting state of crucible. Large undercooling and surface tension difference of the two liquids tend to facilitate further separation and cause severe macrosegregation. On the other hand, rapid cooling and low gravity effectively suppress the coalescence of the minority phase. Severe macrosegregation patterns are formed in the bulk samples processed by both DTA and EML. In contrast, disperse structures with fine spherical Cu-rich spheres homogeneously distributed in the matrix of Co-rich phase have been obtained in drop tube.     

Nonequilibrium dynamical phase transition of 3D kinetic Ising／Heisenberg spin system

We have studied the nonequilibrium dynamic phase transitions of both three-dimensional (3D) kinetic Ising and Heisenberg spin systems in the presence of a perturbative magnetic field by Monte Carlo simulation. The feature of the phase transition is characterized by studying the distribution of the dynamical order parameter. In the case of anisotropic Ising spin system (ISS), the dynamic transition is discontinuous and continuous under low and high temperatures respectively, which indicates the existence of a tri-critical point (TCP) on the phase boundary separating low-temperature order phase and high-temperature disorder phase. The TCP shifts towards the higher temperature region with the decrease of frequency, i.e. T_{TCP}=1.33×exp(-ω/30.7). In the case of the isotropic Heisenberg spin system (HSS), however, the situation on dynamic phase transition of HSS is quite different from that of ISS in that no stable dynamical phase transition was observed in kinetic HSS after a threshold time. The evolution of magnetization in the HSS driven by a symmetrical external field after a certain duration always tends asymptotically to a disorder state no matter what an initial state the system starts with. The threshold time τ depends upon the amplitude H_{0}, reduced temperature T/T_C and the frequency ω as τ=C·ω^α·H_0^{-β}·(T/T_C)^{-γ}.     

基于伊辛模型的Fe_(0.5)Mn_(0.1)Al_(0.4)合金磁化强度和磁熵变蒙特卡洛模拟

基于伊辛模型的单自旋反转蒙特卡洛算法,考虑了粒子间的最近邻以及次近邻相互作用,研究了无序Fe0.5Mn0.1Al0.4合金的磁化强度和磁熵变.首先,强调了粒子间的次近邻相关作用对体系的磁性和热力学性质的影响,明确了次近邻相互作用系数,证实了低温合金阻挫的存在;其次,研究了在相变温度处(不同磁场下)磁化强度随外加磁场(温度)的变化情况以及磁性粒子对磁化强度的贡献,发现反铁磁性粒子Mn在低温区对Fe0.5Mn0.1Al0.4合金的相变起了主要作用,而高温区体系的相变是由铁磁性粒子Fe贡献的;最后,分析了体系在相变温度处磁熵变数值随外加磁场的变化情况以及磁熵变在不同的外磁场下随温度的变化情况,当外加磁场H=0.14(a.u.)时,Mn粒子在冻结温度处的平均磁化强度为零,体系处于最无序的状态,对应的磁熵变ΔS(0.1,0.14)达到了正向最大值,极值的位置对应于体系的相变温度.