BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Friederike Schmid

Friederike Schmid appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.softBiological Physicsphysics.flu-dyncond-mat.mtrl-sciphysics.comp-phcond-mat.stat-mechBiomoleculescond-mat.mes-hallcond-mat.otherphysics.chem-phCell BehaviorSubcellular Processes

Catalog footprint

What is connected

48works
12topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

48 published item(s)

preprint2026arXiv

A general model for frictional contacts in colloidal systems

In simulations of colloidal matter, frictional contacts between particles are often neglected. For spherical colloids, such an approximation can be problematic, since frictional contacts couple translational and rotational degrees of freedom, which may affect the collective behavior of, e.g., colloids under shear and chiral active matter. Deterministic models for frictional contacts have been proposed in the granular matter community. On the colloidal scale, however, thermal fluctuations are important and should be included in a thermodynamically consistent manner. Here, we derive the correct fluctuation-dissipation relation for linear and nonlinear instantaneous frictional contact interactions. Among other, this generates a new generalized class of dissipative particle dynamics (DPD) thermostats with rotation-translation coupling. We demonstrate effects of frictional contact interactions using the examples of Poiseuille flow and motility induced phase separation in active Langevin particles.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Adsorption-active polydisperse brush with tunable molecular mass distribution

Recently a novel class of responsive uncharged polymer brushes has been proposed [Klushin et al, J. Chem. Phys. 154, 074904 (2021)] where the brush-forming chains have an affinity to the substrate. For sufficiently strong surface interactions, a fraction of chains condenses into a near-surface layer, while the remaining ones form the outer brush with a reduced grafting density. The dense layer and the more tenuous outer brush can be seen as coexisting microphases. The effective grafting density of the outer brush is controlled by the adsorption strength and can be changed reversibly as a response to changes in environmental parameters. In this paper we use numerical self-consistent field calculations and theoretical considerations to study this phenomenon in polydisperse brushes. Our results reveal an unexpected effect: Although all chains are chemically identical, shorter chains are adsorbed preferentially. Hence, with the increase in the surface affinity parameter, a reduction in the surface grafting density of the residual brush is accompanied by a change in the shape of its molecular mass distribution. In particular, an originally bidisperse brush can be effectively transformed into a nearly monodisperse one containing only the longer chain fraction.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Dynamic self-consistent field approach for studying kinetic processes in multiblock copolymer melts

The self-consistent field theory is a popular and highly successful theoretical framework for studying equilibrium (co)polymer systems at the mesoscopic level. Dynamic density functionals allow one to use this framework for studying dynamical processes in the diffusive, non-inertial regime. The central quantity in these approaches is the mobility function, which describes the effect of chain connectivity on the nonlocal response of monomers to thermodynamic driving fields. In a recent study [Mantha et al, Macromolecules 53, 3409 (2020)], we have developed a method to systematically construct mobility functions from reference fine-grained simulations. Here we focus on melts of linear chains in the Rouse regime and show how the mobility functions can be calculated semi-analytically for multiblock copolymers with arbitrary sequences without resorting to simulations. In this context, an accurate approximate expression for the single-chain dynamic structure factor is derived. Several limiting regimes are discussed. Then we apply the resulting density functional theory to study ordering processes in a two-length scale block copolymer system after instantaneous quenches into the ordered phase. Different dynamical regimes in the ordering process are identified: At early times, the ordering on short scales dominates; at late times, the ordering on larger scales takes over. For large quench depths, the system does not necessarily relax into the true equilibrium state. Our density functional approach could be used for the computer-assisted design of quenching protocols in order to create novel nonequilibrium materials.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Quorum-sensing active particles with discontinuous motility

We develop a dynamic mean-field theory for polar active particles that interact through a self-generated field, in particular one generated through emitting a chemical signal. While being a form of chemotactic response, it is different from conventional chemotaxis in that particles discontinuously change their motility when the local concentration surpasses a threshold. The resulting coupled equations for density and polarization are linear and can be solved analytically for simple geometries, yielding inhomogeneous density profiles. Specifically, here we consider a planar and circular interface. Our theory thus explains the observed coexistence of dense aggregates with an active gas. There are, however, differences to the more conventional picture of liquid-gas coexistence based on a free energy, most notably the absence of a critical point. We corroborate our analytical predictions by numerical simulations of active particles under confinement and interacting through volume exclusion. Excellent quantitative agreement is reached through an effective translational diffusion coefficient. We finally show that an additional response to the chemical gradient direction is sufficient to induce vortex clusters. Our results pave the way to engineer motility responses in order to achieve aggregation and collective behavior even at unfavorable conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A hybrid particle-continuum resolution method and its application to a homopolymer solution

We discuss in detail a recently proposed hybrid particle-continuum scheme for complex fluids and evaluate it at the example of a confined homopolymer solution in slit geometry. The hybrid scheme treats polymer chains near the impenetrable walls as particles keeping the configuration details, and chains in the bulk region as continuous density fields. Polymers can switch resolutions on the fly, controlled by an inhomogeneous tuning function. By properly choosing the tuning function, the representation of the system can be adjusted to reach an optimal balance between physical accuracy and computational efficiency. The hybrid simulation reproduces the results of a reference particle simulation and is significantly faster (about a factor of 3.5 in our application example).

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Collective behavior of quorum-sensing run-and-tumble particles in confinement

We study a generic model for quorum-sensing bacteria in circular confinement. Every bacterium produces signaling molecules, the local concentration of which triggers a response when a certain threshold is reached. If this response lowers the motility then an aggregation of bacteria occurs, which differs fundamentally from standard motility-induced phase separation due to the long-ranged nature of the concentration of signal molecules. We analyze this phenomenon analytically and by numerical simulations employing two different protocols leading to stationary cluster and ring morphologies, respectively.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Computing bulk and shear viscosities from simulations of fluids with dissipative and stochastic interactions

Exact values for bulk and shear viscosity are important to characterize a fluid and they are a necessary input for a continuum description. Here we present two novel methods to compute bulk viscosities by non-equilibrium molecular dynamics (NEMD) simulations of steady-state systems with periodic boundary conditions -- one based on frequent particle displacements and one based on the application of external bulk forces with an inhomogeneous force profile. In equilibrium simulations, viscosities can be determined from the stress tensor fluctuations via Green-Kubo relations; however, the correct incorporation of random and dissipative forces is not obvious. We discuss different expressions proposed in the literature and test them at the example of a dissipative particle dynamics (DPD) fluid.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Physical mechanisms of micro- and nanodomain formation in multicomponent lipid membranes

This article summarizes a variety of physical mechanisms proposed in the literature, which can generate micro- and nanodomains in multicomponent lipid bilayers and biomembranes. It mainly focusses on lipid-driven mechanisms that do not involve direct protein-protein interactions. Specifically, it considers (i) equilibrium mechanisms based on lipid-lipid phase separation such as critical cluster formation close to critical points, and multiple domain formation in curved geometries, (ii) equilibrium mechanisms that stabilize two-dimensional microemulsions, such as the effect of linactants and the effect of curvature-composition coupling in bilayers and monolayers, and (iii) non-equilibrium mechanisms induced by the interaction of a biomembrane with the cellular environment, such as membrane recycling and the pinning effects of the cytoplasm. Theoretical predictions are discussed together with simulations and experiments. The presentation is guided by the theory of phase transitions and critical phenomena, and the appendix summarizes the mathematical background in a concise way within the framework of the Ginzburg-Landau theory.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Protein corona composition of PEGylated nanoparticles correlates strongly with amino acid composition of protein surface

Extensive molecular dynamics simulations reveal that the interactions between proteins and poly(ethylene glycol)(PEG) can be described in terms of the surface composition of the proteins. PEG molecules accumulate around non-polar residues while avoiding polar ones. A solvent-accessible-surface-area model of protein adsorption on PEGylated nanoparticles accurately fits a large set of data on the composition of the protein corona recently obtained by label-free proteomic mass spectrometry.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

An efficient dissipative particle dynamics-based algorithm for simulating electrolyte solutions

We propose an efficient simulation algorithm based on the dissipative particle dynamics (DPD) method for studying electrohydrodynamic phenomena in electrolyte fluids. The fluid flow is mimicked with DPD particles while the evolution of the concentration of the ionic species is described using Brownian pseudo particles. The method is designed especially for systems with high salt concentrations, as explicit treatment of the salt ions becomes computationally expensive. For illustration, we apply the method to electro-osmotic flow over patterned, superhydrophobic surfaces. The results are in good agreement with recent theoretical predictions.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Computer simulations of single particles in external electric fields

Applying electric fields is an attractive way to control and manipulate single particles or molecules, e.g., in lab-on-a-chip devices. However, the response of nanosize objects in electrolyte solution to external fields is far from trivial. It is the result of a variety of dynamical processes taking place in the ion cloud surrounding charged particles and in the bulk electrolyte, and it is governed by an intricate interplay of electrostatic and hydrodynamic interactions. Already systems composed of one single particle in electrolyte solution exhibit a complex dynamical behaviour. In this review, we discuss recent coarse-grained simulations that have been performed to obtain a molecular-level understanding of the dynamic and dielectric response of single particles and single macromolecules to external electric fields. We address both the response of charged particles to constant fields (DC fields), which can be characterized by an electrophoretic mobility, and the dielectric response of both uncharged and charged particles to alternating fields (AC fields), which is described by a complex polarizability. Furthermore, we give a brief survey of simulation algorithms and highlight some recent developments.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Interplay of curvature-induced micro- and nanodomain structures in multicomponent lipid bilayers

We discuss different mechanisms for curvature-induced domain formation in multicomponent lipid membranes and present a theoretical model that allows us to study the interplay between the domains. The model represents the membrane by two coupled monolayers, which each carry an additional order parameter field describing the local lipid composition. The spontaneous curvature of each monolayer is coupled to the local composition, moreover, the lipid compositions on opposing monolayers are coupled to each other. Using this model, we calculate the phase behavior of the bilayer in mean-field approximation. The resulting phase diagrams are surprisingly complex and reveal a variety of phases and phase transitions, including a decorated microdomain phase where nanodomains are aligned along the microdomain boundaries. Our results suggest that external membrane tension can be used to control the lateral organization of nanodomains (which might be associated with lipid "rafts") in a multicomponent lipid bilayer.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Statistical properties of linear-hyperbranched graft copolymers prepared via "hypergrafting" of AB_m monomers from linear B-functional core chains: A Molecular Dynamics simulation

The reaction of $AB_m$ monomers (m=2,3) with a multifunctional B_f-type polymer chain ("hypergrafting") is studied by coarse-grained molecular dynamics simulations. The $AB_m$ monomers are hypergrafted using the slow monomer addition strategy. Fully dendronized, i.e., perfectly branched polymers are also simulated for comparison. The degree of branching DB of the molecules obtained with the "hypergrafting" process critically depends on the rate with which monomers attach to inner monomers compared to terminal monomers. This ratio is more favorable if the $AB_m$ monomers have lower reactivity, since the free monomers then have time to diffuse inside the chain. Configurational chain properties are also determined, showing that the stretching of the polymer backbone as a consequence of the "hypergrafting" procedure is much less pronounced than for perfectly dendronized chains. Furthermore, we analyze the scaling of various quantities with molecular weight M for large M (M>100). The Wiener index scales as $M^2.3$, which is intermediate between linear chains ($M^3$) and perfectly branched polymers ($M^2 \ln(M)$). The polymer size, characterized by the radius of gyration $R_g$ or the hydrodynamic radius $R_h$, is found to scale as $R_{g,h}~M^ν$ with $ν~0.38$, which lies between the exponent of diffusion limited aggregation ($ν=0.4$) and the mean-field exponent predicted by Konkolewicz and coworkers ($ν=0.33$).

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Computational studies of biomembrane systems: Theoretical considerations, simulation models, and applications

This chapter summarizes several approaches combining theory, simulation and experiment that aim for a better understanding of phenomena in lipid bilayers and membrane protein systems, covering topics such as lipid rafts, membrane mediated interactions, attraction between transmembrane proteins, and aggregation in biomembranes leading to large superstructures such as the light harvesting complex of green plants. After a general overview of theoretical considerations and continuum theory of lipid membranes we introduce different options for simulations of biomembrane systems, addressing questions such as: What can be learned from generic models? When is it expedient to go beyond them? And what are the merits and challenges for systematic coarse graining and quasi-atomistic coarse grained models that ensure a certain chemical specificity?

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Computer Simulations of Charged Colloids in Alternating Electric Fields

We briefly review recent theoretical and simulation studies of charged colloidal dispersions in alternating electric fields (AC fields). The response of single colloid to an external field can be characterized by a complex polarizability, which describes the dielectric properties of the colloid and its surrounding electrical double layer. We present computer simulation studies of single spherical colloid, using a coarse-grained mesoscale approach that accounts in full for hydrodynamic and electrostatic interactions as well as for thermal fluctuations. We investigate systematically a number of controlling parameters, such as the amplitude and frequency of the AC-fields. The results are in good agreement with recent theoretical predictions.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Sharp and fast: Sensors and switches based on polymer brushes with adsorption-active minority chains

We propose a design for polymer-based sensors and switches with sharp switching transition and fast response time. The switching mechanism involves a radical change in the conformations of adsorption-active minority chains in a brush. Such transitions can be induced by a temperature change of only about ten degrees, and the characteristic time of the conformational change is less than a second. We present an analytical theory for these switches and support it by self-consistent field calculations and Brownian dynamics simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Structure of Cholesterol in Lipid Rafts

Rafts, or functional domains, are transient nano- or mesoscopic structures in the plasma membrane and are thought to be essential for many cellular processes such as signal transduction, adhesion, trafficking and lipid/protein sorting. Observations of these membrane heterogeneities have proven challenging, as they are thought to be both small and short-lived. With a combination of coarse-grained molecular dynamics simulations and neutron diffraction using deuterium labeled cholesterol molecules we observe raft-like structures and determine the ordering of the cholesterol molecules in binary cholesterol-containing lipid membranes. From coarse-grained computer simulations, heterogenous membranes structures were observed and characterized as small, ordered domains. Neutron diffraction was used to study the lateral structure of the cholesterol molecules. We find pairs of strongly bound cholesterol molecules in the liquid-disordered phase, in accordance with the umbrella model. Bragg peaks corresponding to ordering of the cholesterol molecules in the raft-like structures were observed and indexed by two different structures: a monoclinic structure of ordered cholesterol pairs of alternating direction in equilibrium with cholesterol plaques, i.e., triclinic cholesterol bilayers.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A Dissipative-Particle-Dynamics Model for Simulating Dynamics of Charged Colloid

A mesoscopic colloid model is developed in which a spherical colloid is represented by many interacting sites on its surface. The hydrodynamic interactions with thermal fluctuations are taken accounts in full using Dissipative Particle Dynamics, and the electrostatic interactions are simulated using Particle-Particle-Particle Mesh method. This new model is applied to investigate the electrophoretic mobility of a charged colloid under an external electric field, and the influence of salt concentration and colloid charge are systematically studied. The simulation results show good agreement with predictions from the electrokinetic theory.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

AC-field-induced Polarization for Uncharged Colloids in Salt Solution: A Dissipative Particle Dynamics Simulation

We study the response of a spherical colloid under alternating electric fields (AC-fields) by mesoscopic simulation method, accounting in full for hydrodynamic and electrostatic interactions. We focus on a special case of uncharged colloids. The main polarization mechanism is the "volume-polarization", where ionic currents are deflected by the core of the uncharged colloid. Specifically, we compute the polarizability of a single colloid and systematically investigate the effect of AC field frequency and salt concentration. The simulation results are compared with predictions from classical Maxwell-Wagner theory and electrokinetic theory.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

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 (DC) and alternating fields (AC) 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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Effective slip-length tensor for a flow over weakly slipping stripes

We discuss the flow past a flat heterogeneous solid surface decorated by slipping stripes. The spatially varying slip length, $b(y)$, is assumed to be small compared to the scale of the heterogeneities, $L$, but finite. For such "weakly" slipping surfaces, earlier analyses have predicted that the effective slip length is simply given by the surface-averaged slip length, which implies that the effective slip-length tensor becomes isotropic. Here we show that a different scenario is expected if the local slip length has step-like jumps at the edges of slipping heterogeneities. In this case, the next-to-leading term in an expansion of the effective slip-length tensor in powers of ${max}\,(b(y)/L)$ becomes comparable to the leading-order term, but anisotropic, even at very small $b(y)/L$. This leads to an anisotropy of the effective slip, and to its significant reduction compared to the surface-averaged value. The asymptotic formulae are tested by numerical solutions and are in agreement with results of dissipative particle dynamics simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Effective slippage on superhydrophobic trapezoidal grooves

We study the effective slippage on superhydrophobic grooves with trapezoidal cross-sections of various geometries (including the limiting cases of triangles and rectangular stripes), by using two complementary approaches. First, dissipative particle dynamics (DPD) simulations of a flow past such surfaces have been performed to validate an expression [E.S.Asmolov and O.I.Vinogradova, J. Fluid Mech. \textbf{706}, 108 (2012)] that relates the eigenvalues of the effective slip-length tensor for one-dimensional textures. Second, we propose theoretical estimates for the effective slip length and calculate it numerically by solving the Stokes equation based on a collocation method. The comparison between the two approaches shows that they are in excellent agreement. Our results demonstrate that the effective slippage depends strongly on the area-averaged slip, the amplitude of the roughness, and on the fraction of solid in contact with the liquid. To interpret these results, we analyze flow singularities near slipping heterogeneities, and demonstrate that they inhibit the effective slip and enhance the anisotropy of the flow. Finally, we propose some guidelines to design optimal one-dimensional superhydrophobic surfaces, motivated by potential applications in microfluidics.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Elastic Properties and Line Tension of Self-Assembled Bilayer Membranes

The elastic properties of a self-assembled bilayer membrane are studied using the self-consistent field theory, applied to a model system composed of flexible amphiphilic chains dissolved in hydrophilic polymeric solvents. Examining the free energy of bilayer membranes with different geometries allows us to calculate their bending modulus, Gaussian modulus, two fourth-order membrane moduli, and the line tension. The dependence of these parameters on the microscopic characteristics of the amphiphilic chain, characterized by the volume fraction of the hydrophilic component, is systematically studied. The theoretical predictions are compared with the results from a simple monolayer model, which approximates a bilayer membrane by two monolayers. Finally the region of validity of the linear elasticity theory is analyzed by examining the higher-order contributions.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Fluctuations in lipid bilayers: Are they understood?

We review recent computer simulation studies of undulating lipid bilayers. Theoretical interpretations of such fluctuating membranes are most commonly based on generalized Helfrich-type elastic models, with additional contributions of local "protrusions" and/or density fluctuations. Such models provide an excellent basis for describing the fluctuations of tensionless bilayers in the fluid phase at a quantitative level. However, this description is found to fail for membranes in the gel phase and for membranes subject to high tensions. The fluctuations of tilted gel membranes show a signature of the modulated ripple structure, which is a nearby phase observed in the pretransition regime between the fluid and tilted gel state. This complicates a quantitative analysis on mesoscopic length scales. In the case of fluid membranes under tension, the large-wavelength fluctuation modes are found to be significantly softer than predicted by theory. In the latter context, we also address the general problem of the relation between frame tension and the fluctuation tension, which has been discussed somewhat controversially in recent years. Simulations of very simple model membranes with fixed area show that the fluctuations should be controlled by the frame tension, and not by the internal tension.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Monolayer curvature stabilizes nanoscale raft domains in mixed lipid bilayers

According to the lipid raft hypothesis, biological lipid membranes are laterally heterogeneous and filled with nanoscale ordered "raft" domains, which are believed to play an important role for the organization of proteins in membranes. However, the mechanisms stabilizing such small rafts are not clear, and even their existence is sometimes questioned. Here we report the observation of raft-like structures in a coarse-grained molecular model for multicomponent lipid bilayers. On small scales, our membranes demix into a liquid ordered (lo) and a liquid disordered (ld) phase. On large scales, phase separation is suppressed and gives way to a microemulsion-type state that contains nanometer size lo domains in a ld environment. Furthermore, we introduce a mechanism that generates rafts of finite size by a coupling between monolayer curvature and local composition. We show that mismatch between the spontaneous curvatures of monolayers in the lo and ld phase induces elastic interactions, which reduce the line tension between the lo and ld phase and can stabilize raft domains with a characteristic size of the order of a few nanometers. Our findings suggest that rafts in multicomponent bilayers might be closely related to the modulated ripple phase in one-component bilayers.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Self-Consistent Field Approach for Crosslinked Copolymer Materials

A generalized self-consistent field approach for polymer networks with fixed topology is developed. It is shown that the theory reproduces the localization of crosslinks which is characteristic for gels. The theory is then used to study the order-disorder transition in regular networks of endlinked diblock copolymers. Compared to diblock copolymer melts, the transition is shifted towards lower values of the incompatibility parameter $χ$ (the Flory- Huggins parameter). Moreover, the transition becomes strongly first order already at the mean-field level. If stress is applied, the transition is further shifted and finally vanishes in a critical point.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Using field theory to construct hybrid particle-continuum simulation schemes with adaptive resolution for soft matter systems

We develop a multiscale hybrid scheme for simulations of soft condensed matter systems, which allows one to treat the system at the particle level in selected regions of space, and at the continuum level elsewhere. It is derived systematically from an underlying particle-based model by field theoretic methods. Particles in different representation regions can switch representations on the fly, controlled by a spatially varying tuning function. As a test case, the hybrid scheme is applied to simulate colloid-polymer composites with high resolution regions close to the colloids. The hybrid simulations are significantly faster than reference simulations of a pure particle-based model, and the results are in good agreement.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Anisotropic flow in striped superhydrophobic channels

We report results of dissipative particle dynamics simulations and develop a semi-analytical theory of an anisotropic flow in a parallel-plate channel with two superhydrophobic striped walls. Our approach is valid for any local slip at the gas sectors and an arbitrary distance between the plates, ranging from a thick to a thin channel. It allows us to optimize area fractions, slip lengths, channel thickness and texture orientation to maximize a transverse flow. Our results may be useful for extracting effective slip tensors from global measurements, such as the permeability of a channel, in experiments or simulations, and may also find applications in passive microfluidic mixing.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Are stress-free membranes really 'tensionless'?

In recent years it has been argued that the tension parameter driving the fluctuations of fluid membranes, differs from the imposed lateral stress, the 'frame tension'. In particular, stress-free membranes were predicted to have a residual fluctuation tension. In the present paper, this argument is reconsidered and shown to be inherently inconsistent -- in the sense that a linearized theory, the Monge model, is used to predict a nonlinear effect. Furthermore, numerical simulations of one-dimensional stiff membranes are presented which clearly demonstrate, first, that the internal 'intrinsic' stress in membranes indeed differs from the frame tension as conjectured, but second, that the fluctuations are nevertheless driven by the frame tension. With this assumption, the predictions of the Monge model agree excellently with the simulation data for stiffness and tension values spanning several orders of magnitude.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Dielectric Response of Nanoscopic Spherical Colloids in Alternating Electric Fields: A Dissipative Particle Dynamics Simulation

We study the response of single nanosized spherical colloids in electrolyte solution to an alternating electric field (AC field) by computer simulations. We use a coarse-grained mesoscopic simulation approach that accounts in full for hydrodynamic and electrostatic interactions as well as for thermal fluctuations. The solvent is modeled as a fluid of single Dissipative Particle Dynamics (DPD) beads, and the colloidal particle is modeled as a rigid body made of DPD beads. We compute the mobility and the polarizability of a single colloid and investigate systematically the effect of amplitude and frequency of the AC-fields. Even though the thickness of the Debye layer is not "thin" compared to the radius of the colloid, and the thermal fluctuations are significant, the results are in good agreement with the theoretical prediction of the Maxwell-Wagner-O'Konski theory, especially for uncharged colloids.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Hybrid Lattice Boltzmann / Dynamic Self-Consistent Field Simulations of Microphase Separation and Vesicle Formation in Block Copolymer Systems

We present a hybrid numerical method to introduce hydrodynamics in dynamic self-consistent field (SCF) studies of inhomogeneous polymer systems. It solves a set of coupled dynamical equations: The Navier-Stokes equations for the fluid flow, and SCF-based convection-diffusion equations for the evolution of the local monomer compositions. The Navier-Stokes equaitons are simulated by the lattice Boltzmann method and the dynamic self-consistent equations are solved by a finite difference scheme. Two applications are presented: First, we study microphase separation in symmetric and asymmetric block copolymer melts with various values of shear and bulk viscosities, comparing the results to those obtaiuned with purely diffusive dynamics. Second, we investigate the effect of hydrodynamics on vesicle formation in amphiphilic block copolymer solutions. In agreement with previous studies, hydrodynamic interactions are found to have little effect on the microphase separation at early times, buty they substantially accelerate the process of structure formation at later times. Furthermore, they also contribute to selecting the pathway of vesicle formation, favoring spherical intermediates over aspherical (disklike) ones.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Membrane-mediated protein-protein interaction: A Monte Carlo study

We investigate membrane-mediated interactions between transmembrane proteins using coarse-grained models. We compare the effective potential of mean force (PMF) between two proteins, which are always aligned parallel to the z-axis of the simulation box, with those PMFs obtained for proteins with fluctuating orientations. The PMFs are dominated by an oscillatory packing-driven contribution and a smooth attractive hydrophobic mismatch contribution, which vanishes if the hydrophobic length of the protein matches the thickness of the membrane. If protein orientations are allowed to fluctuate, the oscillations are greatly reduced compared to proteins with fixed orientation. Furthermore, the hydrophobic mismatch interaction has a smaller range. Finally, we compare the two-dimensional thickness profiles around two proteins with the predictions from the elastic theory of two coupled monolayers, and find them to be in very good agreement.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Mesoscopic Simulations of Electroosmotic Flow and Electrophoresis in Nanochannels

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 $κδ_B$, where $1/κ$ is an effective electrostatic screening length at the channel boundaries.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Separation of chiral particles in micro- or nanofluidic channels

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éclet 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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A method to compute absolute free energies or enthalpies of fluids

We propose a new method to compute the free energy or enthalpy of fluids or disordered solids by computer simulation . The main idea is to construct a reference system by freezing one representative configuration, and then carry out a thermodynamic integration. We present a strategy and an algorithm which allows to sample the thermodynamic integration path even in the case of liquids, despite the fact that the particles can diffuse freely through the system. The method is described in detail and illustrated with applications to hard sphere fluids and solids with mobile defects.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Collapse transition of flexible homopolymers with adhesive contacts

The presence of ordered structures such as helices in collapsed states of polymer chains is still an open question challenging physics and biology. In this work, we present a potential model for polymer chains with monomers that are not strictly attractive, but that can make adhesive contacts with other monomers. We find that the chain develop helical order during the process of collapsing from an initially stretched conformation. It seems in this case that the adhesive contacts help the polymer chain to stay trapped in the helix state.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Coarse-Grained Simulations of Membranes under Tension

We investigate the properties of membranes under tension by Monte-Carlo simulations of a generic coarse-grained model for lipid bilayers. We give a comprising overview of the behavior of several membrane characteristics, such as the area per lipid, the monolayer overlap, the nematic order, and pressure profiles. Both the low-temperature regime, where the membranes are in a gel phase, and the high-temperature regime, where they are in the fluid phase, are considered. In the gel state, the membrane is hardly influenced by tension. In the fluid state, high tensions lead to structural changes in the membrane, which result in different compressibility regimes. The ripple state, which is found at tension zero in the transition regime between the fluid and the gel phase, disappears under tension and gives way to an interdigitated phase. We also study the membrane fluctuations in the fluid phase. In the low tension regime the data can be fitted nicely to a suitably extended elastic theory. At higher tensions the elastic fit consistently underestimates the strength of long-wavelength fluctuations. Finally, we investigate the influence of tension on the effective interaction between simple transmembrane inclusions and show that tension can be used to tune the hydrophobic mismatch interaction between membrane proteins.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Fluctuations and Elastic Properties of Lipid Membranes in the Tilted Gel State: A Coarse-Grained Monte Carlo Study

We study the stress distribution profiles and the height and thickness fluctuations of lipid membranes in the tilted gel state by Monte Carlo simulations of a generic coarse-grained model for lipid membranes, which reproduces many known properties of DPPC bilayers. The results are related to the corresponding properties of fluid membranes, and to theoretical predictions for crystalline and hexatic membranes. One striking observation is that the spontaneous curvature of the monolayers changes sign from the fluid to the gel phase. In the gel phase, the long-wavelength height fluctuations are suppressed, and the fluctuation spectrum is highly anisotropic. In the direction of tilt, it carries the signature of soft modes that are compatible with the wavelength of the ripple phase, which emerges in the transition region between the fluid and the gel state. In the direction perpendicular to the tilt, the thickness fluctuations are almost entirely suppressed, and the height fluctuations seem to be dominated by an interfacial energy, i.e., by out-of-layer fluctuations, up to length scales of tens of nanometers.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Mesoscopic simulations of polyelectrolyte electrophoresis in nanochannels

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Polyelectrolyte Electrophoresis in Nanochannels: A Dissipative Particle Dynamics Simulation

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 $κ\: δ_B$, where $δ_B$ is the slip length, and $κ^{-1}$ an effective electrostatic screening length at the channel boundaries.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Theory and Simulation of Multiphase Polymer Systems

The theory of multiphase polymer systems has a venerable tradition. The 'classical' theory of polymer demixing, the Flory-Huggins theory, was developed already in the forties of the last century. It is still the starting point for most current approaches -- be they improved theories for polymer (im)miscibility that take into account the microscopic structure of blends more accurately, or sophisticated field theories that allow to study inhomogeneous multicomponent systems of polymers with arbitrary architectures in arbitrary geometries. In contrast, simulations of multiphase polymer systems are relatively young. They are still limited by the fact that one must simulate a large number of large molecules in order to obtain meaningful results. Both powerful computers and smart modeling and simulation approaches are necessary to overcome this problem. This article gives an overview over the state-of-the art in both areas, theory and simulation. While the theory has reached a fairly mature stage by now, and many aspects of it are covered in textbooks on polymer physics, the information on simulations is much more scattered. This is why some effort has been invested into putting together a representative list of references in this area (up to the year of 2008) -- which is of course still far from complete.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Computing absolute free energies of disordered structures by molecular simulation

We present a Monte Carlo simulation technique by which the free energy of disordered systems can be computed directly. It is based on thermodynamic integration. The central idea is to construct an analytically solvable reference system from a configuration which is representative for the state of interest. The method can be applied to lattice models (e.g., the Ising model) as well as off-lattice molecular models. We focus mainly on the more challenging off-lattice case. We propose a Monte Carlo algorithm, by which the thermodynamic integration path can be sampled efficiently. At the examples of the hard sphere liquid and a hard disk solid with a defect we discuss several properties of the approach.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Mesoscopic simulations of the counterion-induced electroosmotic flow - a comparative study

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.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Bistable anchoring of nematics on rough substrates

We analyze the interplay between wetting and anchoring of nematic liquid crystals on disordering, e.g., rough substrates in the framework of the Landau-de Gennes theory, in situations of competing homeotropic and planar easy axes on the substrate and the nematic-isotropic (NI) interface. The phase diagram for azimuthally symmetric substrates is calculated. We identify two regimes - a strongly coupled regime, where the wetting transition coincides with an anchoring transition, and a weakly coupled regime, where the two are separated. The anchoring transition is first order and switches between homeotropic and planar anchoring. The two competing orientations are metastable over a broad parameter range. Hence such surfaces can be used to generate bistable surfaces.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Fluctuating Interfaces in Liquid Crystals

We review and compare recent work on the properties of fluctuating interfaces between nematic and isotropic liquid-crystalline phases. Molecular dynamics and Monte Carlo simulations have been carried out for systems of ellipsoids and hard rods with aspect ratio 15:1, and the fluctuation spectrum of interface positions (the capillary wave spectrum) has been analyzed. In addition, the capillary wave spectrum has been calculated analytically within the Landau-de Gennes theory. The theory predicts that the interfacial fluctuations can be described in terms of a wave vector dependent interfacial tension, which is anisotropic at small wavelengths (stiff director regime) and becomes isotropic at large wavelengths (flexible director regime). After determining the elastic constants in the nematic phase, theory and simulation can be compared quantitatively. We obtain good agreement for the stiff director regime. The crossover to the flexible director regime is expected at wavelengths of the order of several thousand particle diameters, which was not accessible to our simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Nematic-Isotropic Interfaces Under Shear: A Molecular Dynamics Simulation

We present a large-scale molecular dynamics study of nematic-paranematic interfaces under shear. We use a model of soft repulsive ellipsoidal particles with well-known equilibrium properties, and consider interfaces which are oriented normal to the direction of the shear gradient (common stress case). The director at the interface is oriented parallel to the interface (planar). A fixed average shear rate is imposed with Lees-Edwards boundary conditions, and the heat is dissipated with a profile-unbiased thermostat. First we study the properties of the interface at one particular shear rate in detail. The local interfacial profiles and the capillary wave fluctuations of the interfaces are calculated and compared with those of the corresponding equilibrium interface. Under shear, the interfacial width broadens and the capillary wave amplitudes at large wavelengths increase. The strain is distributed inhomogeneously in the system (shear banding), the local shear rate in the nematic region being distinctly higher than in the paranematic region. Surprisingly, we also observe (symmetry breaking) flow in the {\em vorticity} direction, with opposite direction in the nematic and the paranematic state. Finally, we investigate the stability of the interface for other shear rates and construct a nonequilibrium phase diagram.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

Elastic constants from direct correlation functions in nematic liquid crystals: a computer simulation study

Density functional theories such as the Poniewierski-Stecki theory relate the elastic properties of nematic liquid crystals with their local liquid structure, i.e., with the direct correlation function (DCF) of the particles. We propose a way to determine the DCF in the nematic state from simulations without any approximations, taking into account the dependence of pair correlations on the orientation of the director explicitly. Using this scheme, we evaluate the Frank elastic constants K1, K2 and K3 in a system of soft ellipsoids. The values are in good agreement with those obtained directly from an analysis of order fluctuations. Our method thus establishes a reliable way to calculate elastic constants from pair distributions in computer simulations.

0 citations0 reviewsNo code linkedOpen access