Monte Carlo simulations diffusion

Monte Carlo simulations require less computer time to execute each iteration than a molecular dynamics simulation on the same system. However, Monte Carlo simulations are more limited in that they cannot yield time-dependent information, such as diffusion coefficients or viscosity. As with molecular dynamics, constant NVT simulations are most common, but constant NPT simulations are possible using a coordinate scaling step. Calculations that are not constant N can be constructed by including probabilities for particle creation and annihilation. These calculations present technical difficulties due to having very low probabilities for creation and annihilation, thus requiring very large collections of molecules and long simulation times. 

Molecular dynamics calculations are more time-consuming than Monte Carlo calculations. This is because energy derivatives must be computed and used to solve the equations of motion. Molecular dynamics simulations are capable of yielding all the same properties as are obtained from Monte Carlo calculations. The advantage of molecular dynamics is that it is capable of modeling time-dependent properties, which can not be computed with Monte Carlo simulations. This is how diffusion coefficients must be computed. It is also possible to use shearing boundaries in order to obtain a viscosity. Molec-... 

Other quantum simulations involve simulations with effective Hamiltonians [261-263] or the simulation of ground state wave properties by Green s function Monte Carlo or diffusion Monte Carlo for reviews and further references on these methods see Refs. 162, 264-268. 

Fig. 7 gives an example of such a comparison between a number of different polymer simulations and an experiment. The data contain a variety of Monte Carlo simulations employing different models, molecular dynamics simulations, as well as experimental results for polyethylene. Within the error bars this universal analysis of the diffusion constant is independent of the chemical species, be they simple computer models or real chemical materials. Thus, on this level, the simplified models are the most suitable models for investigating polymer materials. (For polymers with side branches or more complicated monomers, the situation is not that clear cut.) It also shows that the so-called entanglement length or entanglement molecular mass Mg is the universal scaling variable which allows one to compare different polymeric melts in order to interpret their viscoelastic behavior. 

Those Warren-Cowley parameters have been determined in situ above the order-disorder transition temperature by diffuse neutron scattering. From these experimentally determined static correlations, the first nine effective pair interactions have been deduced using inverse Monte Carlo simulations. 

Trinh, SH Arce, P Locke, BR, Effective Diffusivities of Point-Like Molecules in Isohopic Porous Media by Monte Carlo Simulation, Transport in Porous Media 38, 241, 2000. 

A quantitative analysis [34], based on the adsorption isotherms and the intercrystalline porosity, yielded the remarkable result that a satisfactory fit between the experimental data and the estimates of Aong-range = Pinter Anter following Eqs. (3.1.11) and (3.1.12) only lead to coinciding results for tortuosity factors a differing under the conditions of Knudsen diffusion (low temperatures) and bulk-diffusion (high temperatures) by a factor of at least 3. Similar results have recently been obtained by dynamic Monte Carlo simulations [39—41]. 

In general, percolation is one of the principal tools to analyze disordered media. It has been used extensively to study, for example, random electrical networks, diffusion in disordered media, or phase transitions. Percolation models usually require approximate solution methods such as Monte Carlo simulations, series expansions, and phenomenological renormalization [16]. While some exact results are known (for the Bethe lattice, for instance), they are very rare because of the complexity of the problem. Monte Carlo simulations are very versatile but lack the accuracy of the other methods. The above solution methods were employed in determining the critical exponents given in the following section. 

The independent reaction time (1RT) model was introduced as a shortcut Monte Carlo simulation of pairwise reaction times without explicit reference to diffusive trajectories (Clifford et al, 1982b). At first, the initial positions of the reactive species (any number and kind) are simulated by convolving from a given (usually gaussian) distribution using random numbers. These are examined for immediate reaction—that is, whether any interparticle separation is within the respective reaction radius. If so, such particles are removed and the reactions are recorded as static reactions. 

In Sect. 7.4.6, we discussed various stochastic simulation techniques that include the kinetics of recombination and free-ion yield in multiple ion-pair spurs. No further details will be presented here, but the results will be compared with available experiments. In so doing, we should remember that in the more comprehensive Monte Carlo simulations of Bartczak and Hummel (1986,1987, 1993,1997) Hummel and Bartczak, (1988) the recombination reaction is taken to be fully diffusion-controlled and that the diffusive free path distribution is frequently assumed to be rectangular, consistent with the diffusion coefficient, instead of a more realistic distribution. While the latter assumption can be justified on the basis of the central limit theorem, which guarantees a gaussian distribution for a large number of scatterings, the first assumption is only valid for low-mobility liquids. 

A. Elements of Monte Carlo Simulation of Turbulent Diffusion. 288... 

Clearly this is a very interesting problem and of great practical relevance, very well suited to Monte Carlo simulation. At the same time, simulations of such problems have just only begun. In the context of crystal growth kinetics, models where evaporation-condensation processes compete with surface diffusion processes have occasionally been considered before . But many related processes can be envisaged which have not yet been studied at all. 

---