Search

The growing trend towards engineering interfacial complexion (or phase) transitions has been seen in the grain boundary and solid surface systems.Meanwhile, little attention has been paid to the chemically heterogeneous solid/liquid interfaces. In this work, novel in-plane multi-interfacial states coexist within the Cu(111)/Pb(l) interface at a temperature just above the Pb freezing point is uncovered using atomistic simulations.Four monolayer interfacial states, i.e., two CuPb alloy liquids and two pre-freezing Pb solids, are observed coexisting within two interfacial layers sandwiched between the bulk solid Cu and bulk liquid Pb. Through computing the spatial variations of various properties along the direction normal to the in-plane solid-liquid boundary lines for both interfacial layers, a rich and varied picture depicting the inhomogeneity and anisotropy in the mechanical, thermodynamical, and dynamical properties is presented. The bulk values extracted from the in-plane profiles suggest that each interfacial state examined has distinct equilibrium values from each other and significantly deviates from those of the bulk solid and liquid phases, and indicate that the complexion (or phase) diagrams for the Cu(111)/Pb(l) interface bears a resemblance to that of the eutectic binary alloy systems, instead of the monotectic phase diagram for the bulk CuPb alloy. The reported data could support the development of interfacial complexion (or phase) diagrams and interfacial phase rules and provide a new guide for regulating heterogeneous nucleation and wetting processes., Comment: 32 pages, 7 figures

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. Calculations of dynamics properties reveal three types of dynamical heterogeneity hybridizing within the silica HDL-LDL interface. We also observe that as the interface is traversed from HDL to LDL, the Si/O coordination number ratio jumps to an unexpectedly large value, defining a thin region of the interface where HDL and LDL exhibit significant mixing. In addition, the LLPT phase coexistence is interpreted in the framework of the traditional thermodynamics of alloys and phase equilibria., Comment: 13 pages, 8 figures

We 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 speed up of the collective dynamics in all three BCC crystal-melt interfacial orientations is observed. This trend differs significantly different from the previously observed slowing down of the local collective dynamics at the liquid-vapor interface [Acta Mater 2020;198:281]. Examining the interfacial density relaxation times, we revisit the validity of the recently developed time-dependent Ginzburg-Landau (TDGL) theory for the solidification crystal-melt interface kinetic coefficients, resulting in excellent agreement with both the magnitude and the kinetic anisotropy of the CMI kinetic coefficients measured from the non-equilibrium MD simulations, Comment: 10pages, 6 figures

The kinetic mechanisms underlying bottom-up assembly of colloidal particles have been widely investigated in efforts to control crystallization pathways and to direct growth into targeted superstructures for applications including photonic crystals. Current work builds on recent progress in the development of kinetic theories for crystal growth of body-centered-cubic crystals in systems with short-range inter-particle interactions, accounting for a greater diversity of crystal structures and the role of the longer-ranged interactions and orientational degrees of freedom arising in polar systems. We address the importance of orientational ordering processes in influencing crystal growth in such polar systems, thus advancing the theory beyond the treatment of the translational ordering processes considered in previous investigations. The work employs comprehensive molecular-dynamics simulations that resolve key crystallization processes, and are used in the development of a quantitative theoretical framework based on ideas from time-dependent Ginzburg-Landau theory. The significant impact of orientational ordering on the crystallization kinetics could be potentially leveraged to achieve crystallization kinetics steering through external electric or magnetic fields. Our combined theory/simulation approach provides opportunities for future investigations of more complex crystallization kinetics., Comment: 28 pages, 10 figures, 8 tables

he inhomogeneous structure of a fluid at a wall can be characterized in several ways. Within a thermodynamic description the surface free energy $\gamma$ and the excess adsorption $\Gamma$ are of central importance. For theoretical studies closed expression of $\gamma$ and $\Gamma$ can be very valuable; however, even for a well-studied model system such as a hard-sphere fluid at a planar hard wall, the accuracy of existing expressions for $\gamma$ and $\Gamma$, compared to precise computer simulation data, can still be improved. Here, we compare several known expressions for $\gamma$ and $\Gamma$ to the most precise computer simulation data. While good agreement is generally found at low to intermediate fluid densities, the existing parameterizations show significant deviation at high density. In this work, we propose new parameterizations for $\gamma$ and $\Gamma$ that agree with the simulation data within statistical error over the entire fluid density range.

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 induced by charge fluctuations in the electrolyte. In this work, we evaluate and compare this Fixed Charge Method (FCM) with the more realistic Constant Potential Method (CPM), [Reed, et al., J. Chem. Phys., 126, 084704 (2007)], in which the electrode charges fluctuate in order to maintain constant electric potential in each electrode. For this comparison, we utilize a simplified LiClO$_4$-acetonitrile/graphite EDLC. At low potential difference ($\Delta\Psi\le 2V$), the two methods yield essentially identical results for ion and solvent density profiles; however, significant differences appear at higher $\Delta\Psi$. At $\Delta\Psi\ge 4V$, the CPM ion density profiles show significant enhancement (over FCM) of "partially electrode solvated" Li$^+$ ions very close to the electrode surface. The ability of the CPM electrode to respond to local charge fluctuations in the electrolyte is seen to significantly lower the energy (and barrier) for the approach of Li$^+$ ions to the electrode surface., Comment: Corrected typos

We report the observation of a premelting transition at chemically sharp solid-liquid interfaces using molecular-dynamics simulations. The transition is observed in the solid-Al/liquid-Pb system and involves the formation of a liquid interfacial film of Al with a width that grows logarithmically as the bulk melting temperature is approached from below, consistent with current theories of premelting. The premelting behavior leads to a sharp change in the temperature dependence of the diffusion coefficient in the interfacial region, and could have important consequences for phenomena such as particle coalescence and shape equilibration, which are governed by interfacial kinetic processes., Comment: 6 pages, 4 figures

Using molecular-dynamics simulation, we have calculated the interfacial free energy, \gamma, between a hard-sphere fluid and hard spherical and cylindrical colloidal particles, as functions of the particle radius R and the fluid packing fraction \eta= \rho\sigma^3/6, where \rho and \sigma are the number density and hard-sphere diameter, respectively. These results verify that Hadwiger's theorem from integral geometry, which predicts that \gamma for a fluid at a surface, with certain restrictions, should be a linear combination of the average mean and Gaussian surface curvatures, is valid within the precision of the calculation for spherical and cylindrical surfaces up to \eta about 0.42. In addition, earlier results for \gamma for this system [Bryk, et al., Phys. Rev. E, 68, 031602 (2003)] using a geometrically-based classical Density Functional Theory are in excellent agreement with the current simulation results for packing fractions in the range where Hadwiger's theorem is valid. However, above \eta about 0.42, \gamma(R) shows significant deviations from the Hadwiger form indicating limitations to its use for high-density hard-sphere fluids. Using the results of this study together with Hadwiger's theorem allows one, in principle, to determine $\gamma$ for any sufficiently smooth surface immersed in a hard-sphere fluid.

Extending to continuous potentials a cleaving wall molecular-dynamics simulation method recently developed for the hard-sphere system [Phys.Rev.Lett 85, 4751 (2000)], we calculate the crystal-melt interfacial free energies, $\gamma$, for a Lennard-Jones system as functions of both crystal orientation and temperature. At the triple point, T* = 0.617, the results are consistent with an earlier cleaving potential calculation by Broughton and Gilmer [J. Chem. Phys. {\bf 84}, 5759 (1986)], however, the greater precision of the current calculation allows us to accurately determine the anisotropy of $\gamma$. From our data we find that, at all temperatures studied, $\gamma_{111} < \gamma_{110} < \gamma_{100}$. Comparison is made to the results from our previous hard-sphere calculation and to recent results for Ni by Asta, Hoyt and Karma [Phys. Rev. B, 66 100101(R) (2002)]., Comment: 7 pages, 3 figures, 2 tables

Molecular dynamics simulations are performed to study the [100] and [111] orientations of the crystal-melt interface between an ordered two-component hard sphere with a NaCl structure and its coexisting binary hard-sphere fluid. The diameter ratio of the two types of hard spheres making up the mixture is taken to be 0.414. This work complements our earlier interface simulations [J. Chem. Phys.116, 3410] for the same diameter ratio at lower pressures where the smaller component is immiscible in the solid and the fluid mixture coexists with a pure FCC crystal of large particles. Density profiles and diffusion coefficient profiles are presented for the AB interfacial system. We find that for this system, the transition from crystal-like to fluid-like behavior of both the density and diffusion constant profiles occurs over a narrower region than that seen in our previous studies [J. Chem. Phys. 116, 3410] of the FCC/binary fluid system. But similar to what was found in the FCC/binary fluid interface the transition region for the large particle diffusion constant is shifted about the size of the large particles toward the fluid phase relative to that for the small particles., Comment: 8 pages

We present a general molecular-dynamics simulation scheme, based on the Nose' thermostat, for sampling according to arbitrary phase space distributions. We formulate numerical methods based on both Nose'-Hoover and Nose'-Poincare' thermostats for two specific classes of distributions; namely, those that are functions of the system Hamiltonian and those for which position and momentum are statistically independent. As an example, we propose a generalized variable temperature distribution that designed to accelerate sampling in molecular systems., Comment: 10 pages, 3 figures

We present a set of second-order, time-reversible algorithms for the isothermal (NVT) molecular-dynamics (MD) simulation of systems with mixed hard-core/continuous potentials. The methods are generated by combining real-time Nose' thermostats with our previously developed Collision Verlet algorithm [Mol. Phys. 98, 309 (1999)] for constant energy MD simulation of such systems. In all we present 5 methods, one based on the Nose'-Hoover [Phys. Rev. A 31, 1695 (1985)] equations of motion and four based on the Nose'-Poincare' [J.Comp.Phys., 151 114 (1999)] real-time formulation of Nose' dynamics. The methods are tested using a system of hard spheres with attractive tails and all correctly reproduce a canonical distribution of instantaneous temperature. The Nose'-Hoover based method and two of the Nose'-Poincare' methods are shown to have good energy conservation in long simulations., Comment: 9 pages, 5 figures

The crystal-melt interfaces of a binary hard-sphere fluid mixture in coexistence with a single-component hard-sphere crystal is investigated using molecular-dynamics simulation. In the system under study, the fluid phase consists of a two-component mixture of hard spheres of differing size, with a size ratio $\alpha=0.414$. At low pressures this fluid coexists with a pure fcc crystal of the larger particles in which the small particles are immiscible. For two interfacial orientations, [100] and [111], the structure and dynamics within the interfacial region is studied and compared with previous simulations on single component hard-sphere interfaces. Among a variety of novel properties, it is observed that as the interface is traversed from fluid to crystal the diffusion constant of the larger particle vanishes before that of the small particle defining a region of the interface where the large particles are frozen in their crystal lattice, but the small particles exhibit significant mobility. This behavior was not seen in previous binary hard-sphere interface simulations with less asymmetric diameters., Comment: To appear in J. Chem. Phys. Vol. 116, 22 February 2002

Largely due to its role in nucleation and crystal-growth, the free energy of the crystal-melt interfacial free energy is an object of considerable interest across a number of scientific disciplines, especially in the materials-, colloid- and atmospheric sciences. Over fifty years ago, Turnbull observed that the interfacial free energies (scaled by the mean interfacial area per particle) of a variety of metallic elements exhibit a linear correlation with the enthalpy of fusion. This correlation provides an important empirical "rule-of-thumb" for estimating interfacial free energies, but lacks a compelling physical explanation. In this work we show that the interfacial free energies for close-packed metals are linearly correlated with the melting temperature, and are therefore primarily entropic in origin. We also show that the slope of this linear relationship can be determined with quantitative accuracy using a hard-sphere model, and that the correlation with the enthalpy of fusion reported by Turnbull follows as a consequence of the fact that the entropy of fusion for close-packed metals is relatively constant., Comment: 3 pages, 1 figure, to appear in J. Chem. Phys

Discretizations of the Feynman-Kac path integral representation of the quantum mechanical density matrix are investigated. Each infinite-dimensional path integral is approximated by a Riemann integral over a finite-dimensional function space, by restricting the integration to a subspace of all admissible paths. Using this process, a wide class of methods can be derived, with each method corresponding to a different choice for the approximating subspace. The traditional ``short-time'' approximation and ``Fourier discretization'' can be recovered from this approach, using linear and spectral basis functions respectively. As an illustration, a novel method is formulated using cubic elements and is shown to have improved convergence properties when applied to a simple model problem., Comment: 4 pages, 2 figures, submitted to Physical Review Letters

We present a direct calculation by molecular-dynamics computer simulation of the crystal/melt interfacial free energy, $\gamma$, for a system of hard spheres of diameter $\sigma$. The calculation is performed by thermodynamic integration along a reversible path defined by cleaving, using specially constructed movable hard-sphere walls, separate bulk crystal and fluid systems, which are then merged to form an interface. We find the interfacial free energy to be slightly anisotropic with $\gamma$ = 0.62$\pm 0.01$, 0.64$\pm 0.01$ and 0.58$\pm 0.01 k_BT/\sigma^2$ for the (100), (110) and (111) fcc crystal/fluid interfaces, respectively. These values are consistent with earlier density functional calculations and recent experiments measuring the crystal nucleation rates from colloidal fluids of polystyrene spheres that have been interpreted [Marr and Gast, Langmuir {\bf 10}, 1348 (1994)] to give an estimate of $\gamma$ for the hard-sphere system of $0.55 \pm 0.02 k_BT/\sigma^2$, slightly lower than the directly determined value reported here., Comment: 4 pages, 4 figures, submitted to Physical Review Letters

Model interaction potentials for real materials are generally optimized with respect to only those experimental properties that are easily evaluated as mechanical averages (e.g., elastic constants (at T=0 K), static lattice energies and liquid structure). For such potentials, agreement with experiment for the non-mechanical properties, such as the melting point, is not guaranteed and such values can deviate significantly from experiment. We present a method for re-parameterizing any model interaction potential of a real material to adjust its melting temperature to a value that is closer to its experimental melting temperature. This is done without significantly affecting the mechanical properties for which the potential was modeled. This method is an application of Gibbs-Duhem integration [D. Kofke, Mol. Phys.78, 1331 (1993)]. As a test we apply the method to an embedded atom model of aluminum [J. Mei and J.W. Davenport, Phys. Rev. B 46, 21 (1992)] for which the melting temperature for the thermodynamic limit is 826.4 +/- 1.3K - somewhat below the experimental value of 933K. After re-parameterization, the melting temperature of the modified potential is found to be 931.5K +/- 1.5K., Comment: 9 pages, 5 figures, 4 tables

We present a new molecular-dynamics algorithm for integrating the equations of motion for a system of particles interacting with mixed continuous/impulsive forces. This method, which we call Impulsive Verlet, is constructed using operator splitting techniques similar to those that have been used successfully to generate a variety molecular-dynamics integrators. In numerical experiments, the Impulsive Verlet method is shown to be superior to previous methods with respect to stability and energy conservation in long simulations., Comment: 18 pages, 6 postscript figures, uses rotate.sty

We examine, using molecular dynamics simulation, the structure and thermodynamics of the (100) and (111) disordered face-centered cubic (FCC) crystal/melt interfaces for a binary hard-sphere system. This study is an extension of our previous work, [Phys. Rev. E 54, R5905 (1996)], in which preliminary data for the (100) interface were reported. Density and diffusion profiles on both fine- and course-grained scales are calculated and analyzed leading to the conclusion that equilibrium interfacial segregation is minimal in this system., Comment: 7 pages, 7 figures, to appear in Molecular Physics

In order to construct a general density-functional theory for nonuniform fluid mixtures, we propose an extension to multicomponent systems of the weighted-density approximation (WDA) of Curtin and Ashcroft [Phys. Rev. A 32, 2909 (1985)]. This extension corrects a deficiency in a similar extension proposed earlier by Denton and Ashcroft [Phys. Rev. A 42, 7312 (1990)], in that that functional cannot be applied to the multi-component nonuniform fluid systems with spatially varying composition, such as solid-fluid interfaces. As a test of the accuracy of our new functional, we apply it to the calculation of the freezing phase diagram of a binary hard-sphere fluid, and compare the results to simulation and the Denton-Ashcroft extension., Comment: 4 pages, 4 figures, to appear in Phys. Rev. E as Brief Report

We present a new algorithm for isothermal-isobaric molecular-dynamics simulation. The method uses an extended Hamiltonian with an Andersen piston combined with the Nos'e-Poincar'e thermostat, recently developed by Bond, Leimkuhler and Laird [J. Comp. Phys., 151, (1999)]. This Nos'e-Poincar'e-Andersen (NPA) formulation has advantages over the Nos'e-Hoover-Andersen approach in that the NPA is Hamiltonian and can take advantage of symplectic integration schemes, which lead to enhanced stability for long-time simulations. The equations of motion are integrated using a Generalized Leapfrog Algorithm and the method is easy to implement, symplectic, explicit and time reversible. To demonstrate the stability of the method we show results for test simulations using a model for aluminum., Comment: 7 pages