Simulation Requirements

Calculating lattice energy as discussed in the forgoing sections is embodied in, for example, the perfect lattice section of the lattice-defect program CAS-CADE 2) jjj the codes of the CERIUS and BIOSYM packages. It should 

We now describe two methods of static lattice simulation used to investigate the migration of a dopant. 

---

Abstract. Molecular dynamics (MD) simulations of proteins provide descriptions of atomic motions, which allow to relate observable properties of proteins to microscopic processes. Unfortunately, such MD simulations require an enormous amount of computer time and, therefore, are limited to time scales of nanoseconds. We describe first a fast multiple time step structure adapted multipole method (FA-MUSAMM) to speed up the evaluation of the computationally most demanding Coulomb interactions in solvated protein models, secondly an application of this method aiming at a microscopic understanding of single molecule atomic force microscopy experiments, and, thirdly, a new method to predict slow conformational motions at microsecond time scales. 

Long term simulations require structurally stable integrators. Symplec-tic and symmetric methods nearly perfectly reproduce structural properties of the QCMD equations, as, for example, the conservation of the total energy. We introduced an explicit symplectic method for the QCMD model — the Pickaback scheme— and a symmetric method based on multiple time stepping. 

One of the main advantages of the stochastic dynamics methods is that dramatic tirn savings can he achieved, which enables much longer stimulations to he performed. Fc example, Widmalm and Pastor performed 1 ns molecular dynamics and stochastic dynamic simulations of an ethylene glycol molecule in aqueous solution of the solute and 259 vvatc jnolecules [Widmalm and Pastor 1992]. The molecular dynamics simulation require 300 hours whereas the stochastic dynamics simulation of the solute alone required ju 24 minutes. The dramatic reduction in time for the stochastic dynamics calculation is du not only to the very much smaller number of molecules present hut also to the fact the longer time steps can often he used in stochastic dynamics simulations. 

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. 

Overall, WebLab MedChem Explorer is very easy to use. The stepwise job setup works well, assuming that all users will be following a conventional drug rehnement process. It is not a program that can be used for complex simulations requiring the researcher to manually control many details of the simulation. 

High temperature simulations require special consideration in choosing the sampling interval (see Step size on page 89). 

Eor evaluation of flocculants for pressure belt filters, both laboratory-scale filters and filter simulators are available (52,53) in many cases from the manufacturers of the full-scale equipment. The former can be mn either batchwise or continuously the simulators require less substrate and are mn batchwise. The observed parameters include cake moisture, free drainage, release of the cake from the filter cloth, filter blinding, and retention of the flocculated material during appHcation of pressure. 

Proper condensed phase simulations require that the non-bond interactions between different portions of the system under study be properly balanced. In biomolecular simulations this balance must occur between the solvent-solvent (e.g., water-water), solvent-solute (e.g., water-protein), and solute-solute (e.g., protein intramolecular) interactions [18,21]. Having such a balance is essential for proper partitioning of molecules or parts of molecules in different environments. For example, if the solvent-solute interaction of a glutamine side chain were overestimated, there would be a tendency for the side chain to move into and interact with the solvent. The first step in obtaining this balance is the treatment of the solvent-solvent interactions. The majority of biomolecular simulations are performed using the TIP3P [81] and SPC/E [82] water models. 

This kind of simulation requires massive computer power, and much of it is done on so-called supercomputers . This is a reason why much recent research of this kind has been done at Los Alamos. In a survey of research in the American national laboratories, the then director of the Los Alamos laboratory, Siegfried Hecker (1990) explains that the laboratory has worked closely with all supercomputer vendors over the years, typically receiving the serial No. I machine for each successive model . He goes on to exemplify the kinds of problems in materials science that these extremely powerful machines can handle. 

We can also estimate the required computational work. Since the evolution time of the LG is of order Ljv and the time required to update a single step on the lattice is of order Iq/cs, the LG simulation requires roughly L/[IqM) = TZ/M time steps. Assuming that the computational work for each site is of order 1, we have that the LG simulation therefore requires a memory storage space of order S (TZ/M) and a computational work requirement of order W 7 d+l/ d+2... 

Reactor start-up simulations require initial values of A., T(n,t) and T. be zero. Monomer... 

In the previous chapter we examined cellular automata simulations of first-order reactions. Because these reactions involved just transformations of individual ingredients, the simulations were relatively simple and straightforward to set up. Second-order cellular automata simulations require more instructions than do the first-order models described earlier. First of all, since movement is involved and ingredients can only move into vacant spaces on the grid, one must allow a suitable number of vacant cells on the grid for movement to take place in a sensible manner. For a gas-phase reaction one might wish to allow at least 5-10 vacant cells for each ingredient, so that on a 100 x 100 = 10,000... 

Developed to meet the simulation requirements of the European Space Agency used by such leading companies as British Gas, Lucas Aerospace, BNFL, British Aerospace. 

Implementation of Equation 9.18 in spectral simulators requires some extra precautions (Hagen 1981 Hagen et al. 1985d) (A) The increased periodicity now requires one half of the unit sphere to be scanned. (B) The fact that the term within the absolute-value bars in Equation 9.18 can change sign as a function of molecular orientation implies the possibility that for specific orientations the linewidth becomes equal to zero. To avoid a program crash due to a zero divide, e.g., in the expression for the lineshape in Equation 4.8, a residual linewidth W0 has to be introduced ... 

Line shape analysis was performed for the binding of some dihydroxycholate ions to /1-cyclodextrin.205 The dihydroxycholates show different 18-CH3 signals for the complexed and free dihydroxycholate ions. To extract exchange rate constants from the NMR spectra, a complete line-shape simulation was performed. The simulation requires input of the chemical shift difference between the two sites, the line width in the absence of exchange, the residence time in each site (thg and Tg), and the relative population (fHG and fG) of each site (Equation (11)). The values were varied until the simulated and experimental spectra could be superimposed. The dissociation rate... 

Effective computer codes for the optimization of plants using process simulators require accurate values for first-order partial derivatives. In equation-based codes, getting analytical derivatives is straightforward, but may be complicated and subject to error. Analytic differentiation ameliorates error but yields results that may involve excessive computation time. Finite-difference substitutes for analytical derivatives are simple for the user to implement, but also can involve excessive computation time. 

Since the statistical mechanics simulation requires a knowledge of the features of the solute along the reaction coordinate, the calculated geometric parameters and the energy were fitted into analytical functions of rc. However, the absolute magnitude of the variation of rCH was so small that this parameter was held constant in the simulations. 

It is worth noting than adsorption process design is a mature science and there is seldom a need to employ detailed numerical models to solve the pdes described in Section 9.4. There are some spedahzed circumstances that may fall outside the norms for which many of the today s design tools have been formulated. Only in these circumstances is a more rigorous process simulation required for obtaining a design. 

Before leaving this section we consider a slightly different optimization problem that may also be expensive to solve. In flowsheet optimization, the process simulator is based almost entirely on equilibrium concepts. Separation units are described by equilibrium stage models, and reactors are frequently represented by fixed conversion or equilibrium models. More complex reactor models usually need to be developed and added to the simulator by the engineer. Here the modular nature of the simulator requires the reactor model to be solved for every flowsheet pass, a potentially expensive calculation. For simulation, if the reactor is relatively insensitive to the flowsheet, a simpler model can often be substituted. For process optimization, a simpler, insensitive model will necessarily lead to suboptimal (or even infeasible) results. The reactor and flowsheet models must therefore be considered simultaneously in the optimization. 

---