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

The Charge of Glass and Silica Surfaces

Abstract: We present a method of calculating the electric charge density of glass and silica surfaces in contact with aqueous electrolytes for two cases of practical relevance that are not amenable to standard techniques: surfaces of low specific area at low ionic strength and surfaces interacting strongly with a second anionic surface.

Molecular Weight Dependence of Polymersome Membrane Elasticity and Stability

Abstract: Vesicles prepared in water from a series of diblock copolymers and termed "polymersomes" are physically characterized. With increasing molecular weight $\bar{M}_n$, the hydrophobic core thickness $d$ for the self-assembled bilayers of polyethyleneoxide - polybutadiene (PEO-PBD) increases up to 20 $nm$ - considerably greater than any previously studied lipid system. The mechanical responses of these membranes, specifically, the area elastic modulus $K_a$ and maximal areal strain $\alpha_c$ are measured by micromanipulation. As expected for interface-dominated elasticity, $K_a$ ($\simeq$ 100 $pN/nm$) is found to be independent of $\bar{M}_n$. Related mean-field ideas also predict a limiting value for $\alpha_c$ which is universal and about 10-fold above that typical of lipids. Experiments indeed show $\alpha_c$ generally increases with $\bar{M}_n$, coming close to the theoretical limit before stress relaxation is opposed by what might be chain entanglements at the highest $\bar{M}_n$. The results highlight the interfacial limits of self-assemblies at the nano-scale.

Friction law for dense granular flows: application to the motion of a mass down a rough inclined plane

Abstract: The problem of the spreading of a granular mass released at the top of a rough inclined plane was investigated. We experimentally measure the evolution of the avalanche from the initiation up to the deposit using a Moire image processing technique. The results are quantitatively compared with the prediction of an hydrodynamic model based on depth averaged equations. In the model, the interaction between the flowing layer and the rough bottom is described by a non trivial friction force whose expression is derived from measurements on steady uniform flows. We show that the spreading of the mass is quantitatively predicted by the model when the mass is released on a plane free of particles. When an avalanche is triggered on an initially static layer, the model fails in quantitatively predicting the propagation but qualitatively captures the evolution.

On fluid/wall slippage

Abstract: Certain (non polymeric) fluids show an anomalously low friction when flowing against well chosen solid walls. We discuss here one possible explanation, postulating that a gaseous film of small thickness h is present between fluid and wall. When h is smaller than the mean free path l of the gas (Knudsen regime) the Navier length b is expected to be independent of h and very large (microns).

Evolution of the social network of scientific collaborations

Abstract: The co-authorship network of scientists represents a prototype of complex evolving networks. By mapping the electronic database containing all relevant journals in mathematics and neuro-science for an eight-year period (1991-98), we infer the dynamic and the structural mechanisms that govern the evolution and topology of this complex system. First, empirical measurements allow us to uncover the topological measures that characterize the network at a given moment, as well as the time evolution of these quantities. The results indicate that the network is scale-free, and that the network evolution is governed by preferential attachment, affecting both internal and external links. However, in contrast with most model predictions the average degree increases in time, and the node separation decreases. Second, we propose a simple model that captures the network's time evolution. Third, numerical simulations are used to uncover the behavior of quantities that could not be predicted analytically.

The large-scale structure of semantic networks: statistical analyses and a model for semantic growth

Abstract: We present statistical analyses of the large-scale structure of three types of semantic networks: word associations, WordNet, and Roget's thesaurus. We show that they have a small-world structure, characterized by sparse connectivity, short average path-lengths between words, and strong local clustering. In addition, the distributions of the number of connections follow power laws that indicate a scale-free pattern of connectivity, with most nodes having relatively few connections joined together through a small number of hubs with many connections. These regularities have also been found in certain other complex natural networks, such as the world wide web, but they are not consistent with many conventional models of semantic organization, based on inheritance hierarchies, arbitrarily structured networks, or high-dimensional vector spaces. We propose that these structures reflect the mechanisms by which semantic networks grow. We describe a simple model for semantic growth, in which each new word or concept is connected to an existing network by differentiating the connectivity pattern of an existing node. This model generates appropriate small-world statistics and power-law connectivity distributions, and also suggests one possible mechanistic basis for the effects of learning history variables (age-of-acquisition, usage frequency) on behavioral performance in semantic processing tasks.

Echinocyte Shapes: Bending, Stretching and Shear Determine Spicule Shape and Spacing

Abstract: We study the shapes of human red blood cells using continuum mechanics. In particular, we model the crenated, echinocytic shapes and show how they may arise from a competition between the bending energy of the plasma membrane and the stretching/shear elastic energies of the membrane skeleton. In contrast to earlier work, we calculate spicule shapes exactly by solving the equations of continuum mechanics subject to appropriate boundary conditions. A simple scaling analysis of this competition reveals an elastic length which sets the length scale for the spicules and is, thus, related to the number of spicules experimentally observed on the fully developed echinocyte.

Condensed phases of gases inside nanotube bundles

Abstract: An overview is presented of the various phases predicted to occur when gases are absorbed within a bundle of carbon nanotubes. The behavior may be characterized by an effective dimensionality, which depends on the species and the temperature. Small molecules are strongly attracted to the interstitial channels between tubes. There, they undergo transitions between ordered and disordered quasi-one dimensional (1D) phases. Both small and large molecules display 1D and /or 2D phase behavior when adsorbed within the nanotubes, depending on the species and thermodynamic conditions. Finally, molecules adsorbed on the external surface of the bundle exhibit 1D behavior (striped phases), which crosses over to 2D behavior (monolayer film) and eventually 3D behavior (thick film) as the coverage is increased. The various phases exhibit a wide variety of thermal and other properties that we discuss here.

Assessment of exchange-correlation functionals for the calculation of dynamical properties of small clusters in TDDFT

Abstract: We present a detailed study of different exchange-correlation (xc) functionals in describing the dynamical properties of finite systems. For that purpose, we calculated the static polarizabilities, ionization potentials and optical absorption spectrum of four small clusters, Na_2, Na_4, SiH_4 and Si_2H_6, using a real-space, real-time technique. The computed static polarizabilities and ionization potentials seem to be in rather good agreement with the available experimental data, once the proper asymptotics of the potential are taken into account. The same conclusion holds for the absorption spectra, although the xc kernels in use do not provide a sufficiently strong attractive interaction between electrons and holes, leading to spectra slightly shifted towards higher energies. This deficiency is traced back to the insufficient description of dynamical effects in the correlation functional. Furthermore, it is shown that the xc potential used to obtain the ground state is the key factor to get reasonable spectra, whereas the choice of the xc kernel just amounts to small, although important, quantitative changes.

Bose-Einstein Condensation of Trapped Atomic Gases

Abstract: This article reviews recent investigations on the phenomenon of Bose-Einstein condensation of dilute gases. Since the experimental observation of quantum degeneracy in atomic gases, the research activity in the field of coherent matter-waves literally exploded. The present topical review aims to give an introduction into the thermodynamics of Bose-Einstein condensation, a general overview over experimental techniques and investigations, and a theoretical foundation for the description of bosonic many-body quantum systems.

How does the chain extension of poly (acrylic acid) scale in aqueous solution? A combined study with light scattering and computer simulation

Abstract: This work adresses the question of the scaling behaviour of polyelectrolytes in solution for a realistic prototype: We show results of a combined experimental (light scattering) and theoretical (computer simulations) investigation of structural properties of poly (acrylic acid) (PAA). Experimentally, we determined the molecular weight (M_W) and the hydrodynamic radius (R_H) by static light scattering for six different PAA samples in aqueous NaCl-containing solution (0.1-1 mol/L) of polydispersity D_P between 1.5 and 1.8. On the computational side, three different variants of a newly developed mesoscopic force field for PAA were employed to determine R_H for monodisperse systems of the same M_W as in the experiments. The force field effectively incorporates atomistic information and one coarse-grained bead corresponds to one PAA monomer. We find that R_H matches with the experimental data for all investigated samples. The effective scaling exponent for R_H is found to be around 0.55, which is well below its asymptotic value for good solvents. Additionally, data for the radius of gyration (R_G) are presented.

Overcharging: The Crucial Role of Excluded Volume

Abstract: In this Letter we investigate the mechanism for overcharging of a single spherical colloid in the presence of aqueous salts within the framework of the primitive model by molecular dynamics (MD) simulations as well as integral-equation theory. We find that the occurrence and strength of overcharging strongly depends on the salt-ion size, and the available volume in the fluid. To understand the role of the excluded volume of the microions, we first consider an uncharged system. For a fixed bulk concentration we find that upon increasing the fluid particle size one strongly increases the local concentration nearby the colloidal surface and that the particles become laterally ordered. For a charged system the first surface layer is built up predominantly by strongly correlated counterions. We argue that this a key mechanism to produce overcharging with a low electrostatic coupling, and as a more practical consequence, to account for charge inversion with monovalent aqueous salt ions.

Complexation between oppositely charged polyelectrolytes: beyond the Random Phase Approximation

Abstract: We consider the phase behavior of polymeric systems by calculating the structure factors beyond the Random Phase Approximation. The effect of this correction to the mean-field RPA structure factor is shown to be important in the case of coulombic systems. Two examples are given: simple electrolytes and mixtures of incompatible oppositely charged polyelectrolytes. In this last case, all former studies predicted an enhancement of compatibility for increasing charge densities; we also describe the complexation transition between the polyelectrolytes. We determine a phase diagram of the polyelectrolyte mixture that includes both complexation and incompatibility.

Lees-Edwards boundary conditions for lattice Boltzmann

Abstract: Lees Edwards boundary conditions (LEbc) for Molecular Dynamics simulations are an extension of the well known periodic boundary conditions and allow the simulation of bulk systems in a simple shear flow. We show how the idea of LEbc can be implemented in lattice Boltzmann simulations and how LEbc can be used to overcome the problem of a maxinum shear rate that is limited to less than 1/Ly (with Ly the transverse system size) in traditional lattice Boltzmann implementations of shear flow.

Folding of a Small Helical Protein Using Hydrogen Bonds and Hydrophobicity Forces

Abstract: A reduced protein model with five to six atoms per amino acid and five amino acid types is developed and tested on a three-helix-bundle protein, a 46-amino acid fragment from staphylococcal protein A. The model does not rely on the widely used Go approximation where non-native interactions are ignored. We find that the collapse transition is considerably more abrupt for the protein A sequence than for random sequences with the same composition. The chain collapse is found to be at least as fast as helix formation. Energy minimization restricted to the thermodynamically favored topology gives a structure that has a root-mean-square deviation of 1.8 A from the native structure. The sequence-dependent part of our potential is pairwise additive. Our calculations suggest that fine-tuning this potential by parameter optimization is of limited use.

New mechanism of membrane fusion

Abstract: We have carried out Monte Carlo simulation of the fusion of bilayers of single chain amphiphiles which show phase behavior similar to that of biological lipids. The fusion mechanism we observe is very different from the ``stalk'' hypothesis. Stalks do form on the first stage of fusion, but they do not grow radially to form a hemifused state. Instead, stalk formation destabilizes the membranes and results in hole formation in the vicinity of the stalks. When holes in each bilayer nucleate spontaneously next to the same stalk, an incomplete fusion pore is formed. The fusion process is completed by propagation of the initial connection, the stalk, along the edges of the aligned holes.

Collapse Dynamics of a Homopolymer: Theory and Simulation

Abstract: We present a scaling theory describing the collapse of a homopolymer chain in poor solvent. At time t after the beginning of the collapse, the original Gaussian chain of length N is streamlined to form N/g segments of length R(t), each containing g ~ t monomers. These segments are statistical quantities representing cylinders of length R ~ t^{1/2} and diameter d ~ t^{1/4}, but structured out of stretched arrays of spherical globules. This prescription incorporates the capillary instability. We compare the time-dependent structure factor derived for our theory with that obtained from ultra-large-scale molecular dynamics simulation with explicit solvent. This is the first time such a detailed comparison of theoretical and simulation predictions of collapsing chain structure has been attempted. The favorable agreement between the theoretical and computed structure factors supports the picture of the coarse-graining process during polymer collapse.

Spinor condensates and light scattering from bose-einstein condensates

Abstract: These notes discuss two aspects of the physics of atomic Bose-Einstein condensates: optical properties and spinor condensates. The first topic includes light scattering experiments which probe the excitations of a condensate in both the free-particle and phonon regime. At higher light intensity, a new form of superradiance and phase-coherent matter wave amplification were observed. We also discuss properties of spinor condensates and describe studies of ground-state spin domain structures and dynamical studies which revealed metastable excited states and quantum tunneling.

Density profiles and surface tensions of polymers near colloidal surfaces

Abstract: The surface tension of interacting polymers in a good solvent is calculated theoretically and by computer simulations for a planar wall geometry and for the insertion of a single colloidal hard-sphere. This is achieved for the planar wall and for the larger spheres by an adsorption method, and for smaller spheres by a direct insertion technique. Results for the dilute and semi-dilute regimes are compared to results for ideal polymers, the Asakura-Oosawa penetrable-sphere model, and to integral equations, scaling and renormalization group theories. The largest relative changes with density are found in the dilute regime, so that theories based on non-interacting polymers rapidly break down. A recently developed ``soft colloid'' approach to polymer-colloid mixtures is shown to correctly describe the one-body insertion free-energy and the related surface tension.

Bending Frustration of Lipid-Water Mesophases Based on Cubic Minimal Surfaces

Abstract: Inverse bicontinuous cubic phases are ubiquitous in lipid-water mixtures and consist of a lipid bilayer forming a cubic minimal surface, thereby dividing space into two cubic networks of water channels. For small hydrocarbon chain lengths, the monolayers can be modeled as parallel surfaces to a minimal midsurface. The bending energy of the cubic phases is determined by the distribution of Gaussian curvature over the minimal midsurfaces which we calculate for seven different structures (G, D, P, I-WP, C(P), S and F-RD). We show that the free-energy densities of the structures G, D and P are considerably lower than those of the other investigated structures due to their narrow distribution of Gaussian curvature. The Bonnet transformation between G, D, and P implies that these phases coexist along a triple line, which also includes an excess water phase. Our model includes thermal membrane undulations. Our qualitative predictions remain unchanged when higher order terms in the curvature energy are included. Calculated phase diagrams agree well with the experimental results for 2:1 lauric acid/dilauroyl phosphatidylcholine and water.

Giant Charge Inversion of a Macroion Due to Multivalent Counterions and Monovalent Coions: Molecular Dynamics Study

Abstract: We report molecular dynamics simulation of the (overall neutral) system consisting of an immobile macroion surrounded by the electrolyte of multivalent counterions and monovalent coions. As expected theoretically, counterions adsorb on the macroion surface in the amount much exceeding neutralization requirement, thus effectively inverting the sign of the macroion charge. We find two conditions necessary for charge inversion, namely, counterions must be multivalently charged and Coulomb interactions must be strong enough compared to thermal energy. On the other hand, coion condensation on the multivalent counterions similar to Bjerrum pairing is the major factor restricting the amount of charge inversion. Depending on parameters, we observe inverted charge up to about 200% the original charge of the macroion in absolute value. Our findings are consistent with the theory of "giant charge inversion" [Phys.Rev.Lett., v. 85, 1568 (2000)].

Thermal Folding and Mechanical Unfolding Pathways of Protein Secondary Structures

Abstract: Mechanical stretching of secondary structures is studied through molecular dynamics simulations of a Go-like model. Force vs. displacement curves are studied as a function of the stiffness and velocity of the pulling device. The succession of stretching events, as measured by the order in which contacts are ruptured, is compared to the sequencing of events during thermal folding and unfolding. Opposite cross-correlations are found for an $\alpha$-helix and a $\beta$-hairpin structure. In a tandem of two $\alpha$-helices, the two constituent helices unravel nearly simultaneously. A simple condition for simultaneous vs. sequential unraveling of repeat units is presented.

Cell model and Poisson-Boltzmann theory: A brief introduction

Abstract: The cell-model and its treatment on the Poisson-Boltzmann level are two important concepts in the theoretical description of charged macromolecules. In this brief contribution we provide an introduction to both ideas and summarize a few important results which can be obtained from them. Our article is organized as follows: Section 1 outlines the sequence of approximations which ultimately lead to the cell-model. Section 2 is devoted to two exact results, namely, an expression for the osmotic pressure and a formula for the ion density at the surface of the macromolecule, known as the contact value theorem. Section 3 provides a derivation of the Poisson-Boltzmann equation from a variational principle and the assumption of a product state. Section 4 applies Poisson-Boltzmann theory to the cell-model of linear polyelectrolytes. In particular, the behavior of the exact solution in the limit of zero density is compared to the concept of Manning condensation. Finally, Section 5 shows how a system described by a cell model can be coupled to a salt reservoir, i.e., how the so-called Donnan equilibrium is established.

Phase Transition in Strongly Degenerate Hydrogen Plasma

Abstract: Direct fermionic path-integral Monte-Carlo simulations of strongly coupled hydrogen are presented. Our results show evidence for the hypothetical plasma phase transition. Its most remarkable manifestation is the appearance of metallic droplets which are predicted to be crucial for the electrical conductivity allowing to explain the rapid increase observed in recent shock compression measurments.

Osmotic Stabilisation of Concentrated Emulsions and Foams

Abstract: In the absence of coalescence, coarsening of emulsions (and foams) is controlled by molecular diffusion of dispersed phase between droplets/bubbles. Studies of dilute emulsions have shown how the osmotic pressure of a trapped species within droplets can ``osmotically stabilise'' the emulsion. Webster and Cates (Langmuir, 1998, 14, 2068-2079) gave rigorous criteria for osmotic stabilisation of mono- and polydisperse emulsions, in the dilute regime. We consider here whether analogous criteria exist for the osmotic stabilisation of mono- and polydisperse concentrated emulsions and foams, and suggest that the pressure differences driving coarsening are small compared to the mean Laplace pressure. An exact calculation confirms this for a monodisperse 2D model, finding a bubble's pressure as P_i = P + Pi + P_i^G, with P, Pi the atmospheric and osmotic pressures, and P_i^G a ``geometric pressure'' that reduces to the Laplace pressure only for a spherical bubble. For Princen's 2D emulsion model, P_i^G is only 5% larger in the dry limit than the dilute limit. We conclude that osmotic stabilisation of dense systems typically requires a pressure of trapped molecules in each droplet that is comparable to the Laplace pressures the same droplets would have if spherical, as opposed to the much larger Laplace pressures present in the system. We study coarsening of foams and concentrated emulsions when there is insufficient of the trapped species present. Rate-limiting mechanisms are considered, their applicability and associated droplet growth rates discussed. In a concentrated foam or emulsion, a finite yield threshold for droplet rearrangement may be enough to prevent coarsening of the remainder.

Theory of Sol-Gel Transition in Thermoreversible Gels with due Regard for Fundamental Role of Mesoscopic Cyclization Effects. I. Thermodynamic and Structural Characteristics of Gel Phase

Abstract: The sol-gel transition (SGT), upon which the infinite cluster (IC) of thermoreversibly bonded particles (gel fraction) appears against a background of a set of finite clusters (sol fraction), is first quantitatively considered with due regard for large and complicated (mesoscopic) cycles inevitably present in the IC. To this end we present a new approach based on a concept of the monomer identity breaking and density functional description. We strictly derive, via a proper choice of basic structural units of the gel fraction, all statements usually supposed to be just Flory (Stockmayer) assumptions. A further analysis of the IC structure reveals some new IC structural units (those involved into mesoscopic cycles) overlooked in both Flory and Stockmayer approaches and to be described by a new order parameter characteristic only of the gel phase. As a result, the SGT is found to transform from a geometric phenomenon to a genuine 1st order phase transition always followed by a phase separation into sol and gel phases. The free energy, total conversion, volume fractions of the gel fraction and dangling monomers as well as other structural quantities are calculated as functions of a reduced monomer density and analyzed for all the existing models. The Flory approach is found to be superior to the Stockmayer-Tanaka one and satisfactorily describe some of the dense weak gel properties but fail (even qualitatively) in a quite extended vicinity of the SGT.

Coil-Globule Transition for Regular, Random and Specially Designed Copolymers: Monte Carlo Simulation and Self-Consistent Field Theory

Abstract: The coil-globule transition has been studied for A-B copolymer chains both by means of lattice Monte Carlo simulations using bond fluctuation algorithm and by a numerical self-consistent field method. Copolymer chains of fixed length with A and B monomeric units with regular, random and specially designed (protein-like) primary sequences have been investigated. The dependence of the transition temperature on the AB sequence has been analyzed. A protein-like copolymer is more stable than a copolymer with statistically random sequence. The transition is more sharpe for random copolymers. It is found that there exist a temperature below which the chain appears to be in the lowest energy state (ground state). Both for random and protein-like sequences and for regular copolymers with a relatively long repeating block, a molten globule regime is found between the ground state temperature and the transition temperature. For regular block copolymers the transition temperature increases with block size. Qualitatively, the results from both methods are in agreement. Differences between the methods result from approximations in the SCF theory and equilibration problems in MC simulations. The two methods are thus complementary.

On the Plants Leaves Boundary, "Jupe \`a Godets" and Conformal Embeddings

Abstract: The stable profile of the boundary of a plant's leaf fluctuating in the direction transversal to the leaf's surface is described in the framework of a model called a "surface \`a godets". It is shown that the information on the profile is encoded in the Jacobian of a conformal mapping (the coefficient of deformation) corresponding to an isometric embedding of a uniform Cayley tree into the 3D Euclidean space. The geometric characteristics of the leaf's boundary (like the perimeter and the height) are calculated. In addition a symbolic language allowing to investigate statistical properties of a "surface \`a godets" with annealed random defects of curvature of density $q$ is developed. It is found that at $q=1$ the surface exhibits a phase transition with critical exponent $\alpha=1/2$ from the exponentially growing to the flat structure.

Monte Carlo Simulation of Liquid-Crystal Alignment and of Chiral Symmetry-Breaking

Abstract: We carry out Monte Carlo simulations to investigate the effect of molecular shape on liquid-crystal order. In our approach, each model mesogen consists of several soft spheres bonded rigidly together. The arrangement of the spheres may be straight (to represent uniaxial molecules), Z-shaped (for biaxial molecules), or banana-shaped (for bent-core molecules). Using this approach, we investigate the alignment of the nematic phase by substrates decorated with parallel ridges. We compare results for wide and narrow ridge spacing and examine local order near the substrates, and show that our results are consistent with the predictions of Landau theory. We also investigate chiral symmetry-breaking in systems of bent-core molecules. We find a chiral crystalline phase as well as a nonchiral smectic-A phase, but not a chiral smectic-C phase.

Quantifying dynamics of the financial correlations

Abstract: A novel application of the correlation matrix formalism to study dynamics of the financial evolution is presented. This formalism allows to quantify the memory effects as well as some potential repeatable intradaily structures in the financial time-series. The present study is based on the high-frequency Deutsche Aktienindex (DAX) data over the time-period between November 1997 and December 1999 and demonstrates a power of the method. In this way two significant new aspects of the DAX evolution are identified: (i) the memory effects turn out to be sizably shorter than what the standard autocorrelation function analysis seems to indicate and (ii) there exist short term repeatable structures in fluctuations that are governed by a distinct dynamics. The former of these results may provide an argument in favour of the market efficiency while the later one may indicate origin of the difficulty in reaching a Gaussian limit, expected from the central limit theorem, in the distribution of returns on longer time-horizons.