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

DLS and zeta potential - What they are and what they are not?

A diverse spectrum of technology drivers such as improved thermal barriers, higher efficiency thermoelectric energy conversion, phase-change memory, heat-assisted magnetic recording, thermal management of nanoscale electronics, and nanoparticles for thermal medical therapies are motivating studies of the applied physics of thermal transport at the nanoscale. This review emphasizes developments in experiment, theory, and computation in the past ten years and summarizes the present status of the field. Interfaces become increasingly important on small length scales. Research during the past decade has extended studies of interfaces between simple metals and inorganic crystals to interfaces with molecular materials and liquids with systematic control of interface chemistry and physics. At separations on the order of ∼1 nm, the science of radiative transport through nanoscale gaps overlaps with thermal conduction by the coupling of electronic and vibrational excitations across weakly bonded or rough interface...

/pdf/nanoscale-thermal-transport-ii-2003-2012-4upcs4bq86.pdf

Nanoscale thermal transport. II. 2003–2012

Emerging challenges and materials for thermal management of electronics

Properties of phonons-quanta of the crystal lattice vibrations-in graphene have recently attracted significant attention from the physics and engineering communities. Acoustic phonons are the main heat carriers in graphene near room temperature, while optical phonons are used for counting the number of atomic planes in Raman experiments with few-layer graphene. It was shown both theoretically and experimentally that transport properties of phonons, i.e. energy dispersion and scattering rates, are substantially different in a quasi-two-dimensional system such as graphene compared to the basal planes in graphite or three-dimensional bulk crystals. The unique nature of two-dimensional phonon transport translates into unusual heat conduction in graphene and related materials. In this review, we outline different theoretical approaches developed for phonon transport in graphene, discuss contributions of the in-plane and cross-plane phonon modes, and provide comparison with available experimental thermal conductivity data. Particular attention is given to analysis of recent results for the phonon thermal conductivity of single-layer graphene and few-layer graphene, and the effects of the strain, defects, and isotopes on phonon transport in these systems.

https://iopscience.iop.org/article/10.1088/0953-8984/24/23/233203/pdf

Two-dimensional phonon transport in graphene

In-plane lattice thermal conductivities of multilayer graphene films

Interfacial thermal resistance in multilayer graphene structures

Molecular dynamics simulation is performed to extract the phonon dispersion and phonon lifetime of single layer graphene. The mode dependent thermal conductivity is calculated from the phonon kinetic theory. The predicted thermal conductivity at room temperature exhibits important quantum effects due to the high Debye temperature of graphene. But the quantum effects are reduced significantly when the simulated temperature is as high as 1000 K. Our calculations show that out-of-plane modes contribute about 41.1% to the total thermal conductivity at room temperature. The relative contribution of out-of-plane modes has a little decrease with the increase of temperature. Contact with substrate can reduce both the total thermal conductivity of graphene and the relative contribution of out-of-plane modes, in agreement with previous experiments and theories. Increasing the coupling strength between graphene and substrate can further reduce the relative contribution of out-of-plane modes. The present investigations also show that the relative contribution of different mode phonons is not sensitive to the grain size of graphene. The obtained phonon relaxation time provides useful insight for understanding the phonon mean free path and the size effects in graphene.

Mode dependent lattice thermal conductivity of single layer graphene

https://link.aps.org/accepted/10.1103/PhysRevB.87.125426

Anisotropic Debye model for the thermal boundary conductance

Phonon transport in the c-axis direction of graphite thin films has been studied using non-equilibrium molecular dynamics (MD) simulation. The simulation results show that the c-axis thermal conductivities for films of thickness ranging from 20 to 500 atomic layers are significantly lower than the bulk value. Based on the MD data, a method is developed to construct the c-axis thermal conductivity as an accumulation function of phonon mean free path (MFP), from which we show that phonons with MFPs from 2 to 2000 nm contribute ∼80% of the graphite c-axis thermal conductivity at room temperature, and phonons with MFPs larger than 100 nm contribute over 40% to the c-axis thermal conductivity. These findings indicate that the commonly believed value of just a few nanometers from the simple kinetic theory drastically underestimates the c-axis phonon MFP of graphite.

Zhiyong Wei

Papers

In-plane lattice thermal conductivities of multilayer graphene films

Interfacial thermal resistance in multilayer graphene structures

Mode dependent lattice thermal conductivity of single layer graphene

Anisotropic Debye model for the thermal boundary conductance

Phonon mean free path of graphite along the c-axis