Abstract. In this review we focus on the determination of phase diagrams by computer simulation with particular attention to the fluid-solid and solid-solid equilibria. The methodology to compute the free energy of solid phases will be discussed. In particular, the Einstein crystal and Einstein molecule methodologies are described in a comprehensive way. It is shown that both methodologies yield the same free energies and that free energies of solid phases present noticeable finite size effects. In fact this is the case for hard spheres in the solid phase. Finite size corrections can be introduced, although in an approximate way, to correct for the dependence of the free energy on the size of the system. The computation of free energies of solid phases can be extended to molecular fluids. The procedure to compute free energies of solid phases of water (ices) will be described in detail. The free energies of ices Ih, II, III, IV, V, VI, VII, VIII, IX, XI and XII will be presented for the SPC/E and TIP4P models of water. Initial coexistence points leading to the determination of the phase diagram of water for these two models will be provided. Other methods to estimate the melting point of a solid, as the direct fluid-solid coexistence or simulations of the free surface of the solid will be discussed. It will be shown that the melting points of ice Ih for several water models, obtained from free energy calculations, direct coexistence simulations and free surface simulations, agree within their statistical uncertainty. Phase diagram calculations can indeed help to improve potential models of molecular fluids. For instance, for water, the potential model TIP4P/2005 can be regarded as an improved version of TIP4P. Here we will review some recent work on the phase diagram of the simplest ionic model, the restricted primitive model. Although originally devised to describe ionic liquids, the model is becoming quite popular to describe the behaviour of charged colloids. Besides the possibility of obtaining fluid-solid equilibria for simple protein models will be discussed. In these primitive models, the protein is described by a spherical potential with certain anisotropic bonding sites (patchy sites).

Surfactant-assisted seeded growth of metal nanoparticles (NPs) can be engineered to produce anisotropic gold nanocrystals with high chiroptical activity through the templating effect of chiral micelles formed in the presence of dissymmetric cosurfactants. Mixed micelles adsorb on gold nanorods, forming quasihelical patterns that direct seeded growth into NPs with pronounced morphological and optical handedness. Sharp chiral wrinkles lead to chiral plasmon modes with high dissymmetry factors (~0.20). Through variation of the dimensions of chiral wrinkles, the chiroptical properties can be tuned within the visible and near-infrared electromagnetic spectrum. The micelle-directed mechanism allows extension to other systems, such as the seeded growth of chiral platinum shells on gold nanorods. This approach provides a reproducible, simple, and scalable method toward the fabrication of NPs with high chiral optical activity.

We systematically study the design of simple patchy sphere models that reversibly self-assemble into monodisperse icosahedral clusters. We find that the optimal patch width is a compromise between structural specificity (the patches must be narrow enough to energetically select the desired clusters) and kinetic accessibility (they must be sufficiently wide to avoid kinetic traps). Similarly, for good yields the temperature must be low enough for the clusters to be thermodynamically stable, but the clusters must also have enough thermal energy to allow incorrectly formed bonds to be broken. Ordered clusters can form through a number of different dynamic pathways, including direct nucleation and indirect pathways involving large disordered intermediates. The latter pathway is related to a reentrant liquid-to-gas transition that occurs for intermediate patch widths upon lowering the temperature. We also find that the assembly process is robust to inaccurate patch placement up to a certain threshold and that it is possible to replace the five discrete patches with a single ring patch with no significant loss in yield.

The ability to control the crystallization behaviour (including its absence) of particles, be they biomolecules such as globular proteins, inorganic colloids, nanoparticles, or metal atoms in an alloy, is of both fundamental and technological importance. Much can be learnt from the exquisite control that biological systems exert over the behaviour of proteins, where protein crystallization and aggregation are generally suppressed, but where in particular instances complex crystalline assemblies can be formed that have a functional purpose. We also explore the insights that can be obtained from computational modelling, focussing on the subtle interplay between the interparticle interactions, the preferred local order and the resulting crystallization kinetics. In particular, we highlight the role played by "frustration", where there is an incompatibility between the preferred local order and the global crystalline order, using examples from atomic glass formers and model anisotropic particles.

In this paper a new method to evaluate the free energy of solids is proposed. The method can be regarded as a variant of the method proposed by Frenkel and Ladd [J. Chem. Phys. 81, 3188 (1984)]. The main equations of the method can be derived in a simple way. The method can be easily implemented within a Monte Carlo program. We have applied the method to determine the free energy of hard spheres in the solid phase for several system sizes. The obtained free energies agree within the numerical uncertainty with those obtained by Polson et al. [J. Chem. Phys. 112, 5339 (2000)]. The fluid-solid equilibria has been determined for several system sizes and compared to the values published previously by Wilding and Bruce [Phys. Rev. Lett. 85, 5138 (2000)] using the phase switch methodology. It is shown that both the free energies and the coexistence pressures present a strong size dependence and that the results obtained from free energy calculations agree with those obtained using the phase switch method, which constitutes a cross-check of both methodologies. From the results of this work we estimate the coexistence pressure of the fluid-solid transition of hard spheres in the thermodynamic limit to be p*=11.54(4), which is slightly lower than the classical value of Hoover and Ree (p*=11.70) [J. Chem. Phys. 49, 3609 (1968)]. Taking into account the strong size dependence of the free energy of the solid phase, we propose to introduce finite size corrections, which allow us to estimate approximately the free energy of the solid phase in the thermodynamic limit from the known value of the free energy of the solid phase with N molecules. We have also determined the free energy of a Lennard-Jones solid by using both the methodology of this work and the finite size correction. It is shown how a relatively good estimate of the free energy of the system in the thermodynamic limit is obtained even from the free energy of a relatively small system.