Modeling of many-body interactions between elastic spheres through symmetry functions

Simple models for spherical particles with a soft shell have been shown to self-assemble into numerous crystal phases and even quasicrystals. However, most of these models rely on a simple pairwise interaction, which is usually a valid approximation only in the limit of small deformations, i.e. low densities. In this work, we consider a many-body yet simple model for the evaluation of the elastic energy associated with the deformation of a spherical shell. The resulting energy evaluation, however, is relatively expensive for direct use in simulations. We significantly reduce the associated numerical cost by fitting the potential using a set of symmetry functions. We propose a method for selecting a suitable set of symmetry functions that capture the most relevant features of the particle environment in a systematic manner. The fitted interaction potential is then used in Monte Carlo simulations to draw the phase diagram of the system in two dimensions. The system is found to form both a fluid and a hexagonal crystal phase.Comment: 10 pages, 9 figure

Methane and carbon dioxide adsorption on edge-functionalized graphene: A comparative DFT study

With a view towards optimizing gas storage and separation in crystalline and disordered nanoporous carbon-based materials, we use ab initio density functional theory calculations to explore the effect of chemical functionalization on gas binding to exposed edges within model carbon nanostructures. We test the geometry, energetics, and charge distribution of in-plane and out-of-plane binding of CO2 and CH4 to model zigzag graphene nanoribbons edge-functionalized with COOH, OH, NH2, H2PO3, NO2, and CH3. Although different choices for the exchange-correlation functional lead to a spread of values for the binding energy, trends across the functional groups are largely preserved for each choice, as are the final orientations of the adsorbed gas molecules. We find binding of CO2 to exceed that of CH4 by roughly a factor of two. However, the two gases follow very similar trends with changes in the attached functional group, despite different molecular symmetries. Our results indicate that the presence of NH2, H2PO3, NO2, and COOH functional groups can significantly enhance gas binding with respect to a hydrogen-passivated edge, making the edges potentially viable binding sites in materials with high concentrations of edge carbons. To first order, in-plane binding strength correlates with the larger permanent and induced dipole moments on these groups. Implications for tailoring carbon structures for increased gas uptake and improved CO2/CH4 selectivity are discussed.Comment: 12 pages, 7 figure

Toward a molecular theory of homogeneous bubble nucleation in superheated liquids

Superheated liquids play an important role in various processes in industry, the laboratory and nature. Examples include hazardous vapor explosions, sonochemistry, ascent of tree sap, etc. Superheated liquids have finite lifetimes, after which they phase separate to the vapor phase by the process of bubble nucleation. This thesis describes an investigation of the molecular mechanisms of homogeneous bubble nucleation via the development of a rigorous molecular theory of bubble nucleation. The molecular theory proposed here demonstrates a potential for overcoming the shortcomings of earlier theories such as the classical nucleation theory. Within the framework of the proposed molecular theory, it was found that cavity formation played an important; if not dominant, role in the process of bubble nucleation. A consideration of various ideas set forth within the scaled particle theory of hard particle fluids suggests that critical cavity size exists for negative pressure super-heated liquids, leading to instability inside the liquid for larger cavity sizes. Monte Carlo simulations and density functional theory (DFT) studies of cavities inside a model superheated liquid verify the existence of a critical cavity size for both negative and positive pressure superheated liquids. DFT studies also show that the critical cavity in a true thermodynamic limit of stability. DFT calculations also show that the work of critical cavity formation is a tight upper bound to the work of forming of the critical bubble and the radius of the critical cavity is a lower bound to the radius of the critical bubble. This relationship between the critical cavity and the critical bubble provides a new insight into the molecular mechanism of bubble nucleation. Additional results from DFT calculations include universal scaling behavior of various quantities associated with the critical cavity across temperatures and intermolecular potentials just as was shown for quantities associated with the critical bubble. A new semi-empirical approach based on the scaled particle theory and the free energy perturbation methods that predicts the work of formation of a cavity inside a dense liquid is also presented. Finally, the generalized nucleation theorem is verified for an internally constrained one dimensional hard rod fluid

Note: Effect of fluid phase compositions on the formation of substitutionally ordered solid phases from a binary mixture of oppositely charged colloidal suspensions

A binary mixture of oppositely charged colloidal particles can self-assemble into either a substitutionally ordered or substitutionally disordered crystalline phase depending on the nature and strength of interactions among the particles. An earlier study had mapped out favorable inter-particle interactions for the formation of substitutionally ordered crystalline phases from a fluid phase using Monte Carlo molecular simulations along with the Gibbs-Duhem integration technique. In this paper, those studies are extended to determine the effect of fluid phase composition on formation of substitutionally ordered solid phases

A Gibbs-ensemble based technique for Monte Carlo simulation of electric double layer capacitors (EDLC) at constant voltage

Current methods for molecular simulations of Electric Double Layer Capacitors (EDLC) have both the electrodes and the electrolyte region in a single simulation box. This necessitates simulation of the electrode-electrolyte region interface. Typical capacitors have macroscopic dimensions where the fraction of the molecules at the electrode-electrolyte region interface is very low. Hence, large systems sizes are needed to minimize the electrode-electrolyte region interfacial effects. To overcome these problems, a new technique based on the Gibbs Ensemble is proposed for simulation of an EDLC. In the proposed technique, each electrode is simulated in a separate simulation box. Application of periodic boundary conditions eliminates the interfacial effects. This in addition to the use of constant voltage ensemble allows for a more convenient comparison of simulation results with experimental measurements on typical EDLCs. (C) 2014 AIP Publishing LLC

Determination of favorable inter-particle interactions for formation of substitutionally ordered solid phases from a binary mixture of oppositely charged colloidal suspensions

The solid phase formed by a binary mixture of oppositely charged colloidal particles can be either substitutionally ordered or substitutionally disordered depending on the nature and strength of interactions among the particles. In this work, we use Monte Carlo molecular simulations along with the Gibbs-Duhem integration technique to map out the favorable inter-particle interactions for the formation of substitutionally ordered crystalline phases from a fluid phase. The inter-particle interactions are modeled using the hard core Yukawa potential but the method can be easily extended to other systems of interest. The study obtains a map of interactions depicting regions indicating the type of the crystalline aggregate that forms upon phase transition

Calculation of three-phase methane-ethane binary clathrate hydrate phase equilibrium from Monte Carlo molecular simulations

Methane and ethane are the simplest hydrocarbon molecules that can form clathrate hydrates. Previous studies have reported methods for calculating the three-phase equilibrium using Monte Carlo simulation methods in systems with a single component in the gas phase. Here we extend those methods to a binary gas mixture of methane and ethane. Methane-ethane system is an interesting one in that the pure components form sII clathrate hydrate whereas a binary mixture of the two can form the sII clathrate. The phase equilibria computed from Monte Carlo simulations show a good agreement with experimental data and are also able to predict the sI-sII structural transition in the clathrate hydrate. This is attributed to the quality of the TIP4P/Ice and TRaPPE models used in the simulations. (C) 2014 Elsevier B.V. All rights reserved

Calculation of Chemical Potentials and Occupancies in Clathrate Hydrates through Monte Carlo Molecular Simulations

The flexibility of the water lattice in clathrate hydrates and guest-guest interactions has been shown in previous studies to significantly affect the values of the thermodynamic properties, such as chemical potentials and free energies. Here we describe methods for computing occupancies, chemical potentials, and free energies that account for the flexibility of water lattice and guest-guest interactions in the hydrate phase. The methods are validated for a wide variety of guest molecules, such as methane, ethane, carbon dioxide, and tetrahydrodfuran by comparing the predicted occupancy values of guest molecules with those obtained from isothermal isobaric semigrand Monte Carlo simulations. The proposed methods extend the van der Waals and Platteuw theory for clathrate hydrates, and the Langmuir constant is calculated based on the structure of the empty hydrate lattice. These methods in combination with development of advanced molecular models for water and guest molecules should lead to a more thermodynamically consistent theory for clathrate hydrates

Improving the Rigor and Consistency of the Thermodynamic Theory for Clathrate Hydrates through Incorporation of Movement of Water Molecules of Hydrate Lattice

Current applications of statistical thermodynamic theories for clathrate hydrates do not incorporate the translational and rotational movement of water molecules of the hydrate lattice,in a rigorous manner. Previous studies have shown that the movement of water molecules has a significant effect on the properties of clathrate hydrates. In this Article, a method is presented to incorporate the effect of water movement with as much rigor as possible. This method is then used to calculate the Langmuir constant of the guest species in a clathrate hydrate. Unlike previous studies on modeling of clathrate hydrate thermodynamics, the method presented in this paper does not regress either the intermolecular potentials or the properties of the empty hydrate from clathrate phase equilibria data. Also the properties of empty hydrate used in the theory do not depend on the nature and composition of the guest molecules. The predicted phase equilibria from the resulting theory are shown to be highly accurate and thermodynamically consistent by comparing them with the phase equilibria computed directly from molecular simulations

Free energy barriers for homogeneous crystal nucleation in a eutectic system of binary hard spheres

In this study, the free energy barriers for homogeneous crystal nucleation in a system that exhibits a eutectic point are computed using Monte Carlo simulations. The system studied is a binary hard sphere mixture with a diameter ratio of 0.85 between the smaller and larger hard spheres. The simulations of crystal nucleation are performed for the entire range of fluid compositions. The free energy barrier is found to be the highest near the eutectic point and is nearly five times that for the pure fluid, which slows down the nucleation rate by a factor of 10(-31). These free energy barriers are some of highest ever computed using simulations. For most of the conditions studied, the composition of the critical nucleus corresponds to either one of the two thermodynamically stable solid phases. However, near the eutectic point, the nucleation barrier is lowest for the formation of the metastable random hexagonal closed packed (rhcp) solid phase with composition lying in the two-phase region of the phase diagram. The fluid to solid phase transition is hypothesized to proceed via formation of a metastable rhcp phase followed by a phase separation into respective stable fcc solid phases