Prediction of a structural transition in the hard disk fluid

Starting from the second equilibrium equation in the BBGKY hierarchy under the Kirkwood superposition closure, we implement a new method for studying the asymptotic decay of correlations in the hard disk fluid in the high density regime. From our analysis and complementary numerical studies, we find that exponentially damped oscillations can occur only up to a packing fraction η(?)～0.718, a value that is in substantial agreement with the packing fraction, η～0.723, believed to characterize the transition from the ordered solid phase to a dense fluid phase, as inferred from Mak's Monte Carlo simulations [Phys. Rev. E 73, 065104 (2006)]. Next, we show that the same method of analysis predicts that the exponential damping of oscillations in the hard sphere fluid becomes impossible when λ=4nπσ(3)[1+H(1)]≥34.81, where H(1) is the contact value of the correlation function, n is the number density, and σ is the sphere diameter in exact agreement with the condition, λ≥34.8, which is first reported in a numerical study of the Kirkwood equation by Kirkwood et al. [J. Chem. Phys. 18, 1040 (1950)]. Finally, we show that our method confirms the absence of any structural transition in hard rods for the entire range of densities below close packing.     

Substrate concentration dependence of the diffusion-controlled steady-state rate constant

The Smoluchowski approach to diffusion-controlled reactions is generalized to interacting substrate particles by including the osmotic pressure and hydrodynamic interactions of the nonideal particles in the Smoluchoswki equation within a local-density approximation. By solving the strictly linearized equation for the time-independent case with absorbing boundary conditions, we present an analytic expression for the diffusion-limited steady-state rate constant for small substrate concentrations in terms of an effective second virial coefficient B2*. Comparisons to Brownian dynamics simulations excluding hydrodynamic interactions show excellent agreement up to bulk number densities of B2*rho0 < approximately = 0.4 for hard sphere and repulsive Yukawa-like interactions between the substrates. Our study provides an alternative way to determine the second virial coefficient of interacting macromolecules experimentally by measuring their steady-state rate constant in diffusion-controlled reactions at low densities.     

Brownian dynamics simulations of the self- and collective rotational diffusion coefficients of rigid long thin rods

Recently a microscopic theory for the dynamics of suspensions of long thin rigid rods was presented, confirming and expanding the well-known theory by Doi and Edwards [The Theory of Polymer Dynamics (Clarendon, Oxford, 1986)] and Kuzuu [J. Phys. Soc. Jpn. 52, 3486 (1983)]. Here this theory is put to the test by comparing it against computer simulations. A Brownian dynamics simulation program was developed to follow the dynamics of the rods, with a length over a diameter ratio of 60, on the Smoluchowski time scale. The model accounts for excluded volume interactions between rods, but neglects hydrodynamic interactions. The self-rotational diffusion coefficients D(r)(phi) of the rods were calculated by standard methods and by a new, more efficient method based on calculating average restoring torques. Collective decay of orientational order was calculated by means of equilibrium and nonequilibrium simulations. Our results show that, for the currently accessible volume fractions, the decay times in both cases are virtually identical. Moreover, the observed decay of diffusion coefficients with volume fraction is much quicker than predicted by the theory, which is attributed to an oversimplification of dynamic correlations in the theory.     

Diffusion equations and hard collisions in multiple scattering of charged particles

The processes of angular-spatial evolution of multiple scattering of charged particles are described by the Lewis (special case of Boltzmann) integro-differential equation. The underlying stochastic process for this evolution is the compound Poisson process with transition densities satisfying the Lewis equation. In this paper we derive the Lewis equation from the compound Poisson process and show that the effective method of the solution of this equation can be based on the idea of decomposition of the compound Poisson process into processes of soft and hard collisions. Formulas for transition densities of soft and hard collision processes are provided in this paper together with the formula expressing the general solution of the Lewis equation in terms of those transition densities.     

Diffusion of spheres in crowded suspensions of rods

Translational tracer diffusion of spherical macromolecules in crowded suspensions of rodlike colloids is investigated. Experiments are done using several kinds of spherical tracers in fd-virus suspensions. A wide range of size ratios L/2a of the length L of the rods and the diameter 2a of the tracer sphere is covered by combining several experimental methods: fluorescence correlation spectroscopy for small tracer spheres, dynamic light scattering for intermediate sized spheres, and video microscopy for large spheres. Fluorescence correlation spectroscopy is shown to measure long-time diffusion only for relatively small tracer spheres. Scaling of diffusion coefficients with a/xi, predicted for static networks, is not found for our dynamical network of rods (with xi the mesh size of the network). Self-diffusion of tracer spheres in the dynamical network of freely suspended rods is thus fundamentally different as compared to cross-linked networks. A theory is developed for the rod-concentration dependence of the translational diffusion coefficient at low rod concentrations for freely suspended rods. The proposed theory is based on a variational solution of the appropriate Smoluchowski equation without hydrodynamic interactions. The theory can, in principle, be further developed to describe diffusion through dynamical networks at higher rod concentrations with the inclusion of hydrodynamic interactions. Quantitative agreement with the experiments is found for large tracer spheres, and qualitative agreement for smaller spheres. This is probably due to the increasing importance of hydrodynamic interactions as compared to direct interactions as the size of the tracer sphere decreases.     

Modeling and simulation of equation of state for nonadditive hard disk mixtures

We present exact results for mixtures of nonadditive hard disks and use some of them to derive a consistent model for the equation of state. We also performed molecular dynamics simulation for hard disks over a wide range of size ratios. Comparison of the model to the data shows that the model is accurate for all densities in the case of additive and slightly nonadditive (nonadditivity parameter within ±0.1) mixtures. For large nonadditivity, the model is accurate for low to moderate densities only, and starts to deteriorate at high densities.     

Isotropic-nematic transition of hard rods immersed in random sphere matrices

Using replica density functional theory and Monte Carlo computer simulations we investigate a system of annealed hard spherocylinders adsorbed in a matrix of quenched hard spheres. Theoretical predictions for the partition coefficient, defined as the ratio of density of rods in the matrix and that in a reservoir, agree well with simulation results. Theory predicts the isotropic-nematic transition to remain first order upon increasing sphere packing fraction, and to shift towards lower rod densities. This scenario is consistent with our simulation results that clearly show a jump in the nematic order parameter upon increasing the rod density at constant matrix packing fraction, corresponding to the isotropic-nematic transition, even for sphere matrix packing fractions < or approximately equal to 0.3.     

Pair diffusion, hydrodynamic interactions, and available volume in dense fluids

We calculate the pair diffusion coefficient D(r) as a function of the distance r between two hard sphere particles in a dense monodisperse fluid. The distance-dependent pair diffusion coefficient describes the hydrodynamic interactions between particles in a fluid that are central to theories of polymer and colloid dynamics. We determine D(r) from the propagators (Green's functions) of particle pairs obtained from molecular dynamics simulations. At distances exceeding ~3 molecular diameters, the calculated pair diffusion coefficients are in excellent agreement with predictions from exact macroscopic hydrodynamic theory for large Brownian particles suspended in a solvent bath, as well as the Oseen approximation. However, the asymptotic 1/r distance dependence of D(r) associated with hydrodynamic effects emerges only after the pair distance dynamics has been followed for relatively long times, indicating non-negligible memory effects in the pair diffusion at short times. Deviations of the calculated D(r) from the hydrodynamic models at short distances r reflect the underlying many-body fluid structure, and are found to be correlated to differences in the local available volume. The procedure used here to determine the pair diffusion coefficients can also be used for single-particle diffusion in confinement with spherical symmetry.     

Momentum and stress relaxation in fluids illustrating caging

The self diffusion coefficient, shear viscosity, and velocity time correlation function are calculated for a hard sphere fluid under a severe assumption, namely, the friction arises from uncorrelated binary collisions and from correlated backscattering (caging) collisions as represented in the memory function. Relaxation of the memory function from its zerotime caging value is described as a diffusion process. Derived diffusion coefficients and the shear viscosities, relative to their Enskog values decrease and increase with density, respectively, in a monotonic and gradual fashion in contrast with simulation values that show a precipitous change near the fluid-solid transition. In the present pair diffusion model, the velocity time correlation function vanishes at the proper time but its tail is overly damped relative to the simulation data. A weak breakdown of the Stokes-Einstein relation is also predicted.     

Heterogeneous dynamics in columnar liquid crystals of parallel hard rods

In the wake of previous studies on the rattling-and-jumping diffusion in smectic liquid crystal phases of colloidal rods, we analyze here for the first time the heterogeneous dynamics in columnar phases. More specifically, we perform computer simulations to investigate the relaxation dynamics of a binary mixture of perfectly aligned hard spherocylinders. We detect that the columnar arrangement of the system produces free-energy barriers that the particles should overcome to jump from one column to another, thus determining a hopping-type diffusion. This phenomenon accounts for the non-Gaussian intercolumn diffusion and shows a two-step structural relaxation that is remarkably analogous to that of out-of-equilibrium glass-forming systems and gels. Surprisingly enough, slight deviations from the behavior of simple liquids due to transient cages is also observed in the direction perpendicular to this plane, where the system is usually referred to as liquidlike.     

Strong effect of hydrodynamic coupling on the electric dichroism of bent rods

The effect of hydrodynamic coupling on the spatial orientation of rigid bent rods in electric fields has been analyzed by Brownian dynamics simulations. Bead models for smoothly bent rods were constructed with dimensions of DNA double helices, and established simulation procedures were used to calculate their diffusion tensor, including the translational-rotational coupling tensor. The electric and optical parameters were assigned on the basis of known properties of double helices. Brownian dynamics simulations of the orientation of these models in electric fields showed that both transients and amplitudes of the calculated dichroism are very strongly dependent on translational-rotational coupling over a wide range of electric field strengths. For example, the stationary dichroism of a smoothly bent 179 bp DNA fragment calculated at low field strengths is positive in the presence and negative in the absence of hydrodynamic coupling. The transients are converted from a biphasic to a monophasic shape, when hydrodynamic coupling is turned off. The large changes resulting from hydrodynamic coupling were controlled by calculations based on analytical expressions derived for electrooptical response curves in the limit of low electric field strengths; the results obtained by this independent approach are in very satisfactory agreement with our Brownian dynamics simulations. The effect is strongly dependent on the electric dipole and on its direction. In the absence of any dipole the coupling effect was not observed. The coupling effect increases with the size of the bent rods. Because most macromolecular structures are known to have induced and/or permanent dipole moments, large effects of hydrodynamic coupling on both the amplitudes and the transients of the electric dichroism/birefringence must be expected in general for structures with nonsymmetric shape.     

Nonequilibrium inertial dynamics of colloidal systems

We consider the properties of a one-dimensional fluid of Brownian inertial hard-core particles, whose microscopic dynamics is partially damped by a heat bath. Direct interactions among the particles are represented as binary, instantaneous elastic collisions. Collisions with the heat bath are accounted for by a Fokker-Planck collision operator, whereas direct collisions among the particles are treated by a well known method of kinetic theory, the revised Enskog theory. By means of a time multiple time-scale method we derive the evolution equation for the average density. Remarkably, for large values of the friction parameter and/or of the mass of the particles we obtain the same equation as the one derived within the dynamic density functional theory (DDF). In addition, at moderate values of the friction constant, the present method allows to study the inertial effects not accounted for by DDF method. Finally, a numerical test of these corrections is provided.     

Application of a new equation of state to liquid refrigerant mixtures

A new general equation of state recently reported for pure liquids has been developed to predict the volumetric and thermodynamic properties of six binary and two ternary liquid refrigerant mixtures (including HCs and HFCs mixtures) at different temperatures, pressures, and compositions. The results show this equation of state can be used to reproduce and predict different thermodynamic properties of liquid refrigerant mixtures within experimental errors. The composition dependence of the parameters of this equation of state has been assumed as quadratic functions of mole fraction. Using these mixing rules, the agreement between calculated and experimental densities is better than 0.6% for binary mixtures and 2.3% for ternary mixtures. To compare the performance of this new equation of state against other well-known methods such as the COSTALD method, the density of some refrigerant mixtures, for which the parameters of COSTALD were available, has been computed and compared with those of this new equation of state.     

Structure of polymer fluids in colloidal dispersions

We have calculated some structural properties of a fluid of hard sphere polymer chains about a variable sized central hard sphere with the Monte Carlo method in the canonical ensemble. We have additionally calculated these structural properties with an integral equation based on density functional theory. The integral equation theory gives good agreement with the simulations at all but the highest densities.     

The van Hove distribution function for brownian hard spheres: dynamical test particle theory and computer simulations for bulk dynamics

We describe a test particle approach based on dynamical density functional theory (DDFT) for studying the correlated time evolution of the particles that constitute a fluid. Our theory provides a means of calculating the van Hove distribution function by treating its self and distinct parts as the two components of a binary fluid mixture, with the "self?" component having only one particle, the "distinct" component consisting of all the other particles, and using DDFT to calculate the time evolution of the density profiles for the two components. We apply this approach to a bulk fluid of Brownian hard spheres and compare to results for the van Hove function and the intermediate scattering function from Brownian dynamics computer simulations. We find good agreement at low and intermediate densities using the very simple Ramakrishnan-Yussouff [Phys. Rev. B 19, 2775 (1979)] approximation for the excess free energy functional. Since the DDFT is based on the equilibrium Helmholtz free energy functional, we can probe a free energy landscape that underlies the dynamics. Within the mean-field approximation we find that as the particle density increases, this landscape develops a minimum, while an exact treatment of a model confined situation shows that for an ergodic fluid this landscape should be monotonic. We discuss possible implications for slow, glassy, and arrested dynamics at high densities.     

Structures and dynamics of thermosensitive microgel suspensions studied with three-dimensional cross-correlated light scattering

The structure factors, short- and long-time diffusion coefficients, and hydrodynamic interactions of concentrated poly(N-isopropylacryamide) microgel suspensions were measured with simultaneous static and dynamic three-dimensional cross-correlated light scattering. The data are interpreted through comparison to hard sphere theory. The structure factors are known to be described well by the hard sphere approximation. When the structure factor is fit to an effective hard sphere volume fraction and radius, the diffusion and hydrodynamic interactions are also well described by the hard sphere model. We demonstrate that one single hard sphere volume fraction is sufficient to describe the microgel structures, hydrodynamic interactions, and long- and short-time collective diffusion coefficients. This result is surprising because the particle form of the microgels at these temperatures is not rigid, but rather "fuzzy" spheres with dangling polymer chains.     

椭球形粒子扩散动力学的研究进展

椭球形粒子作为一类重要的各向异性粒子,在生物、化工以及材料等领域有着广泛的应用。对椭球形粒子的扩散动力学进行深入地研究,不仅可以增进人们对各向异性粒子扩散动力学中基本科学问题的理解,而且可以为含各向异性粒子材料的设计与加工提供理论依据。因此,椭球形粒子的扩散问题一直是粒子扩散研究的焦点。虽然椭球形粒子扩散动力学的研究已经有一百多年的历史,其相关理论、实验和模拟研究工作很多,并取得了较大的进展,但是,系统的综述性文章相对较少。本文系统地总结了椭球形粒子扩散动力学的相关工作,包括:不同流体力学边界条件下椭球形粒子的扩散动力学和平动与转动扩散的耦合行为等。此外,本文还分析了椭球形粒子扩散动力学领域存在的关键科学问题,并对该领域研究方向进行了展望。     

Mutual diffusion coefficient in binary mixtures in different aggregation states

An equation for the mutual diffusion coefficient for the components of a binary mixture occurring in various aggregation states has been derived. Situations of this type arise in narrow pores in which the adsorbate density sharply changes across the pore cross-section. The equation was derived with the assumption that molecules of the mixture components have a spherical shape and are similar in size. The constructed equations are based on the lattice-gas model applicable to any aggregation state. The migration of particles is described in terms of the transition state theory. At low mixture densities corresponding to an ideal gas phase, the equations are based on expressions of the rigorous kinetic theory of gases. The theory takes into account the change in the mechanism of particle migration in different phases: from pair collisions for the gas to the overcoming of the activation barrier by thermofluctuation for dense phases. __________ Published in Russian in Izvestiya Akademii Nauk. Seriya Khimicheskaya, No. 8, pp. 1717–1725, August, 2005.