Monte Carlo simulation of hydrogen adsorption on Ni surfaces

In the present paper the adsorption kinetics of the hydrogen molecule on the (111) and (100) surfaces have been studied with the model proposed by Panczyk and the grand canonical Monte Carlo simulation method. The equilibrium adsorption isotherms are calculated at five different temperatures ranging from 314 K to 376 K and compared with the experimental equilibrium adsorption isotherms. The effects of temperature and pressure on coverage are also analyzed.      

Monte Carlo simulation of mictomagnets

Monte Carlo simulations of a binary alloy with impurity concentrations between 20 and 45 at.% have been carried out. The proportion of large clusters relative to that of small clusters increases with the number of MC diffusion steps as well as impurity concentration. Magnetic susceptibility peaks become more prominent and occur at higher temperatures with increasing impurity concentration. The different peaks in the susceptibility and specific heat curves seem to correspond to different sized clusters. A freezing model would explain the observed behaviour with the large clusters freezing first and the small clusters contributing to susceptibility (specific heat) peaks at lower temperatures.Contribution No. 153 from the Solid State and Structural Chemistry Unit     

Quantum Monte Carlo simulation of the high-pressure molecular-atomic crossover in fluid hydrogen

A first-order liquid-liquid phase transition in high-pressure hydrogen between molecular and atomic fluid phases has been predicted in computer simulations using ab initio molecular dynamics approaches. However, experiments indicate that molecular dissociation may occur through a continuous crossover rather than a first-order transition. Here we study the nature of molecular dissociation in fluid hydrogen using an alternative simulation technique in which electronic correlation is computed within quantum Monte Carlo methods, the so-called coupled electron-ion Monte Carlo method. We find no evidence for a first-order liquid-liquid phase transition.     

A Monte Carlo simulation of coagulation

A Monte Carlo simulation technique is described for the study of the coagulation of suspended particles. The method is computationally efficient since the particle trajectories are not used to determine coagulations. Instead, pairs of particles are assigned probabilities to coagulate and the evolution is computed as a stochastic Markov game. We also describe a simple analytic method to obtain the stationary distribution of sizes for the various mechanisms of relative particle motion. It is demonstrated that the simulation yields the correct stationary size distribution independent of initial condition.     

Monte Carlo simulation of aorta autofluorescence

Results of numerical simulation of autofluorescence of the aorta by the method of Monte Carlo are reported. Two states of the aorta, normal and with atherosclerotic lesions, are studied. A model of the studied tissue is developed on the basis of information about optical, morphological, and physico-chemical properties. It is shown that the data obtained by numerical Monte Carlo simulation are in good agreement with experimental results indicating adequacy of the developed model of the aorta autofluorescence.     

Monte Carlo simulation of nanometric cutting

Nanometric cutting of single-crystal materials at conventional cutting speeds (5?m?s^?1) is simulated for the first time using a new Monte Carlo method that is applicable to systems that are neither canonical nor microcanonical. This is accomplished by defining a local temperature in the cutting zone using the thermal analysis developed by Komanduri and Hou for conventional machining. Extension of this method to the nanometric regime permits an accurate estimate of the local temperature in cutting. This temperature is then employed in the Boltzmann probability distribution function that is used to determine the acceptance–rejection of Monte Carlo moves in the simulation. Since cutting speed is closely related to cutting temperature, the cutting speed enters the calculation via the thermal analysis equations. The method is applied to nanometric cutting of single-crystal aluminium with the crystal oriented in the (001) plane and cut in the [100] direction. Three positive rake cutting tools, namely 10°, 30° and 45°, are employed to investigate the effect of the rake angle on the forces, the specific energy and the nature of the chip formation. The method is evaluated by direct comparison with corresponding molecular dynamics simulations conducted under the same conditions.     

Monte Carlo simulation of liquid neopentane

Monte Carlo simulations are reported for liquid neopentane at 25°C and 0·592 g cm^-3 using three different potential models. These results are compared with X-ray diffraction data and simulations for other tetrahedral molecules.     

Monte Carlo simulation of hydrogen physisorption in K-doped single walled carbon nanotube array

Properties of hydrogen physisorption in K-doped single walled carbon nanotube array (SWCNTA) are investigated in detail by grand canonical Monte Carlo simulation. The optimization of hydrogen storage capacity at 293 K and 10 MPa as a function of K-doping schemes, K atoms’ doped-sites, and SWCNTA configuration is discussed.     

Multilayer adsorption by Monte Carlo simulation

Adsorption phenomena are characterized by models that include free parameters trying to reproduce experimental results. In order to understand the relationship between the model parameters and the material properties, the adsorption of small molecules on a crystalline plane surface has been simulated using the bond fluctuation model. A direct comparison between the Guggenheim–Anderson–de Boer (GAB) model for multilayer adsorption and computer simulations allowed us to establish correlations between the adsorption model parameters and the simulated interaction potentials.     

Monte Carlo simulation of magnetic multi-core nanoparticles

In this paper, a Monte Carlo simulation is carried out to evaluate the equilibrium magnetization of magnetic multi-core nanoparticles in a liquid and subjected to a static magnetic field. The particles contain a magnetic multi-core consisting of a cluster of magnetic single-domains of magnetite. We show that the magnetization of multi-core nanoparticles cannot be fully described by a Langevin model. Inter-domain dipolar interactions and domain magnetic anisotropy contribute to decrease the magnetization of the particles, whereas the single-domain size distribution yields an increase in magnetization. Also, we show that the interactions affect the effective magnetic moment of the multi-core nanoparticles.     

Monte Carlo simulation of photoelectron angular distribution

A practical simulation method has been performed for studies of the influence of surface excitations on the angular distributions of photoelectron peak intensities. The surface effects have been incorporated into simulations by using the surface excitation parameters (SEPs) which have been calculated with the extended Drude dielectric function. Also, elastic scattering cross sections are calculated using the finite difference method for a Hartree-Fock-Wigner-Seitz potential in the Dirac equation to take into account the solid-state effect. Results of Monte Carlo simulations reveal that surface effects lead to a reduction of the intensities at small detection angles and a sharp decrease at large angles since the surface excitation is most probable for glancing electrons. The calculated results taking into account surface effects are in better agreement with the experimental data.     

Monte Carlo simulation of the enantioseparation process

By means of Monte Carlo simulation, a study of enantioseparation by capillary electrophoresis has been carried out. A simplified system consisting of two enantiomers S (R) and a selector chiral C, which reacts with the enantiomers to form complexes RC (SC), has been considered. The dependence of Δμ$\Delta \mu$ (enantioseparation) with the concentration of chiral selector and with temperature have been analyzed by simulation. The effect of the binding constant and the charge of the complexes are also analyzed. The results are qualitatively satisfactory, despite the simplicity of the model.     

Monte Carlo simulation of polymers at interfaces

Polymers at interfaces pose challenging problems to statistical physics because their configurations often differ greatly from the bulk. Computer simulation of coarse-grained models then gives valuable insight and allows stringent tests of various theoretical predictions. Three examples are briefly treated: chain configurations of B-chains in the surface-enriched B-rich layer of an (AB) binary polymer mixture; “frustrated” lamellar ordering in ultra-thin block-copolymer films; and the collapse of polymer brushes in bad solvents.