Kinetic Theory of Flocking: Derivation of Hydrodynamic Equations

It is shown how to explicitly coarse-grain the microscopic dynamics of the Vicsek model for self-propelled agents. The macroscopic transport equations are derived by means of an Enskog-type kinetic theory. Expressions for all transport coefficients at large particle speed are given. The phase transition from a disordered to a flocking state is studied numerically and analytically.Comment: 4 pages, 1 figur

Computational study of the thermal conductivity in defective carbon nanostructures

We use non-equilibrium molecular dynamics simulations to study the adverse role of defects including isotopic impurities on the thermal conductivity of carbon nanotubes, graphene and graphene nanoribbons. We find that even in structurally perfect nanotubes and graphene, isotopic impurities reduce thermal conductivity by up to one half by decreasing the phonon mean free path. An even larger thermal conductivity reduction, with the same physical origin, occurs in presence of structural defects including vacancies and edges in narrow graphene nanoribbons. Our calculations reconcile results of former studies, which differed by up to an order of magnitude, by identifying limitations of various computational approaches

Configurational temperatures and interactions in charge-stabilized colloid

We demonstrate that the configurational temperature formalism can be derived from the classical hypervirial theorem, and introduce a hierarchy of hyperconfigurational temperature definitions, which are particularly well suited for experimental studies. We then use these analytical tools to probe the electrostatic interactions in monolayers of charge-stabilized colloidal spheres confined by parallel glass surfaces. The configurational and hyperconfigurational temperatures, together with a novel thermodynamic sum rule, provide previously lacking self-consistency tests for interaction measurements based on digital video microscopy, and thereby cast new light on controversial reports of confinement-induced like-charge attractions. We further introduce a new method for measuring the pair potential directly that uses consistency of the configurational and hyperconfigurational temperatures as a set of constraints for a model-free search.Comment: 15 pages, 12 figures, submitted to J. Chem. Phy

Calculations of canonical averages from the grand canonical ensemble

Grand canonical and canonical ensembles become equivalent in the thermodynamic limit, but when the system size is finite the results obtained in the two ensembles deviate from each other. In many important cases, the canonical ensemble provides an appropriate physical description but it is often much easier to perform the calculations in the corresponding grand canonical ensemble. We present a method to compute averages in canonical ensemble based on calculations of the expectation values in grand canonical ensemble. The number of particles, which is fixed in the canonical ensemble, is not necessarily the same as the average number of particles in the grand canonical ensemble

A Simple Three-Parameter Model Potential For Diatomic Systems: From Weakly and Strongly Bound Molecules to Metastable Molecular Ions

Based on a simplest molecular orbital theory of H$_{2}^{+}$, a three-parameter model potential function is proposed to describe ground-state diatomic systems with closed-shell and/or S-type valence-shell constituents over a significantly wide range of internuclear distances. More than 200 weakly and strongly bound diatomics have been studied, including neutral and singly-charged diatomics (e.g., H$_{2}$, Li$_{2}$, LiH, Cd$_{2}$, Na$_{2}^{+}$, and RbH$^{-}$), long-range bound diatomics (e.g., NaAr, CdNe, He$_{2}$, CaHe, SrHe, and BaHe), metastable molecular dications (e.g., BeH$^{++}$, AlH$^{++}$, Mg$_{2}^{++}$, and LiBa$^{++}$), and molecular trications (e.g., YHe$^{+++}$ and ScHe$^{+++}$).Comment: 5 pages, 4 figures, accepted by Physical Review Letter

Long wavelength structural anomalies in jammed systems

The structural properties of static, jammed packings of monodisperse spheres in the vicinity of the jamming transition are investigated using large-scale computer simulations. At small wavenumber $k$, we argue that the anomalous behavior in the static structure factor, $S(k) \sim k$, is consequential of an excess of low-frequency, collective excitations seen in the vibrational spectrum. This anomalous feature becomes more pronounced closest to the jamming transition, such that $S(0) \to 0$ at the transition point. We introduce an appropriate dispersion relation that accounts for these phenomena that leads us to relate these structural features to characteristic length scales associated with the low-frequency vibrational modes of these systems. When the particles are frictional, this anomalous behavior is suppressed providing yet more evidence that jamming transitions of frictional spheres lie at lower packing fractions that that for frictionless spheres. These results suggest that the mechanical properties of jammed and glassy media may therefore be inferred from measurements of both the static and dynamical structure factors.Comment: 8 pages, 6 figure captions. Completely revised version to appear in Phys. Rev.

Elastic Energy, Fluctuations and Temperature for Granular Materials

We probe, using a model system, elastic and kinetic energies for sheared granular materials. For large enough $P/E_y$ (pressure/Young's modulus) and $P/\rho v^2$ ($P/$kinetic energy density) elastic dominates kinetic energy, and energy fluctuations become primarily elastic in nature. This regime has likely been reached in recent experiments. We consider a generalization of the granular temperature, $T_g$, with both kinetic and elastic terms and that changes smoothly from one regime to the other. This $T_g$ is roughly consistent with a temperature adapted from equilibrium statistical mechanics.Comment: 4 pages, 4 figure

Analytical Rescaling of Polymer Dynamics from Mesoscale Simulations

We present a theoretical approach to scale the artificially fast dynamics of simulated coarse-grained polymer liquids down to its realistic value. As coarse-graining affects entropy and dissipation, two factors enter the rescaling: inclusion of intramolecular vibrational degrees of freedom, and rescaling of the friction coefficient. Because our approach is analytical, it is general and transferable. Translational and rotational diffusion of unentangled and entangled polyethylene melts, predicted from mesoscale simulations of coarse-grained polymer melts using our rescaling procedure, are in quantitative agreement with united atom simulations and with experiments.Comment: 6 pages, 2 figures, 2 table

On Calculation of Thermal Conductivity from Einstein Relation in Equilibrium MD

In equilibrium molecular dynamics, Einstein relation can be used to calculate the thermal conductivity. This method is equivalent to Green-Kubo relation and it does not require a derivation of an analytical form for the heat current. However, it is not commonly used as Green-Kubo relationship. Its wide use is hindered by the lack of a proper definition for integrated heat current (energy moment) under periodic boundary conditions. In this paper, we developed an appropriate definition for integrated heat current to calculate thermal conductivity of solids under periodic conditions. We applied this method to solid argon and silicon based systems; compared and contrasted with the Green-Kubo approach.Comment: We updated this manuscript from second version by changing the title and abstract. This paper is submitted to J. Chem. Phy