Molecular dynamics simulation continuous methods

There are many variants of the predictor-corrector theme of these, we will only mention the algorithm used by Rahman in the first molecular dynamics simulations with continuous potentials [Rahman 1964]. In this method, the first step is to predict new positions as follows ... 

The rapid rise in computer speed over recent years has led to atom-based simulations of liquid crystals becoming an important new area of research. Molecular mechanics and Monte Carlo studies of isolated liquid crystal molecules are now routine. However, care must be taken to model properly the influence of a nematic mean field if information about molecular structure in a mesophase is required. The current state-of-the-art consists of studies of (in the order of) 100 molecules in the bulk, in contact with a surface, or in a bilayer in contact with a solvent. Current simulation times can extend to around 10 ns and are sufficient to observe the growth of mesophases from an isotropic liquid. The results from a number of studies look very promising, and a wealth of structural and dynamic data now exists for bulk phases, monolayers and bilayers. Continued development of force fields for liquid crystals will be particularly important in the next few years, and particular emphasis must be placed on the development of all-atom force fields that are able to reproduce liquid phase densities for small molecules. Without these it will be difficult to obtain accurate phase transition temperatures. It will also be necessary to extend atomistic models to several thousand molecules to remove major system size effects which are present in all current work. This will be greatly facilitated by modern parallel simulation methods that allow molecular dynamics simulations to be carried out in parallel on multi-processor systems [115]. 

The EAM approach appears to provide a formalism within which realistic potentials which describe atomic dynamics can be developed. It should also provide a method for realistically incorporating adsorbates into dynamics simulations. Both of these applications can be considered significant advances, and will help molecular dynamics simulations to continue to contribute to the understanding of technologically important processes. 

Continued Quantum and Molecular Mechanical Simulations, In this technique, a molecular dynamics simulation includes the treatment of some part of the system wilh a quantum mechanical technique. This approach. yMf.MM. is similar to programs that Use quantum mechanical methods to treat the n-systems of the structures in question separately from the sigma framework. The results are combined ai ihe end to render a slructure which is optimized and energy-refined in satisfy both self-consistent field (SCF) and force field energy convergence. 

This process is repeated until convergence in the solute charges is achieved. In general, only a few cycles, 4-5, of quantum calculation/molecular dynamics simulations are needed for convergence. However, it is convenient to continue the procedure for another 10-15 cycles. In this way, the final results can be obtained with the associated statistical error by averaging over the last 5-10 last cycles. The scheme of the method is shown in Figure 4.28. 

Note that many of the molecules produced have few internal polar fimctional groups to which ions may bind. Instead, it is more likely that ion-water-channel interactions escort the ion through the pore. To that end. many of the models can then be viewed as methods to pull water into the lipidic core of a bilayer membrane and thereby stabilize ions in transport. Recent studies of molecular dynamics simulations of ion transportation in human aquaporin-1 and in the bacterial glycerol facilitator GlpF revealed the key role of water in the stabilization of ions in transit and in the molecular selectivity of channels. Synthetic compoxmds form less-defined stmctures than these complex proteins but apparently act as efficiently as more complex natural materials. It is likely that continued study of synthetic systems will continue to reveal the general details underlying all transport processes. 

In Monte-Carlo simulations, the energy of the molecular system is minimized by randomly moving the molecules in accordance with a desired probability distribution. After each move, the energy of each molecule is computed. When the total energy is reduced, the move is accepted and the molecules are redistributed. Moves are continued until equilibrium is achieved. As for molecular dynamics simulations, potential functions are provided. After convergence, the thermophysical properties, at equilibrium, are computed by averaging. Monte-Carlo methods, which are particularly effective for the calculation of thermophysical properties, including phase equilibria, are considered in detail by Rowley (1994). 

As discussed in the previous section a brute-force molecular-dynamics simulation of gas-surface dynamics, although simple in principle, is a large computational problem. Though these direct methods will continue to be of use, particularly in providing numerical benchmarks for the calibration of more approximate methods, it will prove useful to search for more efficient methods. The chief defect of direct methods when applied to gas-solid scattering is that the essentially harmonic chatracter of the lattice is not fully exploited. We expect that the strong, direct interaction with the solid will involve a relatively small number of lattice... 

Molecular Dynamics Simulation Method Multiscale Modeling and Numerical Simulations Non-Continuous Approaches Pressure Driven Single Phase Gas Flows Spectral Methods... 

All the methods presented later assume matter as a continuous medium, which is valid down to dimensions of the order of a few nanometers. For smaller sizes, the discrete nature of matter can be accounted for using molecular dynamics simulations. This issue is discussed in section 9.3.1.3. 

Osguthorpe [141] employed a continuous model and molecular dynamics simulated annealing. In spite of the use of a quite detailed model (main chain united atoms and up to three united atoms per residue), its very flexible chain geometry enabled efficient sampling. The potentials were derived from the statistics of known protein stmctures. The method enabled us to obtain correct predictions of substantial fractions of the stmcture of the attempted targets, and for one of the difficult targets, the prediction resulting from this method was the most accurate. 

As it is assured that the conformations generated with the Marcelja field are consistent with experimental data, the procedure employing the overlay method and the continuous Marcelja model is a powerful tool for preparation of the initial coordinates of biomembrane simulations. By applying a short molecular dynamics simulation to the initial coordinates prepared by this procedure, the long equilibration time necessary to obtain results consistent with experiments is well circumvented. 

---