A Monte Carlo study of crystal structure transformations

The Metropolis algorithm has been generalized to allow for the variation of shape and size of the MC cell. A calculation using different potentials illustrates how the generalized method can be used for the study of crystal structure transformations. A restricted MC integration in the nine dimensional space of the cell components also leads to the stable structure for the Lennard-Jones potential.

Glass-Transition In Isopentane - A Monte-Carlo Study

Monte Carlo simulations with realistic interaction potentials have been carried out on isopentane to investigate the glass transition. Intermolecular pair-correlation functions of the glass show distinct differences from those of the liquid, the CH-CH pair-correlation function being uniquely different from the other pair-correlation functions. The coordination number of the glass is higher than that of the liquid, and the packing in the glass seems to be mainly governed by the geometrical constraints of the molecule. Annealing affects the properties of the glass significantly.

Orientational glassy phases of C60 and neopentane. A Monte Carlo study

Based on an isothermal, isobaric simulation the structure and properties of the plastic crystalline phases of C60 and neopentane have been examined. Instantaneous cooling of the plastic crystalline phases of both C60 and neopentane leads to orientational glassy phases. These are accompanied by significant slowing down of reorientational motion. Constant pressure quench experiments on C60 yield a glass transition temperature of around 80 K.

Role of molecular reorientation in the vitrification of molecular liquids: A Monte Carlo study of the liquid and glassy states of 2,2-dimethylbutane

The liquid and the glassy phases of 2,2-dimethylbutane have been investigated by isothermal isobaric ensemble Monte Carlo simulation. Thermodynamic Properties and radial distribution functions for both the liquid and the glass have been obtained. The radial distribution functions have been classified into three types based on the accessibility of the group. It has been shown that the structure of the Iiquid and the glass can be understood in terms of the above classification of the radial distribution functions. Molecular reorientation plays an important role in the structural rearrangement accompanying glass formation. As much as 35% of the contribution to the increase in the intermolecular interaction energy on vitrification is due to the reorientation of the neighbouring pairs of molecules. The observed changes in the dimerisation energy and the bonding energy distribution function are consistent with the observed structural changes.

Monte Carlo study of the phase boundary of YBa2Cu3O6 +x

A model hamiltonian previously introduced to study the oxygen ordering is considered. The phase boundary is isolated by studying the relaxation behaviour of the order parameters. Our results are consistent with the published Monte Carlo results except at low temperatures.

Nonequilibrium-phase transition and 'specific-heat' singularity in the kinetic Ising model: A Monte Carlo study

The nonequilibrium-phase transition has been studied by Monte Carlo simulation in a ferromagnetically interacting (nearest-neighbour) kinetic Ising model in presence of a sinusoidally oscillating magnetic field. The ('specific-heat') temperature derivative of energies (averaged over a full cycle of the oscillating field) diverge near the dynamic transition point.

Quantum Monte Carlo study of the CO interaction with a dimer model surface for Cr(110)

The chemisorption of CO on a Cr( 110) surface is investigated using the quantum Monte Carlo method in the diffusion Monte Carlo (DMC) variant and a model Cr2CO cluster. The present results are consistent with the earlier ab initio HF study with this model that showed the tilted/ near-parallel orientation as energetically favoured over the perpendicular arrangement. The DMC energy difference between the two orientations is larger (1.9 eV) than that computed in the previous study. The distribution and reorganization of electrons during CO adsorption on the model surface are analysed using the topological electron localization function method that yields electron populations, charge transfer and clear insight on the chemical bonding that occurs with CO adsorption and dissociation on the model surface.

Asymptotic scaling in the two-dimensional O(3) Nonlinear sigma model: a Monte Carlo study on parallel computers

We investigate the 2d O(3) model with the standard action by Monte Carlo simulation at couplings β up to 2.05. We measure the energy density, mass gap and susceptibility of the model, and gather high statistics on lattices of size L ≤ 1024 using the Floating Point Systems T-series vector hypercube and the Thinking Machines Corp.'s Connection Machine 2. Asymptotic scaling does not appear to set in for this action, even at β = 2.10, where the correlation length is 420. We observe a 20% difference between our estimate m/Λ^─_(Ms) = 3.52(6) at this β and the recent exact analytical result . We use the overrelaxation algorithm interleaved with Metropolis updates and show that decorrelation time scales with the correlation length and the number of overrelaxation steps per sweep. We determine its effective dynamical critical exponent to be z' = 1.079(10); thus critical slowing down is reduced significantly for this local algorithm that is vectorizable and parallelizable.

We also use the cluster Monte Carlo algorithms, which are non-local Monte Carlo update schemes which can greatly increase the efficiency of computer simulations of spin models. The major computational task in these algorithms is connected component labeling, to identify clusters of connected sites on a lattice. We have devised some new SIMD component labeling algorithms, and implemented them on the Connection Machine. We investigate their performance when applied to the cluster update of the two dimensional Ising spin model.

Finally we use a Monte Carlo Renormalization Group method to directly measure the couplings of block Hamiltonians at different blocking levels. For the usual averaging block transformation we confirm the renormalized trajectory (RT) observed by Okawa. For another improved probabilistic block transformation we find the RT, showing that it is much closer to the Standard Action. We then use this block transformation to obtain the discrete β-function of the model which we compare to the perturbative result. We do not see convergence, except when using a rescaled coupling β_E to effectively resum the series. For the latter case we see agreement for m/ Λ^─_(Ms) at , β = 2.14, 2.26, 2.38 and 2.50. To three loops m/Λ^─_(Ms) = 3.047(35) at β = 2.50, which is very close to the exact value m/ Λ^─_(Ms) = 2.943. Our last point at β = 2.62 disagrees with this estimate however.

Effect of polydispersity on the self-assembly structure of diblock copolymers under various confined states: A Monte Carlo study

We studied the self-assembly of polydisperse diblock copolymers under various confined states by Monte Carlo simulation. When the copolymers were confined within two parallel walls, it was found that the ordered strip structures appeared alternately with the increase in wall width. Moreover, the wall width at which the ordered structure appeared tended to increase with an increase in the polydispersity index (PDI). On the other hand, the simulation results showed that the copolymers were likely to form ordered concentric strip structures when they were confined within a circle wall.

Self-assembly of diblock copolymer mixtures in confined states: A Monte Carlo study

The self-assembly of diblock copolymer mixtures (A-b-B/A-b-C or A-b-B/B-b-C mixtures) subjected to cylindrical confinement (two-dimensional confinement) was investigated using a Monte Carlo method. In this study, the boundary surfaces were configured to attract blocks A but repel blocks B and C. Relative to the structures of the individual components, the self-assembled structures of mixtures of the diblock copolymers were more complex and interesting. Under cylindrical confinement, with varying cylinder diameters and interaction energies between the boundary surfaces and the blocks, we observed a variety of interesting morphologies. Upon decreasing the cylinder's diameter, the self-assembled structures of the A(15)B(15)/A(15)C(15) mixtures changed from double-helix/cylinder structures (blocks B and C formed double helices, whereas blocks A formed the outer barrel and inner core) to stacked disk/cylinder structures (blocks B and C formed the stacked disk core, blocks A formed the outer cylindrical barrel), whereas the self-assembled structures of the A(15)B(7)/B7C15 mixtures changed from concentric cylindrical barrel structures to screw/cylinder structures (blocks C formed an inside core winding with helical stripes, whereas blocks A and B formed the outer cylindrical barrels) and then finally to the stacked disk/cylinder structures.

Grafting of maleic anhydride onto linear polyethylene: A Monte Carlo study

Monte Carlo simulation was used to study the graft of maleic anhydride (MAH) onto linear polyethylene (PE-g-MAH) initiated by dicumyl peroxide (DCP). Simulation results revealed that major MAH monomers attached onto PE chains as branched graft at higher MAH content. However, at extremely low MAH content, the fraction of bridged graft was very close to that of branched graft. This conclusion was somewhat different from the conventional viewpoint, namely, the fraction of bridged graft was always much lower than that of branched graft under any condition. Moreover, the results indicated that the grafting degree increased almost linearly to MAH and DCP concentrations. On the other hand, it was found that the amount of grafted MAH dropped sharply with increasing the length of grafted MAH, indicating that MAH monomers were mainly attached onto the PE chain as single MAH groups or very short oligomers. With respect to the crosslink of PE, the results showed that the fraction of PE-(MAH)(n)-PE crosslink structure increased continuously, and hence the fraction of PE-PE crosslink decreased with increasing MAH concentration.