Simple algorithm for isothermal-isobaric molecular dynamics

We review principles of non-Hamiltonian statistical mechanics and present a new set of equations and integration algorithm for isothermal-isobaric dynamics. The chief advantage of the present scheme is that it is somewhat simpler than previous methods. We perform numerical simulations to test the accuracy of the algorithm and compare its stability to that of a "gold standard," a symplectic integrator for Hamiltonian dynamics of the same system. The stability of the isothermal-isobaric algorithm is comparable to the stability of the symplectic integrator.     

Fine structures of zeolite-Linde-L (LTL): surface structures, growth unit and defects

High-resolution electron microscopy (HREM) has been used to image the surface structure of nano- and micrometer-sized synthetic crystals of zeolite-Linde-L (LTL). Columnar holes and rotational, nano-sized, wheel-like defects were observed within the crystals, where the hole has a minimum size equal to that of the rotational defect. Predictions of surface structure from atomistic computer simulation concur with the observations from HREM and provide insight into the crystal growth mechanism of perfect and defective LTL. Analysis of the energetics of the formation of rotational defect structures reveals that the driving force for defect creation is thermodynamic and furthermore, the rotational defects could be created in high concentrations. Formation of a columnar hole is found to be slightly energetically unfavourable and therefore we speculate that the incidence of both rotational and nano-sized vacancy defects is strongly dependent on kinetic factors and reaction conditions. The morphology of nano- and microcrystalline LTL is contradistinct and we use insights from simulation to propose an explanation of the disparity in crystal shape.     

Subtleties in data analysis related to the size of critical region

We comment on the analysis of the critical behavior of a layered driven diffusive system recently done by Achahbar and Marro. We discuss why we believe their method of taking the thermodynamic limit and determining the order-parameter exponent leads to unreliable estimates.     

Formation and morphology of latex monolayers. Computer simulation studies

The results of computer simulations of monolayers created from monodisperse latex particles are presented and discussed. Layers are characterized by the normalized coverage,NC (the actual coverage of the surface related to its maximum possible coverage with particles), and by the average number of neighbors,ANN, calculated as the number of particles being in contact with a given one and averaged over all the particles on the surface. Variable parameters used in simulations include: the rate of particles deposition, the probability of lateral movements, the probability of desorption of particles adsorbed on the surface, the probability of covalent immobilization of adsorbed particles, and the on-sphere slip parameter, OSS (characterizing the scattering of a falling particle on the particles being already attached to the surface). Morphology of monolayers is qualitatively characterized by relations betweenANN andNC. It is shown that for all monolayers formed without adhesion (and without repulsion) between the particles adsorbed on the surface the dependence ofANN vs.NC is described by a characteristic master curve (regardless of the values of probabilities of desorption and lateral movements of particles). For the monolayers created including adhesive forces between the adsorbed particles the plots ofANN vs.NC lie above the master curve, while similar plots obtained for the layers made of particles showing various types of repulsive interactions are always placed below it. Thus, the dependencies ofANN vs.NC, derived from computer simulations, can be used for the determination of the character of the interparticle interactions in the real systems.     

Renormalization group for Markov chains and application to metastability

In this paper we introduce a new renormalization group method for the study of the long-time behavior of Markov chains with finite state space and with transition probabilities exponentially small in an external parameter. A general approach of metastability problems emerges from this analysis and is discussed in detail in the case of a two-dimensional Ising system at low temperature,     

Atomistic Simulation Study of Calcium,Phosphate and Fluoride Ion Association to the Teleopeptide‐Tails of Collagen – Initial Steps to Biomineral Formation

The attachment of single ions to putative adsorption sites in the tails of collagen fibers is investigated by means of molecular dynamics simulations and discussed with respect to the very early steps of apatite/collagen biomineral formation. Our studies clearly demonstrate an increased flexibility of the tails of the triple‐helical collagen protein. Apart from the termini of the backbone, several side chains were also observed to be freely accessible to ion attachment from aqueous solution. The teleopeptide was systematically scanned for suitable adsorption sites for calcium, phosphate and fluoride ions. Association of these ions was then explored from potential of mean force calculations. The resulting energy profiles reveal a variety of favorable protein‐ion bonds and hint at the suitability of the collagen tails to promote apatite aggregation.     

Computer experiments on crystalline nylons: structural analysis of nylons with large aliphatic segments

A theoretical study based on force-field calculations has been performed to investigate the structural preferences of crystalline even nylons n with large and very large aliphatic segments. Atomistic energy calculations and Monte Carlo simulations were carried out considering the conventional and forms of nylons 10, 12, 18, 24, and 32. Results indicated that the form is the most favored for nylons 10, 12, 18, and 24. However, the structure was unstable for nylon 32, a polymer in which the density of hydrogen bonds is almost negligible. In this case, the arrangement is energetically more favored than the one.     

Electrophoresis simulations using Chebyshev pseudo-spectral method on a moving mesh

We present the implementation and demonstration of the Chebyshev pseudo-spectral method coupled with an adaptive mesh method for performing fast and highly accurate electrophoresis simulations. The Chebyshev pseudo-spectral method offers higher numerical accuracy than all other finite difference methods and is applicable for simulating all electrophoresis techniques in channels with open or closed boundaries. To improve the computational efficiency, we use a novel moving mesh scheme that clusters the grid points in the regions with poor numerical resolution. We demonstrate the application of the Chebyshev pseudo-spectral method on a moving mesh for simulating nonlinear electrophoretic processes through examples of isotachophoresis (ITP), isoelectric focusing (IEF), and electromigration-dispersion in capillary zone electrophoresis (CZE) at current densities as high as 1000 A/m. We also show the efficacy of our moving mesh method over existing methods that cluster the grid points in the regions with large concentration gradients. We have integrated the adaptive Chebyshev pseudo-spectral method in the open-source SPYCE simulator and verified its implementation with other electrophoresis simulators.     

Rigid-body dynamics in the isothermal-isobaric ensemble: a test on the accuracy and computational efficiency

We have developed a time-reversible rigid-body (rRB) molecular dynamics algorithm in the isothermal-isobaric (NPT) ensemble. The algorithm is an extension of rigid-body dynamics [Matubayasi and Nakahara, J Chem Phys 1999, 110, 3291] to the NPT ensemble on the basis of non-Hamiltonian statistical mechanics [Martyna, G. J. et al., J Chem Phys 1994, 101, 4177]. A series of MD simulations of water as well as fully hydrated lipid bilayer systems have been undertaken to investigate the accuracy and efficiency of the algorithm. The rRB algorithm was shown to be superior to the state-of-the-art constraint-dynamics algorithm SHAKE/RATTLE/ROLL, with respect to computational efficiency. However, it was revealed that both algorithms produced accurate trajectories of molecules in the NPT as well as NVT ensembles, as long as a reasonably short time step was used. A couple of multiple time-step (MTS) integration schemes were also examined. The advantage of the rRB algorithm for computational efficiency increased when the MD simulation was carried out using MTS on parallel processing computer systems; total computer time for MTS-MD of a lipid bilayer using 64 processors was reduced by about 40% using rRB instead of SHAKE/RATTLE/ROLL.