Effective Confinement as Origin of the Equivalence of Kinetic Temperature and Fluctuation-Dissipation Ratio in a Dense Shear Driven Suspension

We study response and velocity autocorrelation functions for a tagged particle in a shear driven suspension governed by underdamped stochastic dynamics. We follow the idea of an effective confinement in dense suspensions and exploit a time-scale separation between particle reorganization and vibrational motion. This allows us to approximately derive the fluctuation-dissipation theorem in a "hybrid" form involving the kinetic temperature as an effective temperature and an additive correction term. We show numerically that even in a moderately dense suspension the latter is negligible. We discuss similarities and differences with a simple toy model, a single trapped particle in shear flow

Renormalized one-loop theory of correlations in polymer blends

The renormalized one-loop theory is a coarse-grained theory of corrections to the self-consistent field theory (SCFT) of polymer liquids, and to the random phase approximation (RPA) theory of composition fluctuations. We present predictions of corrections to the RPA for the structure function $S(k)$ and to the random walk model of single-chain statics in binary homopolymer blends. We consider an apparent interaction parameter $\chi_{a}$ that is defined by applying the RPA to the small $k$ limit of $S(k)$. The predicted deviation of $\chi_{a}$ from its long chain limit is proportional to $N^{-1/2}$, where $N$ is chain length. This deviation is positive (i.e., destabilizing) for weakly non-ideal mixtures, with \chi_{a} N \alt 1, but negative (stabilizing) near the critical point. The positive correction to $\chi_{a}$ for low values of $\chi_{a} N$ is a result of the fact that monomers in mixtures of shorter chains are slightly less strongly shielded from intermolecular contacts. The depression in $\chi_{a}$ near the critical point is a result of long-wavelength composition fluctuations. The one-loop theory predicts a shift in the critical temperature of ${\cal O}(N^{-1/2})$, which is much greater than the predicted ${\cal O}(N^{-1})$ width of the Ginzburg region. Chain dimensions deviate slightly from those of a random walk even in a one-component melt, and contract slightly with increasing $\chi_{e}$. Predictions for $S(k)$ and single-chain properties are compared to published lattice Monte Carlo simulations.Comment: submitted to J. Chem. Phy

Multiscaling for Systems with a Broad Continuum of Characteristic Lengths and Times: Structural Transitions in Nanocomposites

The multiscale approach to N-body systems is generalized to address the broad continuum of long time and length scales associated with collective behaviors. A technique is developed based on the concept of an uncountable set of time variables and of order parameters (OPs) specifying major features of the system. We adopt this perspective as a natural extension of the commonly used discrete set of timescales and OPs which is practical when only a few, widely-separated scales exist. The existence of a gap in the spectrum of timescales for such a system (under quasiequilibrium conditions) is used to introduce a continuous scaling and perform a multiscale analysis of the Liouville equation. A functional-differential Smoluchowski equation is derived for the stochastic dynamics of the continuum of Fourier component order parameters. A continuum of spatially non-local Langevin equations for the OPs is also derived. The theory is demonstrated via the analysis of structural transitions in a composite material, as occurs for viral capsids and molecular circuits.Comment: 28 pages, 1 figur

Mean-Field Treatment of the Many-Body Fokker-Planck Equation

We review some properties of the stationary states of the Fokker - Planck equation for N interacting particles within a mean field approximation, which yields a non-linear integrodifferential equation for the particle density. Analytical results show that for attractive long range potentials the steady state is always a precipitate containing one cluster of small size. For arbitrary potential, linear stability analysis allows to state the conditions under which the uniform equilibrium state is unstable against small perturbations and, via the Einstein relation, to define a critical temperature Tc separating two phases, uniform and precipitate. The corresponding phase diagram turns out to be strongly dependent on the pair-potential. In addition, numerical calculations reveal that the transition is hysteretic. We finally discuss the dynamics of relaxation for the uniform state suddenly cooled below Tc.Comment: 13 pages, 8 figure

Multiscale Modeling of Binary Polymer Mixtures: Scale Bridging in the Athermal and Thermal Regime

Obtaining a rigorous and reliable method for linking computer simulations of polymer blends and composites at different length scales of interest is a highly desirable goal in soft matter physics. In this paper a multiscale modeling procedure is presented for the efficient calculation of the static structural properties of binary homopolymer blends. The procedure combines computer simulations of polymer chains on two different length scales, using a united atom representation for the finer structure and a highly coarse-grained approach on the meso-scale, where chains are represented as soft colloidal particles interacting through an effective potential. A method for combining the structural information by inverse mapping is discussed, allowing for the efficient calculation of partial correlation functions, which are compared with results from full united atom simulations. The structure of several polymer mixtures is obtained in an efficient manner for several mixtures in the homogeneous region of the phase diagram. The method is then extended to incorporate thermal fluctuations through an effective chi parameter. Since the approach is analytical, it is fully transferable to numerous systems.Comment: in press, 13 pages, 7 figures, 6 table

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.

Anomalous interactions in confined charge-stabilized colloid

Charge-stabilized colloidal spheres dispersed in weak 1:1 electrolytes are supposed to repel each other. Consequently, experimental evidence for anomalous long-ranged like-charged attractions induced by geometric confinement inspired a burst of activity. This has largely subsided because of nagging doubts regarding the experiments' reliability and interpretation. We describe a new class of thermodynamically self-consistent colloidal interaction measurements that confirm the appearance of pairwise attractions among colloidal spheres confined by one or two bounding walls. In addition to supporting previous claims for this as-yet unexplained effect, these measurements also cast new light on its mechanism.Comment: 8 pages, 5 figures, RevTeX4. Conference proceedings for CODEF-04, Colloidal Dispersions in External Fields, March 29 - April 1, 200

Theoretical study of the thermal behavior of free and alumina-supported Fe-C nanoparticles

The thermal behavior of free and alumina-supported iron-carbon nanoparticles is investigated via molecular dynamics simulations, in which the effect of the substrate is treated with a simple Morse potential fitted to ab initio data. We observe that the presence of the substrate raises the melting temperature of medium and large $Fe_{1-x}C_x$ nanoparticles ($x$ = 0-0.16, $N$ = 80-1000, non- magic numbers) by 40-60 K; it also plays an important role in defining the ground state of smaller Fe nanoparticles ($N$ = 50-80). The main focus of our study is the investigation of Fe-C phase diagrams as a function of the nanoparticle size. We find that as the cluster size decreases in the 1.1-1.6-nm-diameter range the eutectic point shifts significantly not only toward lower temperatures, as expected from the Gibbs-Thomson law, but also toward lower concentrations of C. The strong dependence of the maximum C solubility on the Fe-C cluster size may have important implications for the catalytic growth of carbon nanotubes by chemical vapor deposition.Comment: 13 pages, 11 figures, higher quality figures can be seen in article 9 at http://alpha.mems.duke.edu/wahyu

Dimensional crossover of a boson gas in multilayers

We obtain the thermodynamic properties for a non-interacting Bose gas constrained on multilayers modeled by a periodic Kronig-Penney delta potential in one direction and allowed to be free in the other two directions. We report Bose-Einstein condensation (BEC) critical temperatures, chemical potential, internal energy, specific heat, and entropy for different values of a dimensionless impenetrability $P\geqslant 0$ between layers. The BEC critical temperature $T_{c}$ coincides with the ideal gas BEC critical temperature $T_{0}$ when $P=0$ and rapidly goes to zero as $P$ increases to infinity for any finite interlayer separation. The specific heat $C_{V}$ \textit{vs} $T$ for finite $P$ and plane separation $a$ exhibits one minimum and one or two maxima in addition to the BEC, for temperatures larger than $T_{c}$ which highlights the effects due to particle confinement. Then we discuss a distinctive dimensional crossover of the system through the specific heat behavior driven by the magnitude of $P$. For $T<T_{c}$ the crossover is revealed by the change in the slope of $\log C_{V}(T)$ and when $T>T_{c}$, it is evidenced by a broad minimum in $C_{V}(T)$.Comment: Ten pages, nine figure