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.

A Generalized Soft-Sphere for Monte Carlo Simulation

A new collision model, called the generalized soft-sphere (GSS) model, is introduced. It has the same total cross section as the generalized hard-sphere model [Phys. Fluids A 5, 738 (1993)], whereas the deflection angle is calculated by the soft-sphere scattering model [Phys. Fluids A 3, 2459 (1991)]. In virtue of a two-term formula given to fit the numerical solutions of the collision integrals for the Lennard-Jones (6-12) potential and for the Stockmayer potential, the parameters involved in the GSS model are determined explicitly that may fully reproduce the transport coefficients under these potentials. Coefficients of viscosity, self-diffusion and diffusion for both polar and nonpolar molecules given by the GSS model and experiment are in excellent agreement over a wide range of temperature from low to high.

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方法计算氢分子H_2基态(Χ~1∑_g~+)势能曲线

本文用量子MontoCarlo方法中优化试探波函数Ψ_T计算氢分子H_2基态(X~1∑_g~+)势能曲线.文中采用相当简单的波函数形式,并用固定样点优化技术优化试探波函数的参数.确定优化试探波函数后,分别用变分Monte Carlo及固定节面M0nte Carlo计算势能曲线各点能值.二种方法先后得95%和100%的相关能.因此,在量子M0nte Carlo方法中,用本文作者提出的试探波函数计算分子势能面,将会获得很好的结果.从而对分子散射和动力学的研究有重要意义.

Quantum and variational monte-carlo interaction potentials for li2 (x1-sigma-g+)

Optimized trial functions are used in quantum Monte Carlo and variational Monte Carlo calculations of the Li₂(X ¹Σ⁺_g) potential curve. The trial functions used are a product of a Slater determinant of molecular orbitals multiplied by correlation functions of electron—nuclear and electron—electron separation. The parameters of the determinant and correlation functions are optimized simultaneously by reducing the deviations of the local energy E_L (E_L  Ψ⁻¹_THΨ_T, where Ψ_T denotes a trial function) over a fixed sample. At the equilibrium separation, the variational Monte Carlo and quantum Monte Carlo methods recover 68% and 98% of the correlation energy, respectively. At other points on the curves, these methods yield similar accuracies.

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.

Monte Carlo Modeling of Micro Scale Gas Flows

An information preservation (IP) method has been used to simulate many micro scale gas flows. It may efficiently reduce the statistical scatter inherent in conventional particle approaches such as the direct simulation Monte Carlo (DSMC) method. This paper reviews applications of IP to some benchmark problems. Comparison of the IP results with those given by experiment, DSMC, and the linearized Boltzmann equation, as well as the Navier-Stokes equations with a slip boundary condition, and the lattice Boltzmann equation, shows that the IP method is applicable to micro scale gas flows over the entire flow regime from continuum to free molecular.

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.

跨大气层飞行器姿态控制相关三维稀薄气体流动的Monte Carlo模拟

本文利用直接模拟Monte Carlo位置元(DSMC-PE)方案模拟了跨大气层飞行器姿态控制相关的两个典型的稀薄气体流动,真空轴对称射流冲击平板和楔形垂直相交平板稀薄气体绕流。真空轴对称射流在附近平板表面的压力、剪应力和热流分布的计算结果,与Legge(1991)和D(?)ring(1990)的测量数据一致;楔形垂直相交平板表面压力分布计算结果,与Allegre和Raffin(1992)实验数据的比较,也令人满意。上述计算表明,DSMC-PE是处理航天领域稀薄气流的有力工具。