Vapor−Liquid and Vapor−Solid Phase Equilibria of Fullerenes: The Role of the Potential Shape on the Triple Point

Extension of the Aggregation-Volume-Bias Monte Carlo Method to the Calculation of Phase Properties of Solid Systems: A Lattice-Based Cluster Approach

The aggregation-volume-bias Monte Carlo method, which has been successful in the calculation of the formation free energies of liquid clusters, is extended to solid systems. This extension is motivated by early studies where disordered clusters are observed when the original method is applied at a temperature even far below the triple point. In order to avoid the formation of disordered aggregates, the insertion of particles is targeted directly toward those crystal lattice sites. Specifically, the insertion volume used to be defined as a spherical volume centered around a given target molecule is now restricted to be around each of the crystal lattice sites near a given target molecule. The free energies obtained for both liquid and solid clusters are then used to extrapolate bulk-phase information such as the chemical potential of the liquid and solid phases at coexistence. Using the temperature and pressure dependencies of the chemical potential information obtained for both liquid and solid phases, the location of the triple point can be determined. For Lennard-Jonesium, the results were found to be in good agreement with previous simulation studies using other approaches.

Towards Molecular Simulations that are Transparent, Reproducible, Usable By Others, and Extensible (TRUE)

Systems composed of soft matter (e.g., liquids, polymers, foams, gels, colloids, and most biological materials) are ubiquitous in science and engineering, but molecular simulations of such systems pose particular computational challenges, requiring time and/or ensemble-averaged data to be collected over long simulation trajectories for property evaluation. Performing a molecular simulation of a soft matter system involves multiple steps, which have traditionally been performed by researchers in a "bespoke" fashion, resulting in many published soft matter simulations not being reproducible based on the information provided in the publications. To address the issue of reproducibility and to provide tools for computational screening, we have been developing the open-source Molecular Simulation and Design Framework (MoSDeF) software suite. In this paper, we propose a set of principles to create Transparent, Reproducible, Usable by others, and Extensible (TRUE) molecular simulations. MoSDeF facilitates the publication and dissemination of TRUE simulations by automating many of the critical steps in molecular simulation, thus enhancing their reproducibility. We provide several examples of TRUE molecular simulations: All of the steps involved in creating, running and extracting properties from the simulations are distributed on open-source platforms (within MoSDeF and on GitHub), thus meeting the definition of TRUE simulations.

Liquid-vapor phase diagram and surface properties in oppositely charged colloids represented by a mixture of attractive and repulsive Yukawa potentials

The liquid-vapor phase diagrams of equal size diameter σ binary mixtures of screened potentials have been reported for several ranges of interaction using Monte Carlo simulation methods [J. B. Caballero, A. M. Puertas, A. Ferńandez-Barbero, F. J. de las Nieves, J. M. Romero-Enrique, and L. F. Rull, J. Chem. Phys. 124, 054909 (2006); A. Fortini, A.-P. Hynninen, and M. Dijkstra, J. Chem. Phys. 125, 094502 (2006)]. Both works report controversial results about the stability of the phase diagram with the inverse Debye screening length κ. Caballero found stability for values of κσ up to 20 while Fortini reported stability for κσ up to 20 while Fortini reported stability for κσ ≤ 4. In this work a spinodal decomposition process where the liquid and vapor phases coexist through an interface in a slab geometry is used to obtain the phase equilibrium and surface properties using a discontinuous molecular dynamics simulations for mixtures of equal size particles carrying opposite charge and interacting with a mixture of attractive and repulsive Yukawa potentials at different values of κσ. An crude estimation of the triple point temperatures is also reported. The isothermal-isobaric method was also used to determine the phase stability using one phase simulations. We found that liquid-vapor coexistence is stable for values of κσ > 20 and that the critical temperatures have a maximum value at around κσ = 10, in agreement with Caballero et al. calculations. There also exists a controversy about the liquid-vapor envelope stability of the pure component attractive Yukawa model which is also discussed in the text. In addition, details about the equivalence between continuous and discontinuous molecular dynamics simulations are given, in the Appendix, for Yukawa and Lennard-Jones potentials.

Fractal aggregates in protein crystal nucleation

Monte Carlo simulations of homogeneous nucleation for a protein model with an exceedingly short-ranged attractive potential yielded a nonconventional crystal nucleation mechanism, which proceeds by the formation of fractal, low-dimensional aggregates followed by a concurrent collapse and increase of the crystallinity of these aggregates to become compact ordered nuclei. This result corroborates a recently proposed two-step mechanism for protein crystal nucleation from solution.

Glass transition in fullerenes: Mode-coupling theory predictions

We report idealized mode-coupling theory results for the glass transition of ensembles of model fullerenes interacting via phenomenological two-body potentials. Transition lines are found for C60, C70, and C96 in the temperature-density plane. We argue that the observed glass transition behavior is indicative of kinetic arrest that is strongly driven by the interparticle attraction in addition to excluded-volume repulsion. In this respect, these systems differ from most standard glass-forming liquids. They feature arrest that occurs at lower densities and that is stronger than would be expected for repulsion-dominated hard-sphere-like or Lennard-Jones-type systems. The influence of attraction increases with increasing the number of carbon atoms per molecule. However, unrealistically large fullerenes would be needed to yield behavior reminiscent of recently investigated model colloids with strong short-ranged attraction (glass-glass transitions and logarithmic decay of time-correlation functions).

Direct calculation of solid-vapor coexistence points by thermodynamic integration: application to single component and binary systems

We present a new thermodynamic integration method that directly connects the vapor and solid phases by a reversible path. The thermodynamic integration in the isothermal-isobaric ensemble yields the Gibbs free energy difference between the two phases, from which the sublimation temperature can be easily calculated. The method extends to the binary mixture without any modification to the integration path simply by employing the isothermal-isobaric semigrand ensemble. The thermodynamic integration, in this case, yields the chemical potential difference between the solid and vapor phases for one of the components, from which the binary sublimation temperature can be calculated. The coexistence temperatures predicted by our method agree well with those in the literature for single component and binary Lennard-Jones systems.

A general perturbation approach for equation of state development: Applications to simple fluids, ab initio potentials, and fullerenes

A new perturbation scheme based on the Barker-Henderson perturbation theory [J. Chem. Phys. 47, 4714 (1967)] is proposed to predict the thermodynamic properties of spherical molecules. Accurate predictions of second virial coefficients and vapor-liquid coexistence properties are obtained for a large variety of potential functions (square well, Yukawa, Sutherland, Lennard-Jones, Buckingham, Girifalco). New Gibbs ensemble Monte Carlo simulations of the generalized exp-m Buckingham potential are reported. An extension of the perturbation approach to mixtures is proposed, and excellent predictions of vapor-liquid equilibria are obtained for Lennard-Jones mixtures. The perturbation scheme can be applied to complex potential functions fitted to ab initio data to predict the properties of real molecules such as neon. The new approach can also be used as an auxiliary tool in molecular simulation studies, to efficiently optimize an intermolecular potential on macroscopic properties or match force fields based on different potential functions.