Researcher profile

Christian Holm

Christian Holm contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - Emerging
43works
0followers
17topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

43 published item(s)

preprint2025arXiv

Generating Quantum Reservoir State Representations with Random Matrices

We demonstrate a novel approach to reservoir computation measurements using random matrices. We do so to motivate how atomic-scale devices could be used for real-world computational applications. Our approach uses random matrices to construct reservoir measurements, introducing a simple, scalable means of generating state representations. In our studies, two reservoirs, a five-atom Heisenberg spin chain and a five-qubit quantum circuit, perform time series prediction and data interpolation. The performance of the measurement technique and current limitations are discussed in detail, along with an exploration of the diversity of measurements provided by the random matrices. In addition, we explore the role of reservoir parameters such as coupling strength and measurement dimension, providing insight into how these learning machines could be automatically tuned for different problems. This research highlights the use of random matrices to measure simple quantum reservoirs for natural learning devices, and outlines a path forward for improving their performance and experimental realization.

preprint2022arXiv

Order and Information in the Patterns of Spinning Magnetic Micro-disks at the Air-water Interface

The application of the Shannon entropy to study the relationship between information and structures has yielded insights into molecular and material systems. However, the difficulty in directly observing and manipulating atoms and molecules hampers the ability of these systems to serve as model systems for further exploring the links between information and structures. Here, we use, as a model experimental system, hundreds of spinning magnetic micro-disks self-organizing at the air-water interface to generate various spatiotemporal patterns with varying degrees of orders. Using the neighbor distance as the information-bearing variable, we demonstrate the links among information, structure, and interactions. Most importantly, we establish a direct link between information and structure without using explicit knowledge of interactions. Finally, we show that the Shannon entropy by neighbor distances is a powerful observable in characterizing structural changes. Our findings are relevant for analyzing natural self-organizing systems and for designing collective robots.

preprint2020arXiv

Modeling the Current Modulation of Bundled DNA Structures in Nanopores

We investigate the salt-dependent current modulation of bundled DNA nanostructures in a nanopore. To this end, we developed four simulation models for a 2x2 origami structure with increasing level of detail: from the mean-field level to an all-atom representation of the DNA structure. We observe a consistent pore conductivity as a function of salt concentration for all four models. However, a comparison of our data to recent experimental investigations on similar systems displays significant deviations. We discuss possible reasons for the discrepancies and propose extensions to our models for future investigations.

preprint2019arXiv

A Computational Model for Bacterial Run-and-Tumble Motion

In our article we present a computational model for the simulation of self-propelled anisotropic bacteria. To this end we use a self-propelled particle model and augment it with a statistical algorithm for the run-and-tumble motion. We derive an equation for the distribution of reorientations of the bacteria that we use to analyze the statistics of the random walk and that allows us to tune the behavior of our model to the characteristics of an E.coli bacterium. We validate our implementation in terms of a single swimmer and demonstrate that our model is capable of reproducing E.coli's run-and-tumble motion with excellent accuracy.

preprint2019arXiv

Machine Learning Inter-Atomic Potentials Generation Driven by Active Learning: A Case Study for Amorphous and Liquid Hafnium dioxide

We propose a novel active learning scheme for automatically sampling a minimum number of uncorrelated configurations for fitting the Gaussian Approximation Potential (GAP). Our active learning scheme consists of an unsupervised machine learning (ML) scheme coupled to Bayesian optimization technique that evaluates the GAP model. We apply this scheme to a Hafnium dioxide (HfO2) dataset generated from a melt-quench ab initio molecular dynamics (AIMD) protocol. Our results show that the active learning scheme, with no prior knowledge of the dataset is able to extract a configuration that reaches the required energy fit tolerance. Further, molecular dynamics (MD) simulations performed using this active learned GAP model on 6144-atom systems of amorphous and liquid state elucidate the structural properties of HfO2 with near ab initio precision and quench rates (i.e. 1.0 K/ps) not accessible via AIMD. The melt and amorphous x-ray structural factors generated from our simulation are in good agreement with experiment. Additionally, the calculated diffusion constants are in good agreement with previous ab initio studies.

preprint2019arXiv

Modeling Gel Swelling Equilibrium in Mean-Field: From explicit Models to Poisson-Boltzmann

We develop a double mean-field theory for charged macrogels immersed in electrolyte solutions in the spirit of the cell model approach. We first demonstrate that the equilibrium sampling of a single explicit coarse-grained charged polymer in a cell yields accurate predictions of the swelling equilibrium if the geometry is suitably chosen and all pressure contributions have been incorporated accurately. We then replace the explicit flexible chain by a suitably modeled penetrable charged rod that allows to compute all pressure terms within the Poisson-Boltzmann approximation. This model, albeit computationally cheap, yields excellent predictions of swelling equilibria under varying chain length, polymer charge fraction, and external reservoir salt concentrations when compared to coarse-grained molecular dynamics simulations of charged macrogels. We present an extension of the model to the experimentally relevant cases of pH-sensitive gels.

preprint2016arXiv

Clustering of Self-Propelled Triangles with Surface Roughness

Self-propelled particles can spontaneously form dense phases from a dilute suspension in a process referred to as motility-induced phase separation. The properties of the out-of-equilibrium structures that are formed are governed by the specifics of the particle interactions and the strength of the activity. Thus far, most studies into the formation of these structures have focused on spherical colloids, dumbbells, and rod-like particles endowed with various interaction potentials. Only a few studies have examined the collective behavior of more complex particle shapes. Here, we increase the geometric complexity and use Molecular Dynamics simulations to consider the structures formed by triangular self-propelled particles with surface roughness. These triangles either move towards their apex or towards their base, i.e., they possess a polarity. We find that apex-directed triangles cluster more readily, more stably, and have a smoother cluster interface than their base-directed counterparts. A difference between the two polarities is in line with the results of [H.H. Wensink, et al., Phys. Rev. E 89, 010302(R) (2014)], however, we obtain the reversed result when it comes to clustering, namely that apex-directed particles cluster more readily. We further show that reducing the surface roughness negatively impacts the stability of the base-directed structures, suggesting that their formation is in large part due to surface roughness. Our results lay a solid foundation for future experimental and computational studies into the effect of roughness on the collective dynamics of swimmers.

preprint2016arXiv

Ionic Screening and Dissociation are Crucial for Understanding Chemical Self-Propulsion in Water

Water is a polar solvent and hence supports the bulk dissociation of itself and its solutes into ions, and the re-association of these ions into neutral molecules in a dynamic equilibrium, e.g., ${\rm H_2O_2}\leftrightharpoons {\rm H^+~+~HO_2^-}$. Using continuum theory, we study the influence of these reactions on the self-propulsion of colloids driven by surface chemical reactions (chemical swimmers) in aqueous solution. The association-dissociation reactions are here shown to have a strong influence on the swimmers' behaviour, and must therefore be included in future modelling. In particular, such bulk reactions permit charged swimmers to propel electrophoretically even if all species involved in the surface reactions are neutral. The bulk reactions also significantly modify the predicted speed of chemical swimmers propelled by ionic currents, by up to an order of magnitude. For swimmers whose surface reactions produce both anions and cations (ionic self-diffusiophoresis), the bulk reactions produce an additional reactive screening length, analogous to the Debye length in electrostatics. This in turn leads to an inverse relationship between swimmer radius and swimming speed, which could provide an alternative explanation for recent experimental observations on Pt-polystyrene Janus swimmers [S. Ebbens et al., Phys. Rev. E 85, 020401 (2012)]. We also use our continuum theory to investigate the effect of the Debye screening length itself, going beyond the infinitely thin limit approximation used by previous analytical theories. We identify significant departures from this limiting behavior for micron-sized swimmers under typical experimental conditions, and find that the limiting behavior fails entirely for nanoscale swimmers.

preprint2016arXiv

Lattice-Boltzmann Hydrodynamics of Anisotropic Active Matter

A plethora of active matter models exist that describe the behavior of self-propelled particles (or swimmers), both with and without hydrodynamics. However, there are few studies that consider shape-anisotropic swimmers and include hydrodynamic interactions. Here, we introduce a simple method to simulate self-propelled colloids interacting hydrodynamically in a viscous medium using the lattice-Boltzmann technique. Our model is based on raspberry-type viscous coupling and a force/counter-force formalism which ensures that the system is force free. We consider several anisotropic shapes and characterize their hydrodynamic multipolar flow field. We demonstrate that shape-anisotropy can lead to the presence of a strong quadrupole and octupole moments, in addition to the principle dipole moment. The ability to simulate and characterize these higher-order moments will prove crucial for understanding the behavior of model swimmers in confining geometries.

preprint2016arXiv

Limiting Spurious Flow in Simulations of Electrokinetic Phenomena

Electrokinetic transport phenomena can strongly influence the behaviour of macromolecules and colloidal particles in solution, with applications in, e.g., DNA translocation through nanopores, electro-osmotic flow in nanocapillaries, and electrophoresis of charged macromolecules. Numerical simulations are an important tool to investigate these electrokinetic phenomena, but are often plagued by spurious fluxes and spurious flows that can easily exceed physical fluxes and flows. Here, we present a method that reduces one of these spurious currents, spurious flow, by several orders of magnitude. We demonstrate the effectiveness and generality of our method for both electrokinetic lattice-Boltzmann and finite-element-method based algorithms by simulating a charged sphere in an electrolyte solution, and flow through a nanopore. We also show that previous attempts to suppress these spurious currents introduce other sources of error.

preprint2016arXiv

Nanoparticle Translocation through Conical Nanopores: A Finite Element Study of Electrokinetic Transport

Recent years have seen a surge of interest in nanopores because such structures show a strong potential for characterizing nanoparticles, proteins, DNA, and even single molecules. These systems have been extensively studied in experiment as well as by all-atom and coarse-grained simulations, with a strong focus on DNA translocation. However, the equally interesting problem of particle characterization using nanopores has received far less attention. Here, we theoretically investigate the translocation of a nanoparticle through a conical nanocapillary. We use a model based on numerically solving the coupled system of electrokinetic continuum equations, which we introduce in detail. Based on our findings, we formulate basic guidelines for obtaining the maximum current signal during the translocation event, which should be transferable to other nanopore geometries. In addition, the dependence of the signal strength on particle properties, such as surface charge and size, is evaluated. Finally, we identify conditions under which the translocation is prevented by the formation of a strong electro-osmotic barrier and show that the particle may even become trapped at the pore orifice, without imposing an external hydrostatic pressure difference.

preprint2016arXiv

Selective Trapping of DNA using Glass Microcapillaries

We show experimentally that a cheap glass microcapillary can accumulate λ-phage DNA at its tip and deliver the DNA into the capillary using a combination of electro-osmotic flow, pressure-driven flow, and electrophoresis. We develop an efficient simulation model for this phenomenon based on the electrokinetic equations and the finite-element method. Using our model, we explore the large parameter space of the trapping mechanism by varying the salt concentration, the capillary surface charge, the applied voltage, the pressure difference, and the mobility of the analyte molecules. Our simulation results show that this system can be tuned to capture a wide range of analyte molecules, such as DNA or proteins, based on their electrophoretic mobility. Our method for separation and pre-concentration of analytes has implications for the development of low-cost lab-on-a-chip devices.

preprint2016arXiv

The Efficiency of Self-Phoretic Propulsion Mechanisms with Surface Reaction Heterogeneity

We consider the efficiency of self-phoretic colloidal particles (swimmers) as a function of the heterogeneity in the surface reaction rate. The set of fluid, species, and electrostatic continuity equations is solved analytically using a linearization and numerically using a finite-element method. To compare spherical swimmers of different size and with heterogeneous catalytic conversion rates, a 'swimmer efficiency' functional $η$ is introduced. It is proven, that in order to obtain maximum swimmer efficiency the reactivity has to be localized at the pole(s). Our results also shed light on the sensitivity of the propulsion speed to details of the surface reactivity, a property that is notoriously hard to measure. This insight can be utilized in the design of new self-phoretic swimmers.

preprint2016arXiv

The Raspberry Model for Hydrodynamic Interactions Revisited. II. The Effect of Confinement

The so-called 'raspberry' model refers to the hybrid lattice-Boltzmann (LB) and Langevin molecular dynamics scheme for simulating the dynamics of suspensions of colloidal particles, originally developed by [V. Lobaskin and B. Dünweg, New J. Phys. 6, 54 (2004)], wherein discrete surface points are used to achieve fluid-particle coupling. In this paper, we present a follow up to our study of the effectiveness of the raspberry model in reproducing hydrodynamic interactions in the Stokes regime for spheres arranged in a simple-cubic crystal [L. Fischer, et al., J. Chem. Phys. 143, 084108 (2015)]. Here, we consider the accuracy with which the raspberry model is able to reproduce such interactions for particles confined between two parallel plates. To this end, we compare our LB simulation results to established theoretical expressions and finite-element calculations. We show that there is a discrepancy between the translational and rotational mobility when only surface coupling points are used, as also found in Part I of our joint publication. We demonstrate that adding internal coupling points to the raspberry, can be used to correct said discrepancy in confining geometries as well. Finally, we show that the raspberry model accurately reproduces hydrodynamic interactions between a spherical colloid and planar walls up to roughly one LB lattice spacing.

preprint2016arXiv

Understanding the Onset of Oscillatory Swimming in Microchannels

Self-propelled colloids (swimmers) in confining geometries follow trajectories determined by hydrodynamic interactions with the bounding surfaces. However, typically these interactions are ignored or truncated to lowest order. We demonstrate that higher-order hydrodynamic moments cause rod-like swimmers to follow oscillatory trajectories in quiescent fluid between two parallel plates, using a combination of lattice-Boltzmann simulations and far-field calculations. This behavior occurs even far from the confining walls and does not require lubrication results. We show that a swimmer's hydrodynamic quadrupole moment is crucial to the onset of the oscillatory trajectories. This insight allows us to develop a simple model for the dynamics near the channel center based on these higher hydrodynamic moments, and suggests opportunities for trajectory-based experimental characterization of swimmers' hydrodynamic properties.

preprint2015arXiv

Buckling of paramagnetic chains in soft gels

We study the magneto-elastic coupling behavior of paramagnetic chains in soft polymer gels exposed to external magnetic fields. To this end, a laser scanning confocal microscope is used to observe the morphology of the paramagnetic chains together with the deformation field of the surrounding gel network. The paramagnetic chains in soft polymer gels show rich morphological shape changes under oblique magnetic fields, in particular a pronounced buckling deformation. The details of the resulting morphological shapes depend on the length of the chain, the strength of the external magnetic field, and the modulus of the gel. Based on the observation that the magnetic chains are strongly coupled to the surrounding polymer network, a simplified model is developed to describe their buckling behavior. A coarse-grained molecular dynamics simulation model featuring an increased matrix stiffness on the surfaces of the particles leads to morphologies in agreement with the experimentally observed buckling effects.

preprint2015arXiv

Diffusiophoretic Self-Propulsion for Partially Catalytic Spherical Colloids

Colloidal spheres with a partial platinum surface coating perform auto-phoretic motion when suspended in hydrogen peroxide solution. We present a theoretical analysis of the self-propulsion velocity of these particles using a continuum multi-component, self-diffusiophoretic model. With this model as a basis, we show how the slip-layer approximation can be derived and in which limits it holds. First, we consider the differences between the full multi-component model and the slip-layer approximation. Then the slip model is used to demonstrate and explore the sensitive nature of the particle's velocity on the details of the molecule-surface interaction. We find a strong asymmetry in the dependence of the colloid's velocity as a function of the level of catalytic coating, when there is a different interaction between the solute and solvent molecules and the inert and catalytic part of the colloid, respectively. The direction of motion can even be reversed by varying the level of the catalytic coating. Finally, we investigate the robustness of these results with respect to variations in the reaction rate near the edge between the catalytic and inert parts of the particle. Our results are of significant interest to the interpretation of experimental results on the motion of self-propelled particles.

preprint2015arXiv

Ferrogels cross-linked by magnetic particles: Field-driven deformation and elasticity studied using computer simulations

Ferrogels, i.e. swollen polymer networks into which magnetic particles are immersed, can be considered as "smart materials" since their shape and elasticity can be controlled by an external magnetic field. Using molecular dynamics simulations on the coarse-grained level we study a ferrogel in which the magnetic particles act as the cross-linkers of the polymer network. In a homogeneous external magnetic field the direct coupling between the orientation of the magnetic moments and the polymers by means of covalent bonds gives rise to a deformation of the gel, independent of the interparticle dipole-dipole interaction. In this paper we report quantitative measurements of this deformation, the gel's elastic moduli and its magnetic response. Our results demonstrate that these properties depend significantly on the topology of the polymer network.

preprint2015arXiv

Importance of Varying Permittivity on the Conductivity of Polyelectrolyte Solutions

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.

preprint2015arXiv

The influence of charged-induced variations in the local permittivity on the static and dynamic properties of polyelectrolyte solutions

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.

preprint2015arXiv

The Raspberry Model for Hydrodynamic Interactions Revisited. I. Periodic Arrays of Spheres and Dumbbells

The so-called 'raspberry' model refers to the hybrid lattice-Boltzmann and Langevin molecular dynamics scheme for simulating the dynamics of suspensions of colloidal particles, originally developed by [V. Lobaskin and B. Dünweg, New J. Phys. 6, 54 (2004)], wherein discrete surface points are used to achieve fluid-particle coupling. This technique has been used in many simulation studies on the behavior of colloids. However, there are fundamental questions with regards to the use of this model. In this paper, we examine the accuracy with which the raspberry method is able to reproduce Stokes-level hydrodynamic interactions when compared to analytic expressions for solid spheres in simple-cubic crystals. To this end, we consider the quality of numerical experiments that are traditionally used to establish these properties and we discuss their shortcomings. We show that there is a discrepancy between the translational and rotational mobility reproduced by the simple raspberry model and present a way to numerically remedy the problem by adding internal coupling points. Finally, we examine a non-convex shape, namely a colloidal dumbbell, and show that the filled raspberry model replicates the desired hydrodynamic behavior in bulk for this more complicated shape. Our investigation is continued in [J. de Graaf, et al., J. Chem. Phys. 143, 084107 (2015)], wherein we consider the raspberry model in the confining geometry of two parallel plates.

preprint2015arXiv

Towards a scale-bridging description of ferrogels and magnetic elastomers

Ferrogels and magnetic elastomers differentiate themselves from other materials by their unique capability of reversibly changing shape and mechanical properties under the influence of an external magnetic field. A crucial issue in the study of these outstanding materials is the interaction between the mesoscopic magnetic particles and the polymer matrix in which they are embedded. Here we analyze interactions between two such particles connected by a polymer chain, a situation representative for particle-crosslinked magnetic gels. To make a first step towards a scale-bridging description of the materials, effective potentials for mesoscopic configurational changes are specified using microscopic input obtained from simulations. Furthermore, the impact of the presence of magnetic interactions on the probability distributions and thermodynamic quantities of the system is considered. The resulting mesoscopic model potentials can be used to economically model the system on the particle length scales. This first coarse-graining step is important to realize simplified but realistic scale-bridging models for these promising materials.

preprint2014arXiv

Computing Coulomb Interaction in Inhomogeneous Dielectric Media via a Local Electrostatics Lattice Algorithm

The local approach to computing electrostatic interactions proposed by Maggs and adapted by Pasichnyk for molecular dynamics simulations is extended to situations where the dielectric background medium is inhomogeneous. We furthermore correct a problem of the original algorithm related to the correct treatment of the global dipole moment, provide an error estimate for the accuracy of the algorithm, and suggest a different form of the treatment of the self-energy problem. Our implementation is highly scalable on many cores, and we have validated and compared its performance against theoretical predictions and simulation data obtained by other algorithmic approaches.

preprint2014arXiv

Effects of the dipolar interaction on the equilibrium morphologies of a single supramolecular magnetic filament in bulk

We study the equilibrium morphologies of a single supramolecular magnetic filament in a three-dimensional system as a function of the effective strength of the magnetic dipolar interactions. The study is performed by means of Langevin dynamics simulations with a bead-spring chain model of freely rotating dipoles. We demonstrate the existence of three structural regimes as the value of the dipolar coupling parameter is increased: a coil compaction regime, a coil expansion regime and a closed chain regime in which the structures tend progressively to an ideal ring configuration. We discuss the governing effects of each regime, the structural transition between open and closed morphologies, and the reasons why we see no multiloop configurations that have been observed in two-dimensional systems under similar conditions.

preprint2014arXiv

Kinetic and ion pairing contributions in the dielectric spectra of electrolyte aqueous solutions

Understanding dielectric spectra can reveal important information about the dynamics of solvents and solutes from the dipolar relaxation times down to electronic ones. In the late 1970s, Hubbard and Onsager predicted that adding salt ions to a polar solution would result in a reduced dielectric permittivity that arises from the unexpected tendency of solvent dipoles to align opposite to the applied field. So far, this effect has escaped an experimental verification, mainly because of the concomitant appearance of dielectric saturation from which the Hubbard-Onsager decrement cannot be easily separated. Here we develop a novel non-equilibrium molecular dynamics simulation approach to determine this decrement accurately for the first time. Using a thermodynamic consistent all-atom force field we show that for an aqueous solution containing sodium chloride around 4.8 Mol/l, this effect accounts for 12\% of the total dielectric permittivity. The dielectric decrement can be strikingly different if a less accurate force field for the ions is used. Using the widespread GROMOS parameters, we observe in fact an {\it increment} of the dielectric permittivity rather than a decrement. We can show that this increment is caused by ion pairing, introduced by a too low dispersion force, and clarify the microscopic connection between long-living ion pairs and the appearance of specific features in the dielectric spectrum of the solution.

preprint2014arXiv

Phase diagram for a single flexible Stockmayer polymer at zero field

The equilibrium conformations of a flexible permanent magnetic chain that consists of a sequence of linked magnetic colloidal nanoparticles with short-ranged Lennard-Jones attractive interactions (Stockmayer polymer) are thoroughly analysed via Langevin dynamics simulations. A tentative phase diagram is presented for a chain of length $N=100$. The phase diagram exhibits several unusual conformational phases when compared with the non-magnetic chains. These phases are characterised by a large degree of conformational anisotropy, and consist of closed chains, helicoidal-like states, partially collapsed states, and very compact disordered states. The phase diagram contains several interesting features like the existence of at least two 'triple points'.

preprint2014arXiv

Simulation of Electric Double Layers around Charged Colloids in Aqueous Solution of Variable Permittivity

The ion distribution around charged colloids in solution has been investigated intensely during the last decade. However, few theoretical approaches have included the influence of variation in the dielectric permittivity within the system, let alone in the surrounding solvent. In this article, we introduce two relatively new methods that can solve the Poisson equation for systems with varying permittivity. The harmonic interpolation method (HIM) approximately solves the Green's function in terms of a spherical harmonics series, and thus provides analytical ion-ion potentials for the Hamiltonian of charged systems. The Maxwell Equations Molecular Dynamics (MEMD) algorithm features a local approach to electrostatics, allowing for arbitrary local changes of the dielectric constant. We show that the results of both methods are in very good agreement. We also found that the renormalized charge of the colloid, and with it the effective far field interaction, significantly changes if the dielectric properties within the vicinity of the colloid are changed.

preprint2013arXiv

Vapor-liquid coexistence of the Stockmayer fluid in nonuniform external fields

We investigate the structure and phase behavior of the Stockmayer fluid in the presence of nonuniform electric fields using molecular simulation. We find that an initially homogeneous vapor phase undergoes a local phase separation in a nonuniform field due to the combined effect of the field gradient and the fluid vapor--liquid equilibrium. This results in a high density fluid condensing in the strong field region. The system polarization exhibits a strong field dependence due to the fluid condensation.

preprint2012arXiv

Investigation of tracer diffusion in crowded cylindrical channel

Based on a coarse-grained model, we carry out molecular dynamics simulations to analyze the diffusion of a small tracer particle inside a cylindrical channel whose inner wall is covered with randomly grafted short polymeric chains. We observe an interesting transient subdiffusive behavior along the cylindrical axis at high attraction between the tracer and the chains, however, the long time diffusion is always normal. This process is found to be enhanced for the case that we immobilize the grafted chains, i.e. the sub-diffusive behavior sets in at an earlier time and spans over a longer time period before becoming diffusive. Even if the grafted chains are replaced with a frozen sea of repulsive, non-connected particles in the background, the transient subdiffusion is observed. The intermediate subdiffusive behavior only disappears when the grafted chains are replaced with a mobile background sea of mutually repulsive particles. Overall, the long time diffusion coefficient of the tracer along the cylinder axis decreases with the increase in system volume fraction, strength of attraction between the tracer and the background and also on freezing the background. We believe that the simple model presented here could be useful for a qualitative understanding of the process of macromolecular diffusion inside the nuclear pore complex.

preprint2011arXiv

Simulational study of anomalous tracer diffusion in hydrogels

In this article, we analyze different factors that affect the diffusion behavior of small tracer particles (as they are used e.g.in fluorescence correlation spectroscopy (FCS)) in the polymer network of a hydrogel and perform simulations of various simplified models. We observe, that under certain circumstances the attraction of a tracer particle to the polymer network strands might cause subdiffusive behavior on intermediate time scales. In theory, this behavior could be employed to examine the network structure and swelling behavior of weakly crosslinked hydrogels with the help of FCS.

preprint2010arXiv

Influence of pore dielectric boundaries on the translocation barrier of DNA

We investigate the impact of dielectric boundary forces on the translocation process of charged rigid DNA segments through solid neutral nanopores. We assess the electrostatic contribution to the translocation free energy barrier of a model DNA segment by evaluating the potential of mean force in absence and presence of polarization effects by means of coarse-grained molecular dynamics simulations. The effect of induced polarization charges has been taken into account by employing ICC*, a recently developed algorithm that can efficiently compute induced polarization charges induced on suitably discretized dielectric boundaries. Since water has a higher dielectric constant than the pore walls, polarization effects repel charged objects in the vicinity of the interface, with the effect of significantly increasing the free energy barrier. Another investigated side effect is the change of the counterion distribution around the charged polymer in presence of the induced pore charges. Furthermore we investigate the influence of adding salt to the solution.

preprint2010arXiv

The ICC* Algorithm: A fast way to include dielectric boundary effects into molecular dynamics simulations

We employ a fast and accurate algorithm to treat dielectric interfaces within molecular dynamics simulations and demonstrate the importance of dielectric boundary forces (DBFs) in two systems of interests in soft-condensed matter science. We investigate a salt solution confined to a slit pore, and a model of a DNA fragment translocating thorugh a narrow pore.

preprint2009arXiv

Mesoscale modelling of polyelectrolyte electrophoresis

The electrophoretic behaviour of flexible polyelectrolyte chains ranging from single monomers up to long fragments of hundred repeat units is studied by a mesoscopic simulation approach. Abstracting from the atomistic details of the polyelectrolyte and the fluid, a coarse-grained molecular dynamics model connected to a mesoscopic fluid described by the Lattice Boltzmann approach is used to investigate free-solution electrophoresis. Our study demonstrates the importance of hydrodynamic interactions for the electrophoretic motion of polyelectrolytes and quantifies the influence of surrounding ions. The length-dependence of the electrophoretic mobility can be understood by evaluating the scaling behavior of the effective charge and the effective friction. The perfect agreement of our results with experimental measurements shows that all chemical details and fluid structure can be safely neglected, and a suitable coarse-grained approach can yield an accurate description of the physics of the problem, provided that electrostatic and hydrodynamic interactions between all entities in the system, i.e., the polyelectrolyte, dissociated counterions, additional salt and the solvent, are properly accounted for. Our model is able to bridge the single molecule regime of a few nm up to macromolecules with contour lengths of more than 100 nm, a length scale that is currently not accessible to atomistic simulations.

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.

preprint2009arXiv

Optimizing end-labeled free-solution electrophoresis by increasing the hydrodynamic friction of the drag-tag

We study the electrophoretic separation of polyelectrolytes of varying lengths by means of end-labeled free-solution electrophoresis (ELFSE). A coarse-grained molecular dynamics simulation model, using full electrostatic interactions and a mesoscopic Lattice Boltzmann fluid to account for hydrodynamic interactions, is used to characterize the drag coefficients of different label types: linear and branched polymeric labels, as well as transiently bound micelles. It is specifically shown that the label's drag coefficient is determined by its hydrodynamic size, and that the drag per label monomer is largest for linear labels. However, the addition of side chains to a linear label offers the possibility to increase the hydrodynamic size, and therefore the label efficiency, without having to increase the linear length of the label, thereby simplifying synthesis. The third class of labels investigated, transiently bound micelles, seems very promising for the usage in ELFSE, as they provide a significant higher hydrodynamic drag than the other label types. The results are compared to theoretical predictions, and we investigate how the efficiency of the ELFSE method can be improved by using smartly designed drag-tags.

preprint2006arXiv

Finite Size Polyelectrolyte Bundles at Thermodynamic Equilibrium

We present the results of extensive computer simulations performed on solutions of monodisperse charged rod-like polyelectrolytes in the presence of trivalent counterions. To overcome energy barriers we used a combination of parallel tempering and hybrid Monte Carlo techniques. Our results show that for small values of the electrostatic interaction the solution mostly consists of dispersed single rods. The potential of mean force between the polyelectrolyte monomers yields an attractive interaction at short distances. For a range of larger values of the Bjerrum length, we find finite size polyelectrolyte bundles at thermodynamic equilibrium. Further increase of the Bjerrum length eventually leads to phase separation and precipitation. We discuss the origin of the observed thermodynamic stability of the finite size aggregates.

preprint2002arXiv

A molecular dynamics study on the equilibrium magnetization properties and structure of ferrofluids

We investigate in detail the initial susceptibility, magnetization curves, and microstructure of ferrofluids in various concentration and particle dipole moment ranges by means of molecular dynamics simulations. We use the Ewald summation for the long-range dipolar interactions, take explicitly into account the translational and rotational degrees of freedom, coupled to a Langevin thermostat. When the dipolar interaction energy is comparable with the thermal energy, the simulation results on the magnetization properties agree with the theoretical predictions very well. For stronger dipolar couplings, however, we find systematic deviations from the theoretical curves. We analyze in detail the observed microstructure of the fluids under different conditions. The formation of clusters is found to enhance the magnetization at weak fields and thus leads to a larger initial susceptibility. The influence of the particle aggregation is isolated by studying ferro-solids, which consist of magnetic dipoles frozen in at random locations but which are free to rotate. Due to the artificial suppression of clusters in ferro-solids the observed susceptibility is considerably lowered when compared to ferrofluids.

preprint2002arXiv

Structure of Polyelectrolytes in Poor Solvent

We present simulations on charged polymers in poor solvent. First we investigate in detail the dilute concentration range with and without imposed extension constraints. The resulting necklace polymer conformations are analyzed in detail. We find strong fluctuations in the number of pearls and their sizes leading only to small signatures in the form factor and the force-extension relation. The scaling of the peak in the structure factor with the monomer density shows a pertinent different behavior from good solvent chains.

preprint2001arXiv

End-effects of strongly charged polyelectrolytes - a molecular dynamics study

We investigate end-effects in the ion distribution around strongly charged, flexible polyelectrolytes with a quenched charge distribution by molecular dynamics simulations of dilute polyelectrolyte solutions. We take the counterions explicitly into account and calculate the full Coulomb interaction via an Ewald summation method. We find that the free counterions of the solution are distributed in such a way that a fraction of the chain charges is effectively neutralized. This in turn leads to an effective charge distribution which is similar to those found for weakly charged titrating polyelectrolytes that have an annealed charge distribution. The delicate interplay between the electrostatic interactions, the chain conformation and the counterion distribution is studied in detail as a function of different system parameters such as the chain length Nm, the charge fraction f, the charged particle density rho, the ionic strength and the solvent quality. Comparisons are made with predictions from a scaling theory.

preprint2001arXiv

Strong electrostatic interactions in spherical colloidal systems

We investigate spherical macroions in the strong Coulomb coupling regime within the primitive model in salt-free environment. We first show that the ground state of an isolated colloid is naturally overcharged by simple electrostatic arguments illustrated by the Gillespie rule. We furthermore demonstrate that in the strong Coulomb coupling this mechanism leads to ionized states and thus to long range attractions between like-charged spheres. We use molecular dynamics simulations to study in detail the counterion distribution for one and two highly charged colloids for the ground state as well as for finite temperatures. We compare our results in terms of a simple version of a Wigner crystal theory and find excellent qualitative and quantitative agreement.

preprint1996arXiv

Fixed versus random triangulations in 2D simplicial Regge calculus

We study 2D quantum gravity on spherical topologies using the Regge calculus approach with the $dl/l$ measure. Instead of a fixed non-regular triangulation which has been used before, we study for each system size four different random triangulations, which are obtained according to the standard Voronoi-Delaunay procedure. We compare both approaches quantitatively and show that the difference in the expectation value of $R^2$ between the fixed and the random triangulation depends on the lattice size and the surface area $A$. We also try again to measure the string susceptibility exponents through a finite-size scaling Ansatz in the expectation value of an added $R^2$ interaction term in an approach where $A$ is held fixed. The string susceptibility exponent $γ_{str}'$ is shown to agree with theoretical predictions for the sphere, whereas the estimate for $γ_{str}$ appears to be too negative.

preprint1995arXiv

Measuring the string susceptibility in 2D simplicial quantum gravity using the Regge approach

We use Monte Carlo simulations to study pure 2D Euclidean quantum gravity with $R^2$-interaction on spherical topologies, employing Regge's formulation. We attempt to measure the string susceptibility exponent $γ_{\rm str}$ by using a finite-size scaling Ansatz in the expectation value of $R^2$, as has been done in a previous study by Bock and Vink ( hep-lat/9406018 ). By considerably extending the range and statistics of their study we find that this Ansatz is plagued by large systematic errors. The $R^2$ specific string susceptibility exponent $\GS'$ is found to agree with theoretical predictions, but its determination also is subject to large systematic errors and the presence of finite-size scaling corrections. To circumvent this obstacle we suggest a new scaling Ansatz which in principle should be able to predict both, $\GS$ and $\GS'$. First results indicate that this requires large system sizes to reduce the uncertainties in the finite-size scaling Ansätze. Nevertheless, our investigation shows that within the achievable accuracy the numerical estimates are still compatible with analytic predictions, contrary to the recent claim by Bock and Vink.

preprint1995arXiv

The Ising transition in 2D simplicial quantum gravity - can Regge calculus be right?

We report a high statistics simulation of Ising spins coupled to 2D quantum gravity in the Regge calculus approach using triangulated tori with up to $512^2$ vertices. For the constant area ensemble and the $dl/l$ functional measure we definitively can exclude the critical exponents of the Ising phase transition as predicted for dynamically triangulated surfaces. We rather find clear evidence that the critical exponents agree with the Onsager values for static regular lattices, independent of the coupling strength of an $R^2$ interaction term. For exploratory simulations using the lattice version of the Misner measure the situation is less clear.