Editing Kern and Frenkel patchy model
Jump to navigation
Jump to search
The edit can be undone. Please check the comparison below to verify that this is what you want to do, and then publish the changes below to finish undoing the edit.
Latest revision | Your text | ||
Line 1: | Line 1: | ||
The '''Kern and Frenkel''' <ref>[http://dx.doi.org/10.1063/1.1569473 Norbert Kern and Daan Frenkel "Fluid–fluid coexistence in colloidal systems with short-ranged strongly directional attraction", Journal of Chemical Physics 118, 9882 (2003)]</ref> [[Patchy particles |patchy model]] published in 2003 is an amalgamation of the [[hard sphere model]] with | The '''Kern and Frenkel''' <ref>[http://dx.doi.org/10.1063/1.1569473 Norbert Kern and Daan Frenkel "Fluid–fluid coexistence in colloidal systems with short-ranged strongly directional attraction", Journal of Chemical Physics 118, 9882 (2003)]</ref> [[Patchy particles |patchy model]] published in 2003 is an amalgamation of the [[hard sphere model]] with | ||
attractive [[Square well model | square well]] patches (HSSW). The model was originally developed by Bol (1982),<ref>[http://dx.doi.org/10.1080/00268978200100461 W. Bol "Monte Carlo simulations of fluid systems of waterlike molecules", Molecular Physics '''45''' pp. 605-616 (1982)]</ref> and later Chapman (1988) <ref | attractive [[Square well model | square well]] patches (HSSW). The model was originally developed by Bol (1982),<ref>[http://dx.doi.org/10.1080/00268978200100461 W. Bol "Monte Carlo simulations of fluid systems of waterlike molecules", Molecular Physics '''45''' pp. 605-616 (1982)]</ref> and later Chapman(1988) <ref>[W.G. Chapman, Doctoral Thesis, Cornell University (1988)]</ref> <ref>[G. Jackson, W.G. Chapman, K.E. Gubbins, Molecular Physics 65, 1-31 (1988)]</ref> reinvented the model as the basis for numerous articles describing properties of associating particles from molecular simulation and theory. The computational advantage of Bol's model is that only a simple dot product is required to determine if a particle is in the bonding orientation. | ||
The potential has an angular aspect, given by (Eq. 1) | The potential has an angular aspect, given by (Eq. 1) | ||
Line 29: | Line 29: | ||
where <math>\delta</math> is the solid angle of a patch (<math>\alpha, \beta, ...</math>) whose axis is <math>\hat{e}</math> (see Fig. 1 of Ref. 1), forming a conical segment. | where <math>\delta</math> is the solid angle of a patch (<math>\alpha, \beta, ...</math>) whose axis is <math>\hat{e}</math> (see Fig. 1 of Ref. 1), forming a conical segment. | ||
==Multiple patches== | ==Multiple patches== | ||
The "two-patch" and "four-patch" Bol (Chapman or Kern and Frenkel) model was extensively studied by Chapman and co-workers for bulk and interfacial systems using hard sphere and Lennard-Jones | The "two-patch" and "four-patch" Bol (Chapman or Kern and Frenkel) model was extensively studied by Chapman and co-workers for bulk and interfacial systems using hard sphere and Lennard-Jones references. Later other groups, including Sciortino and co-workers, considered stronger association energies for the "two-patch" hard sphere reference <ref name="bianchi">[http://dx.doi.org/10.1063/1.2730797 F. Sciortino, E. Bianchi, J. Douglas and P. Tartaglia "Self-assembly of patchy particles into polymer chains: A parameter-free comparison between Wertheim theory and Monte Carlo simulation", Journal of Chemical Physics '''126''' 194903 (2007)]</ref><ref>[http://dx.doi.org/10.1063/1.3415490 Achille Giacometti, Fred Lado, Julio Largo, Giorgio Pastore, and Francesco Sciortino "Effects of patch size and number within a simple model of patchy colloids", Journal of Chemical Physics 132, 174110 (2010)]</ref><ref name="rovigatti">[http://dx.doi.org/10.1063/1.4737930 José Maria Tavares, Lorenzo Rovigatti, and Francesco Sciortino "Quantitative description of the self-assembly of patchy particles into chains and rings", Journal of Chemical Physics '''137''' 044901 (2012)]</ref>. | ||
==Four patches== | ==Four patches== | ||
Line 41: | Line 41: | ||
</math> | </math> | ||
then the patch cannot be involved in more than one bond. Enforcing this condition makes it possible to compare the simulations results with [[Wertheim's first order thermodynamic perturbation theory (TPT1)| Wertheim theory]] | then the patch cannot be involved in more than one bond. Enforcing this condition makes it possible to compare the simulations results with [[Wertheim's first order thermodynamic perturbation theory (TPT1)| Wertheim theory]] <ref name="bianchi"/><ref name="rovigatti"/> | ||
==Hard ellipsoid model== | ==Hard ellipsoid model== |