BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Gary W. Slater

Gary W. Slater contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
cond-mat.softBiological Physicsphysics.comp-phcond-mat.stat-mechphysics.chem-phBiomoleculesphysics.data-an

Trust snapshot

Quick read

Trust 21 - Emerging
13works
0followers
7topics
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

13 published item(s)

preprint2022arXiv

An empirical method to characterize displacement distribution functions for anomalous and transient diffusion

We propose a practical empirical fitting function to characterize the non-Gaussian displacement distribution functions (DispD) often observed for heterogeneous diffusion problems. We first test this fitting function with the problem of a colloidal particle diffusing between two walls using Langevin Dynamics (LD) simulations of a raspberry particle coupled to a lattice Boltzmann (LB) fluid. We also test the function with a simple model of anomalous diffusion on a square lattice with obstacles. In both cases, the fitting parameters provide more physical information than just the Kurtosis (which is often the method used to quantify the degree of anomaly of the dynamics), including a length scale that marks where the tails of the DispD begin. In all cases, the fitting parameters smoothly converge to Gaussian values as the systems become less anomalous.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Capture and translocation of a rod-like molecule by a nanopore: orientation, charge distribution and hydrodynamics

We investigate the translocation of rods with different charge distributions using hybrid Langevin Dynamics and Lattice Boltzmann (LD-LB) simulations. Electrostatic interactions are added to the system using the $P^3M$ algorithm to model the electrohydrodynamic interactions (EHI). We first examine the free-solution electrophoretic properties of rods with various charge distributions. Our translocation simulation results suggest that the order parameter is asymmetric during the capture and escape processes despite the symmetric electric field lines, while the impacts of the charge distribution on rod orientation are more significant during the capture process. The capture/threading/escape times are under the combined effects of charge screening, rod orientation, and charge distributions. We also show that the mean capture time of a rod is shorter when it is launched near the wall because rods tend to align along the wall and hence with the local field lines. Remarkably, the \textit{orientational capture radius} we proposed previously for uniformly charged rods is still valid in the presence of EHI.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Capture of rod-like molecules by a nanopore: defining an "orientational capture radius"

Both the translational diffusion coefficient $D$ and the electrophoretic mobility $μ$ of a short rod-like molecule (such as dsDNA) that is being pulled towards a nanopore by an electric field should depend on its orientation. Since a charged rod-like molecule tends to orient in the presence of an inhomogeneous electric field, $D$ and $μ$ will change as the molecule approaches the nanopore, and this will impact the capture process. We present a simplified study of this problem using theoretical arguments and Langevin Dynamics simulations. In particular, we introduce a new \textit{orientational capture radius} which we compare to the capture radius for the equivalent point-like particle, and we discuss the different physical regimes of orientation during capture and the impact of initial orientations on the capture time.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Voltage-Driven Translocation: Defining a Capture Radius

Analyte translocation involves three phases: (i) diffusion in the loading solution; (ii) capture by the pore; (iii) threading. The capture process remains poorly characterized because it cannot easily be visualized or inferred from indirect measurements. The capture performance of a device is often described by a \textit{capture radius} generally defined as the radial distance $R^*$ at which diffusion-dominated dynamics cross over to field-induced drift. However, this definition is rather ambiguous and the related models are usually over-simplified and studied in the steady-state limit. We investigate different approaches to defining and estimating $R^*$ for a charged particle diffusing in a liquid and attracted to the nanopore by the electric field. We present a theoretical analysis of the Péclet number as well as Monte Carlo simulations with different simulation protocols. Our analysis shows that the boundary conditions, pore size and finite experimental times all matter in the interpretation and calculation of $R^*$.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Electrophoretic Mobility of Polyelectrolytes within a Confining Well

We present a numerical study of polyelectrolytes electrophoresing in free solution while squeezed by an axisymmetric confinement force transverse to their net displacement. Hybrid multi-particle collision dynamics and molecular dynamics simulations with mean-field finite Debye layers show that even though the polyelectrolyte chains remain "free-draining", their electrophoretic mobility increases with confinement in nanoconfining potential wells. The primary mechanism leading to the increase in mobility above the free-solution value, despite long-range hydrodynamic screening by counterion layers, is the orientation of polymer segments within Debye layers. The observed length-dependence of the electrophoretic mobility arises due to secondary effects of counterion condensation related to confinement compactification.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

"Diffusing diffusivity": A model for anomalous and "anomalous yet Brownian" diffusion

Wang et al. [PNAS 106 (2009) 15160] have found that in several systems the linear time dependence of the mean-square displacement (MSD) of diffusing colloidal particles, typical of normal diffusion, is accompanied by a non-Gaussian displacement distribution (DisD), with roughly exponential tails at short times, a situation they termed "anomalous yet Brownian" diffusion. The diversity of systems in which this is observed calls for a generic model. We present such a model where there is "diffusivity memory" but no "direction memory" in the particle trajectory, and we show that it leads to both a linear MSD and a non-Gaussian DisD at short times. In our model, the diffusivity is undergoing a (perhaps biased) random walk, hence the expression "diffusing diffusivity". The DisD is predicted to be exactly exponential at short times if the distribution of diffusivities is itself exponential, but an exponential remains a good fit to the DisD for a variety of diffusivity distributions. Moreover, our generic model can be modified to produce subdiffusion.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Coarse-Grained Molecular Dynamics Simulations of Depletion-Induced Interactions for Soft Matter Systems

Given the ubiquity of depletion effects in biological and other soft matter systems, it is desirable to have coarse-grained Molecular Dynamics simulation approaches appropriate for the study of complex systems. This paper examines the use of two common truncated Lennard-Jones (WCA) potentials to describe a pair of colloidal particles in a thermal bath of depletants. The shifted-WCA model is the steeper of the two repulsive potentials considered, while the combinatorial-WCA model is the softer. It is found that the depletion-induced well depth for the combinatorial-WCA model is significantly deeper than the shifted-WCA model because the resulting overlap of the colloids yields extra accessible volume for depletants. For both shifted- and combinatorial-WCA simulations, the second virial coefficients and pair potentials between colloids are demonstrated to be well approximated by the Morphometric Thermodynamics (MT) model. This agreement suggests that the presence of depletants can be accurately modelled in MD simulations by implicitly including them through simple, analytical MT forms for depletion-induced interactions. Although both WCA potentials are found to be effective generic coarse-grained simulation approaches for studying depletion effects in complicated soft matter systems, combinatorial-WCA is the more efficient approach as depletion effects are enhanced at lower depletant densities. The findings indicate that for soft matter systems that are better modelled by potentials with some compressibility, predictions from hard-sphere systems could greatly underestimate the magnitude of depletion effects at a given depletant density.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Simulating the Entropic Collapse of Coarse-Grained Chromosomes

Depletion forces play a role in the compaction and de-compation of chromosomal material in simple cells but it remains debatable whether they are sufficient to account for chromosomal collapse. We present coarse-grained molecular dynamics simulations, which reveal that depletion-induced attraction is sufficient to cause the collapse of a flexible chain of large structural monomers immersed in a bath of smaller depletants. These simulations use an explicit coarse-grained computational model that treats both the supercoiled DNA structural monomers and the smaller protein crowding agents as combinatorial, truncated Lennard-Jones spheres. By presenting a simple theoretical model, we quantitatively cast the action of depletants on supercoiled bacterial DNA as an effective solvent quality. The rapid collapse of the simulated flexible chromosome at the predicted volume fraction of depletants is a continous phase transition. Additional physical effects to such simple chromosome models, such as enthalpic interactions between structural monomers or chain rigidity, are required if the collapse is to be a first-order phase transition.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Using a Péclet Number for the Translocation of Polymer through a Nanopore to Tune Coarse-Grained Simulations to Experimental Conditions

Coarse-grained simulations are often employed to study the translocation of DNA through a nanopore. The majority of these studies investigate the translocation process in a relatively generic sense and do not endeavour to match any particular set of experimental conditions. In this manuscript, we use the concept of a Péclet number for translocation, $P_t$, to compare the drift-diffusion balance in a typical experiment vs a typical simulation. We find that the standard coarse-grained approach over-estimates diffusion effects by anywhere from a factor of 5 to 50 compared to experimental conditions using dsDNA. By defining a coarse-graining parameter, $λ$, we are able to correct this and tune the simulations to replicate the experimental $P_t$ (for dsDNA and other scenarios). To show the effect that a particular $P_t$ can have on the dynamics of translocation, we perform simulations across a wide range of $P_t$ values for two different types of driving forces: a force applied in the pore and a pulling force applied to the end of the polymer. As $P_t$ brings the system from a diffusion dominated to a drift dominated regime, a variety of effects are observed including a non-monotonic dependence of the translocation time $τ$ on $P_t$ and a steep rise in the probability of translocating. Comparing the two force cases illustrates the impact of the crowding effects that occur on the trans side: a non-monotonic dependence of the width of the $τ$ distributions is obtained for the in-pore force but not for the pulling force.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Monte Carlo Approaches for Simulating a Particle at a Diffusivity Interface and the "Ito--Stratonovich Dilemma"

The possibility of different interpretations of the stochastic term (or calculi) in the overdamped Langevin equation for the motion of a particle in an inhomogeneous medium is often referred to as the "Ito--Stratonovich dilemma," although there is, in fact, a continuum of choices. We introduce two Monte Carlo (MC) simulation approaches for studying such systems, with both approaches giving the choice between different calculi (in particular, Ito, Stratonovich, and "isothermal"). To demonstrate these approaches, we study the diffusion on a 1D interval of a particle released at an interface (in the middle of the system) between two media where this particle has different diffusivities (e.g., two fluids with different viscosities). For reflecting boundary conditions at the ends of the interval, a discontinuity at the interface in the stationary-state particle distribution is found, except for the isothermal case, as expected. We also study the first-passage problem using absorbing boundary conditions. Good agreement is found when comparing the MC approaches to theory as well as Brownian and Langevin dynamics simulations. Also, the results themselves turn out to be interesting. For instance: 1) for some calculi, there are more particles on the low-viscosity side (LVS) at earlier times and then fewer particles at later times; 2) there is no preference to end up on a particular wall for the Ito variant, but a bias towards the wall on the LVS in all other cases; 3) the mean first-passage time to the wall on the LVS grows as the viscosity on the other side is increased, except for the isothermal case where it approaches a constant; 4) when the viscosity ratio is high, the first-passage-time distribution for the wall on the LVS is much broader than for the other wall, with a power-law dependence of the former in a certain time interval whose exponent depends on the calculus.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Translocation of a Polymer through a Nanopore across a Viscosity Gradient

The translocation of a polymer through a pore in a membrane separating fluids of different viscosities is studied via several computational approaches. Starting with the polymer halfway, we find that as a viscosity difference across the pore is introduced, translocation will predominately occur towards one side of the membrane. These results suggest an intrinsic pumping mechanism for translocation across cell walls which could arise whenever the fluid across the membrane is inhomogeneous. Somewhat surprisingly, the sign of the preferred direction of translocation is found to be strongly dependent on the simulation algorithm: for Langevin Dynamics (LD) simulations, a bias towards the low viscosity side is found while for Brownian Dynamics (BD), a bias towards the high viscosity is found. Examining the translocation dynamics in detail across a wide range of viscosity gradients and developing a simple force model to estimate the magnitude of the bias, the LD results are demonstrated to be more physically realistic. The LD results are also compared to those generated from a simple, one dimensional random walk model of translocation to investigate the role of the internal degrees of freedom of the polymer and the entropic barrier. To conclude, the scaling of the results across different polymer lengths demonstrates the saturation of the preferential direction with polymer length and the non-trivial location of the maximum in the exponent corresponding to the scaling of the translocation time with polymer length.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Optimizing the accuracy of Lattice Monte Carlo algorithms for simulating diffusion

The behavior of a Lattice Monte Carlo algorithm (if it is designed correctly) must approach that of the continuum system that it is designed to simulate as the time step and the mesh step tend to zero. However, we show for an algorithm for unbiased particle diffusion that if one of these two parameters remains fixed, the accuracy of the algorithm is optimal for a certain finite value of the other parameter. In one dimension, the optimal algorithm with moves to the two nearest neighbor sites reproduces the correct second and fourth moments (and minimizes the error for the higher moments at large times) of the particle distribution and preserves the first two moments of the first-passage time distributions. In two and three dimensions, the same level of accuracy requires simultaneous moves along two axes ("diagonal" moves). Such moves attempting to cross an impenetrable boundary should be projected along the boundary, rather than simply rejected. We also treat the case of absorbing boundaries.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access