Three parallel algorithms for classical molecular dynamics are presented. The first assigns each processor a fixed subset of atoms; the second assigns each a fixed subset of inter-atomic forces to compute; the third assigns each a fixed spatial region. The algorithms are suitable for molecular dynamics models which can be difficult to parallelize efficiently—those with short-range forces where the neighbors of each atom change rapidly. They can be implemented on any distributed-memory parallel machine which allows for message-passing of data between independently executing processors. The algorithms are tested on a standard Lennard-Jones benchmark problem for system sizes ranging from 500 to 100,000,000 atoms on several parallel supercomputers--the nCUBE 2, Intel iPSC/860 and Paragon, and Cray T3D. Comparing the results to the fastest reported vectorized Cray Y-MP and C90 algorithm shows that the current generation of parallel machines is competitive with conventional vector supercomputers even for small problems. For large problems, the spatial algorithm achieves parallel efficiencies of 90% and a 1840-node Intel Paragon performs up to 165 faster than a single Cray C9O processor. Trade-offs between the three algorithms and guidelines for adapting them to more complex molecular dynamics simulations are also discussed.

Fast parallel algorithms for short-range molecular dynamics

The stability of two-dimensional (2D) layers and membranes is subject of a long standing theoretical debate. According to the so called Mermin-Wagner theorem, long wavelength fluctuations destroy the long-range order for 2D crystals. Similarly, 2D membranes embedded in a 3D space have a tendency to be crumpled. These dangerous fluctuations can, however, be suppressed by anharmonic coupling between bending and stretching modes making that a two-dimensional membrane can exist but should present strong height fluctuations. The discovery of graphene, the first truly 2D crystal and the recent experimental observation of ripples in freely hanging graphene makes these issues especially important. Beside the academic interest, understanding the mechanisms of stability of graphene is crucial for understanding electronic transport in this material that is attracting so much interest for its unusual Dirac spectrum and electronic properties. Here we address the nature of these height fluctuations by means of straightforward atomistic Monte Carlo simulations based on a very accurate many-body interatomic potential for carbon. We find that ripples spontaneously appear due to thermal fluctuations with a size distribution peaked around 70 \AA which is compatible with experimental findings (50-100 \AA) but not with the current understanding of stability of flexible membranes. This unexpected result seems to be due to the multiplicity of chemical bonding in carbon.

Intrinsic ripples in graphene

A Review on the Application of Nonlocal Elastic Models in Modeling of Carbon Nanotubes and Graphenes

(CdSe)ZnS Core−Shell Quantum Dots: Synthesis and Characterization of a Size Series of Highly Luminescent Nanocrystallites

Bending of Euler–Bernoulli beams using Eringen’s integral formulation: A paradox resolved

Wave propagation in single-walled carbon nanotube under longitudinal magnetic field using nonlocal Euler–Bernoulli beam theory

Buckling analysis of micro-/nano-scale plates based on two-variable refined plate theory incorporating nonlocal scale effects

Nonlocal scale effects on wave propagation in multi-walled carbon nanotubes

Prediction of nonlocal scaling parameter for armchair and zigzag single-walled carbon nanotubes based on molecular structural mechanics, nonlocal elasticity and wave propagation

This article presents the buckling analysis of orthotropic nanoplates such as graphene using the two-variable refined plate theory and nonlocal small-scale effects. The two-variable refined plate theory takes account of transverse shear effects and parabolic distribution of the transverse shear strains through the thickness of the plate, hence it is unnecessary to use shear correction factors. Nonlocal governing equations of motion for the monolayer graphene are derived from the principle of virtual displacements. The closed-form solution for buckling load of a simply supported rectangular orthotropic nanoplate subjected to in-plane loading has been obtained by using the Navier’s method. Numerical results obtained by the present theory are compared with first-order shear deformation theory for various shear correction factors. It has been proven that the nondimensional buckling load of the orthotropic nanoplate is always smaller than that of the isotropic nanoplate. It is also shown that small-scale effects contribute significantly to the mechanical behavior of orthotropic graphene sheets and cannot be neglected. Further, buckling load decreases with the increase of the nonlocal scale parameter value. The effects of the mode number, compression ratio and aspect ratio on the buckling load of the orthotropic nanoplate are also captured and discussed in detail. The results presented in this work may provide useful guidance for design and development of orthotropic graphene based nanodevices that make use of the buckling properties of orthotropic nanoplates.

Saggam Narendar

Papers

Wave propagation in single-walled carbon nanotube under longitudinal magnetic field using nonlocal Euler–Bernoulli beam theory

Buckling analysis of micro-/nano-scale plates based on two-variable refined plate theory incorporating nonlocal scale effects

Nonlocal scale effects on wave propagation in multi-walled carbon nanotubes

Prediction of nonlocal scaling parameter for armchair and zigzag single-walled carbon nanotubes based on molecular structural mechanics, nonlocal elasticity and wave propagation

Scale effects on buckling analysis of orthotropic nanoplates based on nonlocal two-variable refined plate theory