Molecular exchange Monte Carlo: A generalized method for identity exchanges in grand canonical Monte Carlo simulations.

A generalized identity exchange algorithm is presented for Monte Carlo simulations in the grand canonical ensemble. The algorithm, referred to as molecular exchange Monte Carlo, may be applied to multicomponent systems of arbitrary molecular topology and provides significant enhancements in the sampling of phase space over a wide range of compositions and temperatures. Three different approaches are presented for the insertion of large molecules, and the pros and cons of each method are discussed. The performance of the algorithms is highlighted through grand canonical Monte Carlo histogram-reweighting simulations performed on a number of systems, which include methane+n-alkanes, butane+perfluorobutane, water+impurity, and 2,2,4-trimethylpentane+neopentane. Relative acceptance efficiencies for molecule transfers of up to 400 times that of standard configurational-bias Monte Carlo are obtained.

Shape-selective deposition and assembly of anisotropic nanoparticles.

We report the large-area assembly of anisotropic gold nanoparticles into lithographically defined templates with control over their angular position using a capillary force-based approach. We elucidate the role of the geometry of the templates in the assembly of anisotropic nanoparticles consisting of different shapes and sizes. These insights allow us to design templates that immobilize individual triangular nanoprisms and concave nanocubes in a shape-selective manner and filter undesired impurity particles from a mixture of triangular prisms and other polyhedra. Furthermore, by studying the assembly of two particles in the same template, we elucidate the importance of interparticle forces in this method. These advances allow for the construction of face-to-face and edge-to-edge nanocube dimers as well as triangular nanoprism bowtie antennas. As an example of the fundamental studies enabled by this assembly method, we investigate the surface-enhanced Raman scattering (SERS) of face-to-face concave cube dimers both experimentally and computationally and reveal a strong polarization dependence of the local field enhancement.

Toward a molecular theory of homogeneous bubble nucleation: I. Equilibrium embryo definition.

We propose an equilibrium embryo definition for homogeneous bubble nucleation within the pure component superheated Lennard-Jones liquid. The suggested embryo definition serves as an improvement to a previous (n,v) embryo model (Uline and Corti 2007). In that model, the constrained equilibrium between the bubble and the surrounding superheated liquid was maintained by placing n particles within a spherical volume v, without concern for the redundant counting of configurations and the relevance of the model to the dynamics of a nucleation process. To eliminate this redundancy, while only considering those embryos that should most likely appear at a transitional state, we now define the volume of the embryo via the use of a shell particle and only include n particles inside the volume that are deemed to be "vapor-like". The underlying free energy surface of formation of the new (n,v) embryo model is generated via Monte Carlo simulation and also approximately by a suggested density functional theory method. The resulting surface implies that small and locally confined regions of near-zero density serve as precursors initiating homogeneous bubble nucleation. Furthermore, the nonredundant counting of the embryo configurations yields a well-defined and sharp conduit in the free energy surface that directs the initially formed embryos toward the critical nucleus. We discuss how the suggested equilibrium embryo model aids in both the identification of the transitional configurations and calculation of the average number density of the critical nuclei within the superheated liquid phase, which is the focus of the companion paper (DOI 10.1021/jp404151h).

Toward a molecular theory of homogeneous bubble nucleation: II. Calculation of the number density of critical nuclei and the rate of nucleation.

In the present paper, we develop a method to calculate the rate of homogeneous bubble nucleation within a superheated L-J liquid based on the (n,v) equilibrium embryo free energy surface introduced in the first paper (DOI: 10.1021/jp404149n). We express the nucleation rate as the product of the concentration of critical nuclei within the metastable liquid phase and the relevant forward rate coefficient. We calculate the forward rate coefficient of the critical nuclei from their average lifetime as determined from MD simulations of a large number of embryo trajectories initiated from the transitional region of the metastable liquid configuration space. Therefore, the proposed rate coefficient does not rely on any predefined reaction coordinate. In our model, the critical nuclei belong to the region of the configuration space where the committor probability is about one-half, guaranteeing the dynamical relevance of the proposed embryos. One novel characteristic of our approach is that we define a limit for the configuration space of the equilibrium metastable phase and do not include the configurations that have zero committor probability in the nucleation free energy surface. Furthermore, in order to take into account the transitional degrees of freedom of the critical nuclei, we develop a simulation-based approach for rigorously mapping the free energy of the (n,v) equilibrium embryos to the concentration of the critical nuclei within the bulk metastable liquid phase.

On the interfacial thermodynamics of nanoscale droplets and bubbles.

We present a new self-consistent thermodynamic formalism for the interfacial properties of nanoscale embryos whose interiors do not exhibit bulklike behavior and are in complete equilibrium with the surrounding mother phase. In contrast to the standard Gibbsian analysis, whereby a bulk reference pressure based on the same temperature and chemical potentials of the mother phase is introduced, our approach naturally incorporates the normal pressure at the center of the embryo as an appropriate reference pressure. While the interfacial properties of small embryos that follow from the use of these two reference pressures are different, both methods yield by construction the same reversible work of embryo formation as well as consistency between their respective thermodynamic and mechanical routes to the surface tension. Hence, there is no a priori reason to select one method over another. Nevertheless, we argue, and demonstrate via a density-functional theory (with the local density approximation) analysis of embryo formation in the pure component Lennard-Jones fluid, that our new method generates more physically appealing trends. For example, within the new approach the surface tension at all locations of the dividing surface vanishes at the spinodal where the density profile spanning the embryo and mother phase becomes completely uniform (only the surface tension at the Gibbs surface of tension vanishes in the Gibbsian method at this same limit). Also, for bubbles, the location of the surface of tension now diverges at the spinodal, similar to the divergent behavior exhibited by the equimolar dividing surface (in the Gibbsian method, the location of the surface of tension vanishes instead). For droplets, the new method allows for the appearance of negative surface tensions (the Gibbsian method always yields positive tensions) when the normal pressures within the interior of the embryo become less than the bulk pressure of the surrounding vapor phase. Such a prediction, which is allowed by thermodynamics, is consistent with the interpretation that the mother phase's attempted compression of the droplet is counterbalanced by the negative surface tension, or free energy cost to decrease the interfacial area. Furthermore, for these same droplets, the surface of tension can no longer be meaningfully defined (the surface of tension always remains well defined in the Gibbsian method). Within the new method, the dividing surface at which the surface tension equals zero emerges as a new lengthscale, which has various thermodynamic analogs to and similar behavior as the surface of tension.

Homogeneous nucleation and growth in simple fluids. I. Fundamental issues and free energy surfaces of bubble and droplet formation.

The free energy of forming a droplet and a bubble with a given particle number n and volume v within the pure-component Lennard-Jones supercooled vapor and superheated liquid, respectively, are further explored using density-functional theory. Similar to what was found previously [M. J. Uline and D. S. Corti, Phys. Rev. Lett. 99, 076102 (2007); M. J. Uline and D. S. Corti, J. Chem. Phys. 129, 234507 (2008)], the limits of stability again appear within both free energy surfaces evaluated at two other metastability conditions, one closer to the binodal and one closer to the spinodal. Furthermore, an ad hoc bond connectivity criterion is also applied in an attempt, however approximately, to eliminate certain configurational redundancies that arise from the chosen droplet and bubble definitions. What results are free energy surfaces describing the formation of equilibrium embryos that should be an improved representation of the fluctuations that are relevant to those nonequilibrium embryos seen in an actual nucleation event. Finally, we discuss in some detail the use of the (n,v) reaction coordinate within the framework of an equilibrium-based theory and its relation to other descriptions of nucleation.

Homogeneous nucleation and growth in simple fluids. II. Scaling behavior, instabilities, and the (n,v) order parameter.

The free energy of forming a droplet and a bubble with a given number of particles n inside a volume v within the pure component Lennard-Jones supercooled vapor and superheated liquid, respectively, is further explored using density-functional theory. Certain key aspects of the free energy surface for bubble formation, such as the radius of the bubble at a stability limit, are found to scale in a nearly temperature independent manner when plotted versus a parameter that quantifies the location of the given state point in the metastable region. The corresponding work at this stability limit exhibits scaling for small values of n, but shows a strong temperature dependence for large n. No aspect of the free energy surface for droplet formation shows scaling over the full range of metastability conditions, including the work of forming the critical droplet and the radius of a droplet at its stability limit. Hence, there is no "universal" surface for embryo formation in metastable fluids. We also generate by thermodynamic arguments alone droplet and bubble trajectories along the corresponding free energy surfaces that avoid by construction the locus of instabilities, which match quite well the results obtained from other approaches. We also discuss in greater detail the use of the (n,v) order parameter within an equilibrium-based description of embryo formation, focusing on why the density profile of the embryo is found to be discontinuous at the embryo surface and why stability limits are expected to develop at certain bubble radii.

Molecular simulation study of cavity-generated instabilities in the superheated Lennard-Jones liquid.

Previous equilibrium-based density-functional theory (DFT) analyses of cavity formation in the pure component superheated Lennard-Jones (LJ) liquid [S. Punnathanam and D. S. Corti, J. Chem. Phys. 119, 10224 (2003); M. J. Uline and D. S. Corti, Phys. Rev. Lett. 99, 076102 (2007)] revealed that a thermodynamic limit of stability appears in which no liquidlike density profile can develop for cavity radii greater than some critical size (being a function of temperature and bulk density). The existence of these stability limits was also verified using isothermal-isobaric Monte Carlo (MC) simulations. To test the possible relevance of these limits of stability to a dynamically evolving system, one that may be important for homogeneous bubble nucleation, we perform isothermal-isobaric molecular dynamics (MD) simulations in which cavities of different sizes are placed within the superheated LJ liquid. When the impermeable boundary utilized to generate a cavity is removed, the MD simulations show that the cavity collapses and the overall density of the system remains liquidlike, i.e., the system is stable, when the initial cavity radius is below some certain value. On the other hand, when the initial radius is large enough, the cavity expands and the overall density of the system rapidly decreases toward vaporlike densities, i.e., the system is unstable. Unlike the DFT predictions, however, the transition between stability and instability is not infinitely sharp. The fraction of initial configurations that generate an instability (or a phase separation) increases from zero to unity as the initial cavity radius increases over a relatively narrow range of values, which spans the predicted stability limit obtained from equilibrium MC simulations. The simulation results presented here provide initial evidence that the equilibrium-based stability limits predicted in the previous DFT and MC simulation studies may play some role, yet to be fully determined, in the homogeneous nucleation and growth of embryos within metastable fluids.