Monte Carlo simulation of phase behavior of polymer blends with special interactions

Full Paper: The phase, behavior of A-B-random copolymer/C-homopolymer, blends with special interaction was studied by a. Monte, Carlo simulation in two dimensions. The interaction between I segment A and segment C was repulsive, whereas it was attractive between segment B and segment C. The simulation results showed that the blend became two large co-continuous phase domains at lower segment-B component compositions, indicating that the blend showed spinodal decomposition. With an increase of the segment-B component, the miscibility between the copolymer,and the polymer was gradually improved up to being miscible. In addition, it was found that segment B tended to move to the surface of the copolymer phase in the case of a lower component of segment B. On the other hand, if was observed that the average, end-to-end distances ((h) over bar) for both copolymer and polymer changed slowly with increasing segment-B component of the copolymer up to 40%, thereafter they increased considerably with increasing segment B component. Moreover, it was found that the (h) over bar of the copolymer was obviously shorter than that of the homopolymer for the segment-B composition, region from 0% to 80%. Finally, a, phase diagram showing I phase and - II phase regions under the condition of constant-temperature is presented.

Monte Carlo simulation of phase separation of A/B/A-B ternary mixtures

Monte Carlo simulations were used to model A/B/A-B ternary mixtures with different AB diblock copolymer volume fractions for which both the dispersed and continuous phase volume fractions were kept constant. For concentrations of the diblock copolymer below a critical value, the domain size increment of the dispersed phase decreases linearly with the copolymer concentration. This is in agreement with the predictions of Noolandi and Hong. The dependence of the domain size as a function of the copolymer volume fraction can also be fitted by the equation of Tang and Huang. Our simulations indicate, for the first time, that the micelles form before saturation of the interface occurs. This means that the formation of the micelles is not a result of the saturation of the interface.

Monte Carlo simulation of cluster growth in surface defects induced by the tip of a scanning tunnelling microscope

First steps are taken to model the electrochemical deposition of metals in nanometer-sized cavities. In the present work, the electrochemical deposition of Cu atoms in nanometer-sized holes dug on Au(111) is investigated through Monte Carlo simulations using the embedded atom method to represent particle interactions. By sweeping the chemical potential of Cu, a cluster is allowed to grow within the hole rising four atomic layers above the surface. Its lateral extension remains confined to the area defined by the borders of the original defect. (C) 2004 Elsevier B.V. All rights reserved.

Monte Carlo simulation of properties of monolayers and metal islands adsorbed on metallic (111) surfaces

To obtain the surface stress changes due to the adsorption of metal monolayers onto metallic surfaces, a new model derived from thermodynamic considerations is presented. Such a model is based on continuum Monte Carlo simulations with embedded atom method potentials in the canonical ensemble, and it is extended to consider the behavior on different islands adsorbed onto (111) substrate surfaces. Homoepitaxial and heteroepitaxial systems are studied. Pseudomorphic growth is not observed for small metal islands with considerable positive misfit with the substrate. Instead, the islands become compressed upon increase of their size. A simple model is proposed to interpolate between the misfits of atoms in small islands and the pseudomorphic behavior of the monolayer.

Atomic ordering in nano-layered FePt: Multiscale Monte Carlo simulation

Nano- and meso-scale simulation of chemical ordering kinetics in nano-layered L1(0)-AB binary intermetallics was performed. In the nano- (atomistic) scale Monte Carlo (MC) technique with vacancy mechanism of atomic migration implemented with diverse models for the system energetics was used. The meso-scale microstructure evolution was, in turn, simulated by means of a MC procedure applied to a system built of meso-scale voxels ordered in particular L1(0) variants. The voxels were free to change the L1(0) variant and interacted with antiphase-boundary energies evaluated within the nano-scale simulations. The study addressed FePt thin layers considered as a material for ultra-high-density magnetic storage media and revealed metastability of the L1(0) c-variant superstructure with monoatomic planes parallel to the (001)-oriented layer surface and off-plane easy magnetization. The layers, originally perfectly ordered in the c-variant, showed discontinuous precipitation of a- and b-L1(0)-variant domains running in parallel with homogeneous disordering (i.e. generation of antisite defects). The domains nucleated heterogeneously on the free monoatomic Fe surface of the layer, grew inwards its volume and relaxed towards an equilibrium microstructure of the system. Two

Monte Carlo simulation of complex cohesive fracture in random heterogeneous quasi-brittle materials

A numerical method is developed to simulate complex two-dimensional crack propagation in quasi-brittle materials considering random heterogeneous fracture properties. Potential cracks are represented by pre-inserted cohesive elements with tension and shear softening constitutive laws modelled by spatially varying Weibull random fields. Monte Carlo simulations of a concrete specimen under uni-axial tension were carried out with extensive investigation of the effects of important numerical algorithms and material properties on numerical efficiency and stability, crack propagation processes and load-carrying capacities. It was found that the homogeneous model led to incorrect crack patterns and load–displacement curves with strong mesh-dependence, whereas the heterogeneous model predicted realistic, complicated fracture processes and load-carrying capacity of little mesh-dependence. Increasing the variance of the tensile strength random fields with increased heterogeneity led to reduction in the mean peak load and increase in the standard deviation. The developed method provides a simple but effective tool for assessment of structural reliability and calculation of characteristic material strength for structural design.

Monte Carlo simulation of interface alloying

An Ising-like model, with interactions ranging up to next-nearest-neighbor pairs, is used to simulate the process of interface alloying. Interactions are chosen to stabilize an intermediate "antiferromagnetic" ordered structure. The dynamics proceeds exclusively by atom-vacancy exchanges. In order to characterize the process, the time evolution of the width of the intermediate ordered region and the diffusion length is studied. Both lengths are found to follow a power-law evolution with exponents depending on the characteristic features of the model.

Modeling premartensitic effects in Ni2MnGa: A mean-field and Monte Carlo simulation study

ic first-order transition line ending in a critical point. This critical point is responsible for the existence of large premartensitic fluctuations which manifest as broad peaks in the specific heat, not always associated with a true phase transition. The main conclusion is that premartensitic effects result from the interplay between the softness of the anomalous phonon driving the modulation and the magnetoelastic coupling. In particular, the premartensitic transition occurs when such coupling is strong enough to freeze the involved mode phonon. The implication of the results in relation to the available experimental data is discussed.

A cluster computing implementation of a Monte Carlo simulation of a particle growth mechanism

Clusters of computers can be used together to provide a powerful computing resource. Large Monte Carlo simulations, such as those used to model particle growth, are computationally intensive and take considerable time to execute on conventional workstations. By spreading the work of the simulation across a cluster of computers, the elapsed execution time can be greatly reduced. Thus a user has apparently the performance of a supercomputer by using the spare cycles on other workstations.