Researcher profile

Daan Frenkel

Daan Frenkel contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
8works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

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

8 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.

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.

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.

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.

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.

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.

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.

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).