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

Electrokinetic Effects in Catalytic Pt-Insulator Janus Swimmers

Abstract: The effect of added salt on the propulsion of Janus platinum-polystyrene colloids in hydrogen peroxide solution is studied experimentally. It is found that micromolar quantities of potassium and silver nitrate salts reduce the swimming velocity by similar amounts, while leading to significantly different effects on the overall rate of catalytic breakdown of hydrogen peroxide. It is argued that the seemingly paradoxical experimental observations could be theoretically explained by using a generalised reaction scheme that involves charged intermediates and has the topology of two nested loops.

Homogeneous ice nucleation at moderate supercooling from molecular simulation

Abstract: Among all the freezing transitions, that of water into ice is probably the most relevant to biology, physics, geology or atmospheric science. In this work we investigate homogeneous ice nucleation by means of computer simulations. We evaluate the size of the critical cluster and the nucleation rate for temperatures ranging between 15K and 35K below melting. We use the TIP4P/2005 and the TIP4P/Ice water models. Both give similar results when compared at the same temperature difference with the model's melting temperature. The size of the critical cluster varies from $\sim$8000 molecules (radius$ = 4$nm) at 15K below melting to $\sim$600 molecules (radius$ = 1.7$nm) at 35K below melting. We use Classical Nucleation Theory (CNT) to estimate the ice-water interfacial free energy and the nucleation free energy barrier. We obtain an interfacial free energy of 29(3)mN/m from an extrapolation of our results to the melting temperature. This value is in good agreement both with experimental measurements and with previous estimates from computer simulations of TIP4P-like models. Moreover, we obtain estimates of the nucleation rate from simulations of the critical cluster at the barrier top. The values we get for both models agree within statistical error with experimental measurements. At temperatures higher than 20K below melting we get nucleation rates slower than the appearance of a critical cluster in all the water of the hydrosphere in the age of the universe. Therefore, our simulations predict that water freezing above this temperature must necessarily be heterogeneous.

Cooperative Motion of Active Brownian Spheres in Three-Dimensional Dense Suspensions

Abstract: The structural and dynamical properties of suspensions of self-propelled Brownian particles of spherical shape are investigated in three spatial dimensions. Our simulations reveal a phase separation into a dilute and a dense phase, above a certain density and strength of self-propulsion. The packing fraction of the dense phase approaches random close packing at high activity, yet the system remains fluid. Although no alignment mechanism exists in this model, we find long-lived cooperative motion of the particles in the dense regime. This behavior is probably due to an interface-induced sorting process. Spatial displacement correlation functions are nearly scale-free for systems with densities close to or above the glass transition density of passive systems.

Photoactivated Colloidal Dockers for Cargo Transportation

Abstract: We introduce a self-propelled colloidal hematite docker that can be steered to a small particle cargo many times its size, dock, transport the cargo to a remote location, and then release it. The self-propulsion and docking are reversible and activated by visible light. The docker can be steered either by a weak uniform magnetic field or by nanoscale tracks in a textured substrate. The light-activated motion and docking originate from osmotic/phoretic particle transport in a concentration gradient of fuel, hydrogen peroxide, induced by the photocatalytic activity of the hematite. The docking mechanism is versatile and can be applied to various materials and shapes. The hematite dockers are simple single-component particles and are synthesized in bulk quantities. This system opens up new possibilities for designing complex micrometer-size factories as well as new biomimetic systems.

Entropically Patchy Particles: Engineering Valence through Shape Entropy

Abstract: Patchy particles are a popular paradigm for the design and synthesis of nanoparticles and colloids for self-assembly. In "traditional" patchy particles, anisotropic interactions arising from patterned coatings, functionalized molecules, DNA, and other enthalpic means create the possibility for directional binding of particles into higher-ordered structures. Although the anisotropic geometry of non-spherical particles contributes to the interaction patchiness through van der Waals, electrostatic, and other interactions, how particle shape contributes entropically to self-assembly is only now beginning to be understood. It has been recently demonstrated that, for hard shapes, entropic forces are directional. A newly proposed theoretical framework that defines and quantifies directional entropic forces demonstrates the anisotropic--that is, patchy--nature of these emergent, attractive forces. Here we introduce the notion of entropically patchy particles as the entropic counterpart to enthalpically patchy particles. Using three example "families" of shapes, we judiciously modify entropic patchiness by introducing geometric features to the particles so as to target specific crystal structures, which then assembled with Monte Carlo simulations. We quantify the emergent entropic valence via a potential of mean force and torque. We generalize these shape operations to shape anisotropy dimensions, in analogy with the anisotropy dimensions introduced for enthalpically patchy particles. Our findings demonstrate that entropic patchiness and emergent valence provide a way of engineering directional bonding into nanoparticle systems, whether in the presence or absence of additional, non-entropic forces.

Energy barriers govern glassy dynamics in tissues

Abstract: Recent observations demonstrate that densely packed tissues exhibit features of glassy dynamics, such as caging behavior and dynamical heterogeneities, although it has remained unclear how single-cell properties control this behavior. Here we develop numerical and theoretical models to calculate energy barriers to cell rearrangements, which help govern cell migration in cell monolayers. In contrast to work on sheared foams, we find that energy barrier heights are exponentially distributed and depend systematically on the cell's number of neighbors. Based on these results, we predict glassy two-time correlation functions for cell motion, with a timescale that increases rapidly as cell activity decreases. These correlation functions are used to construct simple random walks that reproduce the caging behavior observed for cell trajectories in experiments. This work provides a theoretical framework for predicting collective motion of cells in wound-healing, embryogenesis and cancer tumorigenesis.

Perspective: Coulomb fluids -- weak coupling, strong coupling, in between and beyond

Abstract: We present a personal view on the current state of statistical mechanics of Coulomb fluids with special emphasis on the interactions between macromolecular surfaces, concentrating on the weak and the strong coupling limits. Both are introduced for a (primitive) counterion-only system in the presence of macroscopic, uniformly charged boundaries, where they can be derived systematically. Later we show how this formalism can be generalized to the cases with additional characteristic length scales that introduce new coupling parameters into the problem. These cases most notably include asymmetric ionic mixtures with mono- and multivalent ions that couple differently to charged surfaces, ions with internal charge (multipolar) structure and finite static polarizability, where weak and strong coupling limits can be constructed by analogy with the counterion-only case and lead to important new insights into their properties that can not be derived by any other means.

Anomalous transport in the crowded world of biological cells

Abstract: A ubiquitous observation in cell biology is that diffusion of macromolecules and organelles is anomalous, and a description simply based on the conventional diffusion equation with diffusion constants measured in dilute solution fails. This is commonly attributed to macromolecular crowding in the interior of cells and in cellular membranes, summarising their densely packed and heterogeneous structures. The most familiar phenomenon is a power-law increase of the MSD, but there are other manifestations like strongly reduced and time-dependent diffusion coefficients, persistent correlations, non-gaussian distributions of the displacements, heterogeneous diffusion, and immobile particles. After a general introduction to the statistical description of slow, anomalous transport, we summarise some widely used theoretical models: gaussian models like FBM and Langevin equations for visco-elastic media, the CTRW model, and the Lorentz model describing obstructed transport in a heterogeneous environment. Emphasis is put on the spatio-temporal properties of the transport in terms of 2-point correlation functions, dynamic scaling behaviour, and how the models are distinguished by their propagators even for identical MSDs. Then, we review the theory underlying common experimental techniques in the presence of anomalous transport: single-particle tracking, FCS, and FRAP. We report on the large body of recent experimental evidence for anomalous transport in crowded biological media: in cyto- and nucleoplasm as well as in cellular membranes, complemented by in vitro experiments where model systems mimic physiological crowding conditions. Finally, computer simulations play an important role in testing the theoretical models and corroborating the experimental findings. The review is completed by a synthesis of the theoretical and experimental progress identifying open questions for future investigation.

Nanostructured epoxies based on the self-assembly of block copolymers: a new miscible block that can be tailored to different epoxy formulations

Abstract: Nanostructured thermosets may be obtained by self-assembly of amphiphilic block copolymers in a reactive solvent and fixation of the morphologies by the cross-linking reaction. Nanostructuration requires the presence of a bock that remains miscible in the polymer during polymerization. The selection of the miscible block depends on the particular system and in some cases (e.g., for epoxy-amine network based on diglycidyl ether of bisphenol A, and 4,4'- diaminodiphenylsulfone) it is very difficult to find such a block. In this manuscript it is shown that random copolymers of methyl methacrylate (MMA) and N,N-dimethylacrylamide (DMA) containing different molar fractions of DMA, can be used as a miscible block for the nanostructuration of epoxies. The miscibility of the random copolymer during formation of the epoxy network was first analyzed determining cloud-point conversions as a function of the molar fraction of DMA in the copolymer. A thermodynamic model of the phase separation was performed using the Flory-Huggins model and taking the polydispersities of both polymers into account. A single expression of the interaction parameter based on the theory of random copolymers provided a reasonable fitting of the experimental cloud-point curves. The significant increase in the miscibility produced by using small DMA molar fractions in the copolymer was explained by the high negative value of the binary interaction energy between DMA and the epoxy-amine solvent, associated to the positive value of the interaction energy between DMA and MMA units. Block copolymers with poly(n-butyl acrylate) as the immiscible block and the random copolymer P(MMA-co-DMA) as the miscible block were used for the nanostructuration of epoxy networks. The necessary molar fraction of DMA in the miscible block to stabilize a dispersion of nanosize domains depended on the fraction of the immiscible block.

Instabilities and Topological Defects in Active Nematics

Abstract: We study a continuum model of an extensile active nematic to show that mesoscale turbulence develops in two stages: (i) ordered regions undergo an intrinsic hydrodynamic instability generating walls, lines of stong bend deformations, (ii) the walls relax by forming oppositely charged pairs of defects. Both creation and annihilation of defect pairs reinstate nematic regions which undergo further instabilities, leading to a dynamic steady state. We compare this with the development of active turbulence in a contractile active nematic.

Analytical coupled vibroacoustic modeling of membrane-type acoustic metamaterials: membrane model

Abstract: Membrane-type Acoustic Metamaterials (MAMs) have demonstrated unusual capacity in controlling low-frequency sound transmission/reflection. In this paper, an analytical vibroacoustic membrane model is developed to study sound transmission behavior of the MAM under a normal incidence. The MAM is composed of a prestretched elastic membrane with attached rigid masses. To accurately capture finite-dimension rigid mass effects on the membrane deformation, the point matching approach is adopted by applying a set of distributed point forces along the interfacial boundary between masses and the membrane. The accuracy and capability of the theoretical model is verified through the comparison with the finite element method. In particular, microstructure effects such as weight, size and eccentricity of the attached mass, pretension and thickness of the membrane on the resulting transmission peak and dip frequencies of the MAM are quantitatively investigated. New peak and dip frequencies are found for the MAM with one and multiple eccentric attached masses. The developed model can be served as an efficient tool for design of such membrane-type metamaterials.

Glassy dynamics in three-dimensional embryonic tissues

Abstract: Many biological tissues are viscoelastic, behaving as elastic solids on short timescales and fluids on long timescales. This collective mechanical behavior enables and helps to guide pattern formation and tissue layering. Here we investigate the mechanical properties of three-dimensional tissue explants from zebrafish embryos by analyzing individual cell tracks and macroscopic mechanical response. We find that the cell dynamics inside the tissue exhibit features of supercooled fluids, including subdiffusive trajectories and signatures of caging behavior. We develop a minimal, three-parameter mechanical model for these dynamics, which we calibrate using only information about cell tracks. This model generates predictions about the macroscopic bulk response of the tissue (with no fit parameters) that are verified experimentally, providing a strong validation of the model. The best-fit model parameters indicate that although the tissue is fluid-like, it is close to a glass transition, suggesting that small changes to single-cell parameters could generate a significant change in the viscoelastic properties of the tissue. These results provide a robust framework for quantifying and modeling mechanically-driven pattern formation in tissues.

Elastometry of Deflated Capsules: Elastic Moduli from Shape and Wrinkle Analysis

Abstract: Elastic capsules, prepared from droplets or bubbles attached to a capillary (as in a pendant drop tensiometer), can be deflated by suction through the capillary. We study this deflation and show that a combined analysis of the shape and wrinkling characteristics enables us to determine the elastic properties in situ. Shape contours are analyzed and fitted using shape equations derived from nonlinear membrane-shell theory to give the elastic modulus, Poisson ratio and stress distribution of the membrane. We include wrinkles, which generically form upon deflation, within the shape analysis. Measuring the wavelength of wrinkles and using the calculated stress distribution gives the bending stiffness of the membrane. We compare this method with previous approaches using the Laplace-Young equation and illustrate the method on two very different capsule materials: polymerized octadecyltrichlorosilane (OTS) capsules and hydrophobin (HFBII) coated bubbles. Our results are in agreement with the available rheological data. For hydrophobin coated bubbles the method reveals an interesting nonlinear behavior consistent with the hydrophobin molecules having a rigid core surrounded by a softer shell.

Liquid meniscus friction on a wet plate: Bubbles, lamellae and foams

Abstract: Many microfluidics devices, coating processes or diphasic flows involve the motion of a liquid meniscus on a wet wall. This motion induces a specific viscous force, that exhibits a non-linear dependency in the meniscus velocity. We propose a review of the theoretical and experimental work made on this viscous force, for simple interfacial properties. The interface is indeed assumed either perfectly compressible (mobile interface) or perfectly incompressible (rigid interface). We show that, in the second case, the viscous force exerted by the wall on the meniscus is a combination of two power laws, scaling like $Ca^{1/3}$ and $Ca^{2/3}$, with $Ca$ the capillary number. We provide a prediction for the stress exerted on a foam sliding on a wet solid and compare it with experimental data, for the incompressible case.

Specificity, flexibility and valence of DNA bonds guide emulsion architecture

Abstract: The specificity and thermal reversibility of DNA interactions have enabled the self-assembly of crystal structures, self-replicating materials and colloidal molecules. Grafting DNA onto liquid interfaces of emulsions leads to exciting new architectural possibilities due to the mobility of the DNA ligands and the patches they form between bound droplets. Here we show that the size and number of these adhesion patches (valency) can be controlled. Valence 2 leads to flexible polymers of emulsion droplets, while valence above 4 leads to rigid droplet networks. A simple thermodynamic model quantitatively describes the increase in the patch size with droplet radii, DNA concentration and the stiffness of the tether to the sticky-end. The patches are formed between droplets with complementary DNA strands or alternatively with complementary colloidal nanoparticles to mediate DNA binding between droplets. This emulsion system opens the route to directed self-assembly of more complex structures through distinct DNA bonds with varying strengths and controlled valence and flexibility.

Morphogenesis of defects and tactoids during isotropic-nematic phase transition in self-assembled lyotropic chromonic liquid crystals

Abstract: We explore the structure of nuclei and topological defects in the first-order phase transition between the nematic (N) and isotropic (I) phases in lyotropic chromonic liquid crystals (LCLCs). The LCLCs are formed by self-assembled molecular aggregates of various lengths and show a broad biphasic region. The defects emerge as a result of two mechanisms. 1) Surface anisotropy mechanism that endows each N nucleus (tactoid) with topological defects thanks to preferential (tangential) orientation of the director at the closed I-N interface, and 2) Kibble mechanism with defects forming when differently oriented N tactoids merge with each other. Different scenarios of phase transition involve positive (N-in-I) and negative (I-in-N) tactoids with non-trivial topology of the director field and also multiply connected tactoids-in-tactoids configurations. The closed I-N interface limiting a tactoid shows a certain number of cusps; the lips of the interface on the opposite sides of the cusp make an angle different from pi. The N side of each cusp contains a point defect-boojum. The number of cusps shows how many times the director becomes perpendicular to the I-N interface when one circumnavigates the closed boundary of the tactoid. We derive conservation laws that connect the number of cusps c to the topological strength m of defects in the N part of the simply-connected and multiply-connected tactoids. We demonstrate how the elastic anisotropy of the N phase results in non-circular shape of the disclination cores. A generalized Wulff construction is used to derive the shape of I and N tactoids as the function of I-N interfacial tension anisotropy in the frozen director field of various topological charges m. The complex shapes and structures of tactoids and topological defects demonstrate an important role of surface anisotropy in morphogenesis of phase transitions in liquid crystals.

Orientational order in concentrated suspensions of spherical microswimmers

Abstract: We use numerical simulations to probe the dynamics of concentrated suspensions of spherical microswimmers interacting hydrodynamically. Previous work in the dilute limit predicted orientational instabilities of aligned suspensions for both pusher and puller swimmers, which we confirm computationally. Unlike previous work, we show that isotropic suspensions of spherical swimmers are also always unstable. Both types of initial conditions develop long-time polar order, of a nature which depends on the hydrodynamic signature of the swimmer but very weakly on the volume fraction up to very high volume fractions.

Swarming, swirling and stasis in sequestered bristle-bots

Abstract: The collective ability of organisms to move coherently in space and time is ubiquitous in any group of autonomous agents that can move and sense each other and the environment. Here we investigate the origin of collective motion and its loss using macroscopic self-propelled Bristle-Bots, simple automata made from a toothbrush and powered by an onboard cell phone vibrator-motor, that can sense each other through shape-dependent local interactions, and can also sense the environment non-locally via the effects of confinement and substrate topography. We show that when Bristle-Bots are confined to a limited arena with a soft boundary, increasing the density drives a transition from a disordered and uncoordinated motion to organized collective motion either as a swirling cluster or a collective dynamical stasis. This transition is regulated by a single parameter, the relative magnitude of spinning and walking in a single automaton. We explain this using quantitative experiments and simulations that emphasize the role of the agent shape, environment, and confinement via boundaries. Our study shows how the behavioral repertoire of these physically interacting automatons controlled by one parameter translates into the mechanical intelligence of swarms.

Origin of similarity of phase diagrams in amphiphilic and colloidal systems with competing interactions

Abstract: We show that amphiphilic and colloidal systems with competing interactions can be described by the same Landau-Brazovskii functional. The functional is obtained by a systematic coarse-graining procedure applied tosystems with isotropic interaction potentials. Microscopic expressions for the coefficients in the functional are derived. We propose simple criteria to distinguish the effective interparticle potentials that can lead to macro- or microsegregation. Our considerations concern also charged globular proteins in aqueous solutions and other system with effective short-range attraction long-range repulsion interactions.

Power functional theory for Brownian dynamics

Abstract: Classical density functional theory (DFT) provides an exact variational framework for determining the equilibrium properties of inhomogeneous fluids. We report a generalization of DFT to treat the non-equilibrium dynamics of classical many-body systems subject to Brownian dynamics. Our approach is based upon a dynamical functional consisting of reversible free energy changes and irreversible power dissipation. Minimization of this `free power' functional with respect to the microscopic one-body current yields a closed equation of motion. In the equililibrium limit the theory recovers the standard variational principle of DFT. The adiabatic dynamical density functional theory is obtained when approximating the power dissipation functional by that of an ideal gas. Approximations to the excess (over ideal) power dissipation yield numerically tractable equations of motion beyond the adiabatic approximation, opening the door to the systematic study of systems far from equilibrium.

Low-energy non-linear excitations in sphere packings

Abstract: We study theoretically and numerically how hard frictionless particles in random packings can rearrange. We demonstrate the existence of two distinct unstable non-linear modes of rearrangement, both associated with the opening and the closing of contacts. Mode one, whose density is characterized by some exponent {\theta}', corresponds to motions of particles extending throughout the entire system. Mode two, whose density is characterized by an exponent {\theta} != {\theta}', corresponds to the local buckling of a few particles. Mode one is shown to yield at a much higher rate than mode two when a stress is applied. We show that the distribution of contact forces follows P(f) f^{min({\theta}',{\theta})}, and that imposing that the packing cannot be densified further leads to the bounds {\gamma} >= 1/(2+{\theta}') and {\gamma} >= (1-{\theta})/2, where {\gamma} characterizes the singularity of the pair distribution function g(r) at contact. These results extend the theoretical analysis of [M. Wyart, Phys. Rev. Lett 109, 125502 (2012)] where the existence of mode two was not considered. We perform numerics that support that these bounds are saturated with {\gamma} \approx 0.38, {\theta} \approx 0.17 and {\theta}' \approx 0.44. We measure systematically the stability of all such modes in packings, and confirm their marginal stability. The principle of marginal stability thus allows to make clearcut predictions on the ensemble of configurations visited in these out-of-equilibrium systems, and on the contact forces and pair distribution functions. It also reveals the excitations that need to be included in a description of plasticity or flow near jamming, and suggests a new path to study two-level systems and soft spots in simple amorphous solids of repulsive particles.

Physicochemical Effects in Aging Aqueous Laponite Suspensions

Abstract: We study aging behavior of aqueous suspension of Laponite as a function of concentration of Laponite, concentration of salt, time elapsed since preparation of suspension (idle time) and temperature by carrying extensive rheological and conductivity experiments. We observe that temporal evolution of elastic modulus, which describes structural buildup and aging, shifts to low times for experiments carried out for higher concentration of Laponite, higher concentration of salt, greater temperature and on higher idle time while preserving the curvature of evolution in the solid regime. Consequently appropriate shifting of evolution of elastic modulus in the solid regime leads to aging time idle time salt concentration Laponite concentration temperature superposition. Existence of such superposition suggests generic nature of microstructure buildup irrespective of mentioned variables in the explored range. Behavior of shift factors needed to obtain the superposition indicate that energy barrier associated with structural buildup decreases with increase in idle time and temperature; and decreases linearly with increase in concentration of Laponite and that of salt. The conductivity experiments show that ionic conductivity of the suspension increases with increasing Laponite concentration, salt concentration, temperature and very importantly the idle time. We also analyze the interparticle interactions using DLVO theory that suggests increase in idle time, temperature and salt concentration increases height of repulsive energy barrier while decreases width of the same when particles approach each other in a parallel fashion. Analysis of rheological and conductivity experiments suggest strong influence of attractive interactions on the low energy structures in aqueous suspension of Laponite.

Single-Molecule Observation of Long Jumps in Polymer Adsorption

Abstract: Single-molecule fluorescence imaging of adsorption onto initially-bare surfaces shows that polymer chains need not localize immediately after arrival. In a system optimized to present limited adsorption sites (quartz surface to which polyethylene glycol (PEG) is exposed in aqueous solution at pH = 8.2) we find that some chains diffuse back into bulk solution and re-adsorb at some distance away, sometimes multiple times before either they localize at a stable position or else diffuse away into bulk solution. This mechanism of surface diffusion is considerably more rapid than the classical model in which adsorbed polymers crawl on surfaces while the entire molecule remains adsorbed. The trajectories with jumps follow a truncated Levy distribution of step size with limiting slope -2.5, consistent with a well-defined, rapid surface diffusion coefficient over the times we observe.

The Electrochemical Surface Potential Due to Classical Point Charge Models Drives Anion Adsorption to the Air-Water Interface

Abstract: We demonstrate that the driving forces for ion adsorption to the air-water interface for point charge models results from both cavitation and a term that is of the form of a negative electrochemical surface potential. We carefully characterize the role of the free energy due to the electrochemical surface potential computed from simple empirical models and its role in ionic adsorption within the context of dielectric continuum theory. Our research suggests that the electrochemical surface potential due to point charge models provides anions with a significant driving force to the air-water interface. This is contrary to the results of ab initio simulations that indicate that the average electrostatic surface potential should favor the desorption of anions at the air-water interface. The results have profound implications for the studies of ionic distributions in the vicinity of hydrophobic surfaces and proteins.

On the density of shear transformation zones in amorphous solids

Abstract: We study the stability of amorphous solids, focusing on the distribution P(x) of the local stress increase x that would lead to an instability. We argue that this distribution is singular P(x)x^{\theta}, where the exponent {\theta} is non-zero if the elastic interaction between rearranging regions is non-monotonic, and increases with the interaction range. For a class of finite dimensional models we show that stability implies a lower bound on {\theta}, which is found to lie near saturation. For quadrupolar interactions these models yield {\theta} ~ 0.6 for d=2 and \theta ~ 0.4 in d=3 where d is the spatial dimension, accurately capturing previously unresolved observations in atomistic models, both in quasi-static flow and after a fast quench.

On the contact area and mean gap of rough, elastic contacts: Dimensional analysis, numerical corrections and reference data

Abstract: The description of elastic, nonadhesive contacts between solids with self-affine surface roughness seems to necessitate knowledge of a large number of parameters. However, few parameters suffice to determine many important interfacial properties as we show by combining dimensional analysis with numerical simulations. This insight is used to deduce the pressure dependence of the relative contact area and the mean interfacial separation $\Delta \bar{u}$ and to present the results in a compact form. Given a proper unit choice for pressure $p$, i.e., effective modulus $E^*$ times the root-mean-square gradient $\bar{g}$, the relative contact area mainly depends on $p$ but barely on the Hurst exponent $H$ even at large $p$. When using the root-mean-square height $\bar{h}$ as unit of length, $\Delta \bar{u}$ additionally depends on the ratio of the height spectrum cutoffs at short and long wavelengths. In the fractal limit, where that ratio is zero, solely the roughness at short wavelengths is relevant for $\Delta \bar{u}$. This limit, however, should not be relevant for practical applications. Our work contains a brief summary of the employed numerical method Green's function molecular dynamics including an illustration of how to systematically overcome numerical shortcomings through appropriate finite-size, fractal, and discretization corrections. Additionally, we outline the derivation of Persson theory in dimensionless units. Persson theory compares well to the numerical reference data.

Model Vapor-Deposited Glasses: Growth Front and Composition Effects

Abstract: A growing body of experimental work indicates that physical vapor deposition provides an effective route for preparation of stable glasses, whose properties correspond in some cases to those expected for glasses that have been aged for thousands of years. In this work, model binary glasses are prepared in a process inspired by physical vapor deposition, in which particles are sequentially added to the free surface of a growing film. The resulting glasses are shown to be more stable than those prepared by gradual cooling from the liquid phase. However, it is also shown that the composition of the resulting glass, which is difficult to control in physical vapor deposition simulations of thin films, plays a significant role on the physical characteristics of the material. That composition dependence leads to a re- evaluation of previous results from simulations of thinner films than those considered here, where the equivalent age of the corresponding glasses was overestimated. The simulations presented in this work, which correspond to films that are approximately 38 molecular diameters thick, also enable analysis of the devitrification mechanism of vapor-deposited glasses. Consistent with experiments, it is found that this mechanism consists of a mobility front that propagates from the free interface into the interior of the films. Eliminating surface mobility eliminates this route of transformation into the supercooled liquid.

Interplay of inertia and deformability on rheological properties of a suspension of capsules

Abstract: The interplay of inertia and deformability has a substantial impact on the transport of soft particles suspended in a fluid. However, to date a thorough understanding of these systems is still missing and only a limited number of experimental and theoretical studies is available. We combine the finite-element, immersed-boundary and lattice-Boltzmann methods to simulate three-dimensional suspensions of soft particles subjected to planar Poiseuille flow at finite Reynolds numbers. Our findings confirm that the particle deformation and inclination increase when inertia is present. We observe that the Segre-Silberberg effect is suppressed with respect to the particle deformability. Depending on the deformability and strength of inertial effects, inward or outward lateral migration of the particles takes place. In particular, for increasing Reynolds numbers and strongly deformable particles, a hitherto unreported distinct flow focusing effect emerges which is accompanied by a non-monotonic behaviour of the apparent suspension viscosity and thickness of the particle-free layer close to the channel walls. This effect can be explained by the behaviour of a single particle and the change of the particle collision mechanism when both deformability and inertia effects are relevant.

Spontaneous knotting and unknotting of flexible linear polymers: equilibrium and kinetic aspects

Abstract: We report on a computational study of the statics and dynamics of long flexible linear polymers that spontaneously knot and unknot. Specifically, the equilibrium self-entanglement properties, such as the knotting probability, knot length and position, are investigated with extensive Monte Carlo sampling of chains of up to 15,000 beads. Tens of such equilibrated chains of up to 4, 096 beads are next used as starting points for Langevin dynamics simulations. The complex interplay of chain dynamics and self-knotting is addressed by monitoring the time evolution of various metric and entanglement properties. In particular, the extensive duration of the simulations allows for observing the spontaneous formation and disappearance of prime and composite physical knots in linear chains. Notably, a sizeable fraction of self-knotting and unknotting events is found to involve regions that are far away from the chain termini. To the best of our knowledge this represents the first instance where spontaneous changes in knotting for linear homopolymers are systematically characterized using unbiased dynamics simulations.

Finite-size scaling investigation of the liquid-liquid critical point in ST2 water and its stability with respect to crystallization

Abstract: The liquid-liquid critical point scenario of water hypothesizes the existence of two metastable liquid phases---low-density liquid (LDL) and high-density liquid (HDL)---deep within the supercooled region. The hypothesis originates from computer simulations of the ST2 water model, but the stability of the LDL phase with respect to the crystal is still being debated. We simulate supercooled ST2 water at constant pressure, constant temperature and constant number of molecules N for N<=729 and times up to 1000 ns. We observe clear differences between the two liquids, both structural and dynamical. Using several methods, including finite-size scaling, we confirm the presence of a liquid-liquid phase transition ending in a critical point. We find that the LDL is stable with respect to the crystal in 98% of our runs (we perform 372 runs for LDL or LDL-like states), and in 100% of our runs for the two largest system sizes (N=512 and 729, for which we perform 136 runs for LDL or LDL-like states). In all these runs tiny crystallites grow and then melt within 1000 ns. Only for N<=343 we observe six events (over 236 runs for LDL or LDL-like states) of spontaneous crystallization after crystallites reach an estimated critical size of about 70+/-10 molecules.