Molecular Dynamics Simulation of Sympathetic Crystallization of Molecular Ions

It is shown that the translational degrees of freedom of a large variety of molecules, from light diatomic to heavy organic ones, can be cooled sympathetically and brought to rest (crystallized) in a linear Paul trap. The method relies on endowing the molecules with an appropriate positive charge, storage in a linear radiofrequency trap, and sympathetic cooling. Two well--known atomic coolant species, ${}^9{\hbox{Be}}^+$ and ${}^{137}{\hbox{Ba}}^+$, are sufficient for cooling the molecular mass range from 2 to 20,000 amu. The large molecular charge required for simultaneous trapping of heavy molecules and of the coolant ions can easily be produced using electrospray ionization. Crystallized molecular ions offer vast opportunities for novel studies.Comment: Accepted for publication in Phys. Rev.

Molecular Dynamics of Yukawa System using the Fast Multipole Method

In order to perform the large-scale molecular dynamics simulation of the Yukawa system, a mathematical expression for molecular dynamics using the fast multipole method is described. The model simulations are also performed to test the performance of our implementation of the FMM

Molecular mode-coupling theory for supercooled liquids: Application to water

We present mode-coupling equations for the description of the slow dynamics observed in supercooled molecular liquids close to the glass transition. The mode-coupling theory (MCT) originally formulated to study the slow relaxation in simple atomic liquids, and then extended to the analysis of liquids composed by linear molecules, is here generalized to systems of arbitrarily shaped, rigid molecules. We compare the predictions of the theory for the $q$-vector dependence of the molecular nonergodicity parameters, calculated by solving numerically the molecular MCT equations in two different approximation schemes, with ``exact'' results calculated from a molecular dynamics simulation of supercooled water. The agreement between theory and simulation data supports the view that MCT succeeds in describing the dynamics of supercooled molecular liquids, even for network forming ones.Comment: 22 pages 4 figures Late

Molecular Dynamics Simulation of Vascular Network Formation

Endothelial cells are responsible for the formation of the capillary blood vessel network. We describe a system of endothelial cells by means of two-dimensional molecular dynamics simulations of point-like particles. Cells' motion is governed by the gradient of the concentration of a chemical substance that they produce (chemotaxis). The typical time of degradation of the chemical substance introduces a characteristic length in the system. We show that point-like model cells form network resembling structures tuned by this characteristic length, before collapsing altogether. Successively, we improve the non-realistic point-like model cells by introducing an isotropic strong repulsive force between them and a velocity dependent force mimicking the observed peculiarity of endothelial cells to preserve the direction of their motion (persistence). This more realistic model does not show a clear network formation. We ascribe this partial fault in reproducing the experiments to the static geometry of our model cells that, in reality, change their shapes by elongating toward neighboring cells.Comment: 10 pages, 3 figures, 2 of which composite with 8 pictures each. Accepted on J.Stat.Mech. (2009). Appeared at the poster session of StatPhys23, Genoa, Italy, July 13 (2007

Molecular Dynamics Simulation of Polymer-Metal Bonds

Molecular simulation is becoming a very powerful tool for studying dynamic phenomena in materials. The simulation yields information about interaction at length and time scales unattainable by experimental measurements and unpredictable by continuum theories. This is especially meaningful when referring to bonding between a polymer and a metal substrate. A very important characteristic of polymers is that their physical properties do not rely on the detailed chemical structure of the molecular chains but only on their flexibility, and accordingly they will be able to adopt different conformations. In this paper, a molecular simulation of the bonding between vinyl ester polymer and steel is presented. Four different polymers with increasing chain lengths have been studied. Atomic co-ordinates are adjusted in order to reduce the molecular energy. Conformational changes in the macromolecules have been followed to obtain the polymer pair correlation function. Radius of gyration and end-to-end distance distributions of the individual chains have been used as a quantitative measurement of their flexibility. There exists a correlation between flexibility of the molecular chains and the energy of adhesion between the polymer and the metal substrate. Close contacts between the two materials are established at certain points but every atom up to a certain distance from the interface contributes to the total value of the adhesion energy of the system

Color Molecular-Dynamics for High Density Matter

We propose a microscopic simulation for quark many-body system based on molecular dynamics. Using color confinement and one-gluon exchange potentials together with the meson exchange potentials between quarks, we construct nucleons and nuclear/quark matter. Statistical feature and the dynamical change between confinement and deconfinement phases are studied with this molecular dynamics simulation.Comment: 4 pages, 3 figure

Replica-exchange molecular dynamics simulation for supercooled liquids

We investigate to what extend the replica-exchange Monte Carlo method is able to equilibrate a simple liquid in its supercooled state. We find that this method does indeed allow to generate accurately the canonical distribution function even at low temperatures and that its efficiency is about 10-100 times higher than the usual canonical molecular dynamics simulation.Comment: 6 pages, 5 figure