Simulated evolution methods

After the screen, the proteins are ranked according to their measured value of the figure of merit. Typically, the top x percent of the sequences is kept for the next round of mutation. The parameter x is to be adjusted experimentally. In the simulated evolutions, the value of x = 10% was always found to be optimal. Other methods for selecting the proteins to keep for... 

Natural computation denotes the following methods simulated annealing, simulated evolution strategy, genetic algorithms, artificial... 

Equations for the evolution of the updated PDF (i.e., after assimilating the measurements) similar to the ones given in section Markov Vector Methods can be derived. These are called the Kushner-Stratonovich equations. For most practical problems of interest, with nonlinear process and measurement equations, non-Gaussian noises, this equation remains theoretical in nature and suffers from the moment closure problem (Maybeck 1982). Thus updating reliability models is mostly carried out through the simulation-based methods. 

The calculation of the time evolution operator in multidimensional systems is a fomiidable task and some results will be discussed in this section. An alternative approach is the calculation of semi-classical dynamics as demonstrated, among others, by Heller [86, 87 and 88], Marcus [89, 90], Taylor [91, 92], Metiu [93, 94] and coworkers (see also [83] as well as the review by Miller [95] for more general aspects of semiclassical dynamics). This method basically consists of replacing the 5-fimction distribution in the true classical calculation by a Gaussian distribution in coordinate space. It allows for a simulation of the vibrational... 

The method of molecular dynamics (MD), described earlier in this book, is a powerful approach for simulating the dynamics and predicting the rates of chemical reactions. In the MD approach most commonly used, the potential of interaction is specified between atoms participating in the reaction, and the time evolution of their positions is obtained by solving Hamilton s equations for the classical motions of the nuclei. Because MD simulations of etching reactions must include a significant number of atoms from the substrate as well as the gaseous etchant species, the calculations become computationally intensive, and the time scale of the simulation is limited to the... 

But a computer simulation is more than a few clever data structures. We need algorithms to manipulate our system. In some way, we have to invent ways to let the big computer in our hands do things with the model that is useful for our needs. There are a number of ways for such a time evolution of the system the most prominent is the Monte Carlo procedure that follows an appropriate random path through configuration space in order to investigate equilibrium properties. Then there is molecular dynamics, which follows classical mechanical trajectories. There is a variety of dissipative dynamical methods, such as Brownian dynamics. All these techniques operate on the fundamental degrees of freedom of what we define to be our model. This is the common feature of computer simulations as opposed to other numerical approaches. 

Based on the flame-hole dynamics [59], dynamic evolutions of flame holes were simulated to yield the statistical chance to determine the reacting or quenched flame surface under the randomly fluctuating 2D strain-rate field. The flame-hole d5mamics have also been applied to turbulent flame stabilization by considering the realistic turbulence effects by introducing fluctuating 2D strain-rate field [22] and adopting the level-set method [60]. 

Molecular dynamics (MD) permits the nature of contact formation, indentation, and adhesion to be examined on the nanometer scale. These are computer experiments in which the equations of motion of each constituent particle are considered. The evolution of the system of interacting particles can thus be tracked with high spatial and temporal resolution. As computer speeds increase, so do the number of constituent particles that can be considered within realistic time frames. To enable experimental comparison, many MD simulations take the form of a tip-substrate geometry correspoudiug to scauniug probe methods of iuvestigatiug siugle-asperity coutacts (see Sectiou III.A). 

The Monte Carlo method as described so far is useful to evaluate equilibrium properties but says nothing about the time evolution of the system. However, it is in some cases possible to construct a Monte Carlo algorithm that allows the simulated system to evolve like a physical system. This is the case when the dynamics can be described as thermally activated processes, such as adsorption, desorption, and diffusion. Since these processes are particularly well defined in the case of lattice models, these are particularly well suited for this approach. The foundations of dynamical Monte Carlo (DMC) or kinetic Monte Carlo (KMC) simulations have been discussed by Eichthom and Weinberg (1991) in terms of the theory of Poisson processes. The main idea is that the rate of each process that may eventually occur on the surface can be described by an equation of the Arrhenius type ... 

A complete model for the description of plasma deposition of a-Si H should include the kinetic properties of ion, electron, and neutral fluxes towards the substrate and walls. The particle-in-cell/Monte Carlo (PIC/MC) model is known to provide a suitable way to study the electron and ion kinetics. Essentially, the method consists in the simulation of a (limited) number of computer particles, each of which represents a large number of physical particles (ions and electrons). The movement of the particles is simply calculated from Newton s laws of motion. Within the PIC method the movement of the particles and the evolution of the electric field are followed in finite time steps. In each calculation cycle, first the forces on each particle due to the electric field are determined. Then the... 

The simulation of an electronegative gas discharge converges much more slowly than that of an electropositive discharge. This is mainly caused by tbe slow evolution of the negative-ion density, which depends only on attachment (to create negative ions) and ion-ion recombination (to annihilate negative ions), both processes with a very small cross section. In addition to the common procedures adopted in the literature [222, 223, 272, 273], such as the null collision method, and different superparticle sizes and time steps for different types of particle, two other procedures were used to speed up the calculation [224]. 

In addition to the fact that MPC dynamics is both simple and efficient to simulate, one of its main advantages is that the transport properties that characterize the behavior of the macroscopic laws may be computed. Furthermore, the macroscopic evolution equations can be derived from the full phase space Markov chain formulation. Such derivations have been carried out to obtain the full set of hydrodynamic equations for a one-component fluid [15, 18] and the reaction-diffusion equation for a reacting mixture [17]. In order to simplify the presentation and yet illustrate the methods that are used to carry out such derivations, we restrict our considerations to the simpler case of the derivation of the diffusion equation for a test particle in the fluid. The methods used to derive this equation and obtain the autocorrelation function expression for the diffusion coefficient are easily generalized to the full set of hydrodynamic equations. 

This is, in essence, the modern synthesis of Darwin and Mendel achieved in the 1930s by Ronald Fisher and J. B. S. Haldane. Based on a series of relatively straightforward equations, it also took the study of evolution out of meticulously observed natural history and located it within a more abstract mathematised theory. Indeed, evolution itself came to be defined not in terms of organisms and populations, but as the rate of change of gene frequencies within any given population. One consequence has been a tendency for theoretical evolutionists to retreat further and further into abstract hypotheticals based on computer simulations, and to withdraw from that patient observation of the natural world which so characterised Darwin s own method . 

Equations (8.44) and (8.45) guarantee convergence to a canonical distribution only in the case of fixed B. Because B varies (i.e., the method uses information from momenta sampled in the past in determining the vector B), the evolution is not strictly Markovian. As a consequence, the correlations introduced can lead to the accumulation of systematic errors in the determination of configuration averages [77], However, these correlations can be broken if the update of B is not done each step, but with a lower updating frequency. This is analogous to other approximately Markovian procedures employed in MC simulations (e.g., update of the maximum displacements allowed for individual atoms [78]). 

Most free energy and phase-equilibrium calculations by simulation up to the late 1980s were performed with the Widom test particle method [7]. The method is still appealing in its simplicity and generality - for example, it can be applied directly to MD calculations without disturbing the time evolution of a system. The potential distribution theorem on which the test particle method is based as well as its applications are discussed in Chap. 9. 

Sometimes the theoretical or computational approach to description of molecular structure, properties, and reactivity cannot be based on deterministic equations that can be solved by analytical or computational methods. The properties of a molecule or assembly of molecules may be known or describable only in a statistical sense. Molecules and assemblies of molecules exist in distributions of configuration, composition, momentum, and energy. Sometimes, this statistical character is best captured and studied by computer experiments molecular dynamics, Brownian dynamics, Stokesian dynamics, and Monte Carlo methods. Interaction potentials based on quantum mechanics, classical particle mechanics, continuum mechanics, or empiricism are specified and the evolution of the system is then followed in time by simulation of motions resulting from these direct... 

REAPDOR method to measure the C-Al distances by numerical simulations of the time-dependence of the REAPDOR evolution effect. An intemuclear C-Al distance of 3.1 A was determined for a signal at 57.3 ppm in ZSM-5, which is in excellent agreement with quantum chemical calculations of surface methoxy species [232]. Larger C-Al distances (>4 A) were determined for nearby signals at 60.0 and 61.7 ppm. 

---