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

Superfluidity of Bose-Einstein Condensate in An Optical Lattice: Landau-Zener Tunneling and Dynamical Instability

Abstract: Superflow of Bose-Einstein condensate in an optical lattice is represented by a Bloch wave, a plane wave with periodic modulation of the amplitude. We review the theoretical results on the interaction effects in the energy dispersion of the Bloch waves and in the linear stability of such waves. For sufficiently strong repulsion between the atoms, the lowest Bloch band develops a loop at the edge of the Brillouin zone, with the dramatic consequence of a finite probability of Landau-Zener tunneling even in the limit of a vanishing external force. Superfluidity can exist in the central region of the Brillouin zone in the presence of a repulsive interaction, beyond which Landau instability takes place where the system can lower its energy by making transition into states with smaller Bloch wavenumbers. In the outer part of the region of Landau instability, the Bloch waves are also dynamically unstable in the sense that a small initial deviation grows exponentially in time. In the inner region of Landau instability, a Bloch wave is dynamically stable in the absence of persistent external perturbations. Experimental implications of our findings will be discussed.

Thermally Induced Losses in Ultra-Cold Atoms Magnetically Trapped Near Room-Temperature Surfaces

Abstract: We have measured magnetic trap lifetimes of ultra-cold Rb87 atoms at distances of 5-1000 microns from surfaces of conducting metals with varying resistivity. Good agreement is found with a theoretical model for losses arising from near-field magnetic thermal noise, confirming the complications associated with holding trapped atoms close to conducting surfaces. A dielectric surface (silicon) was found in contrast to be so benign that we are able to evaporatively cool atoms to a Bose-Einstein condensate by using the surface to selectively adsorb higher energy atoms.

Kosmotropes and chaotropes: modelling preferential exclusion, binding and aggregate stability

Abstract: Kosmotropic cosolvents added to an aqueous solution promote the aggregation of hydrophobic solute particles, while chaotropic cosolvents act to destabilise such aggregates. We discuss the mechanism for these phenomena within an adapted version of the two-state Muller-Lee-Graziano model for water, which provides a complete description of the ternary water/cosolvent/solute system for small solute particles. This model contains the dominant effect of a kosmotropic substance, which is to enhance the formation of water structure. The consequent preferential exclusion both of cosolvent molecules from the solvation shell of hydrophobic particles and of these particles from the solution leads to a stabilisation of aggregates. By contrast, chaotropic substances disrupt the formation of water structure, are themselves preferentially excluded from the solution, and thereby contribute to solvation of hydrophobic particles. We use Monte Carlo simulations to demonstrate at the molecular level the preferential exclusion or binding of cosolvent molecules in the solvation shell of hydrophobic particles, and the consequent enhancement or suppression of aggregate formation. We illustrate the influence of structure-changing cosolvents on effective hydrophobic interactions by modelling qualitatively the kosmotropic effect of sodium chloride and the chaotropic effect of urea.

An improved Monte Carlo method for direct calculation of the density of states

Abstract: We present an efficient Monte Carlo algorithm for determining the density of states which is based on the statistics of transition probabilities between states. By measuring the infinite temperature transition probabilities--that is, the probabilities associated with move proposal only--we are able to extract excellent estimates of the density of states. When this estimator is used in conjunction with a Wang-Landau sampling scheme [F. Wang and D. P. Landau, Phys. Rev. Lett. 86, 2050 (2001)], we quickly achieve uniform sampling of macrostates (e.g., energies) and systematically refine the calculated density of states. This approach requires only potential energy evaluations, continues to improve the statistical quality of its results as the simulation time is extended, and is applicable to both lattice and continuum systems. We test the algorithm on the Lennard-Jones liquid and demonstrate good statistical convergence properties.

Single chain properties of polyelectrolytes in poor solvent

Abstract: Using molecular dynamics simulations we study the behavior of a dilute solution of strongly charged polyelectrolytes in poor solvents, where we take counterions explicitly into account. We focus on the chain conformational properties under conditions where chain-chain interactions can be neglected, but the counterion concentration remains finite. We investigate the conformations with regard to the parameters chain length, Coulomb interaction strength, and solvent quality, and explore in which regime the competition between short range hydrophobic interactions and long range Coulomb interactions leads to pearl-necklace like structures. We observe that large number and size fluctuations in the pearls and strings lead to only small direct signatures in experimental observables like the single chain form factor. Furthermore we do not observe the predicted first order collapse of the necklace into a globular structure when counterion condensation sets in. We will also show that the pearl-necklace regime is rather small for strongly charged polyelectrolytes at finite densities. Even small changes in the charge fraction of the chain can have a large impact on the conformation due to the delicate interplay between counterion distribution and chain conformation.

Regular spatial structures in arrays of Bose-Einstein condensates induced by modulational instability

Abstract: We show that the phenomenon of modulational instability in arrays of Bose-Einstein condensates confined to optical lattices gives rise to coherent spatial structures of localized excitations. These excitations represent thin disks in 1D, narrow tubes in 2D, and small hollows in 3D arrays, filled in with condensed atoms of much greater density compared to surrounding array sites. Aspects of the developed pattern depend on the initial distribution function of the condensate over the optical lattice, corresponding to particular points of the Brillouin zone. The long-time behavior of the spatial structures emerging due to modulational instability is characterized by the periodic recurrence to the initial low-density state in a finite optical lattice. We propose a simple way to retain the localized spatial structures with high atomic concentration, which may be of interest for applications. Theoretical model, based on the multiple scale expansion, describes the basic features of the phenomenon. Results of numerical simulations confirm the analytical predictions.

Scaling of Star Polymers with one to 80 Arms

Abstract: We present large statistics simulations of 3-dimensional star polymers with up to $f=80$ arms, and with up to 4000 monomers per arm for small values of $f$. They were done for the Domb-Joyce model on the simple cubic lattice. This is a model with soft core exclusion which allows multiple occupancy of sites but punishes each same-site pair of monomers with a Boltzmann factor $v<1$. We use this to allow all arms to be attached at the central site, and we use the `magic' value $v=0.6$ to minimize corrections to scaling. The simulations are made with a very efficient chain growth algorithm with resampling, PERM, modified to allow simultaneous growth of all arms. This allows us to measure not only the swelling (as observed from the center-to-end distances), but also the partition sum. The latter gives very precise estimates of the critical exponents $\gamma_f$. For completeness we made also extensive simulations of linear (unbranched) polymers which give the best estimates for the exponent $\gamma$.

The fluid-fluid interface in a model colloid-polymer mixture: Application of grand canonical Monte Carlo to asymmetric binary mixtures

Abstract: We present a Monte Carlo method to simulate asymmetric binary mixtures in the grand canonical ensemble. The method is used to study the colloid–polymer model of Asakura and Oosawa. We determine the phase diagram of the fluid–fluid unmixing transition and the interfacial tension, both at high polymer density and close to the critical point. We also present density profiles in the two–phase region. The results are compared to predictions of a recent density functional theory.

Alexander's prescription for colloidal charge renormalization

Abstract: The interactions between charged colloidal particles in an electrolyte may be described by usual Debye-H\"uckel theory provided the source of the electric field is suitably renormalized. For spherical colloids, we reconsider and simplify the treatment of the popular proposal put forward by Alexander {\it et al.} [J. Chem. Phys. {\bf 80}, 5776 (1984)], which has proven efficient in predicting renormalized quantities (charge and salt content). We give explicit formulae for the effective charge and describe the most efficient way to apply Alexander's et al. renormalization prescription in practice. Particular attention is paid to the definition of the relevant screening length, an issue that appears confused in the literature.

A new model for simulating colloidal dynamics

Abstract: We present a new hybrid lattice-Boltzmann and Langevin molecular dynamics scheme for simulating the dynamics of suspensions of spherical colloidal particles. The solvent is modeled on the level of the lattice-Boltzmann method while the molecular dynamics is done for the solute. The coupling between the two is implemented through a frictional force acting both on the solvent and on the solute, which depends on the relative velocity. A spherical colloidal particle is represented by interaction sites at its surface. We demonstrate that this scheme quantitatively reproduces the translational and rotational diffusion of a neutral spherical particle in a liquid and show preliminary results for a charged spherical particle. We argue that this method is especially advantageous in the case of charged colloids.

Effective charge versus bare charge for colloids in the infinite dilution limit

Abstract: We propose an analytical approximation for the dependence of the effective charge on the bare charge for spherical and cylindrical macro-ions as a function of the size of the colloid and salt content, for the situation of a unique colloid immersed in a sea of electrolyte (where the definition of an effective charge is non ambiguous). Our approach is based on the Poisson-Boltzmann (PB) mean-field theory. Mathematically speaking, our estimate is asymptotically exact in the limit $\\kappa a\\gg 1$, where $a$ is the radius of the colloid and $\\kappa$ the inverse screening length. In practice, a careful comparison with effective charges parameters obtained by numerically solving the full non-linear PB theory proves that it is good down to $\\kappa a\\sim 1$. This is precisely the limit appropriate to treat colloidal suspensions. A particular emphasis is put on the range of parameters suitable to describe both single and double strand DNA molecules under physiological conditions.

The physics of chromatin

Abstract: Recent progress has been made in the understanding of the physical properties of chromatin -- the dense complex of DNA and histone proteins that occupies the nuclei of plant and animal cells. Here I will focus on the two lowest levels of the hierarchy of DNA folding into the chromatin complex: (i) the nucleosome, the chromatin repeating unit consisting of a globular aggregate of eight histone proteins with the DNA wrapped around: its overcharging, the DNA unwrapping transition, the ''sliding'' of the octamer along the DNA. (ii) The 30nm chromatin fiber, the necklace-like structure of nucleosomes connected via linker DNA: its geometry, its mechanical properties under stretching and its response to changing ionic conditions. I will stress that chromatin combines two seemingly contradictory features: (1) high compaction of DNA within the nuclear envelope and at the same time (2) accessibility to genes, promoter regions and gene regulatory sequences.

Structurally specific thermal fluctuations identify functional sites for DNA transcription

Abstract: We report results showing that thermally-induced openings of double stranded DNA coincide with the location of functionally relevant sites for transcription. Investigating both viral and bacterial DNA gene promoter segments, we found that the most probable opening occurs at the transcription start site. Minor openings appear to be related to other regulatory sites. Our results suggest that coherent thermal fluctuations play an important role in the initiation of transcription. Essential elements of the dynamics, in addition to sequence specificity, are nonlinearity and entropy, provided by local base-pair constraints.

Monitoring of Orientation in Molecular Ensembles by Polarization Sensitive Nonlinear Microscopy

Abstract: We present high resolution two-photon excitation microscopy studies combining two-photon fluorescence (TPF) and second harmonic generation (SHG) in order to probe orientational distributions of molecular ensembles at room temperature. A detailed polarization analysis of TPF and SHG signals is used in order to unravel the parameters of the molecular orientational statistical distribution, using a technique which can be extended and generalized to a broad variety of molecular arrangements. A polymer film containing molecules active for TPF and/or SHG emission is studied as a model system. Polarized TPF is shown to provide information on specific properties pertaining to incoherent emission in molecular media, such as excitation transfer. SHG, being highly sensitive to a slight departure from centrosymmetry such as induced by an external electric field in the medium, complements TPF. The response of each signal to a variable excitation polarization allows investigation of molecular behavior in complex environments which affect their orientations and interactions.

On the crescentic shape of barchan dune

Abstract: Aeolian sand dunes originate from wind flow and sand bed interactions. According to wind properties and sand availability, they can adopt different shapes, ranging from huge motion-less star dunes to small and mobile barchan dunes. The latter are crescentic and emerge under a unidirectional wind, with a low sand supply. Here, a 3d model for barchan based on existing 2d model is proposed. After describing the intrinsic issues of 3d modeling, we show that the deflection of reptating particules due to the shape of the dune leads to a lateral sand flux deflection, which takes the mathematical form of a non-linear diffusive process. This simple and physically meaningful coupling method is used to understand the shape of barchan dunes.

Fractal and Statistical Properties of Large Compact Polymers: A Computational Study

Abstract: We propose a novel combinatorial algorithm for efficient generation of Hamiltonian walks and cycles on a cubic lattice, modeling the conformations of lattice toy proteins. Through extensive tests on small lattices (allowing complete enumeration of Hamiltonian paths), we establish that the new algorithm, although not perfect, is a significant improvement over the earlier approach by Ramakrishnan et. al., as it generates the sample of conformations with dramatically reduced statistical bias. Using this method, we examine the fractal properties of typical compact conformations. In accordance with Flory theorem celebrated in polymer physics, chain pieces are found to follow Gaussian statistics on the scale smaller than the globule size. Cross-over to this Gaussian regime is found to happen at the scales which are numerically somewhat larger than previously believed. We further used Alexander and Vassiliev degrees 2 and 3 topological invariants to identify the trivial knots among the Hamiltonian loops. We found that the probability of being knotted increases with loop length much faster than it was previously thought, and that chain pieces are consistently more compact than Gaussian if the global loop topology is that of a trivial knot.

Thermodynamics of Lattice Heteropolymers

Abstract: We calculate thermodynamic quantities of HP lattice proteins by means of a multicanonical chain growth algorithm that connects the new variants of the Pruned-Enriched Rosenbluth Method (nPERM) and flat histogram sampling of the entire energy space. Since our method directly simulates the density of states, we obtain results for thermodynamic quantities of the system for all temperatures. In particular, this algorithm enables us to accurately simulate the usually difficult accessible low-temperature region. Therefore, it becomes possible to perform detailed analyses of the low-temperature transition between ground states and compact globules.

Model of Controlled Synthesis of Uniform Colloid Particles: Cadmium Sulfide

Abstract: The recently developed two-stage growth model of synthesis of monodispersed polycrystalline colloidal particles is utilized and improved to explain growth of uniform cadmium sulfide spheres. The model accounts for the coupled processes of nucleation, which yields nanocrystalline precursors, and aggregation of these subunits to form the final particles. The key parameters have been identified that control the size selection and uniformity of the CdS spheres, as well as the dynamics of the process. This approach can be used to generally describe the formation of monodispersed colloids by precipitation from homogeneous solutions.

On the formation of caveolae and similar membrane invaginations

Abstract: We study a physical model for the formation of bud-like invaginations on fluid membranes under tension, and apply this model to caveolae formation. We demonstrate that budding can be driven by membrane-bound inclusions (proteins) provided that they exert asymmetric forces on the membrane that give rise to bending moments. In particular, Caveolae formation may not necessarily require forces to be applied by the cytoskeleton. Our theoretical model is able to explain several features observed experimentally in caveolae, where proteins in the caveolin family are known to play a crucial role in the formation of caveolae buds. These include (i) the formation of caveolae buds with sizes in the 100nm range (ii) that a fairly large variation of bud shape is expected (iii) that certain N and C termini deletion mutants result in vesicles that are an order of magnitude larger. Finally, we discuss the possible origin of the morphological striations that are observed on the surfaces of the caveolae.

Translational Diffusion of Polymer Chains with Excluded Volume and Hydrodynamic Interactions by Brownian Dynamics Simulation

Abstract: Within Kirkwood theory, we study the translational diffusion coefficient of a single polymer chain in dilute solution, and focus on the small difference between the short--time Kirkwood value $D^{(K)}$ and the asymptotic long--time value $D$. We calculate this correction term by highly accurate large--scale Brownian Dynamics simulations, and show that it is in perfect agreement with the rigorous variational result $D D^{(K)}$; the older data are shown to have insufficient statistical accuracy to resolve this question. We then discuss the Green--Kubo integrand in some detail. This function behaves very differently for pre--averaged vs. fluctuating hydrodynamics, as shown for the initial value by analytical considerations corroborated by numerical results. We also present further numerical data on the chain's statics and dynamics.

Isotropic-Nematic transition of long thin hard spherocylinders confined in a quasi-two-dimensional planar geometry

Abstract: We present computer simulations of long thin hard spherocylinders in a narrow planar slit. We observe a transition from the isotropic to a nematic phase with quasi-long-range orientational order upon increasing the density. This phase transition is intrinsically two dimensional and of the Kosterlitz-Thouless type. The effective two-dimensional density at which this transition occurs increases with plate separation. We qualitatively compare some of our results with experiments where microtubules are confined in a thin slit, which gave the original inspiration for this work.

Origin of Scaling Behavior of Protein Packing Density: A Sequential Monte Carlo Study of Compact Long Chain Polymers

Abstract: Single domain proteins are thought to be tightly packed. The introduction of voids by mutations is often regarded as destabilizing. In this study we show that packing density for single domain proteins decreases with chain length. We find that the radius of gyration provides poor description of protein packing but the alpha contact number we introduce here characterize proteins well. We further demonstrate that protein-like scaling relationship between packing density and chain length is observed in off-lattice self-avoiding walks. A key problem in studying compact chain polymer is the attrition problem: It is difficult to generate independent samples of compact long self-avoiding walks. We develop an algorithm based on the framework of sequential Monte Carlo and succeed in generating populations of compact long chain off-lattice polymers up to length $N=2,000$. Results based on analysis of these chain polymers suggest that maintaining high packing density is only characteristic of short chain proteins. We found that the scaling behavior of packing density with chain length of proteins is a generic feature of random polymers satisfying loose constraint in compactness. We conclude that proteins are not optimized by evolution to eliminate packing voids.

Discrete element simulations of stress distributions in silos: crossover from two to three dimensions

Abstract: The transition from two-dimensional (2D) to three-dimensional (3D) granular packings is studied using large-scale discrete element computer simulations. We focus on vertical stress profiles and examine how they change with dimensionality from 2D to 3D. We compare results for packings in 2D, quasi-2D packings between flat plates, and 3D packings. Analysis of these packings suggests that the Janssen theory does not fully describe these packings, especially at the top of the piles, where a hydrostatic-like region of vertical stress is visible in all cases. We find that the interior of the packing is far from incipient failure, while in general, the forces at the walls are close to incipient failure.

Silo Music and Silo Quake: Granular Flow Induced Vibration

Abstract: Acceleration and sound measurements during granular discharge from silos are used to show that silo music is a sound resonance produced by silo quake. The latter is produced by stick-slip friction between the wall and the granular material in tall narrow silos. For the discharge rates studied, the occurrence and frequency of flow pulsations are determined primarily by the surface properties of the granular material and the silo wall. The measurements show that the pulsating motion of the granular material drives the oscillatory motion of the silo and the occurrence of silo quake does not require a resonant interaction between the silo and the granular material.

Small and Large Scale Granular Statics

Abstract: Recent experimental results on the static or quasistatic response of granular materials have been interpreted to suggest the inapplicability of the traditional engineering approaches, which are based on elasto-plastic models (which are elliptic in nature). Propagating (hyperbolic) or diffusive (parabolic) models have been proposed to replace the `old' models. Since several recent experiments were performed on small systems, one should not really be surprised that (continuum) elasticity, a macroscopic theory, is not directly applicable, and should be replaced by a grain-scale (``microscopic'') description. Such a description concerns the interparticle forces, while a macroscopic description is given in terms of the stress field. These descriptions are related, but not equivalent, and the distinction is important in interpreting the experimental results. There are indications that at least some large scale properties of granular assemblies can be described by elasticity, although not necessarily its isotropic version. The purely repulsive interparticle forces (in non-cohesive materials) may lead to modifications of the contact network upon the application of external forces, which may strongly affect the anisotropy of the system. This effect is expected to be small (in non-isostatic systems) for small applied forces and for pre-stressed systems (in particular for disordered systems). Otherwise, it may be accounted for using a nonlinear, incrementally elastic model, with stress-history dependent elastic moduli. Although many features of the experiments may be reproduced using models of frictionless particles, results demonstrating the importance of accounting for friction are presented.

Simple Two-State Protein Folding Kinetics Requires Near-Levinthal Thermodynamic Cooperativity

Abstract: Simple two-state folding kinetics of many small single-domain proteins are characterized by chevron plots with linear folding and unfolding arms consistent with a two-state description of equilibrium thermodynamics. This phenomenon is hereby recognized as a nontrivial heteropolymer property capable of providing fundamental insight into protein energetics. Many current protein chain models, including common lattice and continuum Gō models with explicit native biases, fail to reproduce this generic protein property. Here we show that simple two-state kinetics is obtainable from models with a cooperative interplay between core burial and local conformational propensities or an extra strongly favorable energy for the native structure. These predictions suggest that intramolecular recognition in real two-state proteins is more specific than that envisioned by common Gō-like constructs with pairwise additive energies. The many-body interactions in the present kinetically two-state models lead to high thermodynamic cooperativity as measured by their van't Hoff to calorimetric enthalpy ratios, implying that the native and denatured conformational populations are well separated in enthalpy by a high free energy barrier. Consistent with experiments, comparison among the present models with extra native favorabilities indicate that their folding rate at a given interaction strength decreases with increasing native contact order, although the dependence is weaker than the corresponding trend in real proteins. It has been observed experimentally that deviations from Arrhenius behavior are often more severe for folding than for unfolding. This asymmetry may be rationalized by one ...

Shock (Blast) Mitigation by "Soft" Condensed Matter

Abstract: It is a common point that "soft" condensed matter (like granular materials or foams) can reduce damage caused by impact or explosion. It is attributed to their ability to absorb significant energy. This is certainly the case for a quasistatic type of deformation at low velocity of impact where such materials are widely used for packing of fragile devices. At the same time a mitigation of blast phenomena must take into account shock wave properties of "soft" matter which very often exhibit highly nonlinear, highly heterogeneous and dissipative behavior. This paper considers applications of "soft" condensed matter for blast mitigation using simplified approach, presents analysis of some anomalous effects and suggestions for future research in this exciting area.

Polyelectrolyte Multilayering on a Charged Planar Surface

Abstract: The adsorption of highly \textit{oppositely} charged flexible polyelectrolytes (PEs) on a charged planar substrate is investigated by means of Monte Carlo (MC) simulations. We study in detail the equilibrium structure of the first few PE layers. The influence of the chain length and of a (extra) non-electrostatic short range attraction between the polycations and the negatively charged substrate is considered. We show that the stability as well as the microstructure of the PE layers are especially sensitive to the strength of this latter interaction. Qualitative agreement is reached with some recent experiments.

Molecular Dynamics Simulation of Solvent-Polymer Interdiffusion. I. Fickian diffusion

Abstract: The interdiffusion of a solvent into a polymer melt has been studied using large scale molecular dynamics and Monte Carlo simulation techniques. The solvent concentration profile and weight gain by the polymer have been measured as a function of time. The weight gain is found to scale as t^{1/2}, which is expected for Fickian type of diffusion. The concentration profiles are fit very well assuming Fick's second law with a constant diffusivity. The diffusivity found from fitting Fick's second law is found to be independent of time and equal to the self diffusion constant in the dilute solvent limit. We separately calculated the diffusivity as a function of concentration using the Darken equation and found that the diffusivity is essentially constant for the concentration range relevant for interdiffusion.

Slow dynamics in gelation phenomena: From chemical gels to colloidal glasses

Abstract: We here discuss the results of 3d MonteCarlo simulations of a minimal lattice model for gelling systems. We focus on the dynamics, investigated by means of the time autocorrelation function of the density fluctuations and the particle mean square displacement. We start from the case of chemical gelation, i.e. with permanent bonds, we characterize the critical dynamics as determined by the formation of the percolating cluster, as actually observed in polymer gels. By opportunely introducing a finite bond life time $\tau_b$ the dynamics displays relevant changes and eventually the onset of a glassy regime. This has been interpreted in terms of a crossover to dynamics more typical of colloidal systems and a novel connection between classical gelation and recent results on colloidal systems is suggested. By systematically comparing the results in the case of permanent bonds to finite bond lifetime, the crossover and the glassy regime can be understood in terms of effective clusters.