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

Daan Frenkel

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

ResearcherUnclaimed source record
cond-mat.softcond-mat.stat-mechphysics.comp-phphysics.chem-phBiological Physicscond-mat.mtrl-scicond-mat.mes-hallcond-mat.dis-nnBiomoleculesastro-ph.IMphysics.data-anphysics.flu-dynSubcellular Processes

Catalog footprint

What is connected

50works
13topics
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

50 published item(s)

preprint2021arXiv

Temperature Protocols to Guide Selective Self-Assembly of Competing Structures

Multi-component self-assembly mixtures offer the possibility of encoding multiple target structures with the same set of interacting components. Selective retrieval of one of the stored structures has been attempted by preparing an initial state that favours the assembly of the required target, through seeding, concentration patterning or specific choices of interaction strengths. This may not be possible in an experiment where on-the-fly reconfiguration of the building blocks to switch functionality may be required. In this paper, we explore principles of inverse design of a multi-component self-assembly mixture capable of encoding two competing structures that can be selected through simple temperature protocols. We design the target structures to realise the generic situation in which one of targets has the lower nucleation barrier while the other is globally more stable. We observe that to avoid the formation of spurious or chimeric aggregates, the number of neighbouring component pairs that occur in both structures should be minimal. Our design also requires the inclusion of components that are part only of one of the target structures, but we observe that to maximize the selectivity of retrieval, the component library itself should be maximally shared by the two targets. We demonstrate that temperature protocols can be designed which lead to the formation of either one of the target structures with high selectivity. We discuss the important role played by secondary aggregation products, which we term vestigial aggregates.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

'Rocket propulsion' of Janus micro-swimmers

We report simulations of a spherical Janus particle undergoing exothermic surface reactions around one pole only. Our model excludes self-phoretic transport by design. Nevertheless, net motion occurs from direct momentum transfer between solvent and colloid, with speed scaling as the square root of the energy released during the reaction. We find that such propulsion is dominated by the system's short-time response, when neither the time dependence of the flow around the colloid nor the solvent compressibility can be ignored. Our simulations agree reasonably well with previous experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Effect of the interaction strength and anisotropy on the diffusio-phoresis of spherical colloids

Gradients in temperature, concentration or electrostatic potential cannot exert forces on a bulk fluid; they can, however, exert forces on a fluid in a microscopic boundary layer surrounding a (nano)colloidal solute, resulting in so-called phoretic flow. Here we present a simulation study of phoretic flow around a spherical colloid held fixed in a concentration gradient. We show that the resulting flow velocity depends non-monotonically on the strength of the colloid-fluid interaction. The reason for this non-monotonic dependence is that solute particles are effectively trapped in a shell around the colloid and cannot contribute to diffusio-phoresis. We also observe that the flow depends sensitively on the anisotropy of solute-colloid interaction.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Numerical method for computing the free energy of glasses

We propose a numerical technique to compute the equilibrium free energy of glasses that cannot be prepared quasi-reversibly. For such systems, standard techniques for estimating the free energy by extrapolation, cannot be used. Instead, we use a procedure that samples the equilibrium partition function of the basins of attraction of the different inherent structures (local potential energy minima) of the system. If all relevant inherent structures could be adequately sampled in the (supercooled) liquid phase, our approach would be rigorous. In any finite simulation, we will miss the lower-energy inherent structures that become dominant at very low temperatures. We find that our free energy estimates for a Kob-Andersen glass are lower than those obtained by very slow cooling, even at temperatures down to one third of the glass transition temperature. The current approach could be applied to compute the chemical potential of ultra-stable glassy materials, and should enable the estimation of their solubility.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Studying polymer diffusiophoresis with Non-Equilibrium Molecular Dynamics

We report a numerical study of the diffusiophoresis of short polymers using non-equilibrium molecular dynamics simulations. More precisely, we consider polymer chains in a fluid containing a solute which has a concentration gradient, and examine the variation of the induced diffusiophoretic velocity of the polymer chains as the interaction between the monomer and the solute is varied. We find that there is a non-monotonic relation between the diffusiophoretic mobility and the strength of the monomer-solute interaction. In addition we find a weak dependence of the mobility on the length of the polymer chain, which shows clear difference from the diffusiophoresis of a solid particle. Interestingly, the hydrodynamic flow through the polymer is much less screened than for pressure driven flows.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Using Molecular Simulation to Compute Transport Coefficients of Molecular Gases

The existing kinetic theory of gases is based on an analytical approach that becomes intractable for all but the simplest molecules. Here we propose a simple numerical scheme to compute the transport properties of molecular gases in the limit of infinite dilution. The approach that we propose is approximate, but our results for the diffusivity $D$, the viscosity $η$ and the thermal conductivity $λ$ of hard spheres, Lennard-Jones particles and rough hard spheres, agree well with the standard (lowest order) Chapman-Enskog results. We also present results for a Lennard-Jones-dimer model for nitrogen, for which no analytical results are available. In the case of poly-atomic molecules (we consider n-octane), our method remains simple and gives good predictions for the diffusivity and the viscosity. Computing the thermal conductivity of poly-atomic molecules requires an approximate treatment of their quantized internal modes. We show that a well-known approximation that relates $λ$ to $D$ and $η$, yields good results. We note that our approach should yield a lower limit to the exact value of $D$, $η$ and $λ$. Interestingly, the most sophisticated (higher-order) Chapman-Enskog results for rough hard spheres seem to violate this bound.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Information density, structure and entropy in equilibrium and non-equilibrium systems

During a spontaneous change, a macroscopic physical system will evolve towards a macro-state with more realizations. This observation is at the basis of the Statistical Mechanical version of the Second Law of Thermodynamics, and it provides an interpretation of entropy in terms of probabilities. However, we cannot rely on the statistical-mechanical expressions for entropy in systems that are far from equilibrium. In this paper, we compare various extensions of the definition of entropy, which have been proposed for non-equilibrium systems. It has recently been proposed that measures of information density may serve to quantify entropy in both equilibrium and nonequilibrium systems. We propose a new "bit-wise" method to measure the information density for off lattice systems. This method does not rely on coarse-graining of the particle coordinates. We then compare different estimates of the system entropy, based on information density and on the structural properties of the system, and check if the various entropies are mutually consistent and, importantly, whether they can detect non-trivial ordering phenomena. We find that, except for simple (one-dimensional) cases, the different methods yield answers that are at best qualitatively similar, and often not even that, although in several cases, different entropy estimates do detect ordering phenomena qualitatively. Our entropy estimates based on bit-wise data compression contain no adjustable scaling factor, and show large quantitative differences with the thermodynamic entropy obtained from equilibrium simulations. Hence, our results suggest that, at present, there is not yet a single, structure-based entropy definition that has general validity for equilibrium and non equilibrium systems.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

The Lennard-Jones potential: when (not) to use it

The Lennard-Jones 12-6 potential (LJ) is arguably the most widely used pair potential in Molecular Simulations. In fact, it is so popular that the question is rarely asked whether it is fit for purpose. In this paper, we argue that, whilst the LJ potential was designed for noble gases such as argon, it is often used for systems where it is not expected to be particularly realistic. Under those circumstances, the disadvantages of the LJ potential become relevant: most important among these is that in simulations the LJ potential is always modified such that it has a finite range. More seriously, there is by now a whole family of different potentials that are all called Lennard-Jones 12-6, and that are all different - and that may have very different macroscopic properties. In this paper, we consider alternatives to the LJ 12-6 potential that could be employed under conditions where the LJ potential is only used as a typical short-ranged potential with attraction. We construct a class of potentials that are, in many respects LJ-like but that are by construction finite ranged, vanishing quadratically at the cut-off distance, and that are designed to be computationally cheap. Below, we present this potential and report numerical data for its thermodynamic and transport properties, for the most important cases: cut-off distance rc=2 (LJ-like) and rc=1.2 (a typical colloidal potential).

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

DNA brick self-assembly with an off-lattice potential

We report Monte Carlo simulations of a simple off-lattice patchy-particle model for DNA `bricks'. We relate the parameters that characterise this model with the binding free energy of pairs of single-stranded DNA molecules. We verify that an off-lattice potential parameterised in this way reproduces much of the behaviour seen with a simpler lattice model we introduced previously, although the relaxation of the geometric constraints leads to a more error-prone self-assembly pathway. We investigate the self-assembly process as a function of the strength of the non-specific interactions. We show that our off-lattice model for DNA bricks results in robust self-assembly into a variety of target structures.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Hot nano-particles in polar or paramagnetic liquids interact as monopoles

When neutral nano-particles are heated or cooled in a polar liquid, they will interact with each other as if they carry an electrostatic charge that is proportional to the temperature difference between the particle and the surrounding fluid. The same should hold for paramagnetic liquids, in which case the nano-particles should behave as magnetic monopoles. However, the analogy with electrostatics/magnetostatics is not complete: heated/cooled nano-particles do not move under the influence of an applied field. They should, however, interact as monopoles.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Nanoparticle ordering in sandwiched polymer brushes

The organization of nano-particles inside grafted polymer layers is governed by the interplay of polymer-induced entropic interactions and the action of externally applied fields. Earlier work had shown that strong external forces can drive the formation of colloidal structures in polymer brushes. Here we show that external fields are not essential to obtain such colloidal patterns: we report Monte Carlo and Molecular dynamics simulations that demonstrate that ordered structures can be achieved by compressing a `sandwich' of two grafted polymer layers, or by squeezing a coated nanotube, with nano-particles in between. We show that the pattern formation can be efficiently controlled by the applied pressure, while the characteristic length--scale, i.e. the typical width of the patterns, is sensitive to the length of the polymers. Based on the results of the simulations, we derive an approximate equation of state for nano-sandwiches.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Non-equilibrium simulations of thermally induced electric fields in water

Using non-equilibrium molecular dynamics simulations, it has been recently demonstrated that water molecules align in response to an imposed temperature gradient, resulting in an effective electric field. Here, we investigate how thermally induced fields depend on the underlying treatment of long-ranged interactions. For the short-ranged Wolf method and Ewald summation, we find the peak strength of the field to range between $2 \times 10^7$ and $5 \times 10^7~\text{V/m}$ for a temperature gradient of $5.2~\text{K/}\unicode{x212B}$. Our value for the Wolf method is therefore an order of magnitude lower than the literature value [J. Chem. Phys. 139, 014504 (2013) and 143, 036101 (2015)]. We show that this discrepancy can be traced back to the use of an incorrect kernel in the calculation of the electrostatic field. More seriously, we find that the Wolf method fails to predict correct molecular orientations, resulting in dipole densities with opposite sign to those computed using Ewald summation. By considering two different multipole expansions, we show that, for inhomogeneous polarisations, the quadrupole contribution can be significant and even outweigh the dipole contribution to the field. Finally, we propose a more accurate way of calculating the electrostatic potential and the field. In particular, we show that averaging the microscopic field analytically to obtain the macroscopic Maxwell field reduces the error bars by up to an order of magnitude. As a consequence, the simulation times required to reach a given statistical accuracy decrease by up to two orders of magnitude.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Oligomers of heat-shock proteins: Structures that don't imply function

Most proteins must remain soluble in the cytosol in order to perform their biological functions. To protect against undesired protein aggregation, living cells maintain a population of molecular chaperones that ensure the solubility of the proteome. Here we report simulations of a lattice model of interacting proteins to understand how low concentrations of passive molecular chaperones, such as small heat-shock proteins, suppress thermodynamic instabilities in protein solutions. Given fixed concentrations of chaperones and client proteins, the solubility of the proteome can be increased by tuning the chaperone--client binding strength. Surprisingly, we find that the binding strength that optimizes solubility while preventing irreversible chaperone binding also promotes the formation of weakly bound chaperone oligomers, although the presence of these oligomers does not significantly affect the thermodynamic stability of the solution. Such oligomers are commonly observed in experiments on small heat-shock proteins, but their connection to the biological function of these chaperones has remained unclear. Our simulations suggest that this clustering may not have any essential biological function, but rather emerges as a natural side-effect of optimizing the thermodynamic stability of the proteome.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Simple approach for calculating the binding free energy of a multivalent particle

We present a simple yet accurate numerical approach to compute the free energy of binding of multivalent objects on a receptor-coated surface. The method correctly accounts for the fact that one ligand can bind to at most one receptor. The numerical approach is based on a saddle-point approximation to the computation of a complex residue. We compare our theory with the powerful Valence-Limited Interaction Theory (VLIT) (J. Chem. Phys. 137, 094108(2012), J. Chem. Phys. 138, 021102(2013)) and find excellent agreement in the regime where that theory is expected to work. However, the present approach even works for low receptor/ligand densities, where VLIT breaks down.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Structural analysis of high-dimensional basins of attraction

We propose an efficient Monte Carlo method for the computation of the volumes of high-dimensional bodies with arbitrary shape. We start with a region of known volume within the interior of the manifold and then use the multistate Bennett acceptance-ratio method to compute the dimensionless free-energy difference between a series of equilibrium simulations performed within this object. The method produces results that are in excellent agreement with thermodynamic integration, as well as a direct estimate of the associated statistical uncertainties. The histogram method also allows us to directly obtain an estimate of the interior radial probability density profile, thus yielding useful insight into the structural properties of such a high-dimensional body. We illustrate the method by analyzing the effect of structural disorder on the basins of attraction of mechanically stable packings of soft repulsive spheres.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Switch-like surface binding of competing multivalent particles

Multivalent particles competing for binding on the same surface can exhibit switch-like behaviour, depending on the concentration of receptors on the surface. When the receptor concentration is low, energy dominates the free energy of binding, and particles having a small number of strongly-binding ligands preferentially bind to the surface. At higher receptor concentrations, multivalent effects become significant, and entropy dominates the binding free energy; particles having many weakly-binding ligands preferentially bind to the surface. Between these two regimes there is a "switch-point", at which the surface binds the two species of particles equally strongly. We demonstrate that a simple theory can account for this switch-like behaviour and present numerical calculations that support the theoretical predictions. We argue that binding selectivity based on receptor density, rather than identity, may have practical applications.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Theory and Simulation of DNA-Coated Colloids: a Guide for Rational Design

By exploiting the exquisite selectivity of DNA hybridization, DNA-Coated Colloids (DNACCs) can be made to self-assemble in a wide variety of structures. The beauty of this system stems largely from its exceptional versatility and from the fact that a proper choice of the grafted DNA sequences yields fine control over the colloidal interactions. Theory and simulations have an important role to play in the optimal design of self- assembling DNACCs. At present, the powerful model-based design tools are not widely used, because the theoretical literature is fragmented and the connection between different theories is often not evident. In this Perspective, we aim to discuss the similarities and differences between the different models that have been described in the literature, their underlying assumptions, their strengths and their weaknesses. Using the tools described in the present Review, it should be possible to move towards a more rational design of novel self-assembling structures of DNACCs and, more generally, of systems where ligand-receptors bonds are used to control interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Turning intractable counting into sampling: computing the configurational entropy of three-dimensional jammed packings

We report a numerical calculation of the total number of disordered jammed configurations $Ω$ of $N$ repulsive, three-dimensional spheres in a fixed volume $V$. To make these calculations tractable, we increase the computational efficiency of the approach of Xu et al. (Phys. Rev. Lett. 106, 245502 (2011)) and Asenjo et al. (Phys. Rev. Lett. 112, 098002 (2014)) and we extend the method to allow computation of the configurational entropy as a function of pressure. The approach that we use computes the configurational entropy by sampling the absolute volume of basins of attraction of the stable packings in the potential energy landscape. We find a surprisingly strong correlation between the pressure of a configuration and the volume of its basin of attraction in the potential energy landscape. This relation is well described by a power law. Our methodology to compute the number of minima in the potential energy landscape should be applicable to a wide range of other enumeration problems in statistical physics, string theory, cosmology and machine learning, that aim to find the distribution of the extrema of a scalar cost function that depends on many degrees of freedom.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Consistent treatment of hydrophobicity in protein lattice models accounts for cold denaturation

The hydrophobic effect stabilizes the native structure of proteins by minimizing the unfavourable interactions between hydrophobic residues and water through the formation of a hydrophobic core. Here we include the entropic and enthalpic contributions of the hydrophobic effect explicitly in an implicit solvent model. This allows us to capture two important effects: a length-scale dependence and a temperature dependence for the solvation of a hydrophobic particle. This consistent treatment of the hydrophobic effect explains cold denaturation and heat capacity measurements of solvated proteins.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Effects of co-ordination number on the nucleation behaviour in many-component self-assembly

We report canonical and grand-canonical lattice Monte Carlo simulations of the self-assembly of addressable structures comprising hundreds of distinct component types. The nucleation behaviour, in the form of free-energy barriers to nucleation, changes significantly as the co-ordination number of the building blocks is changed from 4 to 8 to 12. Unlike tetrahedral structures - which roughly correspond to DNA bricks that have been studied in experiment - the shapes of the free-energy barriers of higher co-ordination structures depend strongly on the supersaturation, and such structures require a very significant driving force for structure growth before nucleation becomes thermally accessible. Although growth at high supersaturation results in more defects during self-assembly, we show that high co-ordination number structures can still be assembled successfully in computer simulations and that they exhibit self-assembly behaviour analogous to DNA bricks. In particular, the self-assembly remains modular, enabling in principle a wide variety of nanostructures to be assembled, with a greater spatial resolution than is possible in low co-ordination structures.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Evidence for non-ergodicity in quiescent states of periodically sheared suspensions

We present simulations of an equilibrium statistical-mechanics model that uniformly samples the space of quiescent states of a periodically sheared suspension. In our simulations, we compute the structural properties of this model as a function of density. We compare the results of our simulations with the structural data obtained in the corresponding non-equilibrium model of Corté et al. [Nat. Phys. 4, 420 (2008)]. We find that the structural properties of the non-equilibrium model are very different from those of the equilibrium model, even though the two models have exactly the same set of accessible states. This observation shows that the dynamical protocol does not sample all quiescent states with equal probability. In particular, we find that, whilst quiescent states prepared in a non-equilibrium protocol can be hyperuniform [see Phys. Rev. Lett. 114, 110602 (2015), Phys. Rev. Lett. 114, 148301 (2015), and Phys. Rev. Lett. 115, 108301 (2015)], ergodic sampling never leads to hyperuniformity. In addition, we observe ordering phase transitions and a percolation transition in the equilibrium model that do not show up in the non-equilibrium model. Conversely, the quiescent-to-diffusive transition in the dynamical model does not correspond to a phase transition, nor a percolation transition, in the equilibrium model.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Rational design of self-assembly pathways for complex multicomponent structures

The field of complex self-assembly is moving toward the design of multi-particle structures consisting of thousands of distinct building blocks. To exploit the potential benefits of structures with such `addressable complexity,' we need to understand the factors that optimize the yield and the kinetics of self-assembly. Here we use a simple theoretical method to explain the key features responsible for the unexpected success of DNA-brick experiments, which are currently the only demonstration of reliable self-assembly with such a large number of components. Simulations confirm that our theory accurately predicts the narrow temperature window in which error-free assembly can occur. Even more strikingly, our theory predicts that correct assembly of the complete structure may require a time-dependent experimental protocol. Furthermore, we predict that low coordination numbers result in non-classical nucleation behavior, which we find to be essential for achieving optimal nucleation kinetics under mild growth conditions. We also show that, rather surprisingly, the use of heterogeneous bond energies improves the nucleation kinetics and in fact appears to be necessary for assembling certain intricate three-dimensional structures. This observation makes it possible to sculpt nucleation pathways by tuning the distribution of interaction strengths. These insights not only suggest how to improve the design of structures based on DNA bricks, but also point the way toward the creation of a much wider class of chemical or colloidal structures with addressable complexity.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Theoretical prediction of free-energy landscapes for complex self-assembly

We present a technique for calculating free-energy profiles for the nucleation of multicomponent structures that contain as many species as building blocks. We find that a key factor is the topology of the graph describing the connectivity of the target assembly. By considering the designed interactions separately from weaker, incidental interactions, our approach yields predictions for the equilibrium yield and nucleation barriers. These predictions are in good agreement with corresponding Monte Carlo simulations. We show that a few fundamental properties of the connectivity graph determine the most prominent features of the assembly thermodynamics. Surprisingly, we find that polydispersity in the strengths of the designed interactions stabilizes intermediate structures and can be used to sculpt the free-energy landscape for self-assembly. Finally, we demonstrate that weak incidental interactions can preclude assembly at equilibrium due to the combinatorial possibilities for incorrect association.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Accounting for adsorption and desorption in lattice Boltzmann simulations

We report a Lattice-Boltzmann scheme that accounts for adsorption and desorption in the calculation of mesoscale dynamical properties of tracers in media of arbitrary complexity. Lattice Boltzmann simulations made it possible to solve numerically the coupled Navier-Stokes equations of fluid dynamics and Nernst-Planck equations of electrokinetics in complex, heterogeneous media. Associated to the moment propagation scheme, it became possible to extract the effective diffusion and dispersion coefficients of tracers, or solutes, of any charge, e.g. in porous media. Nevertheless, the dynamical properties of tracers depend on the tracer-surface affinity, which is not purely electrostatic, but also includes a species-specific contribution. In order to capture this important feature, we introduce specific adsorption and desorption processes in a Lattice-Boltzmann scheme through a modified moment propagation algorithm, in which tracers may adsorb and desorb from surfaces through kinetic reaction rates. The method is validated on exact results for pure diffusion and diffusion-advection in Poiseuille flows in a simple geometry. We finally illustrate the importance of taking such processes into account on the time-dependent diffusion coefficient in a more complex porous medium.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Crucial role of non-specific interactions in amyloid nucleation

Protein oligomers have been implicated as toxic agents in a wide range of amyloid-related diseases. Yet it has remained unsolved whether the oligomers are a necessary step in the formation of amyloid fibrils, or just a dangerous by-product. Analogously, it has not been resolved if the amyloid nucleation process is a classical one-step nucleation process, or a two-step process involving pre-nucleation clusters. We use coarse-grained computer simulations to study the effect of non-specific attractions between peptides on the primary nucleation process underlying amyloid fibrillization. We find that for peptides that do not attract, the classical one-step nucleation mechanism is possible, but only at non-physiologically high peptide concentrations. At low peptide concentrations, which mimic the physiologically relevant regime, attractive inter-peptide interactions are essential for fibril formation. Nucleation then inevitably takes place through a two-step mechanism involving prefibrillar oligomers. We show that oligomers not only help peptides meet each other, but create an environment that facilitates the conversion of monomers into the $β$-sheet rich form characteristic of fibrils. Nucleation typically does not proceed via the most prevalent oligomers, but via an oligomer size that is only observed in rare fluctuations, which is why such aggregates might be hard to capture experimentally. Finally, we find that the nucleation of amyloid fibrils cannot be described by classical nucleation theory: in the two-step mechanism the critical nucleus size increases both with an increase in concentration and in the inter-peptide interactions, in direct contrast with predictions from classical nucleation theory.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Designing stimulus-sensitive colloidal walkers

Colloidal particles with DNA `legs' that can bind reversibly to receptors on a surface can be made to `walk' if there is a gradient in receptor concentration. We use a combination of theory and Monte Carlo simulations to explore how controllable parameters, e.g. coating density and binding strength, affect the dynamics of such colloids. We find that competition between thermodynamic and kinetic trends imply that there is an optimal value for both, the binding strength and the number of `legs' for which transport is fastest. Using available thermodynamic data on DNA binding, we indicate how directionally reversible, temperature-controlled transport of colloidal walkers can be achieved. In particular, the present results should make it possible to design a chromatographic technique that can be used to separate colloids with different DNA functionalization.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Gibbs, Boltzmann, and negative temperatures

In a recent paper, Dunkel and Hilbert [Nature Physics 10, 67-72 (2014)] use an entropy definition due to Gibbs to provide a 'consistent thermostatistics' which forbids negative absolute temperatures. Here we argue that the Gibbs entropy fails to satisfy a basic requirement of thermodynamics, namely that when two bodies are in thermal equilibrium, they should be at the same temperature. The entropy definition due to Boltzmann does meet this test, and moreover in the thermodynamic limit can be shown to satisfy Dunkel and Hilbert's consistency criterion. Thus, far from being forbidden, negative temperatures are inevitable, in systems with bounded energy spectra.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Lattice simulation method to model diffusion and NMR spectra in porous materials

A coarse-grained simulation method to predict NMR spectra of ions diffusing in porous carbons is proposed. The coarse-grained model uses input from molecular dynamics simulations such as the free-energy profile for ionic adsorption, and density-functional theory calculations are used to predict the NMR chemical shift of the diffusing ions. The approach is used to compute NMR spectra of ions in slit pores with pore widths ranging from 2 to 10 nm. As diffusion inside pores is fast, the NMR spectrum of an ion trapped in a single mesopore will be a sharp peak with a pore size dependent chemical shift. To account for the experimentally observed NMR line shapes, our simulations must model the relatively slow exchange between different pores. We show that the computed NMR line shapes depend on both the pore size distribution and the spatial arrangement of the pores. The technique presented in this work provides a tool to extract information about the spatial distribution of pore sizes from NMR spectra. Such information is diffcult to obtain from other characterisation techniques.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Mobile linkers on DNA-coated colloids: valency without patches

Colloids coated with single-stranded DNA (ssDNA) can bind selectively to other colloids coated with complementary ssDNA. The fact that DNA-coated colloids (DNACCs) can bind to specific partners opens the prospect of making colloidal `molecules'. However, in order to design DNACC-based molecules, we must be able to control the valency of the colloids, i.e. the number of partners to which a given DNACC can bind. One obvious, but not very simple approach is to decorate the colloidal surface with patches of single-stranded DNA that selectively bind those on other colloids. Here we propose a design principle that exploits many-body effects to control the valency of otherwise isotropic colloids. Using a combination of theory and simulation, we show that we can tune the valency of colloids coated with mobile ssDNA, simply by tuning the non-specific repulsion between the particles. Our simulations show that the resulting effective interactions lead to low-valency colloids self-assembling in peculiar open structures, very different from those observed in DNACCs with immobile DNA linkers.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Numerical calculation of granular entropy

We present numerical simulations that allow us to compute the number of ways in which $N$ particles can pack into a given volume $V$. Our technique modifies the method of Xu et al. (Phys. Rev. Lett. 106, 245502 (2011)) and outperforms existing direct enumeration methods by more than 200 orders of magnitude. We use our approach to study the system size dependence of the number of distinct packings of a system of up to 128 poly-disperse soft disks. We show that, even though granular particles are distinguishable, we have to include a factor $1/N!$ to ensure that the entropy does not change when exchanging particles between systems in the same macroscopic state. Our simulations provide strong evidence that the packing entropy, when properly defined, is extensive. As different packings are created with unequal probabilities, it is natural to express the packing entropy as $S=-\sum_i p_i\ln p_i -\ln N!$, where $p_i$ denotes the probability to generate the $i$-th packing. We can compute this quantity reliably and it is also extensive. The granular entropy thus (re)defined, whilst distinct from the one proposed by Edwards assumed (J. Phys.: Condens.Matter 2, SA63(1990)), does have all the properties Edwards assumed.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Numerical evidence for nucleated self-assembly of DNA brick structures

The observation by Ke et al. [Science 338, 1177 (2012)] that large numbers of short, pre-designed DNA strands can assemble into three-dimensional target structures came as a great surprise, as no colloidal self-assembling system has ever achieved the same degree of complexity. That failure seemed easy to rationalise: the larger the number of distinct building blocks, the higher the expected error rate for self-assembly. The experiments of Ke et al. have disproved this argument. Here, we report Monte Carlo simulations of the self-assembly of a DNA brick cube, comprising approximately 1000 types of DNA strand, using a simple model. We model the DNA strands as lattice tetrahedra with attractive patches, the interaction strengths of which are computed using a standard thermodynamic model. We find that, within a narrow temperature window, the target structure assembles with high probability. Our simulations suggest that mis-assembly is disfavoured because of a slow nucleation step. As our model incorporates no aspect of DNA other than its binding properties, these simulations suggest that, with proper design of the building blocks, other systems, such as colloids, may also assemble into truly complex structures.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Superposition Enhanced Nested Sampling

The theoretical analysis of many problems in physics, astronomy and applied mathematics requires an efficient numerical exploration of multimodal parameter spaces that exhibit broken ergodicity. Monte Carlo methods are widely used to deal with these classes of problems, but such simulations suffer from a ubiquitous sampling problem: the probability of sampling a particular state is proportional to its entropic weight. Devising an algorithm capable of sampling efficiently the full phase space is a long-standing problem. Here we report a new hybrid method for the exploration of multimodal parameter spaces exhibiting broken ergodicity. Superposition enhanced nested sampling (SENS) combines the strengths of global optimization with the unbiased/athermal sampling of nested sampling, greatly enhancing its efficiency with no additional parameters. We report extensive tests of this new approach for atomic clusters that are known to have energy landscapes for which conventional sampling schemes suffer from broken ergodicity. We also introduce a novel parallelization algorithm for nested sampling.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A simple analytical formula for the free-energy of ligand-receptor mediated interactions

Recently \1, we presented a general theory for calculat- ing the strength and properties of colloidal interactions mediated by ligand-receptor bonds (such as those that bind DNA-coated colloids). In this communication, we derive a surprisingly simple analytical form for the inter- action free energy, which was previously obtainable only via a costly numerical thermodynamic integration. As a result, the computational effort to obtain potentials of in- teraction is significantly reduced. Moreover, we can gain insight from this analytic expression for the free energy in limiting cases. In particular, the connection of our general theory to other previous specialised approaches is now made transparent. This important simplification will significantly broaden the scope of our theory.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Procedure to construct a multi-scale coarse-grained model of DNA-coated colloids from experimental data

We present a quantitative, multi-scale coarse-grained model of DNA coated colloids. The parameters of this model are transferable and are solely based on experimental data. As a test case, we focus on nano-sized colloids carrying single-stranded DNA strands of length comparable to the colloids' size. We show that in this regime, the common theoretical approach of assuming pairwise additivity of the colloidal pair interactions leads to quantitatively and sometimes even qualitatively wrong predictions of the phase behaviour of DNA-grafted colloids. Comparing to experimental data, we find that our coarse-grained model correctly predicts the equilibrium structure and melting temperature of the formed solids. Due to limited experimental information on the persistence length of single-stranded DNA, some quantitative discrepancies are found in the prediction of spatial quantities. With the availability of better experimental data, the present approach provides a path for the rational design of DNA-functionalised building blocks that can self-assemble in complex, three-dimensional structures.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Visualizing Basins of Attraction for Different Minimization Algorithms

We report a study of the basins of attraction for potential energy minima defined by different minimisation algorithms for an atomic system. We find that whereas some minimisation algorithms produce compact basins, others produce basins with complex boundaries or basins consisting of disconnected parts. Such basins deviate from the "correct" basin of attraction defined by steepest-descent pathways, and the differences can be controlled to some extent by adjustment of the maximum step size. The choice of the most convenient minimisation algorithm depends on the problem in hand. We show that while L-BFGS is the fastest minimiser, the FIRE algorithm is also quite fast, and can lead to less fragmented basins of attraction.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Why Colloidal Systems can be described by Statistical Mechanics: Some not very original comments on the Gibbs paradox

Colloidal particles are distinguishable. Moreover, their thermodynamic properties are extensive. Statistical Mechanics predicts such behaviour if one accepts that the configurational integral of a system of N colloids must be divided by N!. In many textbooks it is argued that the factor N! corrects for the fact that identical particles (in the quantum mechanical sense) are indistinguishable. Clearly, this argument does not apply to colloids. This articles explains why, nevertheless, all is well. The point has been made before, but has not yet sunk in. I also discuss the effect of polydispersity.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A general theory of DNA-mediated and other valence-limited interactions

We present a general theory for predicting the interaction potentials between DNA-coated colloids, and more broadly, any particles that interact via valence-limited ligand-receptor binding. Our theory correctly incorporates the configurational and combinatorial entropic factors that play a key role in valence-limited interactions. By rigorously enforcing self-consistency, it achieves near-quantitative accuracy with respect to detailed Monte Carlo calculations. With suitable approximations and in particular geometries, our theory reduces to previous successful treatments, which are now united in a common and extensible framework. We expect our tools to be useful to other researchers investigating ligand-mediated interactions. A complete and well-documented Python implementation is freely available at http://github.com/patvarilly/DNACC .

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A parameter-free, solid-angle based, nearest-neighbor algorithm

We propose a parameter-free algorithm for the identification of nearest neighbors. The algorithm is very easy to use and has a number of advantages over existing algorithms to identify nearest- neighbors. This solid-angle based nearest-neighbor algorithm (SANN) attributes to each possible neighbor a solid angle and determines the cutoff radius by the requirement that the sum of the solid angles is 4π. The algorithm can be used to analyze 3D images, both from experiments as well as theory, and as the algorithm has a low computational cost, it can also be used "on the fly" in simulations. In this paper, we describe the SANN algorithm, discuss its properties, and compare it to both a fixed-distance cutoff algorithm and to a Voronoi construction by analyzing its behavior in bulk phases of systems of carbon atoms, Lennard-Jones particles and hard spheres as well as in Lennard-Jones systems with liquid-crystal and liquid-vapor interfaces.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Collective sedimentation of colloids in grafted polymer layers

We present Monte Carlo simulations of colloidal particles pulled into grafted polymer layers by external fields. The insertion free energy of a single colloid into the polymer layer is qualitatively different for surfaces with an ordered and a disordered distribution of grafting points. Moreover, the tendency of colloidal particles to traverse the grafting layer is strongly size dependent. In dense colloidal suspensions, under the influence of sufficiently strong external force, a collective instability allows the colloids to penetrate and form internally ordered, columnar structures spanning the polymer layer. Depending on the conditions, these colloidal clusters may be isolated or laterally percolating. The morphology of the observed patterns can be controlled by the external fields, which opens up new routes for the design of thin structured films.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Emergence of Critical Phenomena in Full Configuration Interaction Quantum Monte Carlo

There has been recent literature discussion on the origin and severity of the `sign problem' in full configuration interaction quantum Monte Carlo (FCIQMC) and its `initiator' adaptation (i-FCIQMC), methods of interest and potential because they allow for exact (FCI) ground-state solutions to be obtained often at a much reduced computational cost. In this study we aim to use a simple order parameter, describing the `sign structure' of the stochastic wavefunction representation, to empirically characterise the fundamentally different collective behaviour of the walker population in both methods.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Layering, freezing and re-entrant melting of hard spheres in soft confinement

Confinement can have a dramatic effect on the behavior of all sorts of particulate systems and it therefore is an important phenomenon in many different areas of physics and technology. Here, we investigate the role played by the softness of the confining potential. Using grand canonical Monte Carlo simulations, we determine the phase diagram of three-dimensional hard spheres that in one dimension are constrained to a plane by a harmonic potential. The phase behavior depends strongly on the density and on the stiffness of the harmonic confinement. Whilst we find the familiar sequence of confined hexagonal and square-symmetric packings, we do not observe any of the usual intervening ordered phases. Instead, the system phase separates under strong confinement, or forms a layered re-entrant liquid phase under weaker confinement. It is plausible that this behavior is due to the larger positional freedom in a soft confining potential and to the contribution that the confinement energy makes to the total free energy. The fact that specific structures can be induced or suppressed by simply changing the confinement conditions (e.g. in a dielectrophoretic trap) is important for applications that involve self-assembled structures of colloidal particles.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Probing ergodicity in granular matter

When a granular system is tapped, its volume changes. Here, using a well-defined macroscopic protocol, we prepare an ensemble of granular systems and track the statistics of volume changes as a function of the number of taps. This is in contrast to previous studies, which have focused on single trajectories and assumed ergodicity. We devise a new method to assess the convergence properties of a sequence of ensemble volume histograms and introduce a reasonable approximate version of an invariant histogram. We then compare these invariant histograms with histograms generated by sampling a long trajectory for one system and observe nonergodicity, which we quantify. Finally, we use the overlapping histogram method to assess potential compatibility with Edwards' canonical assumption. Our histograms are incompatible with this assumption.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Quantitative prediction of the phase diagram of DNA-functionalized nano-colloids

We present a coarse-grained model of DNA-functionalized colloids that is computationally tractable. Importantly, the model parameters are solely based on experimental data. Using this highly simplified model, we can predict the phase behavior of DNA-functionalized nano-colloids without assuming pairwise additivity of the inter-colloidal interactions. Our simulations show that for nano-colloids, the assumption of pairwise additivity leads to substantial errors in the estimate of the free energy of the crystal phase. We compare our results with available experimental data and find that the simulations predict the correct structure of the solid phase and yield a very good estimate of the melting temperature. Current experimental estimates for the contour length and persistence length of single-stranded DNA sequences are subject to relatively large uncertainties. Using the best available estimates, we obtain predictions for the crystal lattice constants that are off by a few percent: this indicates that more accurate experimental data on ssDNA are needed to exploit the full power of our coarse-grained approach.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Running faster together: huge speed up of thermal ratchets due to hydrodynamic coupling

We present simulations that reveal a surprisingly large effect of hydrodynamic coupling on the speed of thermal ratchet motors. The model that we use considers particles performing thermal ratchet motion in a hydrodynamic solvent. Using particle-based, mesoscopic simulations that maintain local momentum conservation, we analyze quantitatively how the coupling to the surrounding fluid affects ratchet motion. We find that coupling can increase the mean velocity of the moving particles by almost two orders of magnitude, precisely because ratchet motion has both a diffusive and a deterministic component. The resulting coupling also leads to the formation of aggregates at longer times. The correlated motion that we describe increases the efficiency of motor-delivered cargo transport and we speculate that the mechanism that we have uncovered may play a key role in speeding up molecular motor-driven intracellular transport.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Simulations: the dark side

This paper discusses the Monte Carlo and Molecular Dynamics methods. Both methods are, in principle, simple. However, simple does not mean risk-free. In the literature, many of the pitfalls in the field are mentioned, but usually as a footnote - and these footnotes are scattered over many papers. The present paper focuses on the `dark side' of simulation: it is one big footnote. I should stress that `dark', in this context, has no negative moral implication. It just means: under-exposed.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Direct determination of the size of basins of attraction of jammed solids

We propose a free-energy based Monte-Carlo method to measure the volume of potential-energy basins in configuration space. Using this approach we can estimate the number of distinct potential-energy minima, even when this number is much too large to be sampled directly. We validate our approach by comparing our results with the direct enumeration of distinct jammed states in small packings of frictionless spheres. We find that the entropy of distinct packings is extensive and that the entropy of distinct hard-sphere packings must have a maximum as a function of packing fraction.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Pair interactions between complex mesoscopic particles from Widom's particle-insertion method

We demonstrate that Widom's particle insertion technique provides a convenient and efficient method to determine the effective pair interaction between complex, composite soft-matter particles in the zero-density limit. By means of three different test systems, i.e. amphiphilic dendrimers, electrostatic polymers and colloids coated with electrostatic polymers, we demonstrate the validity and the power of the presented method.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Real-time monitoring of complex moduli from micro-rheology

We describe an approach to online analysis of micro-rheology data using a multi-scale time-correlation method. The method is particularly suited to process high-volume data streams and compress the relevant information in real time. Using this, we can obtain complex moduli of visco-elastic media without suffering from the high-frequency artefacts that are associated with the truncation errors in the most widely used versions of micro-rheology. Moreover, the present approach obviates the need to choose the time interval for data acquisition beforehand. We test our approach first on an artificial data set and then on experimental data obtained both for an optically trapped colloidal probe in water and a similar probe in poly-ethylene glycol solutions at various concentrations. In all cases, we obtain good agreement with the bulk rheology data in the region of overlap. We compare our method with the conventional Kramers-Kronig transform approach and find that the two methods agree over most of the frequency regime. For the same data set, the present approach is superior to Kramers-Kronig at high frequencies and can be made to perform at least comparable at low frequencies.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Waterlike thermodynamic anomalies in a repulsive-step potential system

We report a computer-simulation study of the equilibrium phase diagram of a three-dimensional system of particles with a repulsive step potential. The phase diagram is obtained using free-energy calculations. At low temperatures, we observe a number of distinct crystal phases. We show that at certain values of the potential parameters the system exhibits the water-like thermodynamic anomalies: density anomaly and diffusion anomaly. The anomalies disappear with increasing the repulsive step width: their locations move to the region inside the crystalline phase.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Evidence for out-of-equilibrium crystal nucleation in suspensions of oppositely charged colloids

We report a numerical study of the rate of crystal nucleation in a binary suspension of oppositely charged colloids. Two different crystal structures compete in the thermodynamic conditions under study. We find that the crystal phase that nucleates is metastable and, more surprisingly, its nucleation free energy barrier is not the lowest one. This implies that, during nucleation, there is insufficient time for sub-critical nuclei to relax to their lowest free-energy structure. Such behavior is in direct contradiction with the common assumption that the phase that crystallizes most readily is the one with the lowest free-energy barrier for nucleation. The phenomenon that we describe should be relevant for crystallization experiments where competing solid structures are not connected by an easy transformation.

0 citations0 reviewsNo code linkedOpen access