Editing Kern and Frenkel patchy model

Jump to navigation Jump to search
Warning: You are not logged in. Your IP address will be publicly visible if you make any edits. If you log in or create an account, your edits will be attributed to your username, along with other benefits.

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]] 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 name="Chapman">[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.
attractive [[Square well model | square well]] patches (HSSW). The potential has an angular aspect, given by (Eq. 1)
The potential has an angular aspect, given by (Eq. 1)




:<math>\Phi_{ij}({\mathbf r}_{ij}; \tilde{ {\mathbf \Omega}}_i, \tilde{ {\mathbf \Omega}}_j)  =\Phi_{ij}^{ \mathrm{HSSW}}({\mathbf r}_{ij}) \cdot f(\tilde{ {\mathbf \Omega}}_i, \tilde{ {\mathbf \Omega}}_j) </math>
:<math>\Phi_{ij}({\mathbf r}_{ij}; \tilde{{\mathbf \Omega}}_i, \tilde{{\mathbf \Omega}}_j)  =\Phi_{ij}^{ \mathrm{HSSW}}({\mathbf r}_{ij}) \cdot f(\tilde{{\mathbf \Omega}}_i, \tilde{{\mathbf \Omega}}_j) </math>




Line 20: Line 19:


:<math>
:<math>
f_{ij} \left(\hat{ {\mathbf r}}_{ij}; \tilde{ {\mathbf \Omega}}_i, \tilde{ {\mathbf \Omega}}_j \right) =  
f_{ij} \left(\hat{ {\mathbf r}}_{ij}; \tilde{{\mathbf \Omega}}_i, \tilde{{\mathbf \Omega}}_j \right) =  
\left\{ \begin{array}{clc}
\left\{ \begin{array}{clc}
1        & \mathrm{if}        & \left\{ \begin{array}{ccc}    &  (\hat{e}_\alpha\cdot\hat{r}_{ij} \geq \cos \delta) & \mathrm{for~some~patch~\alpha~on~}i  \\  
1        & \mathrm{if}        & \left\{ \begin{array}{ccc}    &  (\hat{e}_\alpha\cdot\hat{r}_{ij} \leq \cos \delta) & \mathrm{for~some~patch~\alpha~on~}i  \\  
                                                             \mathrm{and} & (\hat{e}_\beta\cdot\hat{r}_{ji} \geq \cos \delta)  & \mathrm{for~some~patch~\beta~on~}j  \end{array} \right. \\
                                                             \mathrm{and} & (\hat{e}_\beta\cdot\hat{r}_{ji} \leq \cos \delta)  & \mathrm{for~some~patch~\beta~on~}j  \end{array} \right. \\
0        & \mathrm{otherwise} &  \end{array} \right.
0        & \mathrm{otherwise} &  \end{array} \right.
</math>
</math>


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==
==Two 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 reference systems.  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>.
The "two-patch" Kern and Frenkel model has been extensively studied by Sciortino and co-workers <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 40:
</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]] <ref name="Chapman"/><ref name="bianchi"/><ref name="rovigatti"/>
then the patch cannot be involved in more than one bond. Enforcing this condition makes it possible to compare the simulations results with Wertheim theory <ref name="bianchi"/><ref name="rovigatti"/>


==Hard ellipsoid model==
The [[hard ellipsoid model]] has also been used as the 'nucleus' of the Kern and Frenkel patchy model <ref>[http://dx.doi.org/10.1063/1.4969074  T. N. Carpency, J. D. Gunton and J. M. Rickman "Phase behavior of patchy spheroidal fluids", Journal of Chemical Physics '''145''' 214904 (2016)]</ref>.
==References==
==References==
<references/>
<references/>
Line 50: Line 47:
*[http://dx.doi.org/10.1063/1.3689308 Christoph Gögelein, Flavio Romano, Francesco Sciortino, and Achille Giacometti "Fluid-fluid and fluid-solid transitions in the Kern-Frenkel model from Barker-Henderson thermodynamic perturbation theory", Journal of Chemical Physics '''136''' 094512 (2012)]
*[http://dx.doi.org/10.1063/1.3689308 Christoph Gögelein, Flavio Romano, Francesco Sciortino, and Achille Giacometti "Fluid-fluid and fluid-solid transitions in the Kern-Frenkel model from Barker-Henderson thermodynamic perturbation theory", Journal of Chemical Physics '''136''' 094512 (2012)]
*[http://dx.doi.org/10.1063/1.4722477 Emanuela Bianchi, Günther Doppelbauer, Laura Filion, Marjolein Dijkstra, and Gerhard Kahl "Predicting patchy particle crystals: Variable box shape simulations and evolutionary algorithms", Journal of Chemical Physics '''136''' 214102 (2012)]
*[http://dx.doi.org/10.1063/1.4722477 Emanuela Bianchi, Günther Doppelbauer, Laura Filion, Marjolein Dijkstra, and Gerhard Kahl "Predicting patchy particle crystals: Variable box shape simulations and evolutionary algorithms", Journal of Chemical Physics '''136''' 214102 (2012)]
*[http://dx.doi.org/10.1063/1.4960423  Z. Preisler, T. Vissers, F. Smallenburg and F. Sciortino "Crystals of Janus colloids at various interaction ranges", Journal of Chemical Physics '''145''' 064513 (2016)]




[[category: models]]
[[category: models]]
Please note that all contributions to SklogWiki are considered to be released under the Creative Commons Attribution Non-Commercial Share Alike (see SklogWiki:Copyrights for details). If you do not want your writing to be edited mercilessly and redistributed at will, then do not submit it here.
You are also promising us that you wrote this yourself, or copied it from a public domain or similar free resource. Do not submit copyrighted work without permission!

To edit this page, please answer the question that appears below (more info):

Cancel Editing help (opens in new window)