A Monte Carlo simulation for the micellization of ABC 3-miktoarm star terpolymers in a selective solvent

We have used Monte Carlo simulation to study the micellization of ABC 3-miktoarm star terpolymers in a selective solvent (good to A segment, bad to B and C segments). The simulation results reveal that the self-assembled morphology is determined by the block length, molecular architecture, terpolymer concentration and insolubility of insoluble block in the solvent. In dilute solution, symmetric terpolymers (N-B = N-C = 30) tend to aggregate into a novel wormlike pearl-necklace structure linked by an alternating arrangement of B and C spheres, whereas the asymmetric terpolymers (NB = 10, NC = 50) are likely to aggregate into spherical or cylindrical micelles (formed by C blocks) connected with some small B spheres, when the concentration of terpolymer is relatively low (chain number is 100). However, when the concentration of terpolymer is relatively high (chain number is 250), the symmetric terpolymers tend to aggregate into a netlike structure linked by an alternation of B and C spheres, whereas the asymmetric terpolymers are likely to aggregate into wormlike micelles (formed by C blocks) connected with some of small spheres (formed by B blocks). Moreover, when the insolubility of insoluble block in the solvent is weak, the insoluble blocks aggregate into some incompact micelles.

Monte Carlo simulation of the aggregation of rod-flexible triblock copolymers in a thin film

The aggregation of rod-flexible ABA and BAB triblock (A was rod block and repulsive with block B) copolymers in a thin film was studied as a function of varying the rigidity (eta) and the length of the rod block by Monte Carlo simulation. The rigidity of block A was defined as eta = R-c/R-max in this study. R-c, was the end-to-end distance below which the conformation of the block was not allowed, whereas R-max, was the longest end-to-end distance that the block could be. If eta = 0 the block was flexible, whereas if eta = 1 the block was a straight rod. The simulation results showed that the ABA triblock copolymer film were likely to form lamella structure with increasing the rigidity (eta) of block A. The lamellas were parallel each other and perpendicular to the film surface. However, the aggregation of BAB triblock copolymers tended to change from lamella to cylinder structure with increasing the rigidity (eta) of block A. Typical lamella and cylinder co-exist structure was obtained at eta = 0.504 for the BAB copolymer film. On the other hand, the simulation results indicated that the film changed from disorder to order, then to disorder structure with increasing the relative length of B block for both ABA and BAB copolymer films.

Monte Carlo simulation of miscibility of polymer blends with repulsive interactions: effect of chain structure

The effects of the chain structure and the intramolecular interaction energy of an A/B copolymer on the miscibility of the binary blends of the copolymer and homopolymer C have been studied by means of a Monte Carlo simulation. In the system, the interactions between segments A, B and C are more repulsive than those between themselves. In order to study the effect of the chain structure of the A/B copolymer on the miscibility, the alternating, random and block copolymers were introduced in the simulations, respectively. The simulation results show that the miscibility of the binary blends strongly depends on the intramolecular interaction energy ((ε) over bar (AB)) between segments A and B within the A/B copolymers. The higher the repulsive interaction energy, the more miscible the A/B copolymer and homopolymer C are. For the diblock copolymer/homopolymer blends, they tend to form micro phase domains. However, the phase domains become so small that the blend can be considered as a homogeneous phase for the alternating copolymer/ homopolymer blends. Furthermore, the investigation of the average end-to-end distance ((h) over bar) in different systems indicates that the copolymer chains tend to coil with the decrease Of (ε) over bar (AB) whereas the (h) over bar of the homopolymer chains depends on the chain structure of the copolymers.

Computer simulation of miscibility and self-assembly structure for polymer-containing systems with special interactions

The miscibility and structure of A-B copolymer/C homopolymer blends with special interactions were studied by a Monte Carlo simulation in two dimensions. The interaction between segment A and segment C was repulsive, whereas it was attractive between segment B and segment C. In order to study the effect of copolymer chain structure on the morphology and structure of A-B copolymer/C homopolymer blends, the alternating, random and block A-B copolymers were introduced into the blends, respectively. The simulation results indicated that the miscibility of A-B block copolymer/C homopolymer blends depended on the chain structure of the A-B copolymer. Compared with alternating or random copolymer, the block copolymer, especially the diblock copolymer, could lead to a poor miscibility of A-B copolymer/C homopolymer blends. Moreover, for diblock A-B copolymer/C homopolymer blends, obvious self-organized core-shell structure was observed in the segment B composition region from 20% to 60%. However, if diblock copolymer composition in the blends is less than 40%, obvious self-organized core-shell structure could be formed in the B-segment component region from 10 to 90%. Furthermore, computer statistical analysis for the simulation results showed that the core sizes tended to increase continuously and their distribution became wider with decreasing B-segment component.

Monte Carlo simulation of the grafting of maleic anhydride onto polypropylene at higher temperature

The graft of maleic anhydride (MAH) onto isotactic polypropylene (iPP) initiated by dicumyl peroxide (DCP) at 190 degreesC was studied by means of the Monte Carlo method. The ceiling temperature theory, i.e., no possibility for the homopolymerization of MA-H to occur at higher temperatures, was used in this study. The simulation results show that most MAH monomers were grafted onto the radical chain ends arising from beta scission at a lower MAH concentration, whereas the amount of MAH monomers attached to the tertiary carbons was much larger than that grafted onto the radical chain ends at a higher MAH concentration for various DCP concentrations. This conclusion gives a good interpretation for the disagreement on the grafting sites along a PP chain. Moreover, it was found that the grafting degree increased considerably up to a peak value; thereafter, it decreased continuously with increasing MA-H concentration. The peak shifted in the lower MAH concentration direction and became lower and lower with increasing DCP concentration. When the DCP concentration was below 0.1 wt %, the peak was hardly observed. Those results are in good agreement with the experiments.

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.

Formation of Fivefold Deformation Twins in Nanocrystalline Face-Centered-Cubic Copper Based on Molecular Dynamics Simulations

Fivefold deformation twins were reported recently to be observed in the experiment of the nanocrystalline face-centered-cubic metals and alloys. However, they were not predicted previously based on the molecular dynamics (MD) simulations and the reason was thought to be a uniaxial tension considered in the simulations. In the present investigation, through introducing pretwins in grain regions, using the MD simulations, the authors predict out the fivefold deformation twins in the grain regions of the nanocrystal grain cell, which undergoes a uniaxial tension. It is shown in their simulation results that series of Shockley partial dislocations emitted from grain boundaries provide sequential twining mechanism, which results in fivefold deformation twins. (c) 2006 American Institute of Physics.

An analysis on overall crack-number-density of short-fatigue-cracks

The evolution of dispersed short-fatigue-cracks is analysed based on the equilibrium of crack-number-density (CND). By separating the mean value and the stochastic fluctuation of local CND, the equilibrium equation of overall CND is derived. Comparing with the mean-field equilibrium equation, the equilibrium equation of overall CND has different forms in the expression of crack-nucleation-rate or crack-growth-rate. The simulation results are compared with experimental measurements showing the stochastic analyses provide consistent tendency with experiments. The discrepancy in simulation results between overall CND and mean-field CND is discussed.

Size Effect and Geometrical Effect of Polycrystals and Thin Film/Substrate System in Micro-indentation Test

Micro-indentation test at scales on the order of sub-micron has shown that the measured hardness increases strongly with decreasing indent depth or indent size, which is frequently referred to as the size effect. Simultaneously, at micron or sub-micron scale, the material microstructure size also has an important influence on the measured hardness. This kind of effect, such as the crystal grain size effect, thin film thickness effect, etc., is called the geometrical effect by here. In the present research, in order to investigate the size effect and the geometrical effect, the micro-indentation experiments are carried out respectively for single crystal copper and aluminum, for polycrystal aluminum, as well as for a thin film/substrate system, Ti/Si3N4. The size effect and geometrical effect are displayed experimentally. Moreover, using strain gradient plasticity theory, the size effect and the geometrical effect are simulated. Through comparing experimental results with simulation results, length-scale parameter appearing in the strain gradient theory for different cases is predicted. Furthermore, the size effect and the geometrical effect are interpreted using the geometrically necessary dislocation concept and the discrete dislocation theory. Member Price: $0; Non-Member Price: $25.00

Study of Gliding Arc Discharge Plasma

本文研究了滑动弧放电过程中电参数的变化，并对滑动弧等离子体中的非平衡度和各参数之间关系进行了讨论。应用了双通道电弧模型 ，对电弧在气流作用下的运动规律进行了数值模拟。模拟的结果有助于分析滑动弧非平衡等离子体的产生机理。 The elelctric parameters change during discharge is studied and the relationship between non-equilibrium degree and parameters is discussed for gliding arc discharge.Using two-channel model, the rules of arc moving due to effect of the airflow is simulated.The numerical simulation results can help analyzing the generation mechanism of gliding arc non-equilibrium plasma.

Research on the Temperature Field in Laser Hardening

Temperature field in the laser hardening process was numerically simulated by MSC.Marc software. The influence of energy density on laser hardening effect is analyzed. Simulation result is verified through the thermocouple temperature transducer measuring the specimen surface temperature under the laser irradiation. Experimental curves of temperature versus time are in agreement with simulation results. The simulation results can be regarded as a basis for choosing laser technological parameters.

Squeeze-film Effects in MEMS Devices with Perforated Plates for Small Amplitude Vibration

Squeeze-film effects of perforated plates for small amplitude vibration are analyzed through modified Reynolds equation (MRE). The analytical analysis reckons in most important influential factors: compressibility of the air, border effects, and the resistance caused by vertical air flow passing through perforated holes. It is found that consideration of air compressibility is necessary for high operating frequency and small ratio of the plate width to the attenuation length. The analytical results presented in this paper agree with ANSYS simulation results better than that under the air incompressibility assumption. The analytical analysis can be used to estimate the squeeze-film effects causing damping and stiffness added to the system. Since the value of Reynolds number involved in this paper is low (< 1), inertial effects are neglected.

Dislocation evolution in epitaxial multilayers and graded composition buffers

This paper presents models to describe the dislocation dynamics of strain relaxation in an epitaxial uniform layer, epitaxial multilayers and graded composition buffers. A set of new evolution equations for nucleation rate and annihilation rate of threading dislocations is developed. The dislocation interactions are incorporated into the kinetics process by introducing a resistance term, which depends only on plastic strain. Both threading dislocation nucleation and threading dislocation annihilation are characterized. The new evolution equations combined with other evolution equations for the plastic strain rate, the mean velocity and the dislocation density rate of the threading dislocations are tested on GexSi1-x/Si(100) heterostructures, including epitaxial multilayers and graded composition buffers. It is shown that the evolution equations successfully predict a wide range of experimental results of strain relaxation and threading dislocation evolution in the materials system. Meanwhile, the simulation results clearly signify that the threading dislocation annihilation plays a vital role in the reduction of threading dislocation density.

Micro-scale Mechanics of the Surface-Nanocrystalline Al-Alloy Material

Based on the microscopic observations and measurements, the mechanical behavior of the surface-nanocrystallized Al-alloy material at microscale is investigated experimentally and theoretically. In the experimental research, the compressive stress-strain curves and the hardness depth curves are measured. In the theoretical simulation, based on the material microstructure characteristics and the experimental features of the compression and indentation, the microstructure cell models are developed and the strain gradient plasticity theory is adopted. The material compressive stress-strain curves and the hardness depth curves-are predicted and simulated. Through comparison of the experimental results with the simulation results, the material and model parameters are determined.

Molecular Dynamics Study Of Mechanical Behaviour Of Screw Dislocation During Cutting With Diamond Tip On Silicon

By means of Tersoff and Morse potentials, a three-dimensional molecular dynamics simulation is performed to study atomic force microscopy cutting on silicon monocrystal surface. The interatomic forces between the workpiece and the pin tool and the atoms of workpiece themselves are calculated. A screw dislocation is introduced into workpiece Si. It is found that motion of dislocations does not occur during the atomic force microscopy cutting processing. Simulation results show that the shear stress acting on dislocation is far below the yield strength of Si.