This study proposes a new method for predicting the crystal-melt interfacial free energy (γ) using the Ginzburg-Landau (GL) model, enhanced by atomistic simulation data for more accurate density wave profiles. The analysis focuses on the soft-sphere system governed by an inverse power potential that stabilizes both BCC and FCC phases. Equilibrium molecular-dynamics simulations are used to obtain density-wave amplitude distributions, which serve as inputs for the GL model to predict γ and its anisotropy.
View Article and Find Full Text PDFThe study of interfacial properties in liquid-liquid and liquid-vapor systems has a history of nearly 200 years, with significant contributions from scientific luminaries such as Thomas Young and Willard Gibbs. However, a similar level of understanding of solid-liquid interfaces has emerged only more recently, largely because of the numerous complications associated with the thermodynamics needed to describe them. The accurate calculation of the interfacial free energy of solid-liquid systems is central to determining which interfaces will be observed and their properties.
View Article and Find Full Text PDFPer- and polyfluoroalkyl substances (PFAS) are emerging anthropogenic pollutants of concern and are associated with potential human and environmental health concerns. PFAS removal can be achieved using adsorbents such as activated carbon and ion exchange resins. Recently, zeolites have been identified as another potential adsorption technology with increased selectivity and product regenerability that is not currently achieved with other adsorbents.
View Article and Find Full Text PDFBy employing non-equilibrium molecular dynamics (NEMD) simulations and time-dependent Ginzburg-Landau (TDGL) theory for solidification kinetics [Cryst. Growth Des. 20, 7862 (2020)], we predict the kinetic coefficients of FCC(100) crystal-melt interface (CMI) of soft-spheres modeled with an inverse-sixth-power repulsive potential.
View Article and Find Full Text PDFThe microscopic stress field inhomogeneity in the interfacial region adjacent to the liquid surface is the fundamental origin of the liquid surface tension, but because of broadening due to capillary fluctuations, a detailed molecular level understanding of the stress field remains elusive. In this work, we deconvolute the capillary fluctuations to reveal the intrinsic stress field and show that the atomic-level contributions to the surface tension are similar in functional form across a variety of monatomic systems. These contributions are confined to an interfacial region approximately 1.
View Article and Find Full Text PDFThe probability of inserting, without overlap, a hard spherical particle of diameter σ in a hard-sphere fluid of diameter σ0 and packing fraction η determines its excess chemical potential at infinite dilution, μex(σ, η). In our previous work [R. L.
View Article and Find Full Text PDFAmorphous silica (a-SiO) surfaces, when grafted with select metals on the active sites of the functionalized surfaces, can act as useful heterogeneous catalysts. From a molecular modeling perspective, one challenge has been generating a-SiO slab models with controllable surface roughness to facilitate the study of the effect of surface morphology on the material properties. Previous computational methods either generate relatively flat surfaces or periodically corrugated surfaces that do not mimic the full range of potential surface roughness of the amorphous silica material.
View Article and Find Full Text PDFJ Chem Phys
October 2022
The equilibrium silica liquid-liquid interface between the high-density liquid (HDL) phase and the low-density liquid (LDL) phase is examined using molecular-dynamics simulation. The structure, thermodynamics, and dynamics within the interfacial region are characterized in detail and compared with previous studies on the liquid-liquid phase transition (LLPT) in bulk silica, as well as traditional crystal-melt interfaces. We find that the silica HDL-LDL interface exhibits a spatial fragile-to-strong transition across the interface.
View Article and Find Full Text PDFWe present a classical molecular-dynamics study of the collective dynamical properties of the coexisting liquid phase at equilibrium body-centered cubic (BCC) Fe crystal-melt interfaces. For the three interfacial orientations (100), (110), and (111), the collective dynamics are characterized through the calculation of the intermediate scattering functions, dynamical structure factors, and density relaxation times in a sequential local region of interest. An anisotropic speedup of the collective dynamics in all three BCC crystal-melt interfacial orientations is observed.
View Article and Find Full Text PDFThe excess chemical potential μ(σ, η) of a test hard spherical particle of diameter σ in a fluid of hard spheres of diameter σ and packing fraction η can be computed with high precision using Widom's particle insertion method [B. Widom, J. Chem.
View Article and Find Full Text PDFTo explore the curvature dependence of solid-fluid interfacial thermodynamics, we calculate, using Grand Canonical Monte Carlo simulation, the surface free energy for a 2d hard-disk fluid confined in a circular hard container of radius R as a function of the bulk packing fraction η and wall curvature C̄=-1/R. (The curvature is negative because the surface is concave.) Combining this with our previous data [Martin et al.
View Article and Find Full Text PDFCarbon-dioxide (CO)-expanded liquids, in which a significant mole fraction of CO is dissolved into an organic solvent, have been of significant interest, especially as catalytic support media. Because the CO mole fraction and density can be controlled over a significant range by changing the CO partial pressure, the transport properties of these solvents are highly tunable. Recently, these liquids have garnered interest as potential electrolyte solutions for catalytic electrochemistry; however, little is currently known about the influence of the electrolyte on CO expansion.
View Article and Find Full Text PDFJ Phys Chem B
September 2020
In this work, we examine the surface thermodynamics of a hard-disk fluid at curved hard walls using Monte Carlo (MC) simulation and a generalized scaled particle theory (gSPT). The curved walls are modeled as hard disks of varying radii, . The surface free energy, γ, and excess surface volume, , for this system are calculated as functions of both the fluid packing fraction and the wall radius.
View Article and Find Full Text PDFCarbon dioxide-expanded liquids (CXLs) represent an important class of reaction media that provide tunability of mass transport, solvation, and solubility. Their properties have been demonstrated to provide advantages over traditional organic solvents. However, the molecular-level effects of the CO expansion on the structure and dynamics of the liquid that lead to this result have not been fully explored.
View Article and Find Full Text PDFWe report molecular-dynamics (MD) simulation results for the surface free energy of a hard-sphere fluid at cylindrical and spherical hard walls of different radii. The precision of the results is much higher than that in our previous study [B. B.
View Article and Find Full Text PDFA generalized scaled-particle theory for the uniform hard-disk mixture is derived in the spirit of the White Bear II free energy of the hard-sphere fluid [H. Hansen-Goos and R. Roth, J.
View Article and Find Full Text PDFWe study the hard-sphere fluid in contact with a planar hard wall. By combining the inhomogeneous virial series with simulation results, we achieve a new benchmark of accuracy for the calculation of surface thermodynamics properties such as surface adsorption Γ and the surface free energy (or surface tension), γ. We briefly introduce the problem of choosing a position for the dividing surface and avoid it by proposing the use of alternative functions to Γ and γ that are independent of the adopted frame of reference.
View Article and Find Full Text PDFIn this work, we examine the effect of surface structure on the heterogeneous nucleation of Pb crystals from the melt at a Cu substrate using molecular-dynamics (MD) simulation. In a previous work [Palafox-Hernandez et al., Acta Mater.
View Article and Find Full Text PDFJ Phys Condens Matter
November 2016
For the molecular simulation of electric double layer capacitors (EDLCs), a number of methods have been proposed and implemented to determine the one-dimensional electric potential profile between the two electrodes at a fixed potential difference. In this work, we compare several of these methods for a model LiClO4-acetonitrile/graphite EDLC simulated using both the traditional fixed-charged method (FCM), in which a fixed charge is assigned a priori to the electrode atoms, or the recently developed constant potential method (CPM) (2007 J. Chem.
View Article and Find Full Text PDFJ Chem Phys
November 2014
A major challenge in the molecular simulation of electric double layer capacitors (EDLCs) is the choice of an appropriate model for the electrode. Typically, in such simulations the electrode surface is modeled using a uniform fixed charge on each of the electrode atoms, which ignores the electrode response to local charge fluctuations in the electrolyte solution. In this work, we evaluate and compare this Fixed Charge Method (FCM) with the more realistic Constant Potential Method (CPM), [S.
View Article and Find Full Text PDFWe examine the thermodynamics and intrinsic structure of the Al-Pb liquid-liquid interface using molecular dynamics simulation and embedded atom method potentials. The instantaneous interfacial positions, from which the intrinsic structure and the capillary fluctuation spectrum are determined, are calculated using a grid-based method. The interfacial free energy extracted from the capillary fluctuation spectrum is shown to be in excellent agreement with that calculated mechanically by integrating the stress profile.
View Article and Find Full Text PDFUsing molecular-dynamics simulation and Gibbs-Cahn Integration, we calculate the interfacial free energy γ of a binary hard-sphere fluid mixture at a structureless, planar hard wall. The calculation is performed as a function of packing fraction (density) for several values of the diameter ratio α = σ2∕σ1, where σ1 and σ2 are the diameters of the two particle types in the mixture. Our results are compared to those obtained from the bulk version of the White Bear Mark II (WBII) classical density-functional theory, which is a modification of the Fundamental-Measure Theory of Rosenfeld.
View Article and Find Full Text PDF