Search

In this paper we present the results of a large-scale numerical investigation of structural properties of a model of cell membrane, simulated as a bilayer of flexible molecules in vacuum. The study was performed by carrying out extensive Molecular Dynamics simulations, in the (NVE) micro-canonical ensemble, of two systems of different sizes (2x32 and 2x256 molecules), over a fairly large set of temperatures and densities, using parallel platforms and more standard serial computers. Depending on the dimension of the system, the dynamics was followed for physical times that go from few hundred of picoseconds for the largest system to 5--10 nanoseconds for the smallest one. We find that the bilayer remains stable even in the absence of water and neglecting Coulomb interactions in the whole range of temperatures and densities we have investigated. The extension of the region of physical parameters that we have explored has allowed us to study significant points in the phase diagram of the bilayer and to expose marked structural changes as density and temperature are varied, which are interpreted as the system passing from a crystal to a gel phase., Comment: 41 pages, 13 figures

In this talk I discuss the general question of the portability of Molecular Dynamics codes for diffusive systems on parallel computers of the APE family. The intrinsic single precision arithmetics of the today available APE platforms does not seem to affect the numerical accuracy of the simulations, while the absence of integer addressing from CPU to individual nodes puts strong constraints on the possible programming strategies. Liquids can be very satisfactorily simulated using the "systolic" method. For more complex systems, like the biological ones at which we are ultimately interested in, the "domain decomposition" approach is best suited to beat the quadratic growth of the inter-molecular computational time with the number of elementary components of the system. The promising perspectives of using this strategy for extensive simulations of lipid bilayers are briefly reviewed., Comment: 4 pages LaTeX, 2 figures included, espcrc2.sty required