Researcher profile

D. Frenkel

D. Frenkel contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

cond-mat.soft cond-mat.mtrl-sci cond-mat.stat-mech physics.chem-ph physics.comp-ph

Trust snapshot

Quick read

Trust 19 - Baseline

5works

0followers

5topics

4close collaborators

Actions

Decide how to stay connected

Follow researcher0

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.

preprint2015arXiv

An enhanced version of the heat exchange algorithm with excellent energy conservation properties

We propose a new algorithm for non-equilibrium molecular dynamics simulations of thermal gradients. The algorithm is an extension of the heat exchange algorithm developed by Hafskjold and co-workers [Mol. Phys. 80, 1389 (1993); Mol. Phys. 81, 251 (1994)], in which a certain amount of heat is added to one region and removed from another by rescaling velocities appropriately. Since the amount of added and removed heat is the same and the dynamics between velocity rescaling steps is Hamiltonian, the heat exchange algorithm is expected to conserve the energy. However, it has been reported previously that the original version of the heat exchange algorithm exhibits a pronounced drift in the total energy, the exact cause of which remained hitherto unclear. Here, we show that the energy drift is due to the truncation error arising from the operator splitting and suggest an additional coordinate integration step as a remedy. The new algorithm retains all the advantages of the original one whilst exhibiting excellent energy conservation as illustrated for a Lennard-Jones liquid and SPC/E water.

preprint2014arXiv

Living clusters and crystals from low density suspensions of active colloids

Recent studies aimed at investigating artificial analogs of bacterial colonies have shown that low-density suspensions of self-propelled particles confined in two dimensions can assemble into finite aggregates that merge and split, but have a typical size that remains constant (living clusters). In this Letter we address the problem of the formation of living clusters and crystals of active particles in three dimensions. We study two systems: self-propelled particles interacting via a generic attractive potential and colloids that can move towards each other as a result of active agents (e.g. by molecular motors). In both cases fluid-like `living' clusters form. We explain this general feature in terms of the balance between active forces and regression to thermodynamic equilibrium. This balance can be quantified in terms of a dimensionless number that allows us to collapse the observed clustering behaviour onto a universal curve. We also discuss how active motion affects the kinetics of crystal formation.

preprint2012arXiv

Controlling the temperature sensitivity of DNA-mediated colloidal interactions through competing linkages

We propose a new strategy to improve the self-assembly properties of DNA-functionalised colloids. The problem that we address is that DNA-functionalised colloids typically crystallize in a narrow temperature window, if at all. The underlying reason is the extreme sensitivity of DNA-mediated interactions to temperature or other physical control parameters. We propose to widen the window for colloidal crystallization by exploiting the competition between DNA linkages with different nucleotide sequences, which results in a temperature-dependent switching of the dominant bond type. Following such a strategy, we can decrease the temperature dependence of DNA-mediated self assembly to make systems that can crystallize in a wider temperature window than is possible with existing systems of DNA functionalised colloids. We report Monte Carlo simulations that show that the proposed strategy can indeed work in practice for real systems and specific, designable DNA sequences. Depending on the length ratio of the different DNA constructs, we find that the bond switching is either energetically driven (equal length or `symmetric' DNA) or controlled by a combinatorial entropy gain (`asymmetric' DNA), which results from the large number of possible binding partners for each DNA strand. We provide specific suggestions for the DNA sequences with which these effects can be achieved experimentally.

preprint2011arXiv

Irreducible finite-size effects in surface free energies from crystal-nucleation data

In this Letter we report a simulation study in which we compare the solid-liquid interfacial free energy of NaCl at coexistence, with the value that follows from the height of the homogeneous nucleation barrier. We find that the two estimates differ by more than 100%. Similar, although smaller discrepancies are found for crystals of hard-sphere colloids and of Lennard-Jones ("argon") particles. We consider a variety of possible causes for this discrepancy and are forced to conclude that it is due to a finite-size effect that cannot be corrected for by any simple thermodynamic procedure. Importantly, we find that the surface free energies that follow from real nucleation experiments should be subject to large finite size effects. Taking this in to account, we obtain quantitative agreement between the simulation data and the surface free energy of NaCl that follows from nucleation experiments. Our finding suggests that most published solid-liquid surface free energies derived from nucleation experiments will have to be revised.

preprint2009arXiv

Local structure of liquid carbon controls diamond nucleation

Diamonds melt at temperatures above 4000 K. There are no measurements of the steady-state rate of the reverse process: diamond nucleation from the melt, because experiments are difficult at these extreme temperatures and pressures. Using numerical simulations, we estimate the diamond nucleation rate and find that it increases by many orders of magnitude when the pressure is increased at constant supersaturation. The reason is that an increase in pressure changes the local coordination of carbon atoms from three-fold to four-fold. It turns out to be much easier to nucleate diamond in a four-fold coordinated liquid than in a liquid with three-fold coordination, because in the latter case the free-energy cost to create a diamond-liquid interface is higher. We speculate that this mechanism for nucleation control is relevant for crystallization in many network-forming liquids. On the basis of our calculations, we conclude that homogeneous diamond nucleation is likely in carbon-rich stars and unlikely in gaseous planets.