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

The long-standing popularity of thermoelectric materials has contributed to the creation of various thermoelectric devices and stimulated the development of strategies to improve their thermoelectric performance. In this review, we aim to comprehensively summarize the state-of-the-art strategies for the realization of high-performance thermoelectric materials and devices by establishing the links between synthesis, structural characteristics, properties, underlying chemistry and physics, including structural design (point defects, dislocations, interfaces, inclusions, and pores), multidimensional design (quantum dots/wires, nanoparticles, nanowires, nano- or microbelts, few-layered nanosheets, nano- or microplates, thin films, single crystals, and polycrystalline bulks), and advanced device design (thermoelectric modules, miniature generators and coolers, and flexible thermoelectric generators). The outline of each strategy starts with a concise presentation of their fundamentals and carefully selected examples. In the end, we point out the controversies, challenges, and outlooks toward the future development of thermoelectric materials and devices. Overall, this review will serve to help materials scientists, chemists, and physicists, particularly students and young researchers, in selecting suitable strategies for the improvement of thermoelectrics and potentially other relevant energy conversion technologies.

Advanced Thermoelectric Design: From Materials and Structures to Devices

A strain sensor has been fabricated from a polymer nanocomposite with multiwalled carbon nanotube (MWNT) fillers. The piezoresistivity
of this nanocomposite strain sensor has been investigated based on an improved three-dimensional (3D) statistical resistor network
model incorporating the tunneling effect between the neighboring carbon nanotubes (CNTs), and a fiber reorientation model. The
numerical results agree very well with the experimental measurements. As compared with traditional strain gauges, much higher sensitivity
can be obtained in the nanocomposite sensors when the volume fraction of CNT is close to the percolation threshold. For a small
CNT volume fraction, weak nonlinear piezoresistivity is observed both experimentally and from numerical simulation. The tunneling
effect is considered to be the principal mechanism of the sensor under small strains.

Yaolu Liu

Papers

Molecular dynamics simulation on interfacial mechanical properties of polymer nanocomposites with wrinkled graphene

Interlaminar mechanical properties of carbon fiber reinforced plastic laminates modified with graphene oxide interleaf

Strengthening effects of twin interface in Cu/Ni multilayer thin films – A molecular dynamics study

Investigation on the interfacial mechanical properties of hybrid graphene-carbon nanotube/polymer nanocomposites

Environmentally-Friendly and Multifunctional Graphene-Silk Fabric Strain Sensor for Human-Motion Detection