Effect of shape anisotropy on the precipitation of dimeric nanoparticles.

We use grand canonical transition matrix Monte Carlo simulations to study the precipitation of dimeric nanoparticles. The dimers are composed of two particles having different chemical features and separated by a fixed distance. The non-attractive and attractive parts of the dimer are modeled using hard-sphere and square-well potentials, respectively. The shape anisotropy is altered by changing the relative sizes of the two particles. We observe that the stability of the nanosuspension increases with the increase in the size of the non-attractive part of the dimer. The precipitates of dimers having larger non-attractive parts have lower packing densities, contain large cavities, and show evidence of self-assembly in the bulk and on the surface. We also use the results from our simulations and the classical nucleation theory to study the kinetics of precipitation. At a given temperature and relative supersaturation, the rate of homogeneous nucleation increases with the increase in the size of the non-attractive parts. Finally, we use an example to show how our results can guide the design of nanosuspensions containing chemically anisotropic dimers that are stable under particular conditions.

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. The interface potentials that emerge from these calculations provide direct access to spreading and drying coefficients. These macroscopic properties are then used to compute the liquid-vapor surface tension and the contact angle of a liquid droplet in contact with the substrate. The isothermal-isobaric ensemble provides a means to change the thickness of the fluid film adjacent to the substrate by modifying the volume of the simulation box. Molecular insertions and removals are not necessary. We introduce a framework for performing local volume change moves wherein one attempts to modify the density of a narrow region of the simulation box. We show that such moves improve the sampling efficiency of inhomogeneous systems. The approach is applied to a model system consisting of a monatomic Lennard-Jones fluid in the vicinity of a structureless substrate. Results are provided for direct spreading and drying interface potential calculations at several temperatures and substrate strengths. Expanded ensemble techniques are used to evaluate interfacial properties over a wide range of temperatures and substrate strengths. The results obtained using the isothermal-isobaric approach are compared with those previously obtained via a grand canonical approach.

Understanding the influence of capillary waves on solvation at the liquid-vapor interface.

This work investigates the question if surface capillary waves (CWs) affect interfacial solvation thermodynamic properties that determine the propensity of small molecules toward the liquid-vapor interface. We focus on (1) the evaluation of these properties from molecular simulations in a practical manner and (2) understanding them from the perspective of theories in solvation thermodynamics, especially solvent reorganization effects. Concerning the former objective, we propose a computational method that exploits the relationship between an external field acting on the liquid-vapor interface and the magnitude of CWs. The system considered contains the solvent, an externally applied field (f) and the solute molecule fixed at a particular location. The magnitude of f is selected to induce changes in CWs. The difference between the solvation free energies computed in the presence and in the absence of f is then shown to quantify the contribution of CWs to interfacial solvation. We describe the implementation of this method in the canonical ensemble by using a Lennard-Jones solvent and a non-ionic solute. Results are shown for three types of solutes that differ in the nature of short-ranged repulsive (hard-core) interactions. Overall, we observe that CWs have a negligible or very small effect on the interfacial solvation free energy of a solute molecule fixed near the liquid-vapor interface for the above systems. We also explain how the effects of pinning or dampening of CWs caused by a fixed solute are effectively compensated and do not contribute to the solvation free energy.

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. We employ Monte Carlo simulations and an interface-potential-based approach to calculate the liquid-vapor and substrate-fluid interfacial properties. Surface tensions for each fluid are evaluated over a range of temperatures that spans from a reduced temperature of approximately 0.6 to the critical point. Contact angles are calculated at select temperatures over a range of substrate-fluid interaction strengths that spans from the near-drying regime to the wetting regime. We observe that an increase in the relative strength of Coulombic interactions between ions leads to increasing deviation from Guggenheim's corresponding states theory. We show how this deviation is related to lower values of liquid-vapor excess entropies observed for strongly ionic fluids. Our results show that the qualitative nature of wetting behavior is significantly influenced by the competition between dispersion and electrostatic interactions. We discuss the influence of electrostatic interactions on the nature of wetting and drying transitions and corresponding states like behavior observed for contact angles. For all of the fluids studied, we observe a relatively narrow range of substrate-fluid interaction strengths wherein the contact angle is nearly independent of temperature. The influence of the ionic nature of the fluid on the temperature dependence of contact angle is also discussed.

Effect of Binding-Affinity and ATPase Activity on the Velocities of Kinesins Using Ratchet Models.

We use two-state ratchet models containing single and coupled Brownian motors to understand the role of motor-microtubule binding, ATPase reaction rate and dimerisation on the translational velocities of Kinesin motors. We use model parameters derived from the experimental measurements on KIF1A, KIF13A, KIF13B, and KIF16B motors to compute velocities in µm/s. We observe that both the models show the same trend in velocities (KIF1A > KIF13A > KIF13B > KIF16B) as the experimental results. However, the models significantly underpredict the velocities when compared with the experiments. The predictions of the coupled-motor model are closer to the experiments than those of the single-motor model. Our results indicate that the variation of ATPase reaction rate governs the trend in velocities for the above four motors. The variation of motor-microtubule binding affinity and the coupling strength between the motor domains may only have a secondary effect. More rigorous models that incorporate the power-stroke mechanism are necessary for better quantitative compliance with the experiments.

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. The former provides a direct measure of the spreading coefficient while the latter provides an analogous drying coefficient. When coupled with an independent measure of the liquid-vapor surface tension, these coefficients enable one to compute the contact angle. We also show how one can combine information gathered from application of the spreading and drying methods at a common state point to obtain direct measures of the contact angle and liquid-vapor surface tension. The computational strategies introduced here are applied to two model systems. One includes a monatomic Lennard-Jones fluid that interacts with a structureless substrate via a long-ranged substrate potential. The second model contains a monatomic Lennard-Jones fluid that interacts with an atomistically detailed substrate via a short-ranged potential. Expanded ensemble techniques are coupled with the interface potential approach to compile the temperature- and substrate strength-dependence of various interfacial properties for these systems. Overall, we find that the approach pursued here provides an efficient and precise means to calculate the wetting and drying properties of model systems.

Using Molecular Simulations to Understand the Effect of Dodecyl Sulfate on the Calcium-Binding Ability of Polystyrene Sulfonate.

We demonstrate the potential to tune the binding of calcium ions with polystyrene sulfonate (PSS) in the presence of dodecyl sulfate (DS). This can aid the design of surfactant-responsive water-softening agents for applications in detergency. We use molecular dynamics simulations to study the effect of the concentration of DS ions and the degree of sulfonation on the propensity of calcium ions toward PSS. We observe that the presence of DS ions increases the propensity of calcium ions toward 100% sulfonated PSS. The above phenomenon is due to the hydrophobic attraction between PSS and DS at low DS concentrations and the formation of calcium ion bridges between sulfonate and sulfate groups at moderate to high DS concentrations. We also observe the formation of calcium ion bridges between the sulfonate groups at high DS concentrations. The presence of DS ions also increases the propensity of calcium ions toward 20% sulfonated PSS. This is mainly due to the hydrophobic attraction between PSS and DS ions. The calcium ion bridges between sulfonate and sulfate groups are less prevalent than those of 100% sulfonated PSS. We do not observe calcium ion bridges between sulfonate groups of 20% sulfonated PSS. We use the above-mentioned insights to suggest potential strategies for designing an anionic polyelectrolyte having a suitable calcium-binding ability at a given concentration of the anionic surfactant. Finally, strong PSS-DS affinity is detrimental to the activity of surfactants because less surfactant ions are available for detergency. Our results also indicate the possibility of altering the PSS-DS affinity by changing the degree of sulfonation.

Multistep Fractionation of Coal and Application for Graphene Synthesis.

Despite its complex structure, coal has shown to be a promising precursor for graphene synthesis by chemical vapor deposition (CVD). However, the presence of heteroatoms and aliphatic chains in coal can lead to defects in the graphene lattice, preventing the formation of pristine graphene layers. Therefore, the goal of this study was to formulate a multistep coal fractionation scheme to extract and characterize the most aromatic fractions and explore their potential as graphene precursors. The scheme consisted of direct coal liquefaction under different conditions, Soxhlet extraction with heptane then toluene, and preparative liquid chromatography on silica gel using heptol solutions with different heptane/toluene ratios. The fractions obtained by this process were analyzed by proton nuclear magnetic resonance, thermogravimetric and elemental analyses, and automated SAR-AD (saturates, aromatics, resins-asphaltene determinator) separations. This characterization allowed the identification of two aromatic fractions with and without heteroatoms, which were subsequently used for graphene synthesis by CVD on nickel and copper foils. Raman spectrometry revealed that both fractions primarily formed defect-free multilayered graphene with approximately 11 layers on nickel due to the high solubility of carbon and the defect-healing effect of nickel. On the other hand, these fractions generated amorphous carbon on copper due to the high solubility of hydrogen in copper, which competed with carbon. Molecules in the more aromatic heteroatom-free fraction still contained alkyl pendant substituents and did not share the same planarity and symmetry to form defect-free graphene on copper. Thus, the quality of graphene was governed by the substrate on nickel and by the precursor quality on copper. When deposited directly on lacey carbon-coated copper grids of a transmission electron microscope, the heteroatom-free fraction gave rise to much larger graphene domains. The presence of heteroatoms promoted the formation of small self-assembled agglomerates of amorphous carbon.

Using Grand Canonical Monte Carlo Simulations to Understand the Role of Interfacial Fluctuations on Solvation at the Water-Vapor Interface.

The present work investigates the effect of interfacial fluctuations (predominantly capillary wave-like fluctuations) on the solvation free energy (Δµ) of a monatomic solute at the water-vapor interface. We introduce a grand-canonical-ensemble-based simulation approach that quantifies the contribution of interfacial fluctuations to Δµ. This approach is used to understand how the above contribution depends on the strength of dispersive and electrostatic solute-water interactions at the temperature of 400 K. At this temperature, we observe that interfacial fluctuations do play a role in the variation of Δµ with the strength of the electrostatic solute-water interaction. We also use grand canonical simulations to further investigate how interfacial fluctuations affect the propensity of the solute toward the water-vapor interface. To this end, we track a quantity called the interface potential (surface excess free energy) with the number of water molecules. With increasing number of water molecules, the liquid-vapor interface moves across a solute, which is kept at a fixed position in the simulation. Hence, the dependence of the interface potential on the number of waters models the process of moving the solute through the water-vapor interface. We analyze the change of the interface potential with the number of water molecules to explain that solute-induced changes in the interfacial fluctuations, like the pinning of capillary-wave-like undulations, do not play any role in the propensity of solutes toward water-vapor interfaces. The above analysis also shows that the dampening of interfacial fluctuations accompanies the adsorption of any solute at the liquid-vapor interface, irrespective of the chemical nature of the solute and solvent. However, such a correlation does not imply that dampening of fluctuations causes adsorption.

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. We also explore how the identity of the anion influences saturation properties. Specifically, we compare results for [C(4)mim][NTf2] with those for 1-butyl-3-methylimidazolium tetrafluoroborate ([C(4)mim][BF4]) and 1-butyl-3-methylimidazolium hexafluorophosphate ([C(4)mim][PF6]). Simulations are completed with a recently developed realistic united-atom force field. A combination of direct grand canonical and isothermal-isobaric temperature expanded ensemble simulations are used to construct phase diagrams. Our results are compared with experimental data and Gibbs ensemble simulation data. Overall, we find good agreement between our results and those measured experimentally. We find that the vapor pressures and enthalpies of vaporization show a strong dependence on the size of the alkyl chain at low temperatures, whereas no particular trend is observed at high temperatures. Finally, we also discuss the effect of temperature on liquid phase nanodomains observed in RTILs with large hydrophobic groups. We do not observe a drastic change in liquid phase structure upon variation of the temperature, which suggests there is not a sharp phase transition between a nanostructured and homogeneous liquid, as has been suggested in earlier studies.

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. Both direct methods, which enable one to locate saturation points at a given temperature, and temperature expanded ensemble methods, which provide a means to trace saturation lines to low temperature, are discussed. We study the liquid-vapor phase behavior of the restricted primitive model (RPM) and a realistic model for 1,3-dimethylimidazolium tetrafluoroborate ([C1mim][BF4]). Results are presented to show the dependence of saturation properties of the RPM and [C1mim][BF4] on the size of the simulation box and the boundary condition used for the Ewald summation. For [C1mim][BF4] we also demonstrate the ability of our strategy to sample ion clusters that form in the vapor phase. Finally, we provide the liquid-vapor saturation properties of these models over a wide range of temperature. Overall, we observe that the choice of system size and boundary condition have a non-negligible effect on the calculated properties, especially at high temperature. Also, we find that the combination of grand canonical MC simulation and isothermal-isobaric temperature expanded ensemble MC simulation provides a computationally efficient means to calculate liquid-vapor saturation properties of ionic liquids.

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. In the former case, we show how reservoir and growth expanded ensemble techniques can be used to facilitate the creation and insertion of complex molecules within a grand canonical simulation. We next focus on grand canonical and isothermal-isobaric temperature expanded ensemble techniques that provide a means to trace saturation lines over a wide range of temperatures. To demonstrate the utility of the strategies introduced here, we present phase coexistence data for a series of molecules, including n-octane, cyclohexane, water, 1-propanol, squalane, and pyrene. Overall, we find the direct grand canonical approach to be the most effective means to directly locate a coexistence point at a given temperature and the isothermal-isobaric temperature expanded ensemble scheme to provide the most effective means to follow a saturation curve to low temperature.