# Reference Interaction Site Model
and Optimized Perturbation theories of

colloidal dumbbells with increasing anisotropy

###### Abstract

We investigate thermodynamic properties of anisotropic colloidal dumbbells in the frameworks provided by the Reference Interaction Site Model (RISM) theory and an Optimized Perturbation Theory (OPT), this latter based on a fourth-order high-temperature perturbative expansion of the free energy, recently generalized to molecular fluids. Our model is constituted by two identical tangent hard spheres surrounded by square-well attractions with same widths and progressively different depths. Gas-liquid coexistence curves are obtained by predicting pressures, free energies, and chemical potentials. In comparison with previous simulation results, RISM and OPT agree in reproducing the progressive reduction of the gas-liquid phase separation as the anisotropy of the interaction potential becomes more pronounced; in particular, the RISM theory provides reasonable predictions for all coexistence curves, bar the strong anisotropy regime, whereas OPT performs generally less well. Both theories predict a linear dependence of the critical temperature on the interaction strength, reproducing in this way the mean-field behavior observed in simulations; the critical density – that drastically drops as the anisotropy increases — turns to be less accurate. Our results appear as a robust benchmark for further theoretical studies, in support to the simulation approach, of self-assembly in model colloidal systems.

^{†}

^{†}thanks: Corresponding author, email: gmunao@unime.it

## I Introduction

In recent years, a considerable number of experimental Bon:14 (); Ma:14 (); Kraft:12 (); Nagao:11 (); Yoon-Chem (); Forster (); Hosein () and numerical Chapela:11 (); Chapela2:11 (); Miller:09 (); Del-gado:11 (); Chong-prl () studies have been devoted to the investigation of phase behavior and self-assembly properties of colloidal dumbbells. Specific interest in this class stems essentially from the possibility to accurately tune the aspect ratio of the constituting spheres, as well as their interaction properties, so to obtain rich and fascinating phase behaviors.Munao:14 (); Munao:15 (); Avvisati-JCP () In particular, if one of the two particles is solvophilic and the other one is solvophobic, colloidal dumbbells represent a simple example of surfactant (Janus dumbbells),Janusdumbbell (); Janusdumbbell1 () constituting a molecular generalization of Janus spherical particles,granick (); Janusweitz (); Adv:10 () largely investigated in the last years because of their peculiar self-assembly properties and phase behaviors.cacciuto (); Chen:11 (); Granick-Lang (); Granick-Nat (); janusprl (); januslong (); Vissers:jcp ()

While a large number of experimental and numerical investigations have involved colloidal dumbbells, considerable less attention has been paid to theoretical studies: the framework of potential energy landscape has been adopted to describe the self-assembly of dumbbells into helices;Fejer:14 () the competition between self-assembly and phase separation has been documented by integral equation theories;Munao:PCCP () the same approach Munao-cpl (); Munao:09 () and fundamental measure theory Marechal:11 () have been adopted to investigate structure and thermodynamics of hard dumbbells. Notwithstanding all such studies, a systematic theoretical description of the phase behavior of colloidal dumbbells is still lacking. In order to fill this gap, here we calculate systematically free energy, pressure, chemical potential and phase diagram of anisotropic colloidal dumbbells, so to investigate how such properties change by tuning the interaction potential. We consider to this task a series of model dumbbells constituted by two identical tangent hard spheres surrounded by square-well attractions with same widths and progressively different depths. We carry out our investigation in the frameworks provided by the Reference Interaction Site Model theory (RISM)chandler:1972 () and a recently proposed Optimized Perturbation Theory (OPT).DelRio:09 (); Gamez:14 () Theoretical predictions for the gas-liquid phase coexistence are compared with previous Monte Carlo results by usMunao:14 () and other authors.Cai:12 () The present work extends our preliminary OPT analysis of the symmetric case Gamez:14 () and our RISM study on a slightly different dumbbell model.Munao:PCCP ()

As for the RISM theory,chandler:1972 () it constitutes a molecular generalization of the Ornstein-Zernike theory of simple fluids.Hansennew () Originally devised for rigid molecules constituted by hard spheres,lowden:5228 () the theory was later extended to more complex systems, like colloidal models; in particular, it was used to investigate the thermodynamic and structural properties of discotic lamellar colloids,Harnau:01 (); Harnau:05 () the self-assembly in diblock copolymers modeled as “ultrasoft” colloids,Hansen:06 () the interaction between colloidal particles and macromolecules,Khalatur:97 () the crystallization and solvation properties of nanoparticles in aqueous solutions,Kung:10 () the liquid structure of tetrahedral colloidal particlesMunao:11 (), and the self-assembly properties of Janus rods.Tripathy-RISM ()

As for the perturbation theory, this scheme played a significant role in describing the thermodynamic properties of a large variety of fluids.Valadez:12 (); Gamez:11 (); Benavides:07 (); Elliott:05 (); Cui:02 (); Chapela:89 () The theory was originally proposed as a high-temperature expansion by several groups and then recast in a very efficient computational tool by Barker and Henderson.BH1 (); BH2 (); BH3 () It has been recently adapted to cope with Janus Giacometti:14 () and tri-block Janus Gogelein:12 () fluids. Here we use an optimized version of this theory based on a fourth-order high-temperature perturbative expansion of the free energy, recently developed to fit a large set of simulation data for square-well fluids.DelRio:09 () Recently, one of the authors generalized this approach to deal with more complex systems, like molecular fluids.Gamez:14 () Within this formalism, the effective energy and range of square-well interactions depend on the molecular anisotropy. In this way, it is possible to calculate analytical equations of state for a large variety of molecular models like for instance dimers, trimers, chains and spherocylinders.Gamez:14 ()

The importance of assessing the accuracy of these theoretical approaches is crucial, as a large number of demanding numerical calculations Munao:15 () is typically required to investigate the experimentally relevant heterogeneous dumbbells case.Kraft:12 () A reliable theoretical tool would then help to focus on those systems of great interest.

The paper is organized as follows: in Section II we provide the details of interaction potential and RISM and OPT schemes. Results are presented and discussed in Section III and conclusions follow in Section IV.

## Ii Models and theories

### ii.1 Interaction potentials

The sequence of models investigated in this work is schematically represented in Fig. 1: they are constituted by two identical tangent hard spheres of diameter , labeled as 1 and 2, interacting with sites 1 and 2 of another dumbbell at distance via square-well (SW) attractions with same width () and different depths :

(1) |

Parameters and provide, respectively, the unit of length and energy, in which terms we define the reduced temperature (with as the Boltzmann constant), number density and pressure . In all calculations we have fixed ; then, we have studied a sequence of cases obtained by progressively reducing starting from the symmetric case, , to , 0.3, 0.1 (in units of ).

### ii.2 Reference interaction site model theory

In the RISM frameworkchandler:1972 () the pair structure of a fluid composed by identical two-site molecules is characterized by a set of four site-site intermolecular pair correlation functions where and are the site-site radial distribution functions. The are related to a set of intermolecular direct correlation functions by a matrix generalization of the Ornstein-Zernike equation for simple fluids, Hansennew () expressed in the -space as:

(2) |

where , , and are symmetric matrices; the elements are the Fourier transforms of the intramolecular correlation functions, written explicitly as:

(3) |

where the bond length is given either by , if , or by , otherwise. In order to solve the RISM equation (2), we have adopted in this work the Kovalenko-Hirata (KH) closure Kovalenko:99 (); Kovalenko:01 () that assumes — beside the exact expression — the direct correlation functions outside the hard-core to be approximated by a combination of the Mean Spherical Approximation (MSA) and HyperNetted Chain (HNC) expression:Hansennew ()

(4) | |||||

where and . We have implemented the numerical solution of the RISM/KH scheme by means of a standard iterative Picard algorithm, on a grid of 8192 points with a mesh .

In order to calculate the thermodynamic properties, we have used standard thermodynamic integration formulæ Hansennew () that already proven to be successful for SW dumbbells.Munao:PCCP () Specifically, we have calculated the excess free energy per particle via the energy route along constant-density paths according to:

(5) |

where is the internal energy of the fluid at a given inverse temperature and corresponds to the excess free energy of tangent hard-sphere dumbbells. For this, the most classical expressions known in the literature are an analytic expression due to Tildesley and Street Tildesley:80 () (and devised to fit their Monte Carlo data) and the Wertheim equation of state for tangent hard dumbbells,WerthPolym () reading:

(6) | |||||

where is the dumbbell packing fraction. We have explicitly verified that such equation agrees with the free energy expression for hard dumbbells provided in Ref. Tildesley:80, . Equation (6) has the additional advantage to avoid free parameters and the fit of simulation data, so that we will use it through all the calculations.

The excess pressure is then obtained through the (numerical) derivative of the l.h.s. of Eq. (5) with respect to the density, taken at constant temperature, while coexistence conditions are determined by searching for equal chemical potentials and pressures (at a given temperature) of the two coexisting phases. Finally, critical temperatures and densities are obtained by fitting the coexistence points through the scaling law for the densities and the law of rectilinear diameters with an effective critical exponent =0.32.frenkelsmit () We refer the reader to our previous work Munao:PCCP () for full details on the above procedure.

Using the symmetric model as a test case, we have also investigated the quality of predictions obtained by different approaches within the RISM framework. Specifically, we have studied the HNC closure, that assumes expressions in the top lines of Eq. (4) holding for all independently on the value of . Moreover, within the KH closure, we have studied a different route from structure to thermodynamics, based on closed formulæ (i.e requiring no thermodynamic integration) to calculate free energy and pressure. Such expressions, derived by Kovalenko and Hirata in Refs. Kovalenko:99, ; Kovalenko:01, , closely follow the corresponding relations deduced long ago Morita:60 () within the Ornstein-Zernike/HNC scheme for simple fluids and later generalized to the RISM/HNC scheme.Singer:85 () We refer the reader to the original works and our recent application Munao:11 () for full details.

### ii.3 Optimized perturbation theory

The starting point of OPT is given by the optimized equation of state recently proposed in Ref. DelRio:09, , where the free energy of a fluid composed by SW particles with well depth and width , at temperature and density is expressed as a perturbation series in the form:

(7) |

In Eq. (7) is the free energy of the reference hard-sphere fluid plus the ideal contribution, and are -th order perturbation terms. In Ref. DelRio:09, the series is truncated, and are calculated up to the fourth-order. The theory includes the Barker-Henderson BH1 () result as a truncation of the expansion to second order. The third and fourth orders, however, are not included as explicit higher order corrections to the Barker-Henderson theory, that would require uncontrolled truncations of the hierarchy of correlation functions, but as an optimized phenomenological calculation of coefficients, hinging on simulations results. DelRio:09 () It turns out that this theory for the SW fluid is accurate over nearly the whole density range and subcritical (low) temperatures. We refer the reader to the original work DelRio:09 () for detailed expressions of , as well as all other necessary machinery.

The OPT generalization at issue, pertinent to a fluid composed by SW dumbbells, is written explicitly as a sum of SW contributions, each taking separately into account the four interactions of Eq. (1), in the form:Gamez:14 ()

(8) |

As visible, each term in the free energy expansion is weighted by a temperature reduced by the corresponding well depth , is the density of sites 1 or 2, and takes the fixed value 0.5. Given the spherical symmetry of each site-site interaction, the molecular anisotropy is taken into account by the reference free energy, provided in our case by the Wertheim equation of state for tangent hard dumbbells of Eq. (6).WerthPolym (). This choice gives better results for non-convex particles in comparison with the possibility to define an effective hard-sphere diameter with the corresponding scaling of density.Gamez:14 (); it turns to be accurate for a moderate number of beads per chain.

Given Eq. (8), one can obtain the thermodynamic properties from usual thermodynamic relations. For instance, the compressibility factor

(9) |

can also be expressed as a high-temperature expansion,

(10) |

where .

As in the RISM framework, we have determined the gas-liquid phase boundaries by the condition of thermal, mechanical, and chemical equilibrium. Numerical methods have been used to this task, based in this case on a combination of the Levenberg-Marquardt algorithm with Gauss and steepest-descent methods.Meeter ()

## Iii Results

In Fig. 2 we compare the energy-route KH and OPT predictions for symmetric colloidal dumbbells, i.e. with in Eq. (1). Results for the pressure (a) and free energy (b) are reported along different isotherms across the critical temperature (estimated in our previous work).Munao:14 () KH and OPT predict the same trends both for the pressure and free energy; in particular, they provide evidence for a supercritical fluid at , whereas a van der Waals loop — heralding a gas-liquid phase separation — develops at and becomes well defined at . The quantitative agreement between the two theories slightly worsens at low temperatures, where OPT curves fall progressively below the corresponding KH ones; this circumstance implies a higher OPT critical temperature and, correspondingly, a wider gas-liquid coexistence curve. In Fig. 2c we compare KH gas-liquid coexistence curves — as obtained by the energy route and the closed formulæ — with Monte Carlo (MC) data by us Munao:14 () and other authors.Cai:12 () OPT predictions, already provided in Ref. Gamez:14, , are also reported for comparison. One can notice that KH reproduces reasonably well all features, whereas OPT, as observed, overestimates the two-phase region. Such a discrepancy could be possibly ascribed to the truncation involved in the OPT procedure (see Eq. (8)). In this case improvements could be achieved by incorporating higher order terms in the free energy expansion, so to take into account more accurately density fluctuations, along the lines described for instance in Ref. Salvino:92, . In general, the OPT scheme for SW dumbbells shares the same level of accuracy of the original approach for the SW atomic fluid, DelRio:09 () the latter performing slightly better in the critical region.

As for HNC, we have verified that at high temperatures no significant differences arise with the pressure calculated within the KH scheme (see the curves corresponding to in the inset of Fig. 2a). On the other hand, upon lowering the temperature, the numerical convergence of HNC becomes progressively more problematic, thus precluding from further exploration a progressively larger interval of conditions. As a consequence, the conditions for gas-liquid coexistence turn to be completely unattainable. This issue is even more pronounced as the molecular anisotropy increases, leaving KH and OPT schemes as the only reliable approaches to study the phase separation properties of our models.

As for KH closed formulæ, pressures and free energies at different temperatures are also reported in Fig. 2. As visible, a van der Waals loop in the pressure is observed in the same range of temperatures of KH energy route and OPT, even if systematically higher and shallower estimates are obtained. The same discrepancy appears to affect the predictions for the free energy. As a result (see Fig. 2c), the liquid branch of phase coexistence appears too low in comparison with simulations, KH energy route and OPT, whereas the critical temperature and the whole gas branch turn to be quite well predicted. Overall, the KH energy route appears more accurate and therefore we shall adopt such a scheme — together with the differently based OPT approach — to investigate the phase behavior of other cases with .

As visible from Fig. 3, the agreement between KH and OPT predictions for the pressure essentially persists as is reduced down to . In Ref. Munao:14, we have documented the progressive shrinkage of the coexistence curves toward low temperatures as is reduced. As visible from Fig. 4, where coexistence curves for , 0.5, and are reported, this feature is reproduced by both theories, essentially at the same level of accuracy observed for the symmetric case: KH still provides reasonable estimates of the gas-liquid coexistence curve, with critical temperatures and densities fairly well predicted, whereas OPT overestimates the amplitude of the two-phase region and, accordingly, the critical temperature.

As documented in our recent simulation study,Munao:14 () increasing the asymmetry of total interaction promotes the development of aggregates in the fluid; specifically at , dumbbells self-assemble into micelles at low densities and planar structures (lamellæ) at intermediate/high densities, provided the temperature is low enough. Moreover the self-assembly process tends to destabilize the gas-liquid phase separation – completely disappearing for Janus dumbbells, i.e. with — with the coexistence curve stretched toward very low gas densities and a critical temperature drastically lower than in the symmetric case. At variance with previous cases, both KH and OPT theories fail to reproduce such a behavior, largely overestimating the coexistence region at , as visible in Fig. 5.

Indeed, we may surmise that liquid state theories, explicitly tailored to describe a homogeneous fluid environment, progressively worsen as local inhomogeneities tend to develop in the system. Also, the progressive displacement of binodal curves toward lower temperatures may negatively affect the OPT predictions based, as described, on a high-temperature expansion of the free energy.

1.0 | MC | 1.566 | 0.146 | -5.680 |
---|---|---|---|---|

KH | 1.616 | 0.135 | -5.973 | |

OPT | 1.842 | 0.106 | -6.549 | |

0.7 | MC | 1.213 | 0.146 | -4.524 |

KH | 1.252 | 0.133 | -4.610 | |

OPT | 1.430 | 0.105 | -5.025 | |

0.5 | MC | 0.980 | 0.143 | -3.725 |

KH | 1.017 | 0.130 | -3.762 | |

OPT | 1.162 | 0.103 | -4.107 | |

0.3 | MC | 0.739 | 0.125 | -2.997 |

KH | 0.794 | 0.126 | -2.964 | |

OPT | 0.905 | 0.097 | -3.248 | |

0.1 | MC | 0.466 | 0.068 | -2.635 |

KH | 0.591 | 0.116 | -2.284 | |

OPT | 0.682 | 0.080 | -2.205 |

Numerical values of KH, OPT and MC critical parameters for all investigated in this work are collectively reported in Tab. 1. Critical temperatures and densities as functions of are shown in Fig. 6. Remarkably, both theories agree with simulations in providing a linear dependence of the critical temperature on . Such a circumstance is compatible with a mean field behavior,Munao:14 () with the only exception (for MC data) of the case. Conversely, the critical density follows quite a different trend, keeping an almost constant value from the symmetric case down to , with a marked drop upon further lowering this parameter. Surprisingly enough, OPT qualitatively reproduces such a behavior, actually predicting a critical density close to simulation data for . On the other hand, the dependence of KH on is less pronounced, failing to reproduce in particular the very low values attained at .

In a recent study,Munao:PCCP () we have analyzed the RISM performances for a close class of SW dumbbells, characterized for instance by an attraction range , shorter than that investigated here (namely ). In that case we have shown that the HNC theory is able to describe the low temperature regime whereupon the gas-liquid coexistence takes place; in particular, the theory predicts the disappearance of a stable phase separation for , i.e. at a value slightly larger than that observed in simulations (). Also for RISM agrees with simulations in predicting the linear scaling of the critical temperature discussed in Fig. 6a, even if with a larger discrepancy with respect to the present case; a better agreement is instead observed for the critical density, essentially because this quantity, at , keeps an almost constant value as a function of , at variance with the sudden drop documented in Fig. 6b.

## Iv Conclusions

We have presented a theoretical study of thermodynamic properties and gas-liquid phase coexistence of anisotropic colloidal dumbbells, where the anisotropy stems from the fact that, although the two colloids forming the dumbbell are identical in sizes, their interactions are not, as one of the two attractive square-well depth is gradually reduced to zero toward the Janus dumbbells limit. We have discussed and contrasted to each other the performances of two venerable and well established tools, the RISM and OPT theories of liquids. As for the RISM, this theory hinges on the solution of the integral equations of liquids, thus providing the intermolecular correlations, whereby all thermodynamic properties can be obtained. Using the symmetric model as a test case, we have shown that the Kovalenko-Hirata closure (combining efficiently MSA and HNC approximations) complemented with the energy route from structure to thermodynamics is superior to the HNC approximation and other procedures to calculate the thermodynamic properties. As for the OPT theory, this method hinges on a fourth-order high-temperature perturbative expansion of the free energy of the square-well fluid, recently generalized to deal with molecular fluids.

The accuracy of both methods has been assessed by a direct comparison with recent numerical simulations. We have found that both approaches provide reasonable predictions, with RISM slightly outperforming OPT in almost all cases with the notable exception of the critical density, where OPT is more accurate. Latter predictions could be improved by finding a way to fully include in the formalism the intramolecular correlations, presently taken into account only at the level of the reference free energy.

We emphasize the importance of having a reliable theoretical tool to provide a quick estimate of thermodynamics of colloidal dumbbells. Recent experiments,kegel () as well as complementary numerical work,Munao:15 () show that in most of the interesting cases the two colloids forming the dumbbells are asymmetric both in size and interaction. This means that a huge number of possible combinations could be in principle studied to probe many topologically different resulting phase diagrams. A reliable theoretical approach would then be invaluable in pinning down the most interesting cases where a more sophisticated, and computationally more demanding, calculation could be performed.

###### Acknowledgements.

GM, DC, CC, FS, and AG gratefully acknowledge support from PRIN-MIUR 2010-2011 Project.## References

- (1) T. S. Skelhon, Y. Chen, and S. A. F. Bon, Soft Matter 10, 7730 (2014)
- (2) F. Ma, S. Wang, H. Zao, D. T. Wu, and N. Wu, Soft Matter 10, 8349 (2014)
- (3) D. J. Kraft, R. Ni, F. Smallenburg, M. Hermes, K. Yoon, D. A. Weitz, A. van Blaaderen, J. Groenewold, M. Dijkstra, and W. Kegel, Proceedings of the National Academy of Science (USA) 109, 10787 (2012)
- (4) D. Nagao, K. Goto, H. Ishii, and M. Konno, Langmuir 27, 13302 (2011)
- (5) K. Yoon, D. Lee, J. W. Kim, J. Kim, and D. A. Weitz, Chem. Comm. 48, 9056 (2012)
- (6) J. D. Forster, J. G. Park, M. Mittal, H. Noh, C. F. Schreck, C. S. O’Hern, H. Cao, E. M. Furst, and E. R. Dufresne, ACS Nano 5, 6695 (2011)
- (7) I. D. Hosein and C. M. Liddell, Langmuir 23, 10479 (2007)
- (8) G. A. Chapela, F. del Río, and J. Alejandre, J. Chem. Phys. 134, 224105 (2011)
- (9) G. A. Chapela and J. Alejandre, J. Chem. Phys. 135, 084126 (2011)
- (10) M. A. Miller, R. Blaak, C. N. Lumb, and J.-P. Hansen, J. Chem. Phys. 130, 114507 (2009)
- (11) P. Ilg and E. Del Gado, Soft Matter 7, 163 (2011)
- (12) S. H. Chong, A. J. Moreno, F. Sciortino, and W. Kob, Phys. Rev. Lett. 94, 215701 (2005)
- (13) G. Munaò, P. O’Toole, T. S. Hudson, D. Costa, C. Caccamo, A. Giacometti, and F. Sciortino, Soft Matter 10, 5269 (2014)
- (14) G. Munaò, P. O’Toole, T. S. Hudson, D. Costa, C. Caccamo, F. Sciortino, and A. Giacometti, J. Phys. Condens. Matter, in press (2015)
- (15) G. Avvisati, T. Vissers, and M. Dijkstra, J. Chem. Phys. 142, 084905 (2015)
- (16) F. Tu, B. J. Park, and D. Lee, Langmuir 29, 12679 (2013)
- (17) B. J. Park and D. Lee, ACS Nano 6, 782 (2012)
- (18) L. Hong, A. Cacciuto, E. Luijten, and S. Granick, Langmuir 24, 621 (2008)
- (19) C.-H. Chen, R. K. Shah, A. R. Abate, and D. A. Weitz, Langmuir 25, 4320 (2009)
- (20) B. S. Jiang, Q. Chen, M. Tripathy, E. Luijten, K. S. Schweizer, and S. Granick, Advanced Materials 22, 1060 (2010)
- (21) L. Hong, A. Cacciuto, E. Luijten, and S. Granick, Nano Letters 6, 2510 (2006)
- (22) Q. Chen, J. K. Whitmer, S. Jiang, S. C. Bae, E. Luijten, and S. Granick, Science 331, 199 (2011)
- (23) Q. Chen, J. Yan, J. Zhang, S. C. Bae, and S. Granick, Langmuir 28, 1355 (2012)
- (24) J. Yan, M. Bloom, S. C. Bae, E. Luijten, and S. Granick, Nature 491, 578 (2012)
- (25) F. Sciortino, A. Giacometti, and G. Pastore, Phys. Rev. Lett. 103, 237801 (Dec. 2009)
- (26) F. Sciortino, A. Giacometti, and G. Pastore, Phys. Chem. Chem. Phys. 12, 11869 (2010)
- (27) T. Vissers, Z. Preisler, F. Smallenburg, M. Dijkstra, and F. Sciortino, J. Chem. Phys. 138, 164505 (2013)
- (28) S. N. Fejer, D. Chakrabarti, H. Kusumaatmaja, and D. J. Wales, Nanoscale 6, 9448 (2014)
- (29) G. Munaò, D. Costa, A. Giacometti, C. Caccamo, and F. Sciortino, Phys. Chem. Chem. Phys. 15, 20590 (2013)
- (30) G. Munaò, D. Costa, and C. Caccamo, Chem. Phys. Lett. 470, 240 (2009)
- (31) G. Munaò, D. Costa, and C. Caccamo, J. Chem. Phys. 130, 144504 (2009)
- (32) M. Marechal, H. H. Goetzke, A. Härtel, and H. Löwen, J. Chem. Phys. 135, 234510 (2011)
- (33) D. Chandler and H. C. Andersen, J. Chem. Phys. 57, 1930 (1972)
- (34) R. Espíndola-Heredia, F. del Río, and A. Malijevský, J. Chem. Phys. 130, 024509 (2009)
- (35) F. Gámez, J. Chem. Phys. 140, 234504 (2014)
- (36) L. Li, K. Tang, L. Wu, W. Zhao, and J. Cai, J. Chem. Phys. 136, 214508 (2012)
- (37) J.-P. Hansen and I. R. McDonald, Theory of simple liquids, 3rd Ed. (Academic Press, New York, 2006)
- (38) L. J. Lowden and D. Chandler, J. Chem. Phys. 61, 5228 (1974)
- (39) L. Harnau, J.-P. Hansen, and D. Costa, Europhys. Lett. 53, 729 (2001)
- (40) D. Costa, J.-P. Hansen, and L. Harnau, Mol. Phys. 103, 1917 (2005)
- (41) J.-P. Hansen and C. Pearson, Mol. Phys. 104, 3389 (2006)
- (42) P. G. Khalatur, L. V. Zherenkova, and A. R. Khokhlov, J. Phys. II (France) 7, 543 (1997)
- (43) W. Kung, P. González-Mozuelos, and M. O. de la Cruz, Soft Matter 6, 331 (2010)
- (44) G. Munaò, D. Costa, F. Sciortino, and C. Caccamo, J. Chem. Phys. 134, 194502 (2011)
- (45) M. Tripathy and K. S. Schweizer, J. Phys. Chem. B 117, 373 (2013)
- (46) N. E. Valadez-Pérez, A. L. Benavides, E. Schöll-Passinger, and R. Castañeda-Priego, J. Chem. Phys. 137, 084905 (2012)
- (47) A. L. Benavides and F. Gámez, J. Chem. Phys. 135, 134511 (2011)
- (48) A. L. Benavides, L. A. Cervantes, and J. Torres, J. Phys. Chem. C 111, 16006 (2007)
- (49) J. R. Elliott and N. H. Gray, J. Chem. Phys. 123, 184902 (2005)
- (50) J. Cui and J. R. Elliott, J. Chem. Phys. 116, 8625 (2002)
- (51) G. A. Chapela, L. E. Scriven, and H. T. David, J. Chem. Phys. 91, 4307 (1989)
- (52) J. A. Barker and D. Henderson, J. Chem. Phys. 47, 2856 (1967)
- (53) J. A. Barker and D. Henderson, J. Chem. Phys. 47, 4714 (1967)
- (54) J. A. Barker and D. Henderson, Phys. Rev. A 1, 1266 (1970)
- (55) A. Giacometti, C. Gögelein, F. Lado, F. Sciortino, S. Ferrari, and G. Pastore, J. Chem. Phys. 140, 094104 (2014)
- (56) C. Gögelein, F. Romano, F. Sciortino, and A. Giacometti, J. Chem. Phys. 136, 094512 (2012)
- (57) A. Kovalenko and F. Hirata, J. Chem. Phys. 110, 10095 (1999)
- (58) A. Kovalenko and F. Hirata, Chem. Phys. Lett. 349, 496 (2001)
- (59) D. J. Tildesley and W. B. Streett, Mol. Phys. 41, 85 (1980)
- (60) M. S. Wertheim, J. Chem. Phys. 87, 7323 (1987)
- (61) D. Frenkel and B. Smit, Understanding molecular simulations (Academic, New York, 1996)
- (62) T. Morita and K. Hiroike, Progr. Theor. Phys. (Japan) 23, 1003 (1960)
- (63) S. J. Singer and D. Chandler, Mol. Phys. 55, 621 (1985)
- (64) L. W. Salvino and J. A. White, J. Chem. Phys. 96, 4559 (1992)
- (65) D. Meeter and P. Wolfe, Non-linear least squares (University of Wisconsin Computing Center, 1965)
- (66) D. J. Kraft, J. Groenewold, and W. K. Kegel, Soft Matter 5, 3823 (2009)