Search

By employing large-scale molecular dynamics simulations of atomistically resolved oligoelectrolytes in aqueous solutions, we study in detail the first four layer-by-layer deposition cycles of an oligoelectrolyte multilayer made of poly(diallyl dimethyl ammonium chloride)/poly(styrene sulfonate sodium salt) (PDADMAC/PSS). The multilayers are grown on a silica substrate in 0.1M {NaCl} electrolyte solutions and the swollen structures are then subsequently exposed to varying added salt concentration. We investigated the microscopic properties of the films, analyzing in detail the differences between three- and four-layers systems. Our simulations provide insights on the early stages of growth of a multilayer, which are particularly challenging for experimental observations. We found a rather strong entanglement of the oligoelectrolytes, with a fuzzy layering of the film structure. The main charge compensation mechanism is for all cases intrinsic, whereas extrinsic compensation is relatively ehanced for the layer of the last deposition cycle. In addition, we quantified other fundamental observables of these systems, as the film thickness, water uptake, and overcharge fractions for each deposition layer.

The reaction ensemble and the constant pH method are well-known chemical equilibrium approaches to simulate protonation and deprotonation reactions in classical molecular dynamics and Monte Carlo simulations. In this article, we show similarity between both methods {under certain conditions}. We perform molecular dynamics simulations of a weak polyelectrolyte in order to compare the titration curves obtained by both approaches. Our findings reveal a good agreement between the methods when the reaction ensemble is used to sweep the reaction constant. Pronounced differences between the reaction ensemble and the constant pH method can be observed for stronger acids and bases in terms of adaptive pH values. These deviations are due to the presence of explicit protons in the reaction ensemble method which induce a screening of electrostatic interactions between the charged titrable groups of the polyelectrolyte. The outcomes of our simulation hint to a better applicability of the reaction ensemble method for systems in confined geometries and titrable groups in polyelectrolytes with different pK$_\text{a}$ values., Comment: 3 figures

Core-level spectra of liquids can be difficult to interpret due to the presence of a range of local environments. We present computational methods for investigating core-level spectra based on the idea that both local structural parameters and the X-ray spectra behave as functions of the local atomic configuration around the absorbing site. We identify correlations between structural parameters and spectral intensities in defined regions of interest, using the oxygen K-edge excitation spectrum of liquid water as a test case. Our results show that this kind of analysis can find the main structure-spectral relationships of ice, liquid water, and supercritical water., Comment: 10 Pages, 8 Figures

Dissolved ions can alter the local permittivity of water, nevertheless most theories and simulations ignore this fact. We present a novel algorithm for treating spatial and temporal variations in the permittivity and use it to measure the equivalent conductivity of a salt-free polyelectrolyte solution. Our new approach quantitatively reproduces experimental results unlike simulations with a constant permittivity that even qualitatively fail to describe the data. We can relate this success to a change in the ion distribution close to the polymer due to the built-up of a permittivity gradient.

There is a large body of literature investigating the static and dynamic properties of polyelectrolytes due both to their widespread application in industrial processes and their ubiquitous presence in biology. Because of their highly charged nature, polyelectrolytes tend to alter the local dielectric permittivity of the solution within a few nanometers of their backbone. This effect has, however, been almost entirely ignored in both simulations and theoretical work. In this article we apply our recently developed electrostatic solver based on Maxwell's equations to examine the effects of the permittivity reduction in the vicinity of the polyelectrolyte. We first verify our new approach by calculating and comparing ion distributions around a linear fixed polyelectrolyte and find both quantitative and qualitative changes in the ion distribution. Further simulations with an applied electric field show that the reduction in the local dielectric constant increases the mobility of the chains by approximately ten percent. More importantly, variations in the local dielectric constant lead to qualitatively different behavior of the conductivity.

For the first time, the combination of experimental preparation and results of fully atomistic simulations of an oligoelectrolyte multilayer (OEM) made of poly(diallyl dimethyl ammonium chloride)/poly(styrene sulfonate sodium salt) (PDADMAC/PSS) is presented. The layer-by-layer growth was carried out by dipping silica substrates in oligoelectrolyte solutions and was modeled by means of atomistic molecular dynamics simulations with a protocol that mimics the experimental procedure up to the assembly of four layers. Measurements of OEM thickness, surface roughness and amount of adsorbed oligoelectrolyte chains obtained from both approaches are compared. A good agreement between simulated and experimental results was found, with some deviations due to intrinsic limitations of both methods. However, the combination of information extracted from simulations to support the analysis of experimental data can overcome such restrictions and improve the interpretation of experimental results. On the other hand, processes dominated by slower kinetics, like the destabilization of adsorbed layers upon equilibration with the surrounding environment, are out of reach for the simulation modeling approach, but they can be investigated by monitoring in situ the oligoelectrolyte adsorption during the assembly process. This demonstrates how the synergistic use of simulation and experiments improves the knowledge of OEM properties down to the molecular scale.

The influence of hydroxyectoine on the properties of the aqueous solution in presence of DPPC lipid bilayers is studied via semi-isotropic constant pressure (NPT) Molecular Dynamics simulations. We investigate the solvent-co-solute behavior in terms of Kirkwood-Buff integrals as well as hydrogen bond life times for an increasing hydroxyectoine concentration up to 0.148 mol/L. The observed preferential exclusion mechanism identifies hydroxyectoine as a kosmotropic osmolyte. Our findings in regards to the DPPC lipid bilayer indicate an increase of the surface pressure as well as the solvent accessible surface area in presence of higher hydroxyectoine concentrations. The results are in agreement to the outcome of recent experiments. With this study, we are able to validate the visibility of co-solute-solute-solvent effects for low and physiologically relevant osmolyte concentrations., Comment: 10 pages, 8 figures

We analyze the results of the German Team Handball Bundesliga for ten seasons in a model-free statistical time series approach. We will show that the home advantage is nearly negligible compared to the total sum of goals. Specific interest has been spent on the time evolution of the team fitness expressed in terms of the goal difference. In contrast to soccer, our results indicate a decay of the team fitness values over a season while the long time correlation behavior over years is nearly comparable. We are able to explain the dominance of a few teams by the large value for the total number of goals in a match. A method for the prediction of match winners is presented in good accuracy with the real results. We analyze the properties of promoted teams and indicate drastic level changes between the Bundesliga and the second league. Our findings reflect in good agreement recent discussions on modern successful attack strategies., Comment: 8 pages, 11 figures

We propose a method to separate enantiomers in microfluidic or nanofluidic channels. It requires flow profiles which break chiral symmetry and have regions with high local shear. Such profiles can be generated in channels confined by walls with different hydrodynamic boundary conditions (e.g. slip lengths). Due to a nonlinear hydrodynamic effect, particles with different chirality migrate at different speed and can be separated. The mechanism is demonstrated by computer simulations. We investigate the influence of thermal fluctuations (i.e. the P\'eclet number) and show that the effect disappears in the linear response regime. The details of the microscopic flow are important and determine which volume forces are necessary to achieve separation., Comment: Accepted for publication in Physical Review Letters

We review recent dissipative particle dynamics (DPD) simulations of electrolyte flow in nanochannels. A method is presented by which the slip length $delta_B$ at the channel boundaries can be tuned systematically from negative to infinity by introducing suitably adjusted wall-fluid friction forces. Using this method, we study electroosmotic flow (EOF) in nanochannels for varying surface slip conditions and fluids of different ionic strength. Analytic expressions for the flow profiles are derived from the Stokes equation, which are in good agreement with the numerical results. Finally, we investigate the influence of EOF on the effective mobility of polyelectrolytes in nanochannels. The relevant quantity characterizing the effect of slippage is found to be the dimensionless quantity $\kappa \delta_B$, where $1/\kappa$ is an effective electrostatic screening length at the channel boundaries.

We present Molecular Dynamics simulations of Chymotrypsin inhibitor II and PEO in presence of compatible solutes. Our results indicate that the the native compact state of the studied macromolecules is stabilized in presence of hydroxyectoine. We are able to explain the corresponding mechanism by a variation of the solvent properties for higher hydroxyectoine concentrations. Our results are validated by detailed free energy calculations., Comment: 8 pages, 9 figures

We investigate the unfolding pathway of a cytosine-rich DNA hairpin structure via Molecular Dynamics simulations. Our results indicate a hidden complexity present in the unfolding process. W e show that this complex behavior arises due to non-stochastic contributions which have an important impact on the description of the dynamics in terms of collective variables., Comment: 5 pages, 5 figures

We present Molecular Dynamics simulations of a single stranded unprotonated DNA i-motif in explicit solvent. Our results indicate that the native structure in non-acidic solution at 300 K is unstable and completely vanishes on a time scale up to 10 ns. Two unfolding mechanisms with decreasing connectivity between the initially interacting nucleobases can be identified where one pathway is characterized as entropically more favorable. The entropic preference can be mainly explained by strong water ordering effects due to hydrogen bonds for several occurring structures along the pathways. Finally we are able to indicate via free energy calculations the most stable configurations belonging to distinct hairpin structures in good agreement to experimental results., Comment: 9 pages, 9 figures

We present the results of high temperature 500 K Molecular Dynamics simulations of the DNA i-motif. The essential dynamics and the main unfolding pathways are compared to a biased metadynamics simulation at 300 K. Our results indicate a remarkable agreement of the concerted motion at both temperatures. The transition can be described by a few number of eigenvectors indicating a simple unfolding process. Two main mechanisms for the unfolding pathway at 500 K can be detected which are in good agreement to the results of the biased simulation at 300 K., Comment: 7 pages, 7 figures

We have performed Molecular Dynamics simulations of ectoine, hydroxyectoine and urea in explicit solvent. Special attention has been spent on the local surrounding structure of water molecules. Our results indicate that ectoine and hydroxyectoine are able to accumulate more water molecules than urea by a pronounced ordering due to hydrogen bonds. We have validated that the charging of the molecules is of main importance resulting in a well defined hydration sphere. The influence of a varying salt concentration is also investigated. Finally we present experimental results of a DPPC monolayer phase transition that validate our numerical findings., Comment: 9 pages, 7 figures, accepted for publication in Biophysical Chemistry

We present the results of mesoscopic dissipative particle dynamics (DPD) simulations of coupled electrohydrodynamic phenomena on the micro- and nanoscale. The effects of electroosmotic flow and slippage combined with polyelectrolyte electrophoresis are investigated in detail, taking full account of hydrodynamic and electrostatic interactions. Our numerical results are in excellent agreement with analytical calculations., Comment: 15 pages, 6 figures

We present an efficient method for the calculation of free energy landscapes. Our approach involves a history dependent bias potential which is evaluated on a grid. The corresponding free energy landscape is constructed via a histogram reweighting procedure a posteriori. Due to the presence of the bias potential, it can be also used to accelerate rare events. In addition, the calculated free energy landscape is not restricted to the actual choice of collective variables and can in principle be extended to auxiliary variables of interest without further numerical effort. The applicability is shown for several examples. We present numerical results for the alanine dipeptide and the Met-Enkephalin in explicit solution to illustrate our approach. Furthermore we derive an empirical formula that allows the prediction of the computational cost for the ordinary metadynamics variant in comparison to our approach which is validated by a dimensionless representation., Comment: 14 pages, 16 figures

We present mesoscopic DPD-simulations of polyelectrolyte electrophoresis in confined nanogeometries, for varying salt concentration and surface slip conditions. Special attention is given to the influence of electroosmotic flow (EOF) on the migration of the polyelectrolyte. The effective polyelectrolyte mobility is found to depend strongly on the boundary properties, i.e., the slip length and the width of the electric double layer. Analytic expressions for the electroosmotic mobility and the total mobility are derived which are in good agreement with the numerical results. The relevant quantity characterizing the effect of slippage is found to be the dimensionless quantity $\kappa \: \delta_B$, where $\delta_B$ is the slip length, and $\kappa^{-1}$ an effective electrostatic screening length at the channel boundaries., Comment: 8 pages, 5 figures

We present mesoscopic simulations of the counterion-induced electroosmotic flow in different electrostatic coupling regimes. Two simulation methods are compared, Dissipative Particle Dynamics (DPD) and coupled Lattice-Boltzmann/Molecular Dynamics (LB/MD). A general mapping scheme to match DPD to LB/MD is developed. For the weak-coupling regime, analytic expressions for the flow profiles in the presence of partial-slip as well as no-slip boundary conditions are derived from the Poisson-Boltzmann and Stokes equations, which are in good agreement with the numerical results. The influence of electrofriction and partial slip on the flow profiles is discussed., Comment: 10 pages, 8 figures, 3 tables, additional references and minor changes in the text

We provide numerical evidence that hydrodynamic interactions are screened for charged polymers in salt solution on time scales below the Zimm time. At very short times, a crossover to hydrodynamic behavior is observed. Our conclusions are drawn from extensive coarse-grained computer simulations of polyelectrolytes in explicit solvent and explicit salt, and discussed in terms of analytical arguments based on the Debye-Hueckel approximation., Comment: 4 pages, 4 figures, submitted to Physical Review Letters

On the micro- and nanoscale, classical hydrodynamic boundary conditions such as the no-slip condition no longer apply. Instead, the flow profiles exhibit ``slip`` at the surface, which is characterized by a finite slip length (partial slip). We present a new, systematic way of implementing partial-slip boundary conditions with arbitrary slip length in coarse-grained computer simulations. The main idea is to represent the complex microscopic interface structure by a spatially varying effective viscous force. An analytical equation for the resulting slip length can be derived for planar and for curved surfaces. The comparison with computer simulations of a DPD (dissipative particle dynamics) fluid shows that this expression is valid from full-slip to no-slip., Comment: submitted to Eur. Phys. J. E (accepted)