ABSTRACT
In the theory of solidification, the kinetic coefficient multiplies the local supercooling to give the solid-liquid interface velocity. The same coefficient should drive interface migration at the coexistence temperature in proportion to a curvature force. This work computes the ice-water kinetic coefficient from molecular simulations starting from a sinusoidal ice-water interface at the coexistence temperature. We apply this method to the basal and prismatic ice planes and compare results to previous estimates from equilibrium correlation functions and simulations at controlled supercooling.
ABSTRACT
Wettability of rock surfaces with respect to oil and water, which is characterized by the contact angle, is an important factor that determines the efficacy of enhanced oil recovery operations. Experimental determination of contact angles for oil-water-rock systems is expensive and time-consuming due to the extremely long times needed for the establishment of adsorption equilibrium at the liquid-solid interface. Hence, molecular simulations form an attractive tool for computing contact angles. In this work, we use the cleaving wall technique that was developed previously in our group [R. K. R. Addula and S. N. Punnathanam, J. Chem. Phys. 153, 154504 (2020)] to compute solid-liquid interfacial free energy, which is then combined with Young's equation to compute the oil-water contact angle on silica surfaces. The silica surface is modeled with the INTERFACE force field that has been developed to accurately reproduce experimental data. We have considered three different surface chemistries of silica, namely, Q2, Q3, and Q4, in this study. Our calculations reveal that while the Q2 and Q3 surfaces are completely wetted by water, the Q4 surface is partially non-wetted by water. All the simulations needed for this calculation can be performed using the Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) molecular package. This should facilitate wider adoption of the Young's equation route to compute contact angles for systems comprised of complex molecules.
ABSTRACT
In this Letter, we present a molecular theory of nucleation from dilute phases such as vapors or dilute solutions. The theory can model the nonclassical two-step crystal nucleation seen in many systems. When applied to study and analyze the crystal nucleation pathways from Lennard-Jones vapor, we find that prior explanations of the two-step mechanism based on lower barrier height for liquid nuclei is incomplete. The analysis from the molecular theory reveal that a complete explanation would also require consideration of anisotropy in the diffusion constants for growth of liquid droplets vis-á-vis the crystal nuclei.
ABSTRACT
In this article, we present two methods based on thermodynamic integration for computing solid-fluid interfacial free energy for a molecular system. As a representative system, we choose two crystal polymorphs of orcinol (5-methylbenzene-1,3-diol) as the solid phase and chloroform and nitromethane as the liquid phase. The computed values of the interfacial free energy are then used in combination with the classical nucleation theory to predict solvent induced polymorph selectivity during crystallization of orcinol from solution.
