Efficiency Comparison of Single- and Multiple-Macrostate Grand Canonical Ensemble Transition-Matrix Monte Carlo Simulations.

Recent interest in parallelizing flat-histogram transition-matrix Monte Carlo simulations in the grand canonical ensemble, due to its demonstrated effectiveness in studying phase behavior, self-assembly and adsorption, has led to the most extreme case of single-macrostate simulations, where each macrostate is simulated independently with ghost particle insertions and deletions. Despite their use in several studies, no efficiency comparisons of these single-macrostate simulations have been made with multiple-macrostate simulations. We show that multiple-macrostate simulations are up to 3 orders of magnitude more efficient than single-macrostate simulations, which demonstrates the remarkable efficiency of flat-histogram biased insertions and deletions, even with low acceptance probabilities.

Coupled Monte Carlo and Molecular Dynamics Simulations on Interfacial Properties of Antifouling Polymer Membranes.

We use molecular simulation to study the wetting behavior of antifouling polymer-tethered membranes. We obtain the interfacial properties (e.g.

Improving the efficiency of Monte Carlo simulations of ions using expanded grand canonical ensembles.

While ionic liquids have promising applications as industrial solvents, predicting their fluid phase properties and coexistence remains a challenge. Grand canonical Monte Carlo simulation is an effective method for such predictions, but equilibration is hampered by the apparent requirement to insert and delete neutral sets of ions simultaneously in order to maintain charge neutrality. For relatively high densities and low temperatures, previously developed methods have been shown to be essential in improving equilibration by gradual insertion and deletion of these neutral sets of ions.

Construction of the interface potential from a series of canonical ensemble simulations.

We introduce a method to construct the interface potential from a series of molecular dynamics simulations conducted within the canonical ensemble. The interface potential provides the surface excess free energy associated with the growth of a fluid film from a surface. We collect the force that the fluid exerts on the surface (disjoining pressure) at a series of film thicknesses.

Application of the interface potential approach for studying wetting behavior within a molecular dynamics framework.

We introduce a means to implement the interface potential approach for computing wetting properties within a molecular dynamics framework. The general approach provides a means to determine the contact angle of a liquid droplet on a solid substrate in a mother vapor. We present a framework for implementing "spreading" and "drying" versions of the method within an isothermal-isobaric ensemble.

Effect of Carboxylic Acid on the Wetting Properties of a Model Water-Octane-Silica System.

Monte Carlo simulations are employed to determine the effects of acetic acid on the wetting properties of a model water-octane-silica system. We first compute the bulk liquid-vapor saturation properties of pure acetic acid and subsequently explore the bulk liquid-liquid saturation properties of the ternary water-octane-acid system. We introduce an expanded ensemble approach to compute the coexistence properties of the ternary system.

Using isothermal-isobaric Monte Carlo simulation to study the wetting behavior of model systems.

We introduce a molecular simulation method to compute the interfacial properties of model systems within the isothermal-isobaric ensemble. We use a free-energy-based approach in which Monte Carlo simulations are employed to obtain an interface potential associated with the growth of a fluid film from a solid substrate. The general method is implemented within "spreading" and "drying" frameworks.

Calculation of the Saturation Properties of a Model Octane-Water System Using Monte Carlo Simulation.

We use Monte Carlo simulation to compute the saturation properties of a model octane-water system. The system phase separates into water-rich liquid and vapor phases, octane-rich liquid and vapor phases, and water-rich liquid and octane-rich liquid phases at various conditions. We outline a strategy for determining the saturation properties of the mixture over a wide range of temperatures, pressures, and compositions.

Multivariable extrapolation of grand canonical free energy landscapes.

We derive an approach for extrapolating the free energy landscape of multicomponent systems in the grand canonical ensemble, obtained from flat-histogram Monte Carlo simulations, from one set of temperature and chemical potentials to another. This is accomplished by expanding the landscape in a Taylor series at each value of the order parameter which defines its macrostate phase space. The coefficients in each Taylor polynomial are known exactly from fluctuation formulas, which may be computed by measuring the appropriate moments of extensive variables that fluctuate in this ensemble.

Connection Between Thermodynamics and Dynamics of Simple Fluids in Pores: Impact of Fluid-Fluid Interaction Range and Fluid-Solid Interaction Strength.

Using molecular simulations, we investigate how the range of fluid-fluid (adsorbate-adsorbate) interactions and the strength of fluid-solid (adsorbate-adsorbent) interactions impact the strong connection between distinct adsorptive regimes and distinct self-diffusivity regimes reported in [Krekelberg, W. P.; Siderius, D.

Position-Dependent Dynamics Explain Pore-Averaged Diffusion in Strongly Attractive Adsorptive Systems.

Using molecular simulations, we investigate the relationship between the pore-averaged and position-dependent self-diffusivity of a fluid adsorbed in a strongly attractive pore as a function of loading. Previous work (Krekelberg, W. P.

Temperature extrapolation of multicomponent grand canonical free energy landscapes.

We derive a method for extrapolating the grand canonical free energy landscape of a multicomponent fluid system from one temperature to another. Previously, we introduced this statistical mechanical framework for the case where kinetic energy contributions to the classical partition function were neglected for simplicity [N. A.

Predicting low-temperature free energy landscapes with flat-histogram Monte Carlo methods.

We present a method for predicting the free energy landscape of fluids at low temperatures from flat-histogram grand canonical Monte Carlo simulations performed at higher ones. We illustrate our approach for both pure and multicomponent systems using two different sampling methods as a demonstration. This allows us to predict the thermodynamic behavior of systems which undergo both first order and continuous phase transitions upon cooling using simulations performed only at higher temperatures.

Understanding the influence of Coulomb and dispersion interactions on the wetting behavior of ionic liquids.

We study the role of dispersion and electrostatic interactions in the wetting behavior of ionic liquids on non-ionic solid substrates. We consider a simple model of an ionic liquid consisting of spherical ions that interact via Lennard-Jones and Coulomb potentials. Bulk and interfacial properties are computed for five fluids distinguished by the strength of the electrostatic interaction relative to the dispersion interaction.

Saturation properties of 1-alkyl-3-methylimidazolium based ionic liquids.

We study the liquid-vapor saturation properties of room temperature ionic liquids (RTILs) belonging to the homologous series 1-alkyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide ([Cnmim][NTf2]) using Monte Carlo simulation. We examine the effect of temperature and cation alkyl chain length n on the saturated densities, vapor pressures, and enthalpies of vaporization. These properties are explicitly calculated for temperatures spanning from 280 to 1000 K for RTILs with n = 2, 4, 6, 8, 10, and 12.

Predicting how nanoconfinement changes the relaxation time of a supercooled liquid.

The properties of nanoconfined fluids can be strikingly different from those of bulk liquids. A basic unanswered question is whether the equilibrium and dynamic consequences of confinement are related to each other in a simple way. We study this question by simulation of a liquid comprising asymmetric dumbbell-shaped molecules, which can be deeply supercooled without crystallizing.

Connection between thermodynamics and dynamics of simple fluids in highly attractive pores.

Using molecular simulations, we investigate the structural and diffusive dynamics properties of a model fluid in highly absorptive cylindrical pores. At subcritical temperatures, self-diffusivity displays three distinct regimes as a function of average pore density ρ: (1) a decrease in self-diffusivity with increasing ρ at low ρ, (2) constant self-diffusivity with respect to varying ρ at moderate density, and (3) a decrease in self-diffusivity with increasing ρ at high density. These regimes are closely linked to the thermodynamic properties of the fluid in the pore, specifically, the adsorption isotherm, isosteric heat of adsorption, and the density profile.

Communication: Phase behavior of materials with isotropic interactions designed by inverse strategies to favor diamond and simple cubic lattice ground states.

We use molecular simulation to construct equilibrium phase diagrams for two recently introduced model materials with isotropic, soft-repulsive pair interactions designed to favor diamond and simple cubic lattice ground states, respectively, over a wide range of densities [Jain et al., Soft Matter 9, 3866 (2013)]. We employ free energy based Monte Carlo simulation techniques to precisely trace the inter-crystal and fluid-crystal coexistence curves.

Impact of small-scale geometric roughness on wetting behavior.

We examine the extent to which small-scale geometric substrate roughness influences the wetting behavior of fluids at solid surfaces. Molecular simulation is used to construct roughness wetting diagrams wherein the progression of the contact angle is traced from the Cassie to Wenzel to impregnation regime with increasing substrate strength for a collection of systems with rectangularly shaped grooves. We focus on the evolution of these diagrams as the length scale of the substrate features approaches the size of a fluid molecule.

Understanding wetting of immiscible liquids near a solid surface using molecular simulation.

We introduce Monte Carlo simulation methods for determining interfacial properties of fluids that exhibit bulk liquid-liquid immiscibility. An interface-potential-based approach, in which the interfacial properties of a system are related to the surface excess free energy of a thin fluid film in contact with a surface, is utilized to deduce the wetting characteristics of these systems. We present a framework for implementing this general method within both the grand canonical and semigrand isobaric-isothermal ensembles.

Monte Carlo Simulation Methods for Computing Liquid-Vapor Saturation Properties of Model Systems.

We discuss molecular simulation methods for computing the phase coexistence properties of complex molecules. The strategies that we pursue are histogram-based approaches in which thermodynamic properties are related to relevant probability distributions. We first outline grand canonical and isothermal-isobaric methods for directly locating a saturation point at a given temperature.

Using Monte Carlo simulation to compute liquid-vapor saturation properties of ionic liquids.

We discuss Monte Carlo (MC) simulation methods for calculating liquid-vapor saturation properties of ionic liquids. We first describe how various simulation tools, including reservoir grand canonical MC, growth-expanded ensemble MC, distance-biasing, and aggregation-volume-biasing, are used to address challenges commonly encountered in simulating realistic models of ionic liquids. We then indicate how these techniques are combined with histogram-based schemes for determining saturation properties.

Monte Carlo simulation strategies to compute interfacial and bulk properties of binary fluid mixtures.

We introduce Monte Carlo simulation methods for determining interfacial properties of binary fluid mixtures. The interface potential approach, in which the interfacial properties of a system are related to the surface excess free energy of a thin fluid film in contact with a surface, is utilized to deduce the wetting characteristics of a fluid mixture. The strategy described here provides an effective means to obtain the evolution of interfacial properties with the chemical composition of the fluid.

Calculation of inhomogeneous-fluid cluster expansions with application to the hard-sphere∕hard-wall system.

We examine the suitability of cluster expansion methods for the description of inhomogeneous fluids. In particular, we apply these methods to characterize the density profile, surface tension, and excess adsorption for a hard-sphere fluid near a hard wall. Coefficients for these series up to seventh order are evaluated by the Mayer-sampling Monte Carlo method.

Monte Carlo simulation methods for computing the wetting and drying properties of model systems.

We introduce general Monte Carlo simulation methods for determining the wetting and drying properties of model systems. We employ an interface-potential-based approach in which the interfacial properties of a system are related to the surface excess free energy of a thin fluid film in contact with a surface. Two versions of this approach are explored: a "spreading" method focused on the growth of a thin liquid film from a surface in a mother vapor and a "drying" method focused on the growth of a thin vapor film from a surface in a mother liquid.