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

BACKGROUND Materials by design is an appealing idea that is very hard to realize in practice. Combining the best of different ingredients in one ultimate material is a task for which we currently have no general solution. However, we do have some successful examples to draw upon: Composite materials and III-V heterostructures have revolutionized many aspects of our lives. Still, we need a general strategy to solve the problem of mixing and matching crystals with different properties, creating combinations with predetermined attributes and functionalities. ADVANCES Two-dimensional (2D) materials offer a platform that allows creation of heterostructures with a variety of properties. One-atom-thick crystals now comprise a large family of these materials, collectively covering a very broad range of properties. The first material to be included was graphene, a zero-overlap semimetal. The family of 2D crystals has grown to includes metals (e.g., NbSe 2 ), semiconductors (e.g., MoS 2 ), and insulators [e.g., hexagonal boron nitride (hBN)]. Many of these materials are stable at ambient conditions, and we have come up with strategies for handling those that are not. Surprisingly, the properties of such 2D materials are often very different from those of their 3D counterparts. Furthermore, even the study of familiar phenomena (like superconductivity or ferromagnetism) in the 2D case, where there is no long-range order, raises many thought-provoking questions. A plethora of opportunities appear when we start to combine several 2D crystals in one vertical stack. Held together by van der Waals forces (the same forces that hold layered materials together), such heterostructures allow a far greater number of combinations than any traditional growth method. As the family of 2D crystals is expanding day by day, so too is the complexity of the heterostructures that could be created with atomic precision. When stacking different crystals together, the synergetic effects become very important. In the first-order approximation, charge redistribution might occur between the neighboring (and even more distant) crystals in the stack. Neighboring crystals can also induce structural changes in each other. Furthermore, such changes can be controlled by adjusting the relative orientation between the individual elements. Such heterostructures have already led to the observation of numerous exciting physical phenomena. Thus, spectrum reconstruction in graphene interacting with hBN allowed several groups to study the Hofstadter butterfly effect and topological currents in such a system. The possibility of positioning crystals in very close (but controlled) proximity to one another allows for the study of tunneling and drag effects. The use of semiconducting monolayers leads to the creation of optically active heterostructures. The extended range of functionalities of such heterostructures yields a range of possible applications. Now the highest-mobility graphene transistors are achieved by encapsulating graphene with hBN. Photovoltaic and light-emitting devices have been demonstrated by combining optically active semiconducting layers and graphene as transparent electrodes. OUTLOOK Currently, most 2D heterostructures are composed by direct stacking of individual monolayer flakes of different materials. Although this method allows ultimate flexibility, it is slow and cumbersome. Thus, techniques involving transfer of large-area crystals grown by chemical vapor deposition (CVD), direct growth of heterostructures by CVD or physical epitaxy, or one-step growth in solution are being developed. Currently, we are at the same level as we were with graphene 10 years ago: plenty of interesting science and unclear prospects for mass production. Given the fast progress of graphene technology over the past few years, we can expect similar advances in the production of the heterostructures, making the science and applications more achievable.

2D materials and van der Waals heterostructures

We present the science and technology roadmap for graphene, related two-dimensional crystals, and hybrid systems, targeting an evolution in technology, that might lead to impacts and benefits reaching into most areas of society. This roadmap was developed within the framework of the European Graphene Flagship and outlines the main targets and research areas as best understood at the start of this ambitious project. We provide an overview of the key aspects of graphene and related materials (GRMs), ranging from fundamental research challenges to a variety of applications in a large number of sectors, highlighting the steps necessary to take GRMs from a state of raw potential to a point where they might revolutionize multiple industries. We also define an extensive list of acronyms in an effort to standardize the nomenclature in this emerging field.

/pdf/science-and-technology-roadmap-for-graphene-related-two-4q9m7czlkc.pdf

Science and technology roadmap for graphene, related two-dimensional crystals, and hybrid systems

Changes in portlandite morphology with solvent composition: Atomistic simulations and experiment

The widespread use of controlled molecular-level motion in key natural processes suggests that great rewards could come from bridging the gap between the present generation of synthetic molecular systems, which by and large rely upon electronic and chemical effects to carry out their functions, and the machines of the macroscopic world, which utilize the synchronized movements of smaller parts to perform specific tasks. This is a scientific area of great contemporary interest and extraordinary recent growth, yet the notion of molecular-level machines dates back to a time when the ideas surrounding the statistical nature of matter and the laws of thermodynamics were first being formulated. Here we outline the exciting successes in taming molecular-level movement thus far, the underlying principles that all experimental designs must follow, and the early progress made towards utilizing synthetic molecular structures to perform tasks using mechanical motion. We also highlight some of the issues and challenges that still need to be overcome.

Synthetic molecular motors and mechanical machines.

The retraction energy of a multiwalled carbon nanotube with an extruded core causes the core to oscillate with respect to its fully retracted position where the van der Waals potential energy is minimized. This phenomenon, together with the available nanotechnology, leads to the creation of nanomechanical systems of operating frequency up to several gigahertz.

Multiwalled carbon nanotubes as gigahertz oscillators

The use of nanomaterials for strain sensors has attracted attention due to their unique electromechanical properties. However, nanomaterials have yet to overcome many technological obstacles and thus are not yet the preferred material for strain sensors. In this work, we investigated graphene woven fabrics (GWFs) for strain sensing. Different than graphene films, GWFs undergo significant changes in their polycrystalline structures along with high-density crack formation and propagation mechanically deformed. The electrical resistance of GWFs increases exponentially with tensile strain with gauge factors of ~103 under 2~6% strains and ~106 under higher strains that are the highest thus far reported, due to its woven mesh configuration and fracture behavior, making it an ideal structure for sensing tensile deformation by changes in strain. The main mechanism is investigated, resulting in a theoretical model that predicts very well the observed behavior.

/pdf/stretchable-and-highly-sensitive-graphene-on-polymer-strain-4ojlvuqrxq.pdf

Stretchable and highly sensitive graphene-on-polymer strain sensors

Upon shearing a microscale lithographically defined graphite mesa, the sheared section retracts spontaneously to minimize interface energy. Here, we demonstrate a sixfold symmetry of the self-retraction and provide a first experimental estimate of the frictional force involved, as direct evidence that the self-retraction is due to superlubricity, where ultralow friction occurs between incommensurate surfaces. The effect is remarkable because it occurs reproducibly under ambient conditions and over a contact area of up to 10×10  μm2, more than 7 orders of magnitude larger than previous scanning-probe-based studies of superlubricity in graphite. By analyzing the sheared interface, we show how the grain structure of highly oriented pyrolitic graphite determines the probability of self-retraction. Our results demonstrate that such self-retraction provides a novel probe of superlubricity, and the robustness of the phenomenon opens the way for practical applications of superlubricity in micromechanical systems.

Observation of Microscale Superlubricity in Graphite

Theory of Representations for Tensor Functions—A Unified Invariant Approach to Constitutive Equations

The underlying mechanisms of stability, metastability, or instability of the Cassie-Baxter and Wenzel wetting modes and their transitions on superhydrophobic surfaces decorated with periodic micropillars are quantitatively studied in this article Hydraulic pressure, which may be generated by the water-air interfacial tension of water droplets or external factors such as raining impact, is shown to be a key to understanding these mechanisms A detailed transition process driven by increasing hydraulic pressure is numerically simulated The maximum sustainable or critical pressure of the Cassie-Baxter wetting state on a pillarlike microstructural surface is formulated for the first time in a simple, unified, and precise form This analytic result reveals the fact that reducing the microstructural scales (eg, the pillars' diameters and spacing) is probably the most efficient measure needed to enlarge the critical pressure significantly We also introduce a dimensionless parameter, the pillar slenderness ratio, to characterize the stability of either the Cassie-Baxter or the Wenzel wetting state and show that the energy barrier for transitioning from the Cassie-Baxter to the Wenzel wetting mode is proportional to both the slenderness ratio and the area fraction Thus, the Cassie-Baxter wetting mode may collapse under a hydraulic pressure lower than the critical one if the slenderness ratio is improperly small This quantitative study explains fairly well some experimental observations of contact angles that can be modeled by neither Wenzel nor Cassie-Baxter contact angles and eventually leads to our proposals for a mixed (or coexisting) wetting mode

Quanshui Zheng

Papers

Multiwalled carbon nanotubes as gigahertz oscillators

Stretchable and highly sensitive graphene-on-polymer strain sensors

Observation of Microscale Superlubricity in Graphite

Theory of Representations for Tensor Functions—A Unified Invariant Approach to Constitutive Equations

Effects of hydraulic pressure on the stability and transition of wetting modes of superhydrophobic surfaces.