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

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

Random walk approach to the d-dimensional disordered Lorentz gas

A correlated random walk approach to diffusion is applied to the disordered nonoverlapping Lorentz gas. By invoking the Lu-Torquato theory for chord-length distributions in random media [J. Chem. Phys. 98, 6472 (1993)], an analytic expression for the diffusion constant in arbitrary number of dimensions d is obtained. The result corresponds to an Enskog-like correction to the Boltzmann prediction, being exact in the dilute limit, and better or nearly exact in comparison to renormalized kinetic theory predictions for all allowed densities in d=2,3. Extensive numerical simulations were also performed to elucidate the role of the approximations involved.Comment: 5 pages, 5 figure

Phonon Life-times from first principles self consistent lattice dynamics

Phonon lifetime calculations from first principles usually rely on time consuming molecular dynamics calculations, or density functional perturbation theory (DFPT) where the zero temperature crystal structure is assumed to be dynamically stable. Here a new and effective method for calculating phonon lifetimes from first principles is presented, not limited to crystal structures stable at 0 K, and potentially much more effective than most corresponding molecular dynamics calculations. The method is based on the recently developed self consistent lattice dynamical method and is here tested by calculating the bcc phase phonon lifetimes of Li, Na, Ti and Zr, as representative examples.Comment: 4 pages, 4 figur

Stability of adhesion clusters under constant force

We solve the stochastic equations for a cluster of parallel bonds with shared constant loading, rebinding and the completely dissociated state as an absorbing boundary. In the small force regime, cluster lifetime grows only logarithmically with bond number for weak rebinding, but exponentially for strong rebinding. Therefore rebinding is essential to ensure physiological lifetimes. The number of bonds decays exponentially with time for most cases, but in the intermediate force regime, a small increase in loading can lead to much faster decay. This effect might be used by cell-matrix adhesions to induce signaling events through cytoskeletal loading.Comment: Revtex, 4 pages, 4 Postscript files include

Phonon Band Structure and Thermal Transport Correlation in a Layered Diatomic Crystal

To elucidate the relationship between a crystal's structure, its thermal conductivity, and its phonon dispersion characteristics, an analysis is conducted on layered diatomic Lennard-Jones crystals with various mass ratios. Lattice dynamics theory and molecular dynamics simulations are used to predict the phonon dispersion curves and the thermal conductivity. The layered structure generates directionally dependent thermal conductivities lower than those predicted by density trends alone. The dispersion characteristics are quantified using a set of novel band diagram metrics, which are used to assess the contributions of acoustic phonons and optical phonons to the thermal conductivity. The thermal conductivity increases as the extent of the acoustic modes increases, and decreases as the extent of the stop bands increases. The sensitivity of the thermal conductivity to the band diagram metrics is highest at low temperatures, where there is less anharmonic scattering, indicating that dispersion plays a more prominent role in thermal transport in that regime. We propose that the dispersion metrics (i) provide an indirect measure of the relative contributions of dispersion and anharmonic scattering to the thermal transport, and (ii) uncouple the standard thermal conductivity structure-property relation to that of structure-dispersion and dispersion-property relations, providing opportunities for better understanding of the underlying physical mechanisms and a potential tool for material design.Comment: 30 pages, 10 figure

Spectra of sparse non-Hermitian random matrices: an analytical solution

We present the exact analytical expression for the spectrum of a sparse non-Hermitian random matrix ensemble, generalizing two classical results in random-matrix theory: this analytical expression forms a non-Hermitian version of the Kesten-Mckay law as well as a sparse realization of Girko's elliptic law. Our exact result opens new perspectives in the study of several physical problems modelled on sparse random graphs. In this context, we show analytically that the convergence rate of a transport process on a very sparse graph depends upon the degree of symmetry of the edges in a non-monotonous way.Comment: 5 pages, 5 figures, 12 pages supplemental materia

Fermion Quasi-Spherical Harmonics

Spherical Harmonics, $Y_\ell^m(\theta,\phi)$, are derived and presented (in a Table) for half-odd-integer values of $\ell$ and $m$. These functions are eigenfunctions of $L^2$ and $L_z$ written as differential operators in the spherical-polar angles, $\theta$ and $\phi$. The Fermion Spherical Harmonics are a new, scalar and angular-coordinate-dependent representation of fermion spin angular momentum. They have $4\pi$ symmetry in the angle $\phi$, and hence are not single-valued functions on the Euclidean unit sphere; they are double-valued functions on the sphere, or alternatively are interpreted as having a double-sphere as their domain.Comment: 16 pages, 2 Tables. Submitted to J.Phys.

Efficiency of quantum and classical transport on graphs

We propose a measure to quantify the efficiency of classical and quantum mechanical transport processes on graphs. The measure only depends on the density of states (DOS), which contains all the necessary information about the graph. For some given (continuous) DOS, the measure shows a power law behavior, where the exponent for the quantum transport is twice the exponent of its classical counterpart. For small-world networks, however, the measure shows rather a stretched exponential law but still the quantum transport outperforms the classical one. Some finite tree-graphs have a few highly degenerate eigenvalues, such that, on the other hand, on them the classical transport may be more efficient than the quantum one.Comment: 5 pages, 3 figure

A First Principle Approach to Rescale the Dynamics of Simulated Coarse-Grained Macromolecular Liquids

We present a detailed derivation and testing of our approach to rescale the dynamics of mesoscale simulations of coarse-grained polymer melts (I. Y. Lyubimov et al. J. Chem. Phys. \textbf{132}, 11876, 2010). Starting from the first-principle Liouville equation and applying the Mori-Zwanzig projection operator technique, we derive the Generalized Langevin Equations (GLE) for the coarse-grained representations of the liquid. The chosen slow variables in the projection operators define the length scale of coarse graining. Each polymer is represented at two levels of coarse-graining: monomeric as a bead-and-spring model and molecular as a soft-colloid. In the long-time regime where the center-of-mass follows Brownian motion and the internal dynamics is completely relaxed, the two descriptions must be equivalent. By enforcing this formal relation we derive from the GLEs the analytical rescaling factors to be applied to dynamical data in the coarse-grained representation to recover the monomeric description. Change in entropy and change in friction are the two corrections to be accounted for to compensate the effects of coarse-graining on the polymer dynamics. The solution of the memory functions in the coarse-grained representations provides the dynamical rescaling of the friction coefficient. The calculation of the internal degrees of freedom provides the correction of the change in entropy due to coarse-graining. The resulting rescaling formalism is a function of the coarse-grained model and thermodynamic parameters of the system simulated. The rescaled dynamics obtained from mesoscale simulations of polyethylene, represented as soft colloidal particles, by applying our rescaling approach shows a good agreement with data of translational diffusion measured experimentally and from simulations. The proposed method is used to predict self-diffusion coefficients of new polyethylene samples.Comment: 21 pages, 6 figures, 6 tables. Submitted to Phys. Rev.