Linear shrinkage of hydrogel coatings exposed to flow: interplay between dissolution of water and advective transport.

We investigate the shrinkage of a surface-grafted water-swollen hydrogel under shear flows of oils by laser scanning confocal microscopy. Interestingly, external shear flows of oil lead to linear dehydration and shrinkage of the hydrogel for all investigated flow conditions irrespective of the chemical nature of the hydrogel. The reason is that the finite solubility of water in oil removes water from the hydrogel continuously by diffusion. The flow advects the water-rich oil, as demonstrated by numerical solutions of the underlying convection-diffusion equation. In line with this hypothesis, shear does not cause gel shrinkage for water-saturated oils or non-solvents. The solubility of water in the oil will tune the dehydration dynamics.

The Cassie-Wenzel transition of fluids on nanostructured substrates: Macroscopic force balance versus microscopic density-functional theory.

Classical density functional theory is applied to investigate the validity of a phenomenological force-balance description of the stability of the Cassie state of liquids on substrates with nanoscale corrugation. A bulk free-energy functional of third order in local density is combined with a square-gradient term, describing the liquid-vapor interface. The bulk free energy is parameterized to reproduce the liquid density and the compressibility of water. The square-gradient term is adjusted to model the width of the water-vapor interface. The substrate is modeled by an external potential, based upon the Lennard-Jones interactions. The three-dimensional calculation focuses on substrates patterned with nanostripes and square-shaped nanopillars. Using both the force-balance relation and density-functional theory, we locate the Cassie-to-Wenzel transition as a function of the corrugation parameters. We demonstrate that the force-balance relation gives a qualitatively reasonable description of the transition even on the nanoscale. The force balance utilizes an effective contact angle between the fluid and the vertical wall of the corrugation to parameterize the impalement pressure. This effective angle is found to have values smaller than the Young contact angle. This observation corresponds to an impalement pressure that is smaller than the value predicted by macroscopic theory. Therefore, this effective angle embodies effects specific to nanoscopically corrugated surfaces, including the finite range of the liquid-solid potential (which has both repulsive and attractive parts), line tension, and the finite interface thickness. Consistently with this picture, both patterns (stripes and pillars) yield the same effective contact angles for large periods of corrugation.

Dynamic and dielectric response of charged colloids in electrolyte solutions to external electric fields.

Computer simulations are used to investigate the response of a charged colloid and its surrounding microion cloud to an external electric field. Both static fields and alternating fields are considered. A mesoscopic simulation method is implemented to account in full for hydrodynamic and electrostatic interactions. The response of the system can be characterized by two quantities: the mobility and the polarizability. Due to the interplay of the electrostatic attraction and hydrodynamic drag, the response of the microions close to the colloid surface is different from that of the microions far away from the colloid. Both the mobility and polarizability exhibit a dependency on the frequency of the external fields, which can be attributed to the concentration polarization, the mobility of the microions, and the inertia of microions. The effects of the colloidal charge, the salt concentration, and the frequency of the external fields are investigated systematically.

Comparing equilibration schemes of high-molecular-weight polymer melts with topological indicators.

Recent theoretical studies have demonstrated that the behaviour of molecular knots is a sensitive indicator of polymer structure. Here, we use knots to verify the ability of two state-of-the-art algorithms-configuration assembly and hierarchical backmapping-to equilibrate high-molecular-weight (MW) polymer melts. Specifically, we consider melts with MWs equivalent to several tens of entanglement lengths and various chain flexibilities, generated with both strategies. We compare their unknotting probability, unknotting length, knot spectra, and knot length distributions. The excellent agreement between the two independent methods with respect to knotting properties provides an additional strong validation of their ability to equilibrate dense high-MW polymeric liquids. By demonstrating this consistency of knotting behaviour, our study opens the way for studying topological properties of polymer melts beyond time and length scales accessible to brute-force molecular dynamics simulations.

Systematic derivation of hydrodynamic equations for viscoelastic phase separation.

We present a detailed derivation of a simple hydrodynamic two-fluid model, which aims at the description of the phase separation of non-entangled polymer solutions, where viscoelastic effects play a role. It is directly based upon the coarse-graining of a well-defined molecular model, such that all degrees of freedom have a clear and unambiguous molecular interpretation. The considerations are based upon a free-energy functional, and the dynamics is split into a conservative and a dissipative part, where the latter satisfies the Onsager relations and the second law of thermodynamics. The model is therefore fully consistent with both equilibrium and non-equilibrium thermodynamics. The derivation proceeds in two steps: firstly, we derive an extended model comprising two scalar and four vector fields, such that inertial dynamics of the macromolecules and of the relative motion of the two fluids is taken into account. In the second step, we eliminate these inertial contributions and, as a replacement, introduce phenomenological dissipative terms, which can be modeled easily by taking into account the principles of non-equilibrium thermodynamics. The final simplified model comprises the momentum conservation equation, which includes both interfacial and elastic stresses, a convection-diffusion equation where interfacial and elastic contributions occur as well, and a suitably convected relaxation equation for the end-to-end vector field. In contrast to the traditional two-scale description that is used to derive rheological equations of motion, we here treat the hydrodynamic and the macromolecular degrees of freedom on the same basis. Nevertheless, the resulting model is fairly similar, though not fully identical, to models that have been discussed previously. Notably, we find a rheological constitutive equation that differs from the standard Oldroyd-B model. Within the framework of kinetic theory, this difference may be traced back to a different underlying statistical-mechanical ensemble that is used for averaging the stress. To what extent the model is able to reproduce the full phenomenology of viscoelastic phase separation is presently an open question, which shall be investigated in the future.

Analysis of a viscoelastic phase separation model.

A new model for viscoelastic phase separation is proposed, based on a systematically derived conservative two-fluid model. Dissipative effects are included by phenomenological viscoelastic terms. By construction, the model is consistent with the second law of thermodynamics. We study well-posedness of the model in two space dimensions, i.e., existence of weak solutions, a weak-strong uniqueness principle, and stability with respect to perturbations, which are proven by means of relative energy estimates. Our numerical simulations based on the new viscoelastic phase separation model are in good agreement with physical experiments. Furthermore, a good qualitative agreement with mesoscopic simulations is observed.

Semiautomatic construction of lattice Boltzmann models.

A crucial step in constructing a lattice Boltzmann model is the definition of a suitable set of lattice velocities and the correct assignment of the associated weights. For high-order models, the solution of this problem requires a nontrivial effort. This paper outlines the functioning of a publicly available Python script which has been written to assist researchers in that task. The speed of sound c_{s} is considered as a parameter, which can, within limits, be chosen at will. Under this premise, the Maxwell-Boltzmann constraint equations are a system of linear equations to determine the weights and hence amenable to numerical solution by standard linear algebra library routines. By suitable contractions, the tensor equations are mapped to a set of equivalent scalar equations, which simplifies the treatment significantly. For a user-supplied set of velocity shells, the software first checks if a solution for the weights exists and returns it if it also happens to be unique. In such a case, the software also calculates the range of c_{s} values that yield positive weights. Standard models like D3Q19 with a well-defined special c_{s} value then result as limiting cases where one of the weights vanishes. In the case of an infinite set of solutions, the user may find one particular solution by supplying a c_{s} value and then minimizing one or several weights within the framework of standard linear programming. Some examples illustrate the feasibility and usefulness of the approach. A number of models that have been discussed in the literature are nicely reproduced, while the software has also been able to find some new models of even higher order.

Implicit and explicit solvent models for the simulation of a single polymer chain in solution: Lattice Boltzmann versus Brownian dynamics.

We present a comparative study of two computer simulation methods to obtain static and dynamic properties of dilute polymer solutions. The first approach is a recently established hybrid algorithm based on dissipative coupling between molecular dynamics and lattice Boltzmann (LB), while the second is standard Brownian dynamics (BD) with fluctuating hydrodynamic interactions. Applying these methods to the same physical system (a single polymer chain in a good solvent in thermal equilibrium) allows us to draw a detailed and quantitative comparison in terms of both accuracy and efficiency. It is found that the static conformations of the LB model are distorted when the box length L is too small compared to the chain size. Furthermore, some dynamic properties of the LB model are subject to an L(-1) finite-size effect, while the BD model directly reproduces the asymptotic L-->infinity behavior. Apart from these finite-size effects, it is also found that in order to obtain the correct dynamic properties for the LB simulations, it is crucial to properly thermalize all the kinetic modes. Only in this case, the results are in excellent agreement with each other, as expected. Moreover, Brownian dynamics is found to be much more efficient than lattice Boltzmann as long as the degree of polymerization is not excessively large.

Statistical mechanics of the fluctuating lattice Boltzmann equation.

We propose a derivation of the fluctuating lattice Boltzmann equation that is consistent with both equilibrium statistical mechanics and fluctuating hydrodynamics. The formalism is based on a generalized lattice-gas model, with each velocity direction occupied by many particles. We show that the most probable state of this model corresponds to the usual equilibrium distribution of the lattice Boltzmann equation. Thermal fluctuations about this equilibrium are controlled by the mean number of particles at a lattice site. Stochastic collision rules are described by a Monte Carlo process satisfying detailed balance. This allows for a straightforward derivation of discrete Langevin equations for the fluctuating modes. It is shown that all nonconserved modes should be thermalized, as first pointed out by Adhikari et al. [Europhys. Lett. 71, 473 (2005)]; any other choice violates the condition of detailed balance. A Chapman-Enskog analysis is used to derive the equations of fluctuating hydrodynamics on large length and time scales; the level of fluctuations is shown to be thermodynamically consistent with the equation of state of an isothermal, ideal gas. We believe this formalism will be useful in developing new algorithms for thermal and multiphase flows.

High-precision estimate of the hydrodynamic radius for self-avoiding walks.

The universal asymptotic amplitude ratio between the gyration radius and the hydrodynamic radius of self-avoiding walks is estimated by high-resolution Monte Carlo simulations. By studying chains of length of up to N=2^{25}≈34×10^{6} monomers, we find that the ratio takes the value R_{G}/R_{H}=1.5803940(45), which is several orders of magnitude more accurate than the previous state of the art. This is facilitated by a sampling scheme which is quite general and which allows for the efficient estimation of averages of a large class of observables. The competing corrections to scaling for the hydrodynamic radius are clearly discernible. We also find improved estimates for other universal properties that measure the chain dimension. In particular, a method of analysis which eliminates the leading correction to scaling results in a highly accurate estimate for the Flory exponent of ν=0.58759700(40).

Dissipative particle dynamics: a useful thermostat for equilibrium and nonequilibrium molecular dynamics simulations.

We discuss dissipative particle dynamics as a thermostat to molecular dynamics, and highlight some of its virtues: (i) universal applicability irrespective of the interatomic potential; (ii) correct and unscreened reproduction of hydrodynamic correlations; (iii) stabilization of the numerical integration of the equations of motion; and (iv) the avoidance of a profile bias in boundary-driven nonequilibrium simulations of shear flow. Numerical results on a repulsive Lennard-Jones fluid illustrate our arguments.

Hydrodynamic boundary condition of water on hydrophobic surfaces.

By combining total internal reflection fluorescence cross-correlation spectroscopy with Brownian dynamics simulations, we were able to measure the hydrodynamic boundary condition of water flowing over a smooth solid surface with exceptional accuracy. We analyzed the flow of aqueous electrolytes over glass coated with a layer of poly(dimethylsiloxane) (advancing contact angle Θ = 108°) or perfluorosilane (Θ = 113°). Within an error of better than 10 nm the slip length was indistinguishable from zero on all surfaces.

Electrophoresis of colloidal dispersions in the low-salt regime.

We study the electrophoretic mobility of spherical charged colloids in a low-salt suspension as a function of the colloidal concentration. Using an effective particle charge and a reduced screening parameter, we map the data for systems with different particle charges and sizes, including numerical simulation data with full electrostatics and hydrodynamics and experimental data for latex dispersions, on a single master curve. We observe two different volume fraction-dependent regimes for the electrophoretic mobility that can be explained in terms of the static properties of the ionic double layer.

Non-Fickian interdiffusion of dynamically asymmetric species: a molecular-dynamics study.

We use molecular dynamics combined with dissipative particle dynamics to construct a model of a binary mixture where the two species differ only in their dynamic properties (friction coefficients). For an asymmetric mixture of slow and fast particles we study the interdiffusion process. The relaxation of the composition profile is investigated in terms of its Fourier coefficients. While for weak asymmetry we observe Fickian behavior, a strongly asymmetric system exhibits clear indications of anomalous diffusion, which occurs in a crossover region between cases I (Fickian) and II (sharp front moving with constant velocity), and is close to the case II limit.