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

Chantal Valeriani

Chantal Valeriani 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-mechcond-mat.mtrl-sciphysics.chem-phBiological Physicsphysics.comp-ph

Catalog footprint

What is connected

24works
6topics
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

24 published item(s)

preprint2022arXiv

Intrinsic structure perspective for MIPS interfaces in two dimensional systems of Active Brownian Particles

Suspensions of Active Brownian Particles (ABP) undergo motility induced phase separation (MIPS) over a wide range of mean density and activity strength [1], even in the absence of an explicit attraction. Negative values of the mechanical surface tension have been reported, from the total forces across the interface, while the stable fluctuations of the interfacial line would be interpreted as a positive capillary surface tension [2], while in equilibrium liquid surfaces these two magnitudes are equal. We present here the analysis of 2D-ABP interfaces in terms of the intrinsic density and force profiles, calculated with the particle distance to the instantaneous interfacial line. Our results provide a new insight in the origin of the MIPS from the local rectification of the random active force on the particles near the interface. As it had been pointed, that effect acts as an external potential [3] that produces a pressure gradient across the interface, so that the mechanical surface tension of the MIPS cannot be described as that of equilibrium coexisting phases; but our analysis shows that most of that effect comes from the tightly caged particles at the dense (inner) side of the MIPS interface, rather than from the free moving particles at the outer side that collide with the dense cluster. Moreover, a clear correlation appears between the decay of the hexatic order parameter at the dense slab and the end of the MIPS as the strength of the active force is lowered. We test that with the strong active forces required for MIPS, the interfacial structure and properties are very similar for ABP with purely repulsive (WCA-LJ model truncated at its minimum) and when the interaction includes a range of the LJ attractive force.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Rheology of Pseudomonas fluorescens biofilms: from experiments to DPD mesoscopic modelling

The presence of bacterial biofilms in clinical and industrial settings is a major issue worldwide. A biofilm is a viscoelastic matrix, composed of bacteria producing a network of Extracellular Polymeric Substances (EPS) to which bacteria crosslink. Modelling of complex biofilms is relevant to provide accurate descriptions and predictions that include parameters as hydrodynamics, dynamics of the bacterial population and solute mass transport. However, up-to-date numerical modelling, even at a coarse-grained level, is not satisfactorily. In this work, we present a numerical coarse-grain model of a bacterial biofilm, consisting of bacteria immersed in an aqueous matrix of a polymer network, that allows to study rheological properties of a biofilm. We study its viscoelastic modulus, varying topology and composition (such as the number of crosslinks between EPS polymers, the number of bacteria and the amount of solvent), and compare the numerical results with experimental rheological data of Pseudomonas fluorescens biofilms grown under static and shaking conditions, as previously described by Jara et al, Frontiers in Microbiology (2021).

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

In-silico modeling of early-stage biofilm formation

Several bacteria and bacteria strands form biofilms in different environmental conditions, e.g. pH, temperature, nutrients, etc. Biofilm growth, therefore, is an extremely robust process. Because of this, while biofilm growth is a complex process affected by several variables, insights into biofilm formation could be obtained studying simple schematic models. In this manuscript, we describe a hybrid molecular dynamics and Monte Carlo model for the simulation of the early stage formation of a biofilm, to explicitly demonstrate that it is possible to account for most of the processes expected to be relevant. The simulations account for the growth and reproduction of the bacteria, for their interaction and motility, for the synthesis of extracellular polymeric substances and Psl trails. We describe the effect of these processes on the early stage formation of biofilms, in two dimensions, and also discuss preliminary three-dimensional results.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Self-adaptation of Pseudomonas fluorescens biofilms to hydrodynamic stress

In some conditions, bacteria self-organise into biofilms, supracellular structures made of a self-produced embedding matrix, mainly composed on polysaccharides, DNA, proteins and lipids. It is known that bacteria change their colony/matrix ratio in the presence of external stimuli such as hydrodynamic stress. However, little is still known about the molecular mechanisms driving this self-adaptation. In this work, we monitor structural features of Pseudomonas fluorescens biofilms grown with and without hydrodynamic stress. Our measurements show that the hydrodynamic stress concomitantly increases the cell density population and the matrix production. At short growth timescales, the matrix mediates a weak cell-cell attractive interaction due to the depletion forces originated by the polymer constituents. Using a population dynamics model, we conclude that hydrodynamic stress causes a faster diffusion of nutrients and a higher incorporation of planktonic bacteria to the already formed microcolonies. This results in the formation of more mechanically stable biofilms due to an increase of the number of crosslinks, as shown by computer simulations. The mechanical stability also lies on a change in the chemical compositions of the matrix, which becomes enriched in carbohydrates, known to display adhering properties. Overall, we demonstrate that bacteria are capable of self-adapting to hostile hydrodynamic stress by tailoring the biofilm chemical composition, thus affecting both the mesoscale structure of the matrix and its viscoelastic properties that ultimately regulate the bacteria-polymer interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Effect of dissolved salt on the anomalies of water at negative pressure

Adding salt to water at ambient pressure affects its thermodynamic properties. At low salt concentration, anomalies such as the density maximum are shifted to lower temperature, while at large enough salt concentration they cannot be observed any more. Here we investigate the effect of salt on an anomaly recently observed in pure water at negative pressure: the existence of a sound velocity minimum along isochores. We compare experiments and simulations for an aqueous solution of sodium chloride with molality around $1.2\,\mathrm{mol\,kg^{-1}}$, reaching pressures beyond $-100\,\mathrm{MPa}$. We also discuss the origin of the minima in the sound velocity and emphasize the importance of the relative position of the temperatures of sound velocity and density anomalies.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Seeding Approach to Bubble Nucleation in Superheated Lennard Jones Fluids

We investigate vapor homogeneous nucleation in a superheated Lennard Jones liquid with computer simulations. Special simulation techniques are required to address this study since the nucleation of a critical vapor bubble --one that has equal chance to grow or shrink-- in a moderately superheated liquid is a rare event. We use the Seeding method, that combines Classical Nucleation Theory with computer simulations of a liquid containing a vapor bubble to provide bubble nucleation rates in a wide temperature range. Seeding has been successfully applied to investigate the nucleation of crystals in supercooled fluids and here we apply it for the fist time --to the best of our knowledge-- to the liquid-to-vapor transition. We find that Seeding provides nucleation rates that are consistent with independent calculations not based on the assumptions of Classical Nucleation Theory. Different criteria to determine the radius of the critical bubble give different rate values. The accuracy of each criterion depends of the degree of superheating. Moreover, Seeding simulations show that the surface tension depends on pressure for a given temperature. Therefore, using Classical Nucleation Theory with the coexistence surface tension does not provide good estimates of the nucleation rate.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Seeding Approach to nucleation in the NVT ensemble: the case of bubble cavitation in overstretched Lennard Jones fluids

Simulations are widely used to study nucleation in first order phase transitions due to the fact that they have access to the relevant length and time scales. However, simulations face the problem that nucleation is an activated process. Therefore, rare event simulation techniques are needed to promote the formation of the critical nucleus. The Seeding method, where the simulations are started with the nucleus already formed, has proven quite useful in efficiently providing estimates of the nucleation rate for a wide range of orders of magnitude. So far, Seeding has been employed in the NPT ensemble, where the nucleus either grows or redissolves. Thus, several trajectories have to be run in order to find the thermodynamic conditions that make the seeded nucleus critical. Moreover, the nucleus lifetime is short and the statistics for obtaining its properties is consequently poor. To deal with these shortcomings we extend the Seeding method to the NVT ensemble. We focus on the problem of bubble nucleation in a mestastable Lennard Jones fluid. We show that, in the NVT ensemble, it is possible to equilibrate and stabilise critical bubbles for a long time. The nucleation rate inferred from NVT-Seeding is fully consistent with that coming from NPT-Seeding. The former is quite suitable to obtain the nucleation rate along isotherms, whereas the latter is preferable if the dependence of the rate with temperature at constant pressure is required. Care should be taken with finite size effects when using NVT-Seeding. Further work is required to extend NVT seeding to other sorts of phase transitions.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A comprehensive scenario of the thermodynamic anomalies of water using the TIP4P/2005 model

The striking behavior of water has deserved it to be referred to as an "anomalous" liquid. The water anomalies are greatly amplified in metastable (supercooled/stretched) regions. This makes difficult a complete experimental description since, beyond certain limits, the metastable phase necessarily transforms into the stable one. Theoretical interpretation of the water anomalies could then be based on simulation results of well validated water models. But the analysis of the simulations has not yet reached a consensus. In particular, one of the most popular theoretical scenarios -involving the existence of a liquid-liquid critical point (LLCP)- is disputed by several authors. In this work we propose to use a number of exact thermodynamic relations which can be tested in a region of the phase diagram outside the LLCP thus avoiding the problems associated to the coexistence region. The central property connected to other water anomalies is the locus of temperatures at which the density along isobars attain a maximum (TMD line) or a minimum (TmD), computed via simulations for a successful water model, TIP4P/2005. Next, we have evaluated the vapor-liquid spinodal in the region of large negative pressures. The shape of these curves and their connection to the extrema of response functions (isothermal compressibility and heat capacity at constant pressure) may help to elucidate the validity of the theoretical proposals. In this way we are able to present for the first time a comprehensive scenario of the thermodynamic water anomalies for TIP4P/2005 and their relation to the vapor-liquid spinodal. The overall picture shows a remarkable similarity with the corresponding one for the ST2 water model, for which the existence of a LLCP has been demonstrated in recent years. It also provides a hint as to where the long-sought for extrema in response functions might become accessible to experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

An Equation of State For The TIP4P/2005 Model of Water Including Negative Pressures

One of the most promising frameworks for understanding the anomalies of cold and especially supercooled water is that of two-structure thermodynamics, in which water is viewed as a non-ideal mixture of two interconvertible local structures. The non-ideality of this mixture may give rise, at very low temperatures, to a liquid-liquid phase transition (LLPT) and a liquid-liquid critical point (LLCP), at which thermodynamic response functions diverge. Various versions of the "two-structure equation of state" (TSEOS) based on this concept have shown remarkable agreement with both experimental data in real water and simulation results. However, recent experiments probing supercooled water at negative pressures reveal the inadequacy of extrapolations of equations of state developed for positive pressures, and have begun to shed additional light on the source of the anomalies of supercooled water. We have analyzed simulation results for the TIP4P/2005 model over a broad range of positive and negative pressures from ambient temperature to deep supercooling. We find that by explicitly incorporating a liquid-vapor spinodal into a two-structure equation of state, we are able to match the simulation data in TIP4P/2005 with striking accuracy. In particular, our equation of state reproduces the observed lines of minima and maxima in the density, isothermal compressibility, and isobaric heat capacity. Contrary to scenarios in which a retracing spinodal accounts for the thermodynamic anomalies of water, we find that the liquid-vapor spinodal in this model continues monotonically to lower pressures as temperature is decreased, influencing but not giving rise to the locus of density maxima and other thermodynamic anomalies. We explain the behavior of TIP4P/2005 in terms of two phenomena: the competition between two local structures and a monotonic liquid-vapor spinodal.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Molecular mechanism for cavitation in water under tension

Despite its relevance in biology and engineering, the molecular mechanism driving cavitation in water remains unknown. Using computer simulations, we investigate the structure and dynamics of vapor bubbles emerging from metastable water at negative pressures. We find that in the early stages of cavitation, bubbles are irregularly shaped and become more spherical as they grow. Nevertheless, the free energy of bubble formation can be perfectly reproduced in the framework of classical nucleation theory (CNT) if the curvature dependence of the surface tension is taken into account. Comparison of the observed bubble dynamics to the predictions of the macroscopic Rayleigh--Plesset (RP) equation, augmented with thermal fluctuations, demonstrates that the growth of nanoscale bubbles is governed by viscous forces. Combining the dynamical prefactor determined from the RP equation with the free energy of CNT yields an analytical expression for the cavitation rate that reproduces the simulation results very well over a wide range of pressures. Furthermore, our theoretical predictions are in excellent agreement with cavitation rates obtained from inclusion experiments. This suggests that homogeneous nucleation is observed in inclusions, whereas only heterogeneous nucleation on impurities or defects occurs in other experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Morphology of clusters of attractive dry and wet self-propelled spherical particle suspensions

In order to asses the effect of hydrodynamics in the assembly of active attractive spheres, we simulate a semi-dilute suspension of attractive self-propelled spherical particles in a quasi two dimensional geometry comparing the case with and without hydrodynamics interactions. To start with, independently on the presence of hydrodynamics, we observe that depending on the ratio between attraction and propulsion, particles either coarsen or aggregate forming finite-size clusters. Focusing on the clustering regime, we characterize two different clusters parameters, i.e. their morphology and orientational order, and compare the case when active particles behave either as pushers or pullers (always in the regime where inter-particles attractions competes with self-propulsion). Studying cluster phases for squirmers with respect to those obtained for active Brownian disks (indicated as ABP), we have shown that hydrodynamics alone can sustain a cluster phase of active swimmers (pullers), while ABP form cluster phases due to the competition between attraction and self propulsion. The structural properties of the cluster phases of squirmers and ABP are similar, although squirmers show sensitivity to active stresses. Active Brownian disks resemble weakly pusher squirmer suspensions in terms of cluster size distribution, structure of the radius of gyration on cluster size and degree of cluster polarity.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Microphase separation in two dimensional suspensions of self-propelled spheres and dumbbells

We use numerical simulations to study the phase behavior of self-propelled spherical and dumbbellar particles interacting via micro-phase separation inducing potentials. Our results indicate that under the appropriate conditions, it is possible to drive the formation of two new active states; a spinning cluster crystal, i.e. an ordered mesoscopic phase having finite size spinning crystallites as lattice sites, and a fluid of living clusters, i.e. a two dimensional fluid where each "particle" is a finite size living cluster. We discuss the dynamics of these phases and suggest ways of extending their stability under a wide range of active forces.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Activity induced collapse and re-expansion of rigid polymers

We study the elastic properties of a rigid filament in a bath of self-propelled particles. We find that while fully flexible filaments swell monotonically upon increasing the strength of the propelling force, rigid filaments soften for moderate activities, collapse into metastable hairpins for intermediate strengths, and eventually re-expand when the strength of the activity of the surrounding fluid is large. This collapse and re-expansion of the filament with the bath activity is reminiscent of the behavior observed in polyelectrolytes in the presence of different concentrations of multivalent salt.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Avalanches mediate crystallization in a hard-sphere glass

By molecular-dynamics simulations, we have studied the devitrification (or crystallization) of aged hard-sphere glasses. First, we find that the dynamics of the particles are intermittent: Quiescent periods, when the particles simply "rattle" in their nearest-neighbor cages, are interrupted by abrupt "avalanches," where a subset of particles undergo large rearrangements. Second, we find that crystallization is associated with these avalanches but that the connection is not straightforward. The amount of crystal in the system increases during an avalanche, but most of the particles that become crystalline are different from those involved in the avalanche. Third, the occurrence of the avalanches is a largely stochastic process. Randomizing the velocities of the particles at any time during the simulation leads to a different subsequent series of avalanches. The spatial distribution of avalanching particles appears random, although correlations are found among avalanche initiation events. By contrast, we find that crystallization tends to take place in regions that already show incipient local order.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Anomalies in bulk supercooled water at negative pressure

Water is the most familiar liquid, and arguably the most complex. Anomalies of supercooled water have been measured during decades, and competing interpretations proposed. Yet, a decisive experiment remains elusive, because of unavoidable crystallization into ice. We investigate the state of water that is both supercooled and under mechanical tension, or negative pressure. Liquids under negative pressure can be found in plants or fluid inclusions in minerals. Using such water inclusions in quartz, we report the first measurements on doubly metastable water down to $-15^\circ\mathrm{C}$ and around $-100\,\mathrm{MPa}$. We observe sound velocity anomalies that can be reproduced quantitatively with molecular dynamics simulations. These results rule out one proposed scenario for water anomalies, and put further constraints on the remaining ones.

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

From compact to fractal crystalline clusters in concentrated systems of monodisperse hard spheres

We address the crystallization of monodisperse hard spheres in terms of the properties of finite- size crystalline clusters. By means of large scale event-driven Molecular Dynamics simulations, we study systems at different packing fractions ϕ ranging from weakly supersaturated state points to glassy ones, covering different nucleation regimes. We find that such regimes also result in different properties of the crystalline clusters: compact clusters are formed in the classical-nucleation-theory regime (ϕ \leq 0.54), while a crossover to fractal, ramified clusters is encountered upon increasing packing fraction (ϕ \geq 0.56), where nucleation is more spinodal-like. We draw an analogy between macroscopic crystallization of our clusters and percolation of attractive systems to provide ideas on how the packing fraction influences the final structure of the macroscopic crystals. In our previous work (Phys. Rev. Lett., 106, 215701, 2011), we have demonstrated how crystallization from a glass (at ϕ > 0.58) happens via a gradual (many-step) mechanism: in this paper we show how the mechanism of gradual growth seems to hold also in super-saturated systems just above freezing showing that static properties of clusters are not much affected by dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Colloids in a bacterial bath: simulations and experiments

We present a joint experimental and computational study of the effect of bacterial motion on micron-scale colloids contained in a two-dimensional suspension of Bacillus subtilis. With respect to previous work using E. coli, here we introduce a novel experimental set-up that allows us to realise a two-dimensional bacterial suspension insensitive to either evaporation or fluid flow. By analysing the mean square displacements of both bacteria and colloids, we confirm the existence of a crossover from super-diffusive behaviour at short time scales to normal diffusion at longer times. We also study the same two-dimensional system by means of numerical simulations, using a suspension of self-propelled dumbbells or the Vicsek model, which has been previously used to study the dynamics of active particles. Our numerical results obtained with both models are in broad agreement with the experimental trends, but only the dumbbell simulations can match the experimental data quantitatively. The level of agreement we find suggest that steric interactions due to collisions are important players in determining collective motion of the bacterial bath, and should complement hydrodynamic interactions in experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Crystallization Mechanism of Hard Sphere Glasses

In supercooled liquids, vitrification generally suppresses crystallization. Yet some glasses can still crystallize despite the arrest of diffusive motion. This ill-understood process may limit the stability of glasses, but its microscopic mechanism is not yet known. Here we present extensive computer simulations addressing the crystallization of monodisperse hard-sphere glasses at constant volume (as in a colloid experiment). Multiple crystalline patches appear without particles having to diffuse more than one diameter. As these patches grow, the mobility in neighbouring areas is enhanced, creating dynamic heterogeneity with positive feedback. The future crystallization pattern cannot be predicted from the coordinates alone: crystallization proceeds by a sequence of stochastic micro-nucleation events, correlated in space by emergent dynamic heterogeneity.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Ion association in low-polarity solvents: comparisons between theory, simulation, and experiment

The association of ions in electrolyte solutions at very low concentration and low temperature is studied using computer simulations and quasi-chemical ion-pairing theory. The specific case of the restricted primitive model (charged hard spheres) is considered. Specialised simulation techniques are employed that lead to efficient sampling of the arrangements and distributions of clusters and free ions, even at conditions corresponding to nanomolar solutions of simple salts in solvents with dielectric constants in the range 5-10, as used in recent experimental work on charged-colloid sus- pensions. A direct comparison is effected between theory and simulation using a variety of clustering criteria and theoretical approximations. It is shown that conventional distance-based cluster criteria can give erroneous results. A reliable set of theoretical and simulation estimators for the degree of association is proposed. The ion-pairing theory is then compared to experimental results for salt solutions in low-polarity solvents. The agreement is excellent, and on this basis some calculations are made for the screening lengths which will figure in the treatment of colloid-colloid interactions in such solutions. The accord with available experimental results is complete.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Nucleation in a sheared Ising model: effects of external field

Simulations using the Forward Flux Sampling method have shown a nonmonotonic de- pendence of the homogeneous nucleation rate on the shear rate for a sheared two dimensional Ising model [R. J. Allen et al, arXiv cond-mat/0805.3029]. For quasi-equilibrium systems (i.e. in the absence of shear), Classical Nucleation Theory (CNT) predicts the dependence of the critical cluster size and the nucleation rate on the external magnetic field. We investigate the behaviour of the sheared Ising model as a function of the external field. At low exter- nal field strength, the same nonmonotonic behaviour holds and the peak in the nucleation rate is remarkably insensitive to the field strength. This suggests that the same external field-dependence holds for the enhancement of nucleation by shear at low shear rates and the suppression of shear at high shear rates. At high field strength, the nucleation behaviour is qualitatively different. We also analyse the size and shape of the largest cluster in the transition state configurations, as a function of the external field. In the sheared system, the transition state cluster becomes larger and more elongated as the field strength decreases. We compare our results for the sheared system to the predictions of the CNT for the quasi- equilibrium case, and find that the CNT cannot easily be used to describe nucleation in the system under shear.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Computer simulations of the restricted primitive model at very low temperature and density

The problem of successfully simulating ionic fluids at low temperature and low density states is well known in the simulation literature: using conventional methods, the system is not able to equilibrate rapidly due to the presence of strongly associated cation-anion pairs. In this manuscript we present a numerical method for speeding up computer simulations of the restricted primitive model (RPM) at low temperatures (around the critical temperature) and at very low densities (down to $10^{-10}σ^{-3}$, where $σ$ is the ion diameter). Experimentally, this regime corresponds to typical concentrations of electrolytes in nonaqueous solvents. As far as we are aware, this is the first time that the RPM has been equilibrated at such extremely low concentrations. More generally, this method could be used to equilibrate other systems that form aggregates at low concentrations.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Forward Flux Sampling for rare event simulations

Rare events are ubiquitous in many different fields, yet they are notoriously difficult to simulate because few, if any, events are observed in a conventiona l simulation run. Over the past several decades, specialised simulation methods have been developed to overcome this problem. We review one recently-developed class of such methods, known as Forward Flux Sampling. Forward Flux Sampling uses a series of interfaces between the initial and final states to calculate rate constants and generate transition paths, for rare events in equilibrium or nonequilibrium systems with stochastic dynamics. This review draws together a number of recent advances, summarizes several applications of the method and highlights challenges that remain to be overcome.

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