Optimization of tried wave-function in quantum monte-carlo method

A method for optimizing tried wave functions in quantum Monte Carlo method has been found and used to calculate the energies of molecules, such as H-2, Li-2, H-3+, H-3 and H-4. Good results were obtained.

Development of a Cell Size Relaxed Scheme for the Direct Simulation Monte Carlo Method

The conventional direct simulation Monte Carlo (DSMC) method has a strong restriction on the cell size because simulated particles are selected randomly within the cell for collisions. Cells with size larger than the molecular mean free path are generally not allowed in correct DSMC simulations. However, the cell-size induced numerical error can be controlled if the gradients of flow properties are properly involved during collisions. In this study, a large cell DSMC scheme is proposed to relax the cell size restriction. The scheme is applied to simulate several test problems and promising results are obtained even when the cell size is greater than 10 mean free paths of gas molecules. However, it is still necessary, of course, that the cell size be small with respect to the flow field structures that must be resolved.

Improved sampling techniques for the direct simulation Monte Carlo method

The direct simulation Monte Carlo (DSMC) method is a widely used approach for flow simulations having rarefied or nonequilibrium effects. It involves heavily to sample instantaneous values from prescribed distributions using random numbers. In this note, we briefly review the sampling techniques typically employed in the DSMC method and present two techniques to speedup related sampling processes. One technique is very efficient for sampling geometric locations of new particles and the other is useful for the Larsen-Borgnakke energy distribution.

A markov chain Monte Carlo method for estimating population mixing using Y-chromosome markers: Mixing of the Han people in China

We present a new approach for estimating mixing between populations based on non-recombining markers, specifically Y-chromosome microsatellites. A Markov chain Monte Carlo (MCMC) Bayesian statistical approach is used to calculate the posterior probability

Monte-Carlo simulation of surface crack growth rate for offshore structural steel E36-Z35

The Monte- Carlo method is used to simulate the surface fatigue crack growth rate for offshore structural steel E36-Z35, and to determine the distributions and relevance of the parameters in the Paris equation. By this method, the time and cost of fatigue crack propagation testing can be reduced. The application of the method is demonstrated by use of four sets of fatigue crack propagation data for offshore structural steel E36-Z35. A comparison of the test data with the theoretical prediction for surface crack growth rate shows the application of the simulation method to the fatigue crack propagation tests is successful.

Permeability of Fractal Porous Media by Monte Carlo Simulations

The permeability of the fractal porous media is simulated by Monte Carlo technique in this work. Based oil the fractal character of pore size distribution in porous media, the probability models for pore diameter and for permeability are derived. Taking the bi-dispersed fractal porous media as examples, the permeability calculations are performed by the present Monte Carlo method. The results show that the present simulations present a good agreement compared with the existing fractal analytical solution in the general interested porosity range. The proposed simulation method may have the potential in prediction of other transport properties (such as thermal conductivity, dispersion conductivity and electrical conductivity) in fractal porous media, both saturated and unsaturated.

Pressure and thermal effects on elastic properties of single crystal alpha-Al2O3 from Monte-Carlo simulation

A rectangular structural unit cell of a-Al2O3 is generated from its hexagonal one. For the rectangular structural crystal with a simple interatomic potential [Matsui, Mineral Mag. 58A, 571 (1994)], the relations of lattice constants to homogeneous pressure and temperature are calculated by using Monte-Carlo method at temperature 298K and 0 GPa, respectively. Both numerical results agree with experimental ones fairly well. By comparing pair distribution function, the crystal structure of a-Al2O3 has no phase transition in the range of systematic parameters. Based on the potential model, pressure dependence of isothermal bulk moduli is predicted. Under variation of general strains, which include of external and internal strains, elastic constants of a-Al2O3 in the different homogeneous load are determined. Along with increase of pressure, axial elastic constants increase appreciably, but nonaxial elastic constants are slowly changed.

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.

Monte Carlo Simulation Of Thermal Fluctuations Below The Onset Of Rayleigh-Benard Convection

The density fluctuations below the onset of convection in the Rayleigh-Benard problem are studied with the direct simulation Monte Carlo method. The particle simulation results clearly show the connection between the static correlation functions of fluctuations below the critical Rayleigh number and the flow patterns above the onset of convection for small Knudsen number flows (Kn=0.01 and Kn=0.005). Furthermore, the physical nature for no convection in the Rayleigh-Benard problem under large Knudsen number conditions (Kn>0.028) is explained based on the dynamics of fluctuations.

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.

Monte Carlo simulation of the compatibility of graft copolymer compatibilized two incompatible homopolymer blends: Effect of graft structure

Compatibility of graft copolymer compatibilized two incompatible homopolymer A and B blends was simulated by using Monte Carlo method in a two-dimensional lattice model. The copolymers with various graft structures were introduced in order to study the effect of graft structure on the compatibility. Simulation results showed that incorporation of both A-g-B (A was backbone) and B-g-A (B was backbone) copolymers could much improve the compatibility of the blends. However, A-g-B copolymer was more effective to compatibilize the blend if homopolymer A formed dispersed phase. Furthermore, simulation results indicated that A-g-B copolymers tended to locate at the interface and anchor two immiscible components when the side chain is relatively long. However, most of A-g-B copolymers were likely to be dispersed into the dispersed homopolymer A phase domains if the side chains were relatively short. On the other hand, B-g-A copolymers tended to be dispersed into the matrix formed by homopolymer B. Moreover, it was found that more and more B-g-A copolymers were likely to form thin layers at the phase interface with decreasing the length of side chain.