Showing 1–29 of 29 results for author: Bolhuis, P G

Molecular free energies, rates, and mechanisms from data-efficient path sampling simulations

Authors: Gianmarco Lazzeri, Hendrik Jung, Peter G. Bolhuis, Roberto Covino

Abstract: Molecular dynamics is a powerful tool for studying the thermodynamics and kinetics of complex molecular events. However, these simulations can rarely sample the required time scales in practice. Transition path sampling overcomes this limitation by collecting unbiased trajectories capturing the relevant events. Moreover, the integration of machine learning can boost the sampling while simultaneous… ▽ More Molecular dynamics is a powerful tool for studying the thermodynamics and kinetics of complex molecular events. However, these simulations can rarely sample the required time scales in practice. Transition path sampling overcomes this limitation by collecting unbiased trajectories capturing the relevant events. Moreover, the integration of machine learning can boost the sampling while simultaneously learning a quantitative representation of the mechanism. Still, the resulting trajectories are by construction non-Boltzmann-distributed, preventing the calculation of free energies and rates. We developed an algorithm to approximate the equilibrium path ensemble from machine learning-guided path sampling data. At the same time, our algorithm provides efficient sampling, the mechanism, free energy, and rates of rare molecular events at a very moderate computational cost. We tested the method on the folding of the mini-protein chignolin. Our algorithm is straightforward and data-efficient, opening the door to applications on many challenging molecular systems. △ Less

Submitted 28 July, 2023; v1 submitted 20 July, 2023; originally announced July 2023.

Comments: Corrected some typos and updated the bibliography

Force field optimization by imposing kinetic constraints with path reweighting

Authors: P. G. Bolhuis, Z. F. Brotzakis, B. G. Keller

Abstract: Empirical force fields employed in molecular dynamics simulations of complex systems can be optimised to reproduce experimentally determined structural and thermodynamic properties. In contrast, experimental knowledge about the rates of interconversion between metastable states in such systems, is hardly ever incorporated in a force field, due to a lack of an efficient approach. Here, we introduce… ▽ More Empirical force fields employed in molecular dynamics simulations of complex systems can be optimised to reproduce experimentally determined structural and thermodynamic properties. In contrast, experimental knowledge about the rates of interconversion between metastable states in such systems, is hardly ever incorporated in a force field, due to a lack of an efficient approach. Here, we introduce such a framework, based on the relationship between dynamical observables such as rate constants, and the underlying force field parameters, using the statistical mechanics of trajectories. Given a prior ensemble of molecular trajectories produced with imperfect force field parameters, the approach allows the optimal adaption of these parameters, such that the imposed constraint of equal predicted and experimental rate constant is obeyed. To do so, the method combines the continuum path ensemble Maximum Caliber approach with path reweighting methods for stochastic dynamics. When multiple solutions are found, the method selects automatically the combination that corresponds to the smallest perturbation of the entire path ensemble, as required by the Maximum Entropy principle. To show the validity of the approach we illustrate the method on simple test systems undergoing rare event dynamics. Next to simple 2D potentials we explore particle models representing molecular isomerisation reactions as well as protein-ligand unbinding. Besides optimal interaction parameters the methodology gives physical insight into what parts of the model are most sensitive to the kinetics. We discuss the generality and broad implications of the methodology. △ Less

Submitted 10 July, 2022; originally announced July 2022.

Comments: 18 pages, 12 figures

Extended Wertheim theory predicts the anomalous chain length distributions of divalent patchy particles under extreme confinement

Authors: H. J. Jonas, P. Schall, P. G. Bolhuis

Abstract: Colloidal patchy particles with divalent attractive interaction can self-assemble into linear polymer chains. Their equilibrium properties in 2D and 3D are well described by Wertheim's thermodynamic perturbation theory which predicts a well-defined exponentially decaying equilibrium chain length distribution. In experimental realizations, due to gravity, particles sediment to the bottom of the sus… ▽ More Colloidal patchy particles with divalent attractive interaction can self-assemble into linear polymer chains. Their equilibrium properties in 2D and 3D are well described by Wertheim's thermodynamic perturbation theory which predicts a well-defined exponentially decaying equilibrium chain length distribution. In experimental realizations, due to gravity, particles sediment to the bottom of the suspension forming a monolayer of particles with a gravitational height smaller than the particle diameter. In accordance with experiments, an anomalously high monomer concentration is observed in simulations which is not well understood. To account for this observation, we interpret the polymerization as taking place in a highly confined quasi-2D plane and extend the Wertheim thermodynamic perturbation theory by defining addition reactions constants as functions of the chain length. We derive the theory, test it on simple square well potentials, and apply it to the experimental case of synthetic colloidal patchy particles immersed in a binary liquid mixture that are described by an accurate effective critical Casimir patchy particle potential. The important interaction parameters entering the theory are explicitly computed using the integral method in combination with Monte Carlo sampling. Without any adjustable parameter, the predictions of the chain length distribution are in excellent agreement with explicit simulations of self-assembling particles. We discuss generality of the approach, and its application range. △ Less

Submitted 22 June, 2022; originally announced June 2022.

Comments: The following article has been submitted to The Journal of Chemical Physics

A method of incorporating rate constants as kinetic constraints in molecular dynamics simulations

Authors: Z. Faidon Brotzakis, Michele Vendruscolo, Peter. G. Bolhuis

Abstract: From the point of view of statistical mechanics, a full characterisation of a molecular system requires the experimental determination of its possible states, their populations and the respective interconversion rates. Well-established methods can incorporate in molecular dynamics simulations experimental information about states using structural restraints, and about populations using thermodynam… ▽ More From the point of view of statistical mechanics, a full characterisation of a molecular system requires the experimental determination of its possible states, their populations and the respective interconversion rates. Well-established methods can incorporate in molecular dynamics simulations experimental information about states using structural restraints, and about populations using thermodynamic restraints. However, it is still unclear how to include experimental knowledge of interconversion rates. Here we introduce a method of imposing known rate constants as constraints in molecular dynamics simulations, which is based on a combination of the maximum entropy and maximum caliber principles. Starting from an existing ensemble of trajectories, obtained from either molecular dynamics or enhanced trajectory sampling, this method provides a minimally perturbed path distribution consistent with the kinetic constraints, as well as a modified free energy and committor landscape. We illustrate the application of the method to simple toy systems, as well as to all atom molecular simulations of peptide association and folding. We find that by combining experimental rate coefficient data and molecular dynamics simulations we are able to determine new transition states, reaction mechanisms and free energies. For instance, in the case of chignolin protein folding we find that imposing a slower folding rate shifts the transition state to more native like conformations, while it increases the stability of the unfolded region. △ Less

Submitted 1 June, 2020; originally announced June 2020.

Approximating Free Energy and Committor Landscapes in Standard Transition Path Sampling using Virtual Interface Exchange

Authors: Z. Faidon Brotzakis, Peter G. Bolhuis

Abstract: Transition path sampling (TPS) is a powerful technique for investigating rare transitions, especially when the mechanism is unknown and one does not have access to the reaction coordinate. Straightforward application of TPS does not directly provide the free energy landscape nor the kinetics, which motivated the development of path sampling extensions, such as transition interface sampling (TIS),… ▽ More Transition path sampling (TPS) is a powerful technique for investigating rare transitions, especially when the mechanism is unknown and one does not have access to the reaction coordinate. Straightforward application of TPS does not directly provide the free energy landscape nor the kinetics, which motivated the development of path sampling extensions, such as transition interface sampling (TIS), and the reweighted paths ensemble (RPE), that are able to simultaneously access both kinetics and thermodynamics. However, performing TIS is more involved than TPS, and still requires (some) insight in the reaction to define interfaces. While packages that can efficiently compute path ensembles for TIS are now available, it would be useful to directly compute the free energy from a single TPS simulation. To achieve this, we developed an approximate method, denoted Virtual Interface Exchange, that makes use of the rejected pathways in a form of waste recycling. The method yields an approximate reweighted path ensemble that allows an immediate view of the free energy landscape from a single TPS, as well as enables a full committor analysis. △ Less

Submitted 9 July, 2019; originally announced July 2019.

Comments: 10 figures

Ratchet-induced variations in bulk states of an active ideal gas

Authors: Jeroen Rodenburg, Siddharth Paliwal, Marjolein de Jager, Peter G. Bolhuis, Marjolein Dijkstra, René van Roij

Abstract: We study the distribution of active, noninteracting particles over two bulk states separated by a ratchet potential. By solving the steady-state Smoluchowski equations in a flux-free setting, we show that the ratchet potential affects the distribution of particles over the bulks, and thus exerts an influence of infinitely long range. As we show, crucial for having such a long-range influence is an… ▽ More We study the distribution of active, noninteracting particles over two bulk states separated by a ratchet potential. By solving the steady-state Smoluchowski equations in a flux-free setting, we show that the ratchet potential affects the distribution of particles over the bulks, and thus exerts an influence of infinitely long range. As we show, crucial for having such a long-range influence is an external potential that is nonlinear. We characterize how the difference in bulk densities depends on activity and on the ratchet potential, and we identify power law dependencies on system parameters in several limiting cases. While weakly active systems are often understood in terms of an effective temperature, we present an analytical solution that explicitly shows that this is not possible in the current setting. Instead, we rationalize our results by a simple transition state model, that presumes particles to cross the potential barrier by Arrhenius rates modified for activity. While this model does not quantitatively describe the difference in bulk densities for feasible parameter values, it does reproduce - in its regime of applicability - the complete power law behavior correctly. △ Less

Submitted 12 November, 2018; v1 submitted 24 July, 2018; originally announced July 2018.

Comments: 11 pages, 6 figures

Journal ref: J. Chem. Phys. 149, 174910 (2018)

Multiscale simulations of anisotropic particles combining Brownian Dynamics and Green's Function Reaction Dynamics

Authors: Adithya Vijaykumar, Thomas E. Ouldridge, Pieter Rein ten Wolde, Peter G. Bolhuis

Abstract: The modeling of complex reaction-diffusion processes in, for instance, cellular biochemical networks or self-assembling soft matter can be tremendously sped up by employing a multiscale algorithm which combines the mesoscopic Green's Function Reaction Dynamics (GFRD) method with explicit stochastic Brownian, Langevin, or deterministic Molecular Dynamics to treat reactants at the microscopic scale… ▽ More The modeling of complex reaction-diffusion processes in, for instance, cellular biochemical networks or self-assembling soft matter can be tremendously sped up by employing a multiscale algorithm which combines the mesoscopic Green's Function Reaction Dynamics (GFRD) method with explicit stochastic Brownian, Langevin, or deterministic Molecular Dynamics to treat reactants at the microscopic scale [A. Vijaykumar, P.G. Bolhuis and P.R. ten Wolde, J. Chem. Phys. {\bf 43}, 21: 214102 (2015)]. Here we extend this multiscale BD-GFRD approach to include the orientational dynamics that is crucial to describe the anisotropic interactions often prevalent in biomolecular systems. We illustrate the novel algorithm using a simple patchy particle model. After validation of the algorithm we discuss its performance. The rotational BD-GFRD multiscale method will open up the possibility for large scale simulations of e.g. protein signalling networks. △ Less

Submitted 28 November, 2016; originally announced November 2016.

Comments: 14 pages, 11 figures

Self-assembly of Microcapsules via Colloidal Bond Hybridization and Anisotropy

Authors: Chris H. J. Evers, Jurriaan A. Luiken, Peter G. Bolhuis, Willem K. Kegel

Abstract: Particles with directional interactions are promising building blocks for new functional materials and may serve as models for biological structures. Mutually attractive nanoparticles that are deformable due to flexible surface groups, for example, may spontaneously order themselves into strings, sheets and large vesicles. Furthermore, anisotropic colloids with attractive patches can self-assemble… ▽ More Particles with directional interactions are promising building blocks for new functional materials and may serve as models for biological structures. Mutually attractive nanoparticles that are deformable due to flexible surface groups, for example, may spontaneously order themselves into strings, sheets and large vesicles. Furthermore, anisotropic colloids with attractive patches can self-assemble into open lattices and colloidal equivalents of molecules and micelles. However, model systems that combine mutual attraction, anisotropy, and deformability have---to the best of our knowledge---not been realized. Here, we synthesize colloidal particles that combine these three characteristics and obtain self-assembled microcapsules. We propose that mutual attraction and deformability induce directional interactions via colloidal bond hybridization. Our particles contain both mutually attractive and repulsive surface groups that are flexible. Analogous to the simplest chemical bond, where two isotropic orbitals hybridize into the molecular orbital of H2, these flexible groups redistribute upon binding. Via colloidal bond hybridization, isotropic spheres self-assemble into planar monolayers, while anisotropic snowman-like particles self-assemble into hollow monolayer microcapsules. A modest change of the building blocks thus results in a significant leap in the complexity of the self-assembled structures. In other words, these relatively simple building blocks self-assemble into dramatically more complex structures than similar particles that are isotropic or non-deformable. △ Less

Submitted 18 July, 2016; v1 submitted 2 March, 2015; originally announced March 2015.

Journal ref: Nature 534 364-368 (2016)

The competition between surface adsorption and folding of fibril-forming polypeptides

Authors: Ran Ni, J. Mieke Kleijn, Sanne Abeln, Martien A. Cohen Stuart, Peter G. Bolhuis

Abstract: Self-assembly of polypeptides into fibrillar structures can be initiated by planar surfaces that interact favorably with certain residues. Using a coarse grained model, we systematically studied the folding and adsorption behavior of a $β$-roll forming polypeptide. We find that there are two different folding pathways depending on the temperature: (i) at low temperature, the polypeptide folds in s… ▽ More Self-assembly of polypeptides into fibrillar structures can be initiated by planar surfaces that interact favorably with certain residues. Using a coarse grained model, we systematically studied the folding and adsorption behavior of a $β$-roll forming polypeptide. We find that there are two different folding pathways depending on the temperature: (i) at low temperature, the polypeptide folds in solution into a $β$-roll before adsorbing onto the attractive surface, (ii) at higher temperature, the polypeptide first adsorbs in a disordered state, and folds while on the surface. The folding temperature increases with increasing attraction, as the folded $β$-roll is stabilized by the surface. Surprisingly, further increasing the attraction lowers the folding temperature again, as strong attraction also stabilizes the adsorbed disordered state, which competes with folding of the polypeptide. Our results suggest that to enhance the folding, one should use a weakly attractive surface. They also explain the recent experimental observation of the non-monotonic effect of charge on the fibril formation on an oppositely charged surface [C. Charbonneau, et al., ACS Nano, 8, 2328 (2014)]. △ Less

Submitted 6 February, 2015; originally announced February 2015.

Comments: Accepted in Phys. Rev. E

Journal ref: Phys. Rev. E 91, 022711 (2015)

Tunable long range forces mediated by self-propelled colloidal hard spheres

Authors: Ran Ni, Martien A. Cohen Stuart, Peter G. Bolhuis

Abstract: Using Brownian dynamics simulations, we systematically study the effective interaction between two parallel hard walls in a 2D suspension of self-propelled (active) colloidal hard spheres, and we find that the effective force between two hard walls can be tuned from a long range repulsion into a long range attraction by changing the density of active particles. At relatively high densities, the ac… ▽ More Using Brownian dynamics simulations, we systematically study the effective interaction between two parallel hard walls in a 2D suspension of self-propelled (active) colloidal hard spheres, and we find that the effective force between two hard walls can be tuned from a long range repulsion into a long range attraction by changing the density of active particles. At relatively high densities, the active hard spheres can form a dynamic crystalline bridge, which induces a strong oscillating long range dynamic wetting repulsion between the walls. With decreasing density, the dynamic bridge gradually breaks, and an intriguing long range dynamic depletion attraction arises. A similar effect occurs in a quasi-2D suspension of self-propelled colloidal hard spheres by changing the height of the confinement. Our results open up new possibilities to manipulate the motion and assembly of microscopic objects by using active matter. △ Less

Submitted 11 December, 2014; v1 submitted 6 March, 2014; originally announced March 2014.

Comments: Accepted in Physical Review Letters

Journal ref: Phys. Rev. Lett. 114, 018302 (2015)

Crystallizing hard-sphere glasses by doping with active particles

Authors: Ran Ni, Martien A. Cohen Stuart, Marjolein Dijkstra, Peter G. Bolhuis

Abstract: Crystallization and vitrification are two different routes to form a solid. Normally these two processes suppress each other, with the glass transition preventing crystallization at high density (or low temperature). This is even true for systems of colloidal hard spheres, which are commonly used as building blocks for novel functional materials with potential applications, e.g. photonic crystals.… ▽ More Crystallization and vitrification are two different routes to form a solid. Normally these two processes suppress each other, with the glass transition preventing crystallization at high density (or low temperature). This is even true for systems of colloidal hard spheres, which are commonly used as building blocks for novel functional materials with potential applications, e.g. photonic crystals. By performing Brownian dynamics simulations of glassy systems consisting of mixtures of active and passive hard spheres, we show that the crystallization of such hard-sphere glasses can be dramatically promoted by doping the system with small amounts of active particles. Surprisingly, even hard-sphere glasses of packing fraction up to $φ= 0.635$ crystallize, which is around $0.5\%$ below the random close packing at $φ\simeq 0.64$. Our results suggest a novel way of fabricating crystalline materials from (colloidal) glasses. This is particularly important for materials that get easily kinetically trapped in glassy states, and crystal nucleation hardly occurs. △ Less

Submitted 25 August, 2014; v1 submitted 28 October, 2013; originally announced October 2013.

Journal ref: Soft Matter, 10 (35), 6609 (2014)

Interplay between folding and assembly of fibril-forming polypeptides

Authors: Ran Ni, Sanne Abeln, Marieke Schor, Martien A. Cohen Stuart, Peter G. Bolhuis

Abstract: Polypeptides can self-assemble into hierarchically organized fibrils consisting of a stack of individually folded polypeptides driven together by hydrophobic interaction. Using a coarse grained model, we systematically studied this self-assembly as a function of temperature and hydrophobicity of the residues on the outside of the building block. We find the self-assembly can occur via two differen… ▽ More Polypeptides can self-assemble into hierarchically organized fibrils consisting of a stack of individually folded polypeptides driven together by hydrophobic interaction. Using a coarse grained model, we systematically studied this self-assembly as a function of temperature and hydrophobicity of the residues on the outside of the building block. We find the self-assembly can occur via two different pathways - a random aggregation-folding route, and a templated-folding process - thus indicating a strong coupling between folding and assembly. The simulation results can explain experimental evidence that assembly through stacking of folded building blocks is rarely observed, at the experimental concentrations. The model thus provides a generic picture of hierarchical fibril formation. △ Less

Submitted 11 July, 2013; originally announced July 2013.

Comments: Accepted in Physical Review Letters

Journal ref: Phys. Rev. Lett. 111, 058101 (2013)

Phase transition to bundles of flexible supramolecular polymers

Authors: B. A. H. Huisman, P. G. Bolhuis, A. Fasolino

Abstract: We report Monte Carlo simulations of the self-assembly of supramolecular polymers based on a model of patchy particles. We find a first-order phase transition, characterized by hysteresis and nucleation, toward a solid bundle of polymers, of length much greater than the average gas phase length. We argue that the bundling transition is the supramolecular equivalent of the sublimation transition,… ▽ More We report Monte Carlo simulations of the self-assembly of supramolecular polymers based on a model of patchy particles. We find a first-order phase transition, characterized by hysteresis and nucleation, toward a solid bundle of polymers, of length much greater than the average gas phase length. We argue that the bundling transition is the supramolecular equivalent of the sublimation transition, that results from a weak chain-chain interaction. We provide a qualitative equation of state that gives physical insight beyond the specific values of the parameters used in our simulations. △ Less

Submitted 29 November, 2007; originally announced November 2007.

Effect of excluded volume interactions on the interfacial properties of colloid-polymer mixtures

Authors: Andrea Fortini, Peter G. Bolhuis, Marjolein Dijkstra

Abstract: We report a numerical study of equilibrium phase-diagrams and interfacial properties of bulk and confined colloid-polymer mixtures using grand canonical Monte Carlo simulations. Colloidal particles are treated as hard spheres, while the polymer chains are described as soft repulsive spheres. The polymer-polymer, colloid-polymer, and wall-polymer interactions are described by density-dependent po… ▽ More We report a numerical study of equilibrium phase-diagrams and interfacial properties of bulk and confined colloid-polymer mixtures using grand canonical Monte Carlo simulations. Colloidal particles are treated as hard spheres, while the polymer chains are described as soft repulsive spheres. The polymer-polymer, colloid-polymer, and wall-polymer interactions are described by density-dependent potentials derived by Bolhuis and Louis [Macromolecules, 35 (2002), p.1860]. We compared our results with those of the Asakura-Oosawa-Vrij model, that treats the polymers as ideal particles. We find that the number of polymers needed to drive the demixing transition is larger for the interacting polymers, and that the gas-liquid interfacial tension is smaller. When the system is confined between two parallel hard plates, we find capillary condensation. Compared with the AOV model, we find that the excluded volume interactions between the polymers suppress capillary condensation. In order to induce capillary condensation, smaller undersaturations and smaller plate separations are needed in comparison with ideal polymers. △ Less

Submitted 7 November, 2007; originally announced November 2007.

Comments: 9 pages, 10 figures, accepted for publication in the J. Chem. Phys

Journal ref: J. Chem. Phys. 128, 024904 (2008)

Elaborating Transition Interface Sampling Methods

Authors: Titus S. van Erp, Peter G. Bolhuis

Abstract: We review two recently developed efficient methods for calculating rate constants of processes dominated by rare events in high-dimensional complex systems. The first is transition interface sampling (TIS), based on the measurement of effective fluxes through hypersurfaces in phase space. TIS improves efficiency with respect to standard transition path sampling (TPS) rate constant techniques, be… ▽ More We review two recently developed efficient methods for calculating rate constants of processes dominated by rare events in high-dimensional complex systems. The first is transition interface sampling (TIS), based on the measurement of effective fluxes through hypersurfaces in phase space. TIS improves efficiency with respect to standard transition path sampling (TPS) rate constant techniques, because it allows a variable path length and is less sensitive to recrossings. The second method is the partial path version of TIS. Developed for diffusive processes, it exploits the loss of long time correlation. We discuss the relation between the new techniques and the standard reactive flux methods in detail. Path sampling algorithms can suffer from ergodicity problems, and we introduce several new techniques to alleviate these problems, notably path swapping, stochastic configurational bias Monte Carlo shooting moves and order-parameter free path sampling. In addition, we give algorithms to calculate other interesting properties from path ensembles besides rate constants, such as activation energies and reaction mechanisms. △ Less

Submitted 6 May, 2004; originally announced May 2004.

Comments: 36 pages, 5 figures

Investigating Rare Events by Transition Interface Sampling

Authors: Daniele Moroni, Titus S. van Erp, Peter G. Bolhuis

Abstract: We briefly review simulation schemes for the investigation of rare transitions and we resume the recently introduced Transition Interface Sampling, a method in which the computation of rate constants is recast into the computation of fluxes through interfaces dividing the reactant and product state. We briefly review simulation schemes for the investigation of rare transitions and we resume the recently introduced Transition Interface Sampling, a method in which the computation of rate constants is recast into the computation of fluxes through interfaces dividing the reactant and product state. △ Less

Submitted 25 November, 2003; originally announced November 2003.

Comments: 12 pages, 1 figure, contributed paper to the proceedings of NEXT 2003, Second Sardinian International Conference on News and Expectations in Thermostatistics, 21-28 Sep 2003, Cagliari (Italy)

Rate constants for diffusive processes by partial path sampling

Authors: Daniele Moroni, Peter G. Bolhuis, Titus S. van Erp

Abstract: We introduce a path sampling method for the computation of rate constants for systems with a highly diffusive character. Based on the recently developed algorithm of transition interface sampling (TIS) this procedure increases the efficiency by sampling only parts of complete transition trajectories confined within a certain region. The algorithm assumes the loss of memory for highly diffusive p… ▽ More We introduce a path sampling method for the computation of rate constants for systems with a highly diffusive character. Based on the recently developed algorithm of transition interface sampling (TIS) this procedure increases the efficiency by sampling only parts of complete transition trajectories confined within a certain region. The algorithm assumes the loss of memory for highly diffusive progression along the reaction coordinate. We compare the new technique to the TIS method for a simple diatomic system and show that the computation time of the new method scales linearly, instead of quadraticaly, with the length of the diffusive barrier. The validity of the memory loss assumption is also discussed. △ Less

Submitted 20 October, 2003; originally announced October 2003.

Comments: 12 pages, including 8 figures, RevTeX4

Journal ref: J. Chem. Phys. 120, 4055 (2004)

The Asakura-Oosawa model in the protein limit: the role of many-body interactions

Authors: A. Moncho-Jorda, A. A. Louis, P. G. Bolhuis, R. Roth

Abstract: We study the Asakura-Oosawa model in the "protein limit", where the penetrable sphere radius $R_{AO}$ is much greater than the hard sphere radius $R_c$. The phase behaviour and structure calculated with a full many-body treatment show important qualitative differences when compared to a description based on pair potentials alone. The overall effect of the many-body interactions is repulsive. We study the Asakura-Oosawa model in the "protein limit", where the penetrable sphere radius $R_{AO}$ is much greater than the hard sphere radius $R_c$. The phase behaviour and structure calculated with a full many-body treatment show important qualitative differences when compared to a description based on pair potentials alone. The overall effect of the many-body interactions is repulsive. △ Less

Submitted 11 September, 2003; originally announced September 2003.

Comments: 9 pages and 11 figures, submitted to J. Phys.: Condensed Matter, special issue "Effective many-body interactions and correlations in soft matter"

A Novel Path Sampling Method for the Calculation of Rate Constants

Authors: Titus S. van Erp, Daniele Moroni, Peter G. Bolhuis

Abstract: We derive a novel efficient scheme to measure the rate constant of transitions between stable states separated by high free energy barriers in a complex environment within the framework of transition path sampling. The method is based on directly and simultaneously measuring the fluxes through many phase space interfaces and increases the efficiency with at least a factor of two with respect to… ▽ More We derive a novel efficient scheme to measure the rate constant of transitions between stable states separated by high free energy barriers in a complex environment within the framework of transition path sampling. The method is based on directly and simultaneously measuring the fluxes through many phase space interfaces and increases the efficiency with at least a factor of two with respect to existing transition path sampling rate constant algorithms. The new algorithm is illustrated on the isomerization of a diatomic molecule immersed in a simple fluid. △ Less

Submitted 28 October, 2002; originally announced October 2002.

Comments: 14 pages, including 13 figures, RevTeX4

Journal ref: J. Chem. Phys. 118, 7762 (2003)

Colloid-polymer mixtures in the protein limit

Authors: Peter G. Bolhuis, Evert Jan Meijer, Ard A. Louis

Abstract: We computed the phase-separation behavior and effective interactions of colloid-polymer mixtures in the "protein limit", where the polymer radius of gyration is much larger than the colloid radius. For ideal polymers, the critical colloidal packing fraction tends to zero, whereas for interacting polymers in a good solvent the behavior is governed by a universal binodal, implying a constant criti… ▽ More We computed the phase-separation behavior and effective interactions of colloid-polymer mixtures in the "protein limit", where the polymer radius of gyration is much larger than the colloid radius. For ideal polymers, the critical colloidal packing fraction tends to zero, whereas for interacting polymers in a good solvent the behavior is governed by a universal binodal, implying a constant critical colloid packing fraction. In both systems the depletion interaction is not well described by effective pair potentials but requires the incorporation of many-body contributions. △ Less

Submitted 23 October, 2002; originally announced October 2002.

Comments: 4 pages, 3 figures, submitted to Physical Review Letters

Influence of polymer excluded volume on the phase behavior of colloid-polymer mixtures

Authors: P. G. Bolhuis, A. A. Louis, J-P Hansen

Abstract: We determine the depletion-induced phase-behavior of hard sphere colloids and interacting polymers by large-scale Monte Carlo simulations using very accurate coarse-graining techniques. A comparison with standard Asakura-Oosawa model theories and simulations shows that including excluded volume interactions between polymers leads to qualitative differences in the phase diagrams. These effects be… ▽ More We determine the depletion-induced phase-behavior of hard sphere colloids and interacting polymers by large-scale Monte Carlo simulations using very accurate coarse-graining techniques. A comparison with standard Asakura-Oosawa model theories and simulations shows that including excluded volume interactions between polymers leads to qualitative differences in the phase diagrams. These effects become increasingly important for larger relative polymer size. Our simulations results agree quantitatively with recent experiments. △ Less

Submitted 5 June, 2002; originally announced June 2002.

Comments: 5 pages, 4 figures submitted to Physical Review Letters

How to derive and parameterize effective potentials in colloid-polymer mixtures

Authors: P. G. Bolhuis, A. A. Louis

Abstract: Polymer chains in colloid-polymer mixtures can be coarse-grained by replacing them with single soft particles interacting via effective polymer-polymer and polymer-colloid pair potentials. Here we describe in detail how Ornstein-Zernike inversion techniques, originally developed for atomic and molecular fluids, can be generalized to complex fluids and used to derive effective potentials from com… ▽ More Polymer chains in colloid-polymer mixtures can be coarse-grained by replacing them with single soft particles interacting via effective polymer-polymer and polymer-colloid pair potentials. Here we describe in detail how Ornstein-Zernike inversion techniques, originally developed for atomic and molecular fluids, can be generalized to complex fluids and used to derive effective potentials from computer simulations on a microscopic level. In particular, we consider polymer solutions for which we derive effective potentials between the centers of mass, and also between mid-points or end-points from simulations of self-avoiding walk polymers. In addition, we derive effective potentials for polymers near a hard wall or a hard sphere. We emphasize the importance of including both structural and thermodynamic information (through sum-rules) from the underlying simulations. In addition we develop a simple numerical scheme to optimize the parameterization of the density dependent polymer-polymer, polymer-wall and polymer-sphere potentials for dilute and semi-dilute polymer densities, thus opening up the possibility of performing large-scale simulations of colloid-polymer mixtures. The methods developed here should be applicable to a much wider range effective potentials in complex fluids. △ Less

Submitted 21 May, 2002; originally announced May 2002.

Comments: uses revtex4.cls; submitted for archival purposes

Journal ref: Macromolecules 35, 1860 (2002)

Polymer induced depletion potentials in polymer-colloid mixtures

Authors: A. A. Louis, P. G. Bolhuis, E. J. Meijer, J. P. Hansen

Abstract: The depletion interactions between two colloidal plates or between two colloidal spheres, induced by interacting polymers in a good solvent, are calculated theoretically and by computer simulations. A simple analytical theory is shown to be quantitatively accurate for case of two plates. A related depletion potential is derived for two spheres; it also agrees very well with direct computer simul… ▽ More The depletion interactions between two colloidal plates or between two colloidal spheres, induced by interacting polymers in a good solvent, are calculated theoretically and by computer simulations. A simple analytical theory is shown to be quantitatively accurate for case of two plates. A related depletion potential is derived for two spheres; it also agrees very well with direct computer simulations. Theories based on ideal polymers show important deviations with increasing polymer concentration: They overestimate the range of the depletion potential between two plates or two spheres at all densities, with the largest relative change occurring in the dilute regime. They underestimate the well depth at contact for the case of two plates, but overestimate it for two spheres. Depletion potentials are also calculated using a coarse graining approach which represents the polymers as ``soft colloids'': good agreement is found in the dilute regime. Finally, the effect of the polymers on colloid-colloid osmotic virial coefficients is related to phase behavior of polymer-colloid mixtures. △ Less

Submitted 6 March, 2002; originally announced March 2002.

Comments: revtex

Journal ref: J. Chem. Phys. 117, 1893 (2002)

Coarse-graining polymers as soft colloids

Authors: A. A. Louis, P. G. Bolhuis, R. Finken, V. Krakoviack, E. J. Meijer, J. P. Hansen

Abstract: We show how to coarse grain polymers in a good solvent as single particles, interacting with density-independent or density-dependent interactions. These interactions can be between the centres of mass, the mid-points or end-points of the polymers. We also show how to extend these methods to polymers in poor solvents and mixtures of polymers. Treating polymers as soft colloids can greatly speed… ▽ More We show how to coarse grain polymers in a good solvent as single particles, interacting with density-independent or density-dependent interactions. These interactions can be between the centres of mass, the mid-points or end-points of the polymers. We also show how to extend these methods to polymers in poor solvents and mixtures of polymers. Treating polymers as soft colloids can greatly speed up the simulation of complex many-polymer systems, including polymer-colloid mixtures. △ Less

Submitted 6 December, 2001; originally announced December 2001.

Comments: to appear in Physica A, special STATPHYS 2001 edition. Content of invited talk by AAL

Journal ref: Physica A 306, 251 (2002)

Density profiles and surface tensions of polymers near colloidal surfaces

Authors: A. A. Louis, P. G. Bolhuis, E. J. Meijer, J. P. Hansen

Abstract: The surface tension of interacting polymers in a good solvent is calculated theoretically and by computer simulations for a planar wall geometry and for the insertion of a single colloidal hard-sphere. This is achieved for the planar wall and for the larger spheres by an adsorption method, and for smaller spheres by a direct insertion technique. Results for the dilute and semi-dilute regimes are… ▽ More The surface tension of interacting polymers in a good solvent is calculated theoretically and by computer simulations for a planar wall geometry and for the insertion of a single colloidal hard-sphere. This is achieved for the planar wall and for the larger spheres by an adsorption method, and for smaller spheres by a direct insertion technique. Results for the dilute and semi-dilute regimes are compared to results for ideal polymers, the Asakura-Oosawa penetrable-sphere model, and to integral equations, scaling and renormalization group theories. The largest relative changes with density are found in the dilute regime, so that theories based on non-interacting polymers rapidly break down. A recently developed ``soft colloid'' approach to polymer-colloid mixtures is shown to correctly describe the one-body insertion free-energy and the related surface tension. △ Less

Submitted 29 November, 2001; v1 submitted 27 November, 2001; originally announced November 2001.

Journal ref: J. Chem. Phys. 116, 10547 (2002)

Many-body interactions and correlations in coarse-grained descriptions of polymer solutions

Authors: P. G. Bolhuis, A. A. Louis, J. P. Hansen

Abstract: We calculate the two, three, four, and five-body (state independent) effective potentials between the centers of mass (CM) of self avoiding walk polymers by Monte-Carlo simulations. For full overlap, these coarse-grained n-body interactions oscillate in sign as (-1)^n, and decrease in absolute magnitude with increasing n. We find semi-quantitative agreement with a scaling theory, and use this to… ▽ More We calculate the two, three, four, and five-body (state independent) effective potentials between the centers of mass (CM) of self avoiding walk polymers by Monte-Carlo simulations. For full overlap, these coarse-grained n-body interactions oscillate in sign as (-1)^n, and decrease in absolute magnitude with increasing n. We find semi-quantitative agreement with a scaling theory, and use this to discuss how the coarse-grained free energy converges when expanded to arbitrary order in the many-body potentials. We also derive effective {\em density dependent} 2-body potentials which exactly reproduce the pair-correlations between the CM of the self avoiding walk polymers. The density dependence of these pair potentials can be largely understood from the effects of the {\em density independent} 3-body potential. Triplet correlations between the CM of the polymers are surprisingly well, but not exactly, described by our coarse-grained effective pair potential picture. In fact, we demonstrate that a pair-potential cannot simultaneously reproduce the two and three body correlations in a system with many-body interactions. However, the deviations that do occur in our system are very small, and can be explained by the direct influence of 3-body potentials. △ Less

Submitted 9 March, 2001; originally announced March 2001.

Comments: 11 pages, 1 table, 9 figures, RevTeX (revtex.cls)

Journal ref: Phys. Rev. E. 64, 021801 (2001)

Accurate effective pair potentials for polymer solutions

Authors: P. G. Bolhuis, A. A. Louis, J. P. Hansen, E. J. Meijer

Abstract: Dilute or semi-dilute solutions of non-intersecting self-avoiding walk (SAW) polymer chains are mapped onto a fluid of ``soft'' particles interacting via an effective pair potential between their centers of mass. This mapping is achieved by inverting the pair distribution function of the centers of mass of the original polymer chains, using integral equation techniques from the theory of simple… ▽ More Dilute or semi-dilute solutions of non-intersecting self-avoiding walk (SAW) polymer chains are mapped onto a fluid of ``soft'' particles interacting via an effective pair potential between their centers of mass. This mapping is achieved by inverting the pair distribution function of the centers of mass of the original polymer chains, using integral equation techniques from the theory of simple fluids. The resulting effective pair potential is finite at all distances, has a range of the order of the radius of gyration, and turns out to be only moderately concentration-dependent. The dependence of the effective potential on polymer length is analyzed in an effort to extract the scaling limit. The effective potential is used to derive the osmotic equation of state, which is compared to simulation data for the full SAW segment model, and to the predictions of renormalization group calculations. A similar inversion procedure is used to derive an effective wall-polymer potential from the center of mass density profiles near the wall, obtained from simulations of the full polymer segment model. The resulting wall-polymer potential turns out to depend strongly on bulk polymer concentration when polymer-polymer correlations are taken into account, leading to a considerable enhancement of the effective repulsion with increasing concentration. The effective polymer-polymer and wall-polymer potentials are combined to calculate the depletion interaction induced by SAW polymers between two walls. The calculated depletion interaction agrees well with the ``exact'' results from much more computer-intensive direct simulation of the full polymer-segment model, and clearly illustrates the inadequacy -- in the semi-dilute regime -- of the standard Asakura-Oosawa approximation based on the assumption of non-interacting polymer coils. △ Less

Submitted 6 September, 2000; originally announced September 2000.

Comments: 18 pages, 24 figures, ReVTeX, submitted to J. Chem. Phys

Journal ref: J. Chem. Phys. 114, 4296 (2001)

Mean Field Fluid Behavior of the Gaussian Core Model

Authors: A. A. Louis, P. G. Bolhuis, J. P. Hansen

Abstract: We show that the Gaussian core model of particles interacting via a penetrable repulsive Gaussian potential, first considered by Stillinger (J. Chem. Phys. 65, 3968 (1976)), behaves like a weakly correlated ``mean field fluid'' over a surprisingly wide density and temperature range. In the bulk the structure of the fluid phase is accurately described by the random phase approximation for the dir… ▽ More We show that the Gaussian core model of particles interacting via a penetrable repulsive Gaussian potential, first considered by Stillinger (J. Chem. Phys. 65, 3968 (1976)), behaves like a weakly correlated ``mean field fluid'' over a surprisingly wide density and temperature range. In the bulk the structure of the fluid phase is accurately described by the random phase approximation for the direct correlation function, and by the more sophisticated HNC integral equation. The resulting pressure deviates very little from a simple, mean-field like, quadratic form in the density, while the low density virial expansion turns out to have an extremely small radius of convergence. Density profiles near a hard wall are also very accurately described by the corresponding mean-field free-energy functional. The binary version of the model exhibits a spinodal instability against de-mixing at high densities. Possible implications for semi-dilute polymer solutions are discussed. △ Less

Submitted 6 July, 2000; v1 submitted 4 July, 2000; originally announced July 2000.

Comments: 13 pages, 2 columns, ReVTeX epsfig,multicol,amssym, 15 figures; submitted to Phys. Rev. E (change: important reference added)

Journal ref: Phys. Rev. E. 62, 7961 (2000)

Can Polymer Coils be modeled as "Soft Colloids"?

Authors: A. A. Louis, P. G. Bolhuis, J. P. Hansen, E. J. Meijer

Abstract: We map dilute or semi-dilute solutions of non-intersecting polymer chains onto a fluid of ``soft'' particles interacting via a concentration dependent effective pair potential, by inverting the pair distribution function of the centers of mass of the initial polymer chains. A similar inversion is used to derive an effective wall-polymer potential; these potentials are combined to successfully re… ▽ More We map dilute or semi-dilute solutions of non-intersecting polymer chains onto a fluid of ``soft'' particles interacting via a concentration dependent effective pair potential, by inverting the pair distribution function of the centers of mass of the initial polymer chains. A similar inversion is used to derive an effective wall-polymer potential; these potentials are combined to successfully reproduce the calculated exact depletion interaction induced by non-intersecting polymers between two walls. The mapping opens up the possibility of large-scale simulations of polymer solutions in complex geometries. △ Less

Submitted 8 September, 2000; v1 submitted 12 June, 2000; originally announced June 2000.

Comments: 4 pages, 3 figures ReVTeX[epsfig,multicol,amssymb] references updated

Journal ref: Phys. Rev. Lett. 85, 2522 (2000)