Monte-Carlo technique

Ferrenberg A M and Swendsen R H 1988 New Monte-Carlo technique for studying phase-transitions Phys. Rev.L 61 2635-8... 

Senderowltz H, F Guarnieri and W C Still 1995. A Smart Monte Carlo Technique for Free Enerj Simulations of Multicanonical Molecules. Direct Calculations of the Conformational Populatioi of Organic Molecules. Journal of the American Chemical Society 117 8211-8219. 

In analytical chemistry, a number of identical measurements are taken and then an error is estimated by computing the standard deviation. With computational experiments, repeating the same step should always give exactly the same result, with the exception of Monte Carlo techniques. An error is estimated by comparing a number of similar computations to the experimental answers or much more rigorous computations. 

Another method of simulating chemical reactions is to separate the reaction and particle displacement steps. This kind of algorithm has been considered in Refs. 90, 153-156. In particular. Smith and Triska [153] have initiated a new route to simulate chemical equilibria in bulk systems. Their method, being in fact a generalization of the Gibbs ensemble Monte Carlo technique [157], has also been used to study chemical reactions at solid surfaces [90]. However, due to space limitations of the chapter, we have decided not to present these results. 

Within this context, the following sections are devoted to the description of the state of the art in the modeling and simulation of surface chemical reactions of simple systems using Monte Carlo techniques. 

These equations identify the dominant source and loss processes for HO and H02 when NMHC reactions are unimportant. Imprecisions inherent in the laboratory measured rate coefficients used in atmospheric mechanisms (for instance, the rate constants in Equation E6) can, themselves, add considerable uncertainty to computed concentrations of atmospheric constituents. A Monte-Carlo technique was used to propagate rate coefficient uncertainties to calculated concentrations (179,180). For hydroxyl radical, uncertainties in published rate constants propagate to modelled [HO ] uncertainties that range from 25% under low-latitude marine conditions to 72% under urban mid-latitude conditions. A large part of this uncertainty is due to the uncertainty (la=40%) in the photolysis rate of 0(3) to form O D, /j. 

Figure 1.24. Rejection of suspected outliers. A series of normally distributed values was generated by the Monte Carlo technique the mean and the standard deviation were calculated the largest normalized absolute deviate (residual) z = xi - /.i is plotted versus n (black...

In the following, an example from Chapter 4 will be used to demonstrate the concept of statistical ruggedness, by applying the chosen fitting model to data purposely corrupted by the Monte Carlo technique. The data are normalized TLC peak heights from densitometer scans. (See Section 4.2) ... 

Figure 3.11. Smoothing a noisy signal. The synthetic, noise-free signal is given at the top. After the addition of noise by means of the Monte Carlo technique, the panels in the second row are obtained (little noise, left, five times as much noise, right). A seven-point Savitzky-Golay filter of order 2 (third row) and a seven-point moving average (bottom row) filter are...

SIMILAR (SIMulate statistically simILAR data sets) Section 3.5 Monte Carlo Technique... 

Generation of random numbers belonging to the ND(0, 1) (Monte Carlo Technique). 

Giiell, O. A., and Holcombe, J. A., Analytical Applications of Monte Carlo Techniques, ArtflZ. Chem. 62, 1990, 529A-542A. 

We did not extensively discuss the consequences of lateral interactions of surface species adsorbed in adsorption overlayers. They lead to changes in the effective activation energies mainly because of consequences to the interaction energies in coadsorbed pretransition states. At lower temperatures, it can also lead to surface overlayer pattern formation due to phase separation. Such effects cannot be captured by mean-field statistical methods such as the microkinetics approaches but require treatment by dynamic Monte Carlo techniques as discussed in [25]. 

The approach to the evaluation of vibrational spectra described above is based on classical simulations for which quantum corrections are possible. The incorporation of quantum effects directly in simulations of large molecular systems is one of the most challenging areas in theoretical chemistry today. The development of quantum simulation methods is particularly important in the area of molecular spectroscopy for which quantum effects can be important and where the goal is to use simulations to help understand the structural and dynamical origins of changes in spectral lineshapes with environmental variables such as the temperature. The direct evaluation of quantum time- correlation functions for anharmonic systems is extremely difficult. Our initial approach to the evaluation of finite temperature anharmonic effects on vibrational lineshapes is derived from the fact that the moments of the vibrational lineshape spectrum can be expressed as functions of expectation values of positional and momentum operators. These expectation values can be evaluated using extremely efficient quantum Monte-Carlo techniques. The main points are summarized below. 

Such Bayesian models could be couched in terms of parametric distributions, but the mathematics for real problems becomes intractable, so discrete distributions, estimated with the aid of computers, are used instead. The calculation of probability of outcomes from assumptions (inference) can be performed through exhaustive multiplication of conditional probabilities, or with large problems estimates can be obtained through stochastic methods (Monte Carlo techniques) that sample over possible futures. 

Figure 22.4 Monte Carlo techniques were used to simulate different hypothetical individuals for different instances of the trial design, using variability and uncertainty distributions from the model analysis. The result is a collection of predicted outcomes, shown as a binned histogram (top figure). Success was defined as a difference in end point measurement of X or smaller between drug and comparator. Likelihood of success (shown in the bottom figure as a cumulative probability) for this example (low/medium drug dose and high comparator dose) is seen to be low, about 33%.

Rowe RC, York P, Colbourn EA, Roskilly S. The influence of pellet shape, size and distribution on capsule filling-a preliminary evaluation of three dimensional computer simulation using a Monte Carlo technique. Int J Pharm 2005 300 32-7. 

The mechanism of the NO -1- CO reaction at realistic pressures is thus very complicated. In addition to the reaction steps considered above, one also has to take into account that intermediates on the surface may organize into islands or periodically ordered structures. Monte Carlo techniques are needed to account for these effects. Consequently, we are still far from a complete kinetic description of the CO -1- NO reaction. For an interesting review of the mechanism and kinetics of this reaction we refer to Zhdanov and Kasemo [V.P. Zhdanov and B. Kasemo, Suif. Sci. Rep. 29 (1997) 31],... 

Tunon et al.194 studied the water molecule in liquid water. The sample of conformations by the microscopic environment (water in this case) was obtained using Monte Carlo technique. The energy was calculated as in the approach of Stanton et al.189 i.e., using Eqs. 4.25 and 4.26. The solvent induced increase of the dipole moment amounted to 0.61 Debye in line with the results by Wei and Salahub and close to the experimental value of 0.75 Debye. The solvation enthalpy amounted —12.6 kcal/mol, while the value calculated by Salahub and Wei and the experimental ones were —10.4 kcal/mol and —9.9 kcal/mol, respectively. 

Almarza, N. G. Lomba, E., Determination of the interaction potential from the pair distribution function an inverse Monte Carlo technique, Phys. Rev. E 2003, 68, 011202... 

As it will be explained in section 6, the usual way to evaluate the potential energy of a system simulated by Monte Carlo techniques, makes use of the pair potential approximation (although, as it will also be reviewed, several works have already appeared where nonadditivity corrections to the interaction potential have been included). In the pair potential approximation only two body interactions are taken into account. We will briefly explain here how to apply this approximation for the calculation of the potential energy, to the periodic system just described. The interaction potential energy under the pair potential approximation can be written as ... 

---