Showing papers in "arXiv: Soft Condensed Matter in 2007"

Fingerprinting Soft Materials: A Framework for Characterizing Nonlinear Viscoelasticity

Abstract: We introduce a comprehensive scheme to physically quantify both viscous and elastic rheological nonlinearities simultaneously, using an imposed large amplitude oscillatory shear (LAOS) strain. The new framework naturally lends a physical interpretation to commonly reported Fourier coefficients of the nonlinear stress response. Additionally, we address the ambiguities inherent in the standard definitions of viscoelastic moduli when extended into the nonlinear regime, and define new measures which reveal behavior that is obscured by conventional techniques.

Rheology of giant micelles

Abstract: Giant micelles are elongated, polymer-like objects created by the self-assembly of amphiphilic molecules (such as detergents) in solution. Giant micelles are typically flexible, and can become highly entangled even at modest concentrations. The resulting viscoelastic solutions show fascinating flow behaviour (rheology) which we address theoretically in this article at two levels. First, we summarise advances in understanding linear viscoelastic spectra and steady-state nonlinear flows, based on microscopic constitutive models that combine the physics of polymer entanglement with the reversible kinetics of self-assembly. Such models were first introduced two decades ago, and since then have been shown to explain robustly several distinctive features of the rheology in the strongly entangled regime, including extreme shear-thinning. We then turn to more complex rheological phenomena, particularly involving spatial heterogeneity, spontaneous oscillation, instability, and chaos. Recent understanding of these complex flows is based largely on grossly simplified models which capture in outline just a few pertinent microscopic features, such as coupling between stresses and other order parameters such as concentration. The role of `structural memory' (the dependence of structural parameters such as the micellar length distribution on the flow history) in explaining these highly nonlinear phenomena is addressed. Structural memory also plays an intriguing role in the little-understood shear-thickening regime, which occurs in a concentration regime close to but below the onset of strong entanglement, and which is marked by a shear-induced transformation from an inviscid to a gelatinous state.

Screening-Limited Response of NanoBiosensors

Abstract: Despite tremendous potential of highly sensitive electronic detection of bio-molecules by nanoscale biosensors for genomics and proteomic applications, many aspects of experimentally observed sensor response (S) are unexplained within consistent theoretical frameworks of kinetic response or electrical screening. In this paper, we combine analytic solutions of Poisson-Boltzmann and reaction-diffusion equations to show that the electrical response of nanobiosensor varies logarithmically with the concentration of target molecules, time, the salt concentration, and inversely with the fractal dimension of sensor surface. Our analysis provides a coherent theoretical interpretation of wide variety of puzzling experimental data that have so far defied intuitive explanation.

Precise Particle Tracking Against a Complicated Background: Polynomial Fitting with Gaussian Weight

Abstract: We present a new particle tracking software algorithm designed to accurately track the motion of low-contrast particles against a background with large variations in light levels. The method is based on a polynomial fit of the intensity around each feature point, weighted by a Gaussian function of the distance from the centre, and is especially suitable for tracking endogeneous particles in the cell, imaged with bright field, phase contrast or fluorescence optical microscopy. Furthermore, the method can simultaneously track particles of all different sizes, and allows significant freedom in their shape. The algorithm is evaluated using the quantitative measures of accuracy and precision of previous authors, using simulated images at variable signal-to-noise ratios. To these we add a new test of the error due to a non-uniform background. Finally the tracking of particles in real cell images is demonstrated. The method is made freely available for non-commencial use as a software package with a graphical user-inferface, which can be run within the Matlab programming environment.

Solitary wave trains in granular chains: Experiments, theory and simulations

Abstract: The features of solitary waves observed in horizontal monodisperse chain of barely touching beads not only depend on geometrical and material properties of the beads but also on the initial perturbation provided at the edge of the chain An impact of a large striker on a monodisperse chain, and similarly a sharp decrease of bead radius in a stepped chain, generates a solitary wave train containing many single solitary waves ordered by decreasing amplitudes We find, by simple analytical arguments, that the unloading of compression force at the chain edge has a nearly exponential decrease The characteristic time is mainly a function involving the grains' masses and the striker mass Numerical calculations and experiments corroborate these findings

First-principles Study of Physisorption of Nucleic Acid Bases on Small-Diameter Carbon Nanotube

Abstract: We report the results of our first-principles study based on density functional theory on the interaction of the nucleic acid base molecules adenine (A), cytosine (C), guanine (G), thymine (T), and uracil (U), with a single-walled carbon nanotube (CNT). Specifically, the focus is on the physisorption of base molecules on the outer wall of a (5,0) metallic CNT possessing one of the smallest diameters possible. Compared to CNTs with large diameters, the physisorption energy is found to be reduced in the high-curvature case. The base molecules exhibit significantly different interaction strengths, and the calculated binding energies follow the hierarchy G > A > T > C > U, which appears to be independent of the tube curvature. The stabilizing factor in the interaction between the base molecule and CNT is dominated by the molecular polarizability that allows a weakly attractive dispersion force to be induced between them. The present study provides an improved understanding of the role of the base sequence in deoxyribonucleic acid (DNA) or ribonucleic acid (RNA) on their interactions with carbon nanotubes of varying diameters.

Why do ultrasoft repulsive particles cluster and crystallize? Analytical results from density functional theory

Abstract: We demonstrate the accuracy of the hypernetted chain closure and of the mean-field approximation for the calculation of the fluid-state properties of systems interacting by means of bounded and positive-definite pair potentials with oscillating Fourier transforms. Subsequently, we prove the validity of a bilinear, random-phase density functional for arbitrary inhomogeneous phases of the same systems. On the basis of this functional, we calculate analytically the freezing parameters of the latter. We demonstrate explicitly that the stable crystals feature a lattice constant that is independent of density and whose value is dictated by the position of the negative minimum of the Fourier transform of the pair potential. This property is equivalent with the existence of clusters, whose population scales proportionally to the density. We establish that regardless of the form of the interaction potential and of the location on the freezing line, all cluster crystals have a universal Lindemann ratio L = 0.189 at freezing. We further make an explicit link between the aforementioned density functional and the harmonic theory of crystals. This allows us to establish an equivalence between the emergence of clusters and the existence of negative Fourier components of the interaction potential. Finally, we make a connection between the class of models at hand and the system of infinite-dimensional hard spheres, when the limits of interaction steepness and space dimension are both taken to infinity in a particularly described fashion.

Quasi-binary amorphous phase in a 3D system of particles with repulsive-shoulder interactions

Abstract: We report a computer-simulation study of the equilibrium phase diagram of a three-dimensional system of particles with a repulsive step potential. Using free-energy calculations, we have determined the equilibrium phase diagram of this system. At low temperatures, we observe a number of distinct crystal phases. However, under certain conditions the system undergoes a glass transition in a regime where the liquid appears thermodynamically stable. We argue that the appearance of this amorphous low-temperature phase can be understood by viewing this one-component system as a pseudo-binary mixture.

Swinging and tumbling of elastic capsules in shear flow

Abstract: The deformation of an elastic micro-capsule in an infinite shear flow is studied numerically using a spectral method. The shape of the capsule and the hydrodynamic flow field are expanded into smooth basis functions. Analytic expressions for the derivative of the basis functions permit the evaluation of elastic and hydrodynamic stresses and bending forces at specified grid points in the membrane. Compared to methods employing a triangulation scheme, this method has the advantage that the resulting capsule shapes are automatically smooth, and few modes are needed to describe the deformation accurately. Computations are performed for capsules both with spherical and ellipsoidal unstressed reference shape. Results for small deformations of initially spherical capsules coincide with analytic predictions. For initially ellipsoidal capsules, recent approximative theories predict stable oscillations of the tank-treading inclination angle, and a transition to tumbling at low shear rate. Both phenomena have also been observed experimentally. Using our numerical approach we could reproduce both the oscillations and the transition to tumbling. The full phase diagram for varying shear rate and viscosity ratio is explored. While the numerically obtained phase diagram qualitatively agrees with the theory, intermittent behaviour could not be observed within our simulation time. Our results suggest that initial tumbling motion is only transient in this region of the phase diagram.

Ultra-sensitive surface absorption spectroscopy using sub-wavelength diameter optical fibers

Abstract: The guided modes of sub-wavelength diameter air-clad optical fibers exhibit a pronounced evanescent field. The absorption of particles on the fiber surface is therefore readily detected via the fiber transmission. We show that the resulting absorption for a given surface coverage can be orders of magnitude higher than for conventional surface spectroscopy. As a demonstration, we present measurements on sub-monolayers of 3,4,9,10-perylene-tetracarboxylic dianhydride (PTCDA) molecules at ambient conditions, revealing the agglomeration dynamics on a second to minutes timescale.

Simulation of Heme using DFT+U: a step toward accurate spin-state energetics

Abstract: We investigate the DFT+U approach as a viable solution to describe the low-lying states of ligated and unligated iron heme complexes. Besides their central role in organometallic chemistry, these compounds represent a paradigmatic case where LDA, GGA, and common hybrid functionals fail to reproduce the experimental magnetic splittings. In particular, the imidazole pentacoordinated heme is incorrectly described as a triplet by all usual DFT flavors. In this study we show that a U parameter close to 4 eV leads to spin transitions and molecular geometries in quantitative agreement with experiments, and that DFT+U represents an appealing tool in the description of iron porphyrin complexes, at a much reduced cost compared to correlated quantum-chemistry methods. The possibility of obtaining the U parameter from first-principles is explored through a self-consistent linear-response formulation. We find that this approach, which proved to be successful in other iron systems, produces in this case some overestimation with respect to the optimal values of U.

Rheological Behavior of Aqueous Suspensions of Laponite: New Insights into the Ageing Phenomena

Abstract: In this paper, ageing behavior of suspensions of laponite with varying salt concentration is investigated using rheological tools. It is observed that the ageing is accompanied by an increase in the complex viscosity. The succeeding creep experiments performed at various ages showed damped oscillations in the strain. The characteristic time-scale of the damped oscillations, retardation time, showed a prominent decrease with the age of the system. However, this dependence weakens with an increase in the salt concentration, which is known to change microstructure of the system from glass-like to gel-like. We postulate that a decrease in the retardation time can be represented as a decrease in the viscosity (friction) of the dissipative environment surrounding the arrested entities that opposes elastic deformation of the system. We believe that ageing in colloidal glass leads to a greater ordering that enhances relative spacing between the constituents thereby reducing the frictional resistance. However, since a gel state is inherently different in structure (fractal network) than that of a glass (disordered), ageing in the same does not induce ordering. Consequently, we observe inverse dependence of retardation time on age becoming weaker with an increase in the salt concentration. We analyze these results from a perspective of ageing dynamics of both glass state and gel state of laponite suspensions.

Medium dependence of asphaltene agglomeration inhibitor efficiency

Abstract: Applying chemical additives (molecule inhibitors or dispersants) is one of the common ways to control asphaltene agglomeration and precipitation. However, it is not clear why at some conditions the synthetic flocculation inhibitors as well as resins not only do not inhibit the asphaltene agglomeration,, they may also promote it, and why the increasing of the additive concentration may lead to the diminishing of their efficacy. To clarify this issue, in the present work we have performed a set of vapor preassure osmometry experiments investigating the asphaltene agglomeration inhibition by commercial and new inhibitor molecules in toluene and o-diclorobenzene. Monte Carlo computer modeling has been applied to interpret some unexpected trends of molar mass of the Puerto Ceiba asphaltene clusters at different concentrations of inhibitor, assuming that inhibitors efficiency is directly related to their adsorption on the surface of asphaltene or its complexes. It has been found that a self-assembly of inhibitor molecules, induced by relative lyophilic or lyophobic interactions, may be a reason of the inhibitor efficacy declining.

Pore-blockade Times for Field-Driven Polymer Translocation

Abstract: We study pore blockade times for a translocating polymer of length $N$, driven by a field $E$ across the pore in three dimensions. The polymer performs Rouse dynamics, i.e., we consider polymer dynamics in the absence of hydrodynamical interactions. We find that the typical time the pore remains blocked during a translocation event scales as $\sim N^{(1+2 u)/(1+ u)}/E$, where $ u\simeq0.588$ is the Flory exponent for the polymer. In line with our previous work, we show that this scaling behaviour stems from the polymer dynamics at the immediate vicinity of the pore -- in particular, the memory effects in the polymer chain tension imbalance across the pore. This result, along with the numerical results by several other groups, violates the lower bound $\sim N^{1+ u}/E$ suggested earlier in the literature. We discuss why this lower bound is incorrect and show, based on conservation of energy, that the correct lower bound for the pore-blockade time for field-driven translocation is given by $\eta N^{2 u}/E$, where $\eta$ is the viscosity of the medium surrounding the polymer.

Toward the zero surface tension limit: The granular fingering instability

Abstract: The finger-like branching pattern that occurs when a less viscous fluid displaces a more viscous one confined between two parallel plates has been widely studied as a classical example of a mathematically-tractable hydrodynamic instability since the time of Saffman and Taylor. Fingering in such Hele-Shaw geometries have been generated not only with Newtonian fluids but also with various non-Newtonian fluids including fine granular material displaced by gas, liquid, or large grains. Here we study a granular Hele-Shaw system to explore whether the absence of cohesive forces in dry granular material can produce an ideal venue for studying the hitherto-unrealizable singular hydrodynamics predicted in the zero-surface-tension limit. We demonstrate that the grain-gas interface does indeed exhibit fractal structure and sharp cusps associated with finite-time singularities. Above the yield stress, the scaling for the finger width is distinct from that for ordinary fluids, reflecting unique granular properties such as friction-induced dissipation as opposed to viscous damping. Despite such differences, the dimension of the global fractal structure and the shape of the singular cusps on the interface agree with the theories based on simple Laplacian growth of conventional fluid fingering in the zero-surface-tension limit. Our study provides new insights not only on the dynamics of two-phase dense granular flows and granular pattern formation, but also on the fluid dynamics in the zero-surface-tension limit.

On the possibility of extending the Nore-Frenkel generalized law of correspondent states to non-isotropic patchy interactions

Abstract: Colloidal systems (and protein solutions) are often characterized by attractive interactions whose range is much smaller than the particle size. When this is the case and the interaction is spherical, systems obey a generalized law of correspondent states (GLCS), first proposed by Noro and Frenkel [ J.Chem.Phys. 113, 2941 (2000) ]. The thermodynamic properties become insensitive to the details of the potential, depending only on the value of the second virial coefficient B_2 and the density $\rho$. The GLCS does not generically hold for the case of non-spherical potentials. In this Letter we suggest that when particles interact via short-ranged small-angular amplitude patchy interactions (so that the condition of only one bond per patch is fulfilled) it is still possible to generalize the GLCS close to the liquid-gas critical point. Keywords: Colloids, Second Virial Coefficient, Proteins interactions, Short-ranged attractive attractions.

Aging, Fragility and Reversibility Window in Bulk Alloy Glasses

Abstract: Non-reversing relaxation enthalpies (DHnr) at glass transitions Tg(x) in the PxGexSe1-2x ternary display a wide, sharp and deep global minimum (~0) in the 0.09 < x < 0.145 range, within which Tg becomes thermally reversing. In the reversibility window these glasses are found not to age, in contrast to aging observed for fragile glass compositions outside the window. Thermal reversibility and lack of aging are paradigms that molecular glasses in the window share with proteins in transition states, which result from structural self-organization in both systems. In proteins the self-organized structures appear to be at places where life sustaining repeating foldings and unfoldings occur.

Model of nanocrystal formation in solution by burst nucleation and diffusional growth

Abstract: The phenomenon of burst nucleation in solution, in which a period of apparent chemical inactivity is followed by a sudden and explosive growth of nucleated particles from a solute species, has been given a widely accepted qualitative explanation by LaMer and co-workers. Here, we present a model with the assumptions of instantaneous thermalization below the critical nucleus size and irreversible diffusive growth above the critical size, which for the first time formulates LaMer's explanation of burst nucleation in a manner allowing quantitative calculations. The behavior of the model at large times is derived, with the result that the average cluster size, as measured by the number of atoms, grows linearly with time, while the width of the cluster distribution grows as the square root of time. We develop an effective numerical scheme to integrate the equations of the model and compare the asymptotic expressions to results from numerical simulation. Finally, we discuss the physical effects which cause real nucleation processes in solution to deviate from the behavior of the model.

Kohn anomalies and non-adiabaticity in doped carbon nanotubes

Abstract: The high-frequency Raman-active phonon modes of metallic single-walled carbon nanotubes (SWNTs) are thought to be characterized by Kohn anomalies (KAs), which are expected to be modified by the doping-induced tuning of the Fermi energy level $\epsilon_F$, obtained through the intercalation of SWNTs with alkali atoms or by the application of a gate potential. We present a Density-Functional Theory (DFT) study of the phonon properties of a (9,9) metallic SWNT as a function of electronic doping. For such study, we use, as in standard DFT calculations of vibrational properties, the Born-Oppenheimer (BO) approximation. We also develop an analytical model capable of reproducing and interpreting our DFT results. Both DFT calculations and this model predict, for increasing doping levels, a series of EPC-induced KAs in the vibrational mode parallel to the tube axis at the $\mathbf\Gamma$ point of the Brillouin zone, usually indicated in Raman spectroscopy as the $G^-$ peak. Such KAs would arise each time a new conduction band is populated. However, we show that they are an artifact of the BO approximation. The inclusion of non-adiabatic (NA) effects dramatically affects the results, predicting KAs at $\mathbf\Gamma$ only when $\epsilon_F$ is close to a band crossing $E_{X}$. For each band crossing a double KA occurs for $\epsilon_F=E_{X}\pm \hbar\omega/2$, where $\hbar\omega$ is the phonon energy. In particular, for a 1.2 $nm$ metallic nanotube, we predict a KA to occur in the so-called $G^-$ peak at a doping level of about $N_{el}/C=\pm 0.0015$ atom ($\epsilon_F\approx \pm 0.1 ~eV$). Furthermore, we predict that the Raman linewidth of the $G^-$ peak significantly decreases for $|\epsilon_F| \geq \hbar\omega/2$.

Designing phoretic micro- and nano-swimmers

Abstract: Small objects can swim by generating around them fields or gradients which in turn induce fluid motion past their surface by phoretic surface effects. We quantify for arbitrary swimmer shapes and surface patterns, how efficient swimming requires both surface ``activity'' to generate the fields, and surface ``phoretic mobility.'' We show in particular that (i) swimming requires symmetry breaking in either or both of the patterns of "activity" and ``mobility,'' and (ii) for a given geometrical shape and surface pattern, the swimming velocity is size-independent. In addition, for given available surface properties, our calculation framework provides a guide for optimizing the design of swimmers.

Protein-DNA interactions: reaching and recognizing the targets

Abstract: This paper has been withdrawn because of the requirement from the journal where the modified version of the paper is to be published.

Phase diagram of model anisotropic particles with octahedral symmetry

Abstract: We computed the phase diagram for a system of model anisotropic particles with six attractive patches in an octahedral arrangement. We chose to study this model for a relatively narrow value of the patch width where the lowest-energy configuration of the system is a simple cubic crystal. At this value of the patch width, there is no stable vapour-liquid phase separation, and there are three other crystalline phases in addition to the simple cubic crystal that is most stable at low pressure. Firstly, at moderate pressures, it is more favourable to form a body-centred cubic crystal, which can be viewed as two interpenetrating, and almost non-interacting, simple cubic lattices.Secondly, at high pressures and low temperatures, an orientationally ordered face-centred cubic structure becomes favourable. Finally, at high temperatures a face-centred cubic plastic crystal is the most stable solid phase.

Non-linear mechanical response of the Red Blood Cell

Abstract: We measure the dynamical mechanical properties of human red blood cells. Single cell response is measured with optical tweezers. We investigate both the stress relaxation following a fast deformation, and the effect of varying the strain rate. We find a power law decay of the stress as a function of time, down to a plateau stress, and a power law increase of the cell's elasticity as a function of the strain rate. Interestingly, the exponents of these quantities violate the linear superposition principle, indicating a nonlinear response. We propose that this is due to breaking of a fraction of the crosslinks during the deformation process. The Soft Glassy Rheology Model accounts for the relation between the exponents we observe experimentally. This picture is consistent with recent models of bond remodeling in the red blood cell's molecular structure. Our results imply that the blood cell's mechanical behavior depends critically on the deformation process.

A multi-blob representation of semi-dilute polymer solutions

Abstract: A coarse-grained multi-blob description of polymer solutions is presented, based on soft, transferable effective interactions between bonded and non-bonded blobs. The number of blobs is chosen such that the blob density does not exceed their overlap threshold, allowing polymer concentrations to be explored deep into the semi-dilute regime. This quantitative multi-blob description is shown to preserve known scaling laws of polymer solutions and provides accurate estimates of amplitudes, while leading to orders of magnitude increase of simulation efficiency and allowing analytic calculations of structural and thermodynamic properties. Many conformational, structural and thermodynamic properties of semi-dilute polymer solutions, both in the bulk and under confinement, can be qualitatively understood in terms of scaling arguments based on the de Gennes-Pincus “blob” picture [1]. This picture is of value whenever the characteristic length scale of the polymer solution (e.g. the correlation length ξ), or of the confinement is significantly shorter than the radius of gyration Rg of an isolated polymer chain. The blob picture suggests a systematic coarse-graining procedure whereby each polymer chain is divided into a number n of blobs, each containing the same number of monomers or segments of the initial coil, such that blobs of the same or different chains do not, on average, overlap. In this letter we present and validate a quantitative formulation of such a multi-blob representation, which allows a popular single blob coarse graining procedure [2, 3, 4, 5, 6] to be extended to highly concentrated solutions. Consider a solution of N self avoiding polymer chains of L segments (each of size b) in a volume V ; the polymer number density is ρ = N/V , and if Rg ∼ bL ν (with ν ≃ 0.588 the Flory exponent) is the radius of gyration of an isolated chain, the polymer overlap density is ρ � = 3/ 4πR 3 � In the dilute regime, ρ � < ρ, where polymers do not, on average, overlap, they may be represented by a single blob of radius Rg; the effective interaction potential v(r) between the centres of mass (CM) of two blobs can be calculated by averaging over monomer conformations for a given distance r between their CM’s, e.g. by Monte Carlo (MC) simulations of an isolated pair of polymers [2, 3, 4, 5, 6]. The resulting v(r) depends weakly on polymer length L, and in the scaling limit (L → ∞), it is accurately represented by a Gaussian of width ∼ Rg [6]: v(r) kBT ≃ Aexp � −α(r/R g) 2

The Glassy Wormlike Chain

Abstract: We introduce a new model for the dynamics of a wormlike chain in an environment that gives rise to a rough free energy landscape, which we baptise the glassy wormlike chain. It is obtained from the common wormlike chain by an exponential stretching of the relaxation spectrum of its long-wavelength eigenmodes, controlled by a single stretching parameter. Predictions for pertinent observables such as the dynamic structure factor and the microrheological susceptibility exhibit the characteristics of soft glassy rheology and compare favourably with experimental data for reconstituted cytoskeletal networks and live cells. We speculate about the possible microscopic origin of the stretching, implications for the nonlinear rheology, and the potential physiological significance of our results.

CARS polarized microscopy of three-dimensional director structures in liquid crystals

Abstract: We demonstrate three-dimensional vibrational imaging of director structures in liquid crystals using coherent anti-Stokes Raman scattering (CARS) polarized microscopy. Spatial mapping of the structures is based on sensitivity of a polarized CARS signal to orientation of anisotropic molecules in liquid crystals. As an example, we study structures in a smectic material and demonstrate that single-scan CARS and two-photon fluorescence images of molecular orientation patterns are consistent with each other and with the structure model.

Gel-forming patchy colloids and network glass formers: Thermodynamic and Dynamic analogies

Abstract: This article discusses recent attempts to provide a deeper understanding of the thermoreversible "gel" state of colloidal matter and to unravel the analogies between gels at the colloidal level and gels at the molecular level, commonly known as network-forming strong liquids. The connection between gel-forming patchy colloids and strong liquids is provided by the limited valence of the inter-particle interactions, i.e. by the presence of a limit in the number of bonded nearest neighbors.

Linear response and stability of ordered phases of diblock copolymer melts

Abstract: An efficient pseudo-spectral numerical method is introduced for calculating a self-consistent field (SCF) approximation for the linear susceptibility of ordered phases in block copolymer melts (sometimes referred to as the random phase approximation). Our method is significantly more efficient than that used in the first calculations of this quantity by Shi, Laradji and coworkers, allowing for the study of more strongly segregated structures. We have re-examined the stability of several phases of diblock copolymer melts, and find that some conclusions of Laradji et al. regarding the stability of the Gyroid phase were the result of insufficient spatial resolution. We find that an epitaxial (k=0) instability of the Gyroid phase with respect to the hexagonal phase that was considered previously by Matsen competes extremely closely with an instability that occurs at a nonzero crystal wavevector k.

A dynamic density functional theory for particles in a flowing solvent

Abstract: We present a dynamic density functional theory (dDFT) which takes into accou nt the advection of the particles by a flowing solvent. For potential flows we can use the same closure as in the absence of solvent flow. The structure of the resulting advected dDFT suggests that it could be used for non-potential flows as well. We apply this dDFT to Brownian particles (e.g., polymer coils) in a solvent flowing around a spherical obstacle (e.g., a colloid) and compare the results with direct simulations of the underlying Brownian dynamics. Although numerical limitations do not allow for an accurate quantitative check of the advected dDFT both show the same qualitative features. In contrast to previous works which neglected the deformation of the flow by the obstacle, we find that the bow-wave in the density distribution of particles in front of the obstacle as well as the wake behind it are reduced dramatically. As a consequence the friction force exerted by the (polymer) particles on the colloid can be reduced drastically.

Random close packing of granular matter

Abstract: We propose an interpretation of the random close packing of granular materials as a phase transition, and discuss the possibility of experimental verification.